2022年Fortran常见问题汇总 .pdf

《2022年Fortran常见问题汇总 .pdf》由会员分享，可在线阅读，更多相关《2022年Fortran常见问题汇总 .pdf（10页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、Fortran常见问题汇总。 （转自编程爱好者Fortran论坛的 FAQ）以下操作，如无特别说明，都是以CVF6.6 为例。1. 如何加大 Stack size? 选 Project = Settings = Link = Category: Output = Stack allocations Reserve: 这里填新值（默认为1M，若需要 10M，则填 10000000）2. 如何用 Fortran 批量生成文件？设要生成 4000 个文件，文件名为AA1-AA4000 ，如何写循环生成文件 ,而不用写 4000 次 write 命令呢 ? 用内部文件：character(len=80

2、) : 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 case (100:999) write(form,(i3) i 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 10 页 - - - - - - - - - case (1000:9999) write(form,(i4) i end s

3、elect 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) : form real : data(100) integer : i,num data = (/ (i,i=1,100) /)/10.0 read(*,*) num write(form,*

4、) (,num,f10.3) write(*,form) data stop 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 10 页 - - - - - - - - - end 4. MS POWERSTATION4.0 是不是很垃圾？POWERSTATION4.0 是垃圾，其中Bug 太多，多到不可用的地步！在这个主题里，换了CVF后问题就没了的人已有相当的数目。如果你用 POWERSTATION4.0，遇到莫名其妙的错误，建议换CVF6.6，这是一个比较成熟的编译器

5、。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 名师资料总结 - - -精品资料欢迎下载 - -

6、 - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 10 页 - - - - - - - - - end program random 6. 函数/子程序超载的例子设要编一个两个变量值互换的子程序swap(a,b) ，哑元 a,b 可能是实型数，整型数，数组，矩阵，字符串，派生类型等等。但是希望只用一个子程序接口swap(a,b) 来实现。 F90 可以用类属接口来实现这种子程序超载：module Utilities implicit none private I_Swap,R_Swap,RVec_Swap,RMat_Sw

7、ap,Type_Swap public : Swap interface Swap module procedure I_Swap,R_Swap,RVec_Swap,RMat_Swap,Type_Swap end interface contains subroutine i_swap (a,b) ! 整型置换 integer (ikind),intent(in out) : a,b integer (ikind) : t 。 ! 略end subroutine i_swap subroutine r_swap (a,b) ! 实型置换名师资料总结 - - -精品资料欢迎下载 - - - -

8、- - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 10 页 - - - - - - - - - real (rkind), intent(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_S

9、wap (a(i),b(i) end do return end subroutine RVec_swap subroutine RMat_swap (a,b) ! 实型矩阵置换。 ! 略end subroutine RMat_swap subroutine Type_swap (a,b) ! 派生类型置换。 ! 略end subroutine Type_swap end module Utilities 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 10 页 - -

10、- - - - - - - 7. 如何在 CVF中为多行代码加注释? 见： http:/ 8. 推荐好的代码风格根据 F90 子集语言 ELF90 和 F 的要求整理（部分）。“ 强迫用 ” 的语言特性：+ F90 的自由格式的源代码。+ implicit none。+ 子过程的哑元都要有intent 属性。+ 函数子程序的哑元必须指定为intent(in) 。+ 所有子程序和函数都放在模块（module ）中，然后引用（ use）该模块；或者放在program 中。+ 数组哑元要求是假定形状的，或者有固定的维数和大小。字符哑元要求是假定长度的。+ 对于 recursive function

11、（递归函数）语句，必须有result 子句。+ 在所有派生类型 （type ）的定义语句中， 必须用双冒号分隔符（: ）。+ 主程序要求有program 语句。+ 在程序单元的end 语句中要求后跟程序单元的类型和名称。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 10 页 - - - - - - - - - + 在 end type 语句中要求后跟类型的名称。+ end program 前必须有 stop 语句以表示停止执行。+ 子过程中必须有return 语句，以表

12、示返回。+ subroutine s( )并且 call s( )，即必须有括号。“ 不得用 ” 的语言特性：- allocatable 、intent 、pointer 、save、dimension 、parameter和 target 语句形式。 (用属性形式代替。 ) - external 语句形式。 (用显式的接口代替。) - assign、赋值 go to、交错 return 、continue 、entry 、和计算 go to 语句。- include 文件。 (用模块代替。 ) - data 和 block data。(在类型声明语句中进行初始化或赋值。) - common

13、块。 (将全局数据放在模块中，用模块代替。) - equivalence 。(被认为是许多难以查找的编程错误的来源。) - double precision 语句。 (用 real 语句声明双精度的实型数。) - 语句函数。 (用内部函数代替。) - 专用固有函数。 (用类属函数代替。 ) - 假定大小数组。 (用假定形状数组代替。) - do n （其中 n 为语句标号）。 (用 do 和 end do 代替。 ) - 非整数 do 变量和表达式。- 同一行上多条语句。- 逻辑型 case 表达式。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - -

14、- - - - - 名师精心整理 - - - - - - - 第 7 页，共 10 页 - - - - - - - - - - 从 if 块外面分支到end if。- where 语句形式。 (用 where 结构形式。 ) - 在 open 和 inquire 语句中的 blank= 说明符。- 双字关键词之间要求有空格：in out，go to。不能写为 inout ，goto 。9. 将字符串改为大写的子程序subroutine UpCase (str) != ! change to upper case != character(len=*),intent(in out) : str

15、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 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - -

16、 - - - - 第 8 页，共 10 页 - - - - - - - - - return end subroutine UpCase 10. CVF中源代码信息浏览默认情况下，代码信息浏览是失效的（Go to Definition/Reference 都不起作用），你可以将其启用：1. Project/Settings/Fortran, 选中 Generate Source Browse Information 。2. 选 BrowseInfo 卡页，选中 Build Browse info file，点击 OK。3. Build 程序。代码浏览器仅当成功Build 程序后才可应用。4.

17、Tools/Source Browser 。5. 在 Browse 对话框底部，找到Case sensitive 复选框。 Fortran 是大小写不敏感的语言，去除Case sensitive 的选中。点击OK。6. 现在可以用 Browse 对话框来浏览函数调用关系以及变量声明和引用关系。7. 光标置于一个变量名或函数名中，鼠标右键：Go to Definition 使光标跳到该变量声明或函数定义语句；Go to Reference使光标跳到该变量的引用或函数引用语句。注： (Addison 提供 ) IVF 中现在不适用，因为此技术是属于MS 的，不是 INTEL 的，所以，intel 正在与 MS就有关问题进行磋商，目前的进展比较顺利，可能要在9.1 以后的版本中支持，目前还没有时间表。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 10 页 - - - - - - - - - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 10 页 - - - - - - - - -