您所在位置：主题：[原创]Fortran讨论区 FAQ (不定期更新)
mltx
[专家分：20880] 发布于 2006-07-30 10:22:00
为了减少重复回答问题，特编此帖，并不定期添加和更新内容。
错误难免，欢迎讨论，仅供参考。
很多人问哪里可以找到Fortran编译器，有不少热心学友提供网址，特汇集在这里。虽然俺检验过这些链接，但是它们不一定总有效。
[color=FF0000]Fortran 编译器下载：[/color]
CVF 6.6
http://wwater.a365.net/SoftDown.asp?ID=11937
FTN95 (License: Free for personal use)
http://www.download.com/Silverfrost-FTN95/3000-2069_4-10491440.html?tag=lst-0-10
以下操作，如无特别说明，都是以CVF6.6为例。
[color=FF0000]1. 如何加大Stack size?[/color]
选Project => Settings => Link => Category: Output =>
Stack allocations
Reserve: 这里填新值（默认为1M，若需要10M，则填10000000）
[color=FF0000]2. 如何用Fortran批量生成文件？[/color]
设要生成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
[color=FF0000]3. 如何用Fortran动态生成输出格式？[/color]
设有一个数组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
[color=FF0000]4. MS POWERSTATION4.0是不是很垃圾？[/color]
POWERSTATION4.0是垃圾，其中Bug太多，多到不可用的地步！
在这个主题里，换了CVF后问题就没了的人已有相当的数目。
如果你用POWERSTATION4.0，遇到莫名其妙的错误，建议换CVF6.6，这是一个比较成熟的编译器。
[color=FF0000]5. 如何用F90/95生成随机数？[/color]
注意：
现在计算机产生的随机数都是伪随机数。
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
[color=FF0000]6. 函数/子程序超载的例子[/color]
设要编一个两个变量值互换的子程序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
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
[color=FF0000]7. 如何在CVF中为多行代码加注释?[/color]
见：http://www.programfan.com/club/showbbs.asp?id=124318
[color=FF0000]8. 推荐好的代码风格[/color]
根据F90子集语言ELF90和F的要求整理（部分）。
[color=0000FF]“强迫用”[/color]的语言特性：
+ 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( )，即必须有括号。
[color=0000FF]“不得用”[/color]的语言特性：
- 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。
[color=FF0000]9. 将字符串改为大写的子程序[/color]
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
[color=FF0000]10. CVF中源代码信息浏览[/color]
默认情况下，代码信息浏览是失效的（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就有关问题进行磋商，目前的进展比较顺利，可能要在9.1以后的版本中支持，目前还没有时间表。
[color=FF0000]11. 如何在CVF中检验代码是否符合F90/95标准？[/color]
为了检验你的代码是否符合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。
臭石头雪球提供：
[color=0000FF]VS2005 + IVF9.1 组合的方法差不多：[/color]
项目属性，Fortran，Diagnostics，Warn For Non-standard Fortran
选择为Fortran 90 或 Fortran 95
[color=FF0000]12. 如何输出tab字符？[/color]
char(9)。
[color=FF0000]13. 如何在CVF中设置把Tab转为空格[/color]
Tab字符不在Fortran标准中，在不同环境下会有不同的显示，不推荐使用。
采用下列方法，可将Tab键转为3个空格：
依次选菜单Tools/Option，
选Tab卡项，
File type选Fortran，
选Insert spaces，
Tab size选3，
Indent size选3。
[color=FF0000]14. 如何实现“代码自动补全”[/color]
[color=0000FF]在CVF中：[/color]
http://www.programfan.com/club/post-238710.html
[color=0000FF]在IVF中：[/color]
http://www.programfan.com/club/post-241489.html
错误难免，欢迎讨论，仅供参考。
很多人问哪里可以找到Fortran编译器，有不少热心学友提供网址，特汇集在这里。虽然俺检验过这些链接，但是它们不一定总有效。
[color=FF0000]Fortran 编译器下载：[/color]
CVF 6.6
http://wwater.a365.net/SoftDown.asp?ID=11937
FTN95 (License: Free for personal use)
http://www.download.com/Silverfrost-FTN95/3000-2069_4-10491440.html?tag=lst-0-10
以下操作，如无特别说明，都是以CVF6.6为例。
[color=FF0000]1. 如何加大Stack size?[/color]
选Project => Settings => Link => Category: Output =>
Stack allocations
Reserve: 这里填新值（默认为1M，若需要10M，则填10000000）
[color=FF0000]2. 如何用Fortran批量生成文件？[/color]
设要生成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
[color=FF0000]3. 如何用Fortran动态生成输出格式？[/color]
设有一个数组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
[color=FF0000]4. MS POWERSTATION4.0是不是很垃圾？[/color]
POWERSTATION4.0是垃圾，其中Bug太多，多到不可用的地步！
在这个主题里，换了CVF后问题就没了的人已有相当的数目。
如果你用POWERSTATION4.0，遇到莫名其妙的错误，建议换CVF6.6，这是一个比较成熟的编译器。
[color=FF0000]5. 如何用F90/95生成随机数？[/color]
注意：
现在计算机产生的随机数都是伪随机数。
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
[color=FF0000]6. 函数/子程序超载的例子[/color]
设要编一个两个变量值互换的子程序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
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
[color=FF0000]7. 如何在CVF中为多行代码加注释?[/color]
见：http://www.programfan.com/club/showbbs.asp?id=124318
[color=FF0000]8. 推荐好的代码风格[/color]
根据F90子集语言ELF90和F的要求整理（部分）。
[color=0000FF]“强迫用”[/color]的语言特性：
+ 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( )，即必须有括号。
[color=0000FF]“不得用”[/color]的语言特性：
- 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。
[color=FF0000]9. 将字符串改为大写的子程序[/color]
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
[color=FF0000]10. CVF中源代码信息浏览[/color]
默认情况下，代码信息浏览是失效的（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就有关问题进行磋商，目前的进展比较顺利，可能要在9.1以后的版本中支持，目前还没有时间表。
[color=FF0000]11. 如何在CVF中检验代码是否符合F90/95标准？[/color]
为了检验你的代码是否符合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。
臭石头雪球提供：
[color=0000FF]VS2005 + IVF9.1 组合的方法差不多：[/color]
项目属性，Fortran，Diagnostics，Warn For Non-standard Fortran
选择为Fortran 90 或 Fortran 95
[color=FF0000]12. 如何输出tab字符？[/color]
char(9)。
[color=FF0000]13. 如何在CVF中设置把Tab转为空格[/color]
Tab字符不在Fortran标准中，在不同环境下会有不同的显示，不推荐使用。
采用下列方法，可将Tab键转为3个空格：
依次选菜单Tools/Option，
选Tab卡项，
File type选Fortran，
选Insert spaces，
Tab size选3，
Indent size选3。
[color=FF0000]14. 如何实现“代码自动补全”[/color]
[color=0000FF]在CVF中：[/color]
http://www.programfan.com/club/post-238710.html
[color=0000FF]在IVF中：[/color]
http://www.programfan.com/club/post-241489.html
