Fortran常见问题汇总(版).pdfVIP

Fortran 常见问题汇总 1. 如何加大Stack size? 选Project = Settings = Link = Category: Output = Stack allocations Reserve: 这里填新值（默认为 1M，若需要10M，则 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 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) :: form real :: data(100) integer :: i,num data = (/ (i,i=1,100) /)/10.0 read(*,*) num write(form,*) (,num,f10.3) write(*,form) data stop end 4. MS POWERSTATION4.0 是不是很垃圾？ POWERSTATION4.0 是垃圾，其中Bug 太多，多到不可用的地步！ 在这个主题里，换了CVF 后问题就没了的人已有相当的数目。 如果你用POWERSTATION4.0，遇到莫名其妙的错误，建议换CVF6.6，这是一个比较成 熟的编译器。 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 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_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) ! 实型置换 real (rkind), intent(in out) :: a,b real (rkind) :: t