FORTRAN 90 程序编程规范

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

FORTRAN 90 程序编程规范

Fortran 90 编程规范,使程序代码高度组织化,更加易读、易懂、易于维护,程序更加高效。使编出的程序更易懂、易于维护。

1 语言选择

数值预报创新系统软件开发应避免使用Fortran77 的某些过时特征以Fortran 90不一致的特征。选择Fortran 90 作为开发语言,并采用Fortran 90 的新功能,如动态内存的分配(dynamic memory allocation)、递归(recursion ), 模块(modules)、POINTER 、长变量名、自由格式等。

Fortran 77其中某些只是一些冗余的功能,这些功能已经过时,另外,还有一些在Fortran90 中被证明是不好的用法,建议不要使用。

2 Fortran 90 的新特性

2.1.1 建议使用的Fortran 90 新特性

建议使用Fortran 90 提供的模块(module ),并用Use ONLY 指定module 中哪些变量或派生类型定义可用于调用程序。

尽量使用数组下标三元组,这样可优化并减少所需的代码行数。为提高可读性,要在括号内表明数组的维数,例如:

1dArrayA(:) = 1dArrayB(:) + 1dArrayC(:)

2dArray(: , :) = scalar * Another2dArray(: , :) 当访问数组的子集时,例如在有限差分等式中,可以通过使用下标三元组实现。例如:2dArray(: , 2:len2) = scalar *( &

Another2dArray(:, 1:len2 -1) &

- Another2dArray(:, 2:len2) &

)

对程序单元(program units )命名,并使用End program ,End subroutine ,End interface ,End module 等结构再次指定“program unit ”的名称。

在逻辑表达式中使用>、 >=、 ==、 <、 <=、 /=,它们分别代

替.gt.、.ge.、.eq.、.lt.、.le.、.ne. 。新的表示方法更接近标准的数学符号

在变量定义中始终使用“::”;始终用“DIMENSION ”定义数组形状;始终用(len=)的语法格式声明字符变量的长度。

从效率的角度考虑,除递归定义的数据结构外,应尽量避免使用递归程序(它在使用内存和cpu 上会降低效率)。

建议定义新的操作符,而不要重载已有的操作符,这样可以清晰地表明功能,避免降低代码的可读性和可维护性。

2.1.2 关于Interface

使用Interface 接口块。在FORTRAN 语言中,主调程序与被调程序是分别编译的。由于Fortran 90 对过程的许多功能做了扩充,有些功能单靠简单的调用语句已无法反映,因而系统也就无法进行正确的编译。这时需要在主调程序中加入interface 接口块,通过它为主调程序与被调程序指明一个显式的接口。如果被调用程序中哑元含有假定形状(assumed-shape )数组,或含可选变元,或含键盘输入的参数,就需要interface 接口块说明。一般来讲,在Fortran 90 程序之间需要提供的interface 块,建议将interface 接口块明确写入调用程序,并复制被调用程序的参数列表。这种方法简单易用,但也相应增加了维护代码的工作量,因为只要被调用程序的参数列表发生变化就必须相应改变interface 接口块和调用(call )语句。

2.1.3 关于动态内存的使用

并行程序中存在着动态内存使用不充分的情况,合理使用动态内存可提高效率。在Fortran90 中常采用下列方法获取动态内存。

自动数组(automatic arrays):这种数组的大小在子程序中指定,其大小取决于运行时的变量,例如,可通过参数列表向子程序传递实元变量。使用自动数组优先于其它形式的动态内存分配方法。

可分配数组(allocatable arrays):声明为ALLOCATABLE 属性的数组变量可在运行时通过ALLOCATE 命令分配空间。与指针不同的是ALLOCATABLE 属性不能用于派生数据类型中。这种方法申请的空间要使用DEALLOCATE 命令释放。

在一个程序中,不要重复进行“ALLOCATE—DEALLOCATE-然后再ALLOCATE 一个更大的空间”这样的操作,这样会导致大量不可利用的碎片空间产生。应在程序中对要开辟的动态数组空间进行连续分配,或连续释放。

在使用ALLOCATE 和DEALLOCATE 命令之后,要利用ALLOCATE 和DEALLOCATE 命令参数中的返回值来判断操作是否成功;不建议使用指针数组。

2.2 Fortran 90 中不建议使用的过时的特性

不要使用固定书写格式,采用Fortran 90 的自由书写格式。

不要使用旧式的DO 循环和CONTINUE 语句。在以前的DO 循环中,关键字之后大多有一个标号,循环在该标号处结束,如:

DO 10 I = 1, 10

DO 20 j= 1, 20

A(i,j)=0.

20 CONTINUE˜

10 CONTINUE˜

建议采用统一的DO 结构,如:

DO i = 1, 10

DO j = 1, 20

A(i, j) = 0.

END DO

END DO

不能使用实循环变量˜

在DO 循环中若用实数作循环变量,由于截断误差的存在,会导致得不到预期的循环次数,如应将下面的DO 循环

DO r = 0.1 ,0.9 ,0.1

...

END DO

改写为:

DO i = 1 , 9

r = i * 0.1

...

END DO

不要使用DO WHILE 语句,可用不带循环变量的DO 语句配合exit 替代,如:不要使用选择性返回。例如:

CALL Foo( x, y, *100, *200, *300 )

在执行完该子程序后程序根据Foo 的返回语句后的表达式值选择从标号为

100 、200 、

300 处继续执行。可将上述子程改写为结构化CASE 语句:

CALL Foo(x, y, return_code)

SELECT CASE(return_code)

...

不要使用算法IF 语句:使用块IF 结构来代替;

不要使用计算GO TO 语句,如GO TO (10,20,30) i ;

不要使用标号赋值语句,如

ASSIGN 300 i

相关文档
最新文档