《2022年fortran常见问题解决.pdf》由会员分享，可在线阅读，更多相关《2022年fortran常见问题解决.pdf（9页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、为了减少重复回答问题，特编此帖，并不定期添加和更新内容。错误难免，欢迎讨论，仅供参考。很多人问哪里可以找到Fortran编译器， 有不少热心学友提供网址，特汇集在这里。 虽然俺检验过这些链接，但是它们不一定总有效。Fortran 编译器下载：CVF FTN95 (License: Free for personal use) 以下操作，如无特别说明，都是以为例。1. 如何加大 Stack size?选 Project = Settings = Link = Category: Output = Stack allocations Reserve: 这里填新值（默认为1M ，若需要10M ，则填
2、）2. 如何用 Fortran批量生成文件？设要生成4000 个文件，文件名为AA1-AA4000 ， 如何写循环生成文件, 而不用写4000 次 write命令呢 ? 用内部文件：character(len=80) : filename,form integer : i do i=1,4000 select case (i) case (1:9) write(form,(i1) i case (10:99) write(form,(i2) i 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 1 页，共
3、9 页 - - - - - - - - - - case (100:999) write(form,(i3) i case (1000:9999) write(form,(i4) i end select write(filename,*) AA,trim(form),.TXT open(10,file=filename) write(10,*) i close(10) end do stop end 3. 如何用 Fortran动态生成输出格式？设有一个数组data(100)， 输出时， 希望每行输出num个数，而 num由用户输入， 如何实现？用内部文件：character(len=80)
4、 : form real : data(100) integer : i,num data = (/ (i,i=1,100) /)/ read(*,*) num write(form,*) (,num, write(*,form) data stop end 4. MS 是不是很垃圾？是垃圾，其中Bug 太多，多到不可用的地步！在这个主题里，换了CVF后问题就没了的人已有相当的数目。如果你用，遇到莫名其妙的错误，建议换，这是一个比较成熟的编译器。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 2 页，共
5、9 页 - - - - - - - - - - 5. 如何用 F90/95 生成随机数？注意：现在计算机产生的随机数都是伪随机数。random_number(x) 产生一个0到 1 之间的随机数（x 可以是向量），但是每次总是那几个数。用了 random_seed () 后，系统根据日期和时间随机地提供种子，使得随机数更随机了。program random implicit none real : x call random_seed () ! 系统根据日期和时间随机地提供种子call random_number (x) ! 每次的随机数就都不一样了write(*,*) x stop end
6、program random 6. 函数 /子程序超载的例子设要编一个两个变量值互换的子程序swap(a,b) ，哑元 a,b 可能是实型数，整型数，数组，矩阵，字符串，派生类型等等。但是希望只用一个子程序接口swap(a,b) 来实现。 F90 可以用类属接口来实现这种子程序超载：module Utilities implicit none private I_Swap,R_Swap,RVec_Swap,RMat_Swap,Type_Swap public : Swap interface Swap module procedure I_Swap,R_Swap,RVec_Swap,RMat_
7、Swap,Type_Swap end interface contains subroutine i_swap (a,b) ! 整型置换精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 3 页，共 9 页 - - - - - - - - - - integer (ikind),intent(in out) : a,b integer (ikind) : t 。! 略end subroutine i_swap subroutine r_swap (a,b) ! 实型置换real (rkind), intent
8、(in out) : a,b real (rkind) : t t = a a = b b = t return end subroutine r_swap subroutine RVec_swap (a,b) ! 实型向量置换real (rkind), intent(in out) : a(:),b(:) integer (ikind) : i do i=1, size(a) call R_Swap (a(i),b(i) end do return end subroutine RVec_swap subroutine RMat_swap (a,b) ! 实型矩阵置换。! 略end subr
9、outine RMat_swap subroutine Type_swap (a,b) ! 派生类型置换。! 略end subroutine Type_swap end module Utilities 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 4 页，共 9 页 - - - - - - - - - - 7. 如何在 CVF中为多行代码加注释?见：8. 推荐好的代码风格根据 F90 子集语言ELF90和 F 的要求整理（部分）。“强迫用” 的语言特性：+ F90 的自由格式的源代码。+ implic
10、it none。+ 子过程的哑元都要有intent属性。+ 函数子程序的哑元必须指定为intent(in)。+ 所有子程序和函数都放在模块（module）中，然后引用 （use）该模块； 或者放在program中。+ 数组哑元要求是假定形状的，或者有固定的维数和大小。字符哑元要求是假定长度的。+ 对于 recursive function（递归函数）语句，必须有result子句。+ 在所有派生类型（type ）的定义语句中，必须用双冒号分隔符（: ）。+ 主程序要求有program 语句。+ 在程序单元的end 语句中要求后跟程序单元的类型和名称。+ 在 end type 语句中要求后跟类型的
11、名称。+ end program 前必须有stop 语句以表示停止执行。+ 子过程中必须有return语句，以表示返回。+ subroutine s( ) 并且 call s( ) ，即必须有括号。“不得用” 的语言特性：- allocatable、intent、pointer、save、dimension 、parameter 和 target语句形式。 ( 用属性形式代替。) - external语句形式。 ( 用显式的接口代替。) - assign 、赋值 go to 、交错 return 、continue 、entry 、和计算 go to 语句。- include文件。 (用模块代
12、替。) - data 和 block data 。 ( 在类型声明语句中进行初始化或赋值。) - common块。 ( 将全局数据放在模块中，用模块代替。) - equivalence。 (被认为是许多难以查找的编程错误的来源。) - double precision语句。 ( 用 real 语句声明双精度的实型数。) - 语句函数。 ( 用内部函数代替。) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 5 页，共 9 页 - - - - - - - - - - - 专用固有函数。( 用类属函数代替。)
13、 - 假定大小数组。( 用假定形状数组代替。) - do n （其中 n 为语句标号）。( 用 do 和 end do 代替。 ) - 非整数 do 变量和表达式。- 同一行上多条语句。- 逻辑型 case 表达式。- 从 if块外面分支到end if 。- where 语句形式。 ( 用 where 结构形式。 ) - 在 open 和 inquire语句中的blank= 说明符。- 双字关键词之间要求有空格：in out ，go to 。不能写为inout ， goto 。9. 将字符串改为大写的子程序subroutine UpCase (str) != ! change to upper
14、 case != character(len=*),intent(in out) : str integer(4) : icha,LL,icval integer(4),parameter : diff = ichar(a) - ichar(A) LL = len_trim(str) do icha=1,LL icval = ichar(str(icha:icha) if (icval=ichar(a) .and. icval=ichar(z) then str(icha:icha) = char(icval-diff) end if end do return end subroutine
15、UpCase 10. CVF中源代码信息浏览默认情况下，代码信息浏览是失效的（Go to Definition/Reference都不起作用），你可以将其启用：精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 6 页，共 9 页 - - - - - - - - - - 1. Project/Settings/Fortran, 选中 Generate Source Browse Information。2. 选 BrowseInfo卡页，选中Build Browse info file，点击 OK 。3. B
16、uild程序。代码浏览器仅当成功Build程序后才可应用。4. Tools/Source Browser 。5. 在 Browse 对话框底部，找到Case sensitive复选框。 Fortran是大小写不敏感的语言，去除Case sensitive的选中。点击OK 。6. 现在可以用Browse 对话框来浏览函数调用关系以及变量声明和引用关系。7. 光标置于一个变量名或函数名中，鼠标右键：Go to Definition使光标跳到该变量声明或函数定义语句；Go to Reference 使光标跳到该变量的引用或函数引用语句。注： (Addison 提供 ) IVF 中现在不适用，因为此技
17、术是属于MS的，不是INTEL 的，所以， intel正在与 MS就有关问题进行磋商，目前的进展比较顺利，可能要在以后的版本中支持，目前还没有时间表。11. 如何在 CVF中检验代码是否符合F90/95 标准？为了检验你的代码是否符合F90 或 F95 标准，在CVF中选：Project/Settings. 选 Fortran页卡Category 下拉框选： Compilation Diagnostics Fortran Standard Checking 下拉框选： Fortran 90 或Fortran 95 这样，凡是CVF自家扩展的语法就都会有提示了。例：type test integ
18、er, allocatable : a(:) integer, pointer : b(:) end type test 若选了标准检验，会提示：Warning: Allocatable fields of derived types are non-standard A 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 7 页，共 9 页 - - - - - - - - - - 也就是说按照F90/95 标准， type 中不允许用allocatable。臭石头雪球提供：VS2005 + 组合的方法差不多
19、：项目属性， Fortran ，Diagnostics，Warn For Non-standard Fortran 选择为 Fortran 90 或Fortran 95 12. 如何输出tab 字符？char(9) 。13. 如何在 CVF中设置把Tab 转为空格Tab 字符不在Fortran标准中，在不同环境下会有不同的显示，不推荐使用。采用下列方法，可将Tab 键转为 3 个空格：依次选菜单Tools/Option，选 Tab 卡项，File type 选 Fortran ，选 Insert spaces ，Tab size 选 3，Indent size 选 3。14. 如何实现“代码自
20、动补全”在 CVF中：在 IVF 中：问题各位大侠，非常感谢你们的回复。置顶帖中关于批处理读取文件的帖子我已经看了。帖子里是批处理文件，文件名是1-4000 的数字序列。而我要读取的是时间序列的文件，即1 月 1 日-12 月 31 日， year+month+day 的形式。因此置顶帖中的例子不适用我这个问题。文件名为（2006 年 365 天的数据，共365 个文件）pre.dat pre.dat pre.dat 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 8 页，共 9 页 - - - - -
21、- - - - - . . pre.dat pre.dat pre.dat . . pre.dat 单个文件读取已经完成，但是如何批处理还没有解决。解决 1：既然是有规律你就可以做循环来进行读取, 每次循环读取一个文件. 例如做 3 层循环 , 分别是年循环( 如果就一年那就不用循环了), 月循环 , 日循环 . 不难吧 .解决 2我的方法：先把文件名：pre.dat,pre.dat,.输入到一个文件中，然后依次读取，赋给字符型变量数组filenames(i),然后就可以循环打开所有文件。程序如下：program test integer:a(12)=31,28,31,30,31,30,31,
22、31,30,31,30,31 character(len=12),allocatable:filenames(:) m=sum(a(1:12) allocate(filenames(m) open(10,file=) do i=1,12 do j=1,a(i) write(10,(a8,a4) ,i,j,.dat enddo enddo rewind 10 do i=1,m read(10,*) filenames(i) c open(i,file=filenames(i) ！我这没有那些文件，所以注释掉了。enddo close(10) end program test精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 9 页，共 9 页 - - - - - - - - - -