很多人问哪里可以找到Fortran编译器,有不少热心学友提供网址,特汇集在这里。虽然俺检验过这些链接,但是它们不一定总有效。 Fortran 编译器下载: CVF
FTN95 (License: Free for personal use)
以下操作,如无特别说明,都是以为例。
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\ 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) /)/ read(*,*) num
write(form,*) \"(\write(*,form) data stop end
4. MS 是不是很垃圾?
是垃圾,其中Bug太多,多到不可用的地步!
在这个主题里,换了CVF后问题就没了的人已有相当的数目。
如果你用,遇到莫名其妙的错误,建议换,这是一个比较成熟的编译器。
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 = 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 subroutine RMat_swap
subroutine Type_swap (a,b) ! 派生类型置换 。。。 ! 略 end subroutine Type_swap
end module Utilities
t 7. 如何在CVF中为多行代码加注释? 见:
8. 推荐好的代码风格
根据F90子集语言ELF90和F的要求整理(部分)。 “强迫用”的语言特性: + F90的自由格式的源代码。 + implicit none。
+ 子过程的哑元都要有intent属性。 + 函数子程序的哑元必须指定为intent(in)。
+ 所有子程序和函数都放在模块(module)中,然后引用(use)该模块;或者放在program中。
+ 数组哑元要求是假定形状的,或者有固定的维数和大小。字符哑元要求是假定长度的。 + 对于recursive function(递归函数)语句,必须有result子句。 + 在所有派生类型(type)的定义语句中,必须用双冒号分隔符(::)。 + 主程序要求有program语句。
+ 在程序单元的end语句中要求后跟程序单元的类型和名称。 + 在end type语句中要求后跟类型的名称。 + 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文件。(用模块代替。)
- data和block data。(在类型声明语句中进行初始化或赋值。) - common块。(将全局数据放在模块中,用模块代替。) - equivalence。(被认为是许多难以查找的编程错误的来源。) - double precision语句。(用real语句声明双精度的实型数。) - 语句函数。(用内部函数代替。)
- 专用固有函数。(用类属函数代替。) - 假定大小数组。(用假定形状数组代替。)
- 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 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 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. 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就有关问题进行磋商,目前的进展比较顺利,可能要在以后的版本中支持,目前还没有时间表。
11. 如何在CVF中检验代码是否符合F90/95标准? 为了检验你的代码是否符合F90或F95标准,在CVF中选:
Project/Settings... 选Fortran页卡
Category下拉框选:Compilation Diagnostics
Fortran Standard Checking下拉框选:Fortran 90 或 Fortran 95
这样,凡是CVF自家扩展的语法就都会有提示了。 例:
type test
integer, allocatable :: a(:) integer, pointer :: b(:) end type test
若选了标准检验,会提示:
Warning: Allocatable fields of derived types are non-standard [A]
也就是说按照F90/95标准,type中不允许用allocatable。
臭石头雪球提供:
VS2005 + 组合的方法差不多:
项目属性,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. 如何实现“代码自动补全” 在CVF中:
在IVF中: 问题
各位大侠,非常感谢你们的回复。置顶帖中关于批处理读取文件的帖子我已经看了。 帖子里是批处理文件,文件名是1-4000的数字序列。
而我要读取的是时间序列的文件,即1月1日-12月31日,year+month+day 的形式。因此置顶帖中的例子不适用我这个问题。
文件名为 (2006年365天的数据,共365个文件) pre..dat pre..dat pre..dat
... ... 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,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
因篇幅问题不能全部显示,请点此查看更多更全内容