Fortran的主要版本与差别

1、、Fortran的主要版本及差别按其发展历史，Fortran编译器的版本其实很多。现在在广泛使用的是Fortran77和Fortran90oFortran90在Fortran77基础上添加了不少使用的功能，并且改良了77编程的版面格式，所以编程时推荐使用90。鉴于很多现成的程序只有77版本，有必要知道77的一些基本常识，至少保证能够看懂77程序。以下是77和90的一些格式上的区别。Fortran77:固定格式(fixedformat)，程序代码扩展名：.f或.for(1) 若某行以C,c或*开头，则该行被当成注释；每行前六个字符不能写程序代码，可空着，或者15字符以数字表明行代码(用作格式化输

2、入输出等)；772为程序代码编写区；73往后被忽略；太长的话可以续行，所续行的第六个字符必须是"0"以外的任何字符Fortran90:自由格式(freeformat),扩展名：.f90(1) 以"!"引导注释；(2) 每行可132字符，行代码放在每行最前面；(3) 以&J行，放在该行末或下行初。以下都是讨论Fortran90。3、Fortran的一些特点，和C的一些不同其实很多，在下面涉及具体方面时可以看到。这里只是大致提一些。(1) 不分大小写(2) 每句末尾不必要写分号(3) 程序代码命令问的空格没有意义(4) 不像C,Fortran不使用数

3、据类型多出了复数和逻辑判断类型。比如复数类型complex:a!声明复数的方法。复数显然方便了科学计算，满足了工程方面的需求a=(1.0,2.0)!a=1+i多出了乘籍运算(*)。乘籍除了整数还可以是实数形式。如开方，开立方a=4.0*0.5,a=8.0*(1.0/3.0)。(5) 数组有一些整体操作的功能；可以方便的对部分元素进行操作有些情况下可以声明大小待定的数组，很实用的功能4、Fortran的基本程序结构先看一看所谓的"HelloFortran"程序。programmain!程序开始,main是program的名字，完全自定义write(*,*)"Hell

4、o”!主程序stop!终止程序endprogrammain!end用于封装代码，表示代码编写完毕。中的容可省略，下同。再看一段实用一些的程序，好有点感性认识。程序用于计算圆柱的表面积，要求输入底面半径和高。其中展示了Fortran的一些特色用法。程序摘自维基。其实是一个叫www.answers.的网上引用的维基的网页。推荐去看看！能查到不少有意思的东西。programcylinder!给主函数起个名字!Calculatetheareaofacylinder.!Declarevariablesandconstants.!constants=pi!variables=radiussquaredan

5、dheightimplicitnone!Requireallvariablestobeexplicitlydeclared这个一般都是要写上的。下面会进一步说明integer:ierrcharacter:ynreal:radius,height,areareal,parameter:pi=3.1415926536!这是常量的声明方法interactive_loop:do!do循环，Fortran中的循环可以加标签，如d面的interactive_loop就是标Prompttheuserforradiusandheightandreadthem.write(*,*)'Enterradiu

6、sandheight.read(*,*,iostat=ierr)radius,height屏幕输出键盘输入。isotat的值用于判断!输入成功否。!Ifradiusandheightcouldnotbereadfrominput,!thencyclethroughtheloop.if(ierr/=0)thenwrite(*,*)'Error,invalidinput.'cycleinteractive_loop!cycle相当于C里的continueendif!Computearea.The*means"raisetoapower."area=2*pi*(r

7、adius*2+radius*height)!便!Writetheinputvariables(radius,height)!andoutput(area)tothescreen.write(*,'(1x,a7,f6.2,5x,a7,f6.2,5x,a5,f6.2)')&!"&"这里指数运算比C方表小续行O还显示了格式化输出radius=',radius,'height=',height,'area=',areayn=''yn_loop:do!嵌的另一个do循环write(*,*)

8、9;Performanothercalculation?yn'read(*,'(a1)')ynif(yn='y'.or.yn='Y')exityn_loopif(yn='n'.or.yn='N'.or.yn='')exitinteractive_loopenddoyn_loop!结束嵌do循环enddointeractive_loop股还会有些module的部分在主函数前,endprogramcylinderFortran程序的主要结构就是这样了函数在主函数后。三、数据类型及基本输入输出1、

9、数据类型，声明及赋初值(1) integer:短整型kind=2,长整型kind=4integer(kind=2):a=3如果声明成integer:a,则默认为长整型。!":"在声明并同时赋初值时必须要写上；类型名后面有形容词时也必须有；其他情况可略去!所谓形容词，可以看一下这个。比如声明常数real,parameter:pi=3.1415926。parameter就是形容词。(2) real:单精度kind=4(默认)，双精度kind=8real(kind=8):a=3.0还有指数的形式，如1E10为单精度，1D10为双精度(3) complex单精度和双精度comple

10、x(kind=4)b(4) charactercharacter(len=10)c!len为最大长度(5) logicallogical*2:d=.ture.(等价于logical(2):d=.ture.)自定义类型type:类似于C中的structFortran77中给变量赋初值常用DATAOr令，可同时给多个变量赋初值dataa,b,string/1,2.0,'fortran'/与C不同的是，Fortran中变量不声明也能使用，即有默认类型(跟implicit命令有关)。按照默认的规定，以i,j,k,l,m,n开头的变量被定义为integer,其余为real。取消该设置需在

11、程序声明部分之前加implicitnone。国伦建议一般都使用该语句。另一点关于声明的不同是Fortran有”等价声明”：integera,bequivalence(a,b)使得a,b使用同一块存。这样可以节省存；有时可精简代码。如：equivalence(很长名字的变量比如三维数组的某个元素，a),之后使用a来编写程序就简洁多了2、基本输入输出输入：read(*,*)a!从键盘读入输出：write(*,*)"text"!在屏幕上输出。Fortran77用'text'。Fortran90!一般""和''都可print*,&

12、quot;text”!只能屏幕输出(*,*)完整写为(unit=*,fmt=*)。其中unit为输入/输出位置，如屏幕，文件等；fmt为格式。如果这两项都写成*，则按默认的方式进行，即上面描述的。print后面的*表示按默认格式输四、流程控制1、运算符（1）逻辑运算符=/=>>=<<=!Fortran90用法.EQ.NE.GT.GE.LT.LE.!Fortran77用法（2）涉及相互关系的集合运算符.AND.OR.NOT.EQV.NEQV.!仅.NOT.连接一个表达式，其余左右两边都要有表达式（可以是logical类型的变量）!.EQV.:当两边逻辑运算值相同时为真，.

13、NEQV.:当两边逻辑运算值不同时为真2、IF（1）基本：if（逻辑判断式）thenendif如果then后面只有一句，可写为if（逻辑判断式）！then和endif可省略多重判断：if（逻辑判断式）thenelseifelseifelseendif嵌套：if（逻辑判断式）thenif（逻辑判断式）thenif（逻辑判断式）thenelseif（逻辑判断式）thenelseendif!没必要跟C语句多时用(了，因为有endifendifendif算术判断：programexampleimplicitnonerealcwrite(*,*)"inputanumber"read(

14、*,*)cif(c)10,20,30!10,20和30为行代码，根据c小于/等于/大于0,执行10/20/30行的程序10write(*,*)"A"goto40!goto可实现跳到任意前面或后面的行代码处，但用多了破坏程序结构20write(*,*)"B"goto4030write(*,*)"C"goto4040stopend3、SELECTCASE类似于C的switch语句selectcase（变量）case（数值1）!比如case（1:5）代表1<=变量<=5会执行该模块!case（1,3,5）代表变量等于1或3或5会

15、执行该模块case（数值2）!括号中数值只能是integer,character或logical型常量，不台匕北是real型.casedefaultendcase4、PAUSE,CONTINUEpause暂停程序执行，按enter可继续执行continue貌似没什么用处，可用作封装程序的标志五、循环1、DOdocounter=初值，终值，增/减量！counter的值从初值到终值按增/减量变化，!counter每取一个值对应着一次循环。增/减量不写则默认为1.!循环主体也没有必要用enddoFortran77中不是用enddo来终止，而是下面这样子：do循环最后一行的行代码counter=初值，

16、终值，增/减量行代码!这是do的最后一行2、DOWHILEdowhile（逻辑运算）enddo类似于C中的while（逻辑运算）。一开始那个计算圆柱表面积的程序中，应该也算是这一类。不过它是通过部的if语句来控制循环。看来也是可以的，不过在这本书上没看到这样写。其实应该也可以归于下面这种。3、没看到和C里面的dowhile（逻辑运算）；相对应的循环语句，不过可以这样，保证至少做一次循环：dowhile（.ture.）if（逻辑运算）exit!exit就好比C里面的break。C里的continue在Fortran里是cycleenddo4、Fortran的一个特色：带署名的循环可以这样，不易出

17、错：outer:doi=1,3inner:doj=1,3enddoinnerenddoouter还可以这样，很方便：loop1:doi=1,3loop2:doj=1,3终止整个循环loop1跳出loop2的本次循环，进行loop2的if(i=3)exitloop1!exitif(j=2)cycleloop2!cycle一次循环write(*,*)i,jenddoloop2enddoloopl还有一些循环主要用于Fortran中的数组运算，为Fortran特有，很实用六、数组1、数组的声明和C不同的是,Fortran中的数组元素的索引值写在(),且高维的也只用一个(),如下integera(5)

18、!声明一个整型一维数组real:b(3,6)!声明一个实型二维数组类型可以是integer,real,character,logical或type。最高可以至U7维。数组大小必须为常数。但是和C语言不同,Fortran也有办法使用大小可变的数组，方法如下：integer,allocatable:a(:)!声明大小可变经过某个途径得知所需数组大小size之后，用下面的语句：allocate(a(size)!配置存空间之后该数组和通过一般方法声明的数组完全相同。与C不同，Fortran索引值默认为从1开始，而且可以在声明时改变该规则:integera(-3:1)!索引值为-3,-2,-1,0,1i

19、ntegerb(2:3,-1:3)!b(23,-13)为可使用的元素2、数组在存中的存放和C不同，Fortran中的数组比如a(2,2)在存中存放顺序为a(1,1),a(2,1),a(1,2),a(2,2)0原则是先放低维的元素，再放高维的元素。此规则称为columnmajor。3、赋初值(1) 最普通的做法：integera(5)dataa/1,2,3,4,5/或integer:a(5)=(/1,2,3,4,5/)若integer:a(5)=5，贝U5个元素均为5对于integer:a(2,2)=(/1,2,3,4/)根据数组元素在存中存放的方式，等价于赋值a(1,1)=1,a(2,1)=2

20、,a(1,2)=3,a(2,2)=4(2) 利用Fortran的特色：隐含式循环integera(5)integeridata(a(i),i=2,4)/2,3,4/!(a(i),i=2,4)认值1看例子就明白了。表示i从2到4循环，增量为默还可以这样：五个元素分别赋值为1,2,2,2,5五个元素分别赋值为1,2,3,4,integeriinteger:a(5)=(/1,(2,i=2,4),5/)!integer:b(5)=(/i,i=1,5/)!5还可以嵌套data(a(i,j),i=1,2),j=1,2)=/1,2,3,4/!a(1,1)=1,1(2,1)=2,a(1,2)=3,a(2,2)

21、=44、操作整个数组设a,b为相同类型、维数和大小的数组a=5!所有元素赋值为5a=(/1,2,3/)!这里假设a为一维，a(1)=1,a(2)=2,a(3)=3a=b!对应元素赋值，要求a,b,c维数和大小相同，下同a=b+ca=b-ca=b*ca=b/ca=sin(b)!部函数都可以这样用5、操作部分数组元素a为一维数组a(3:5)=(/3,4,5/)=3,a(4)=4,a(5)=5a(1:5:2)=3!a(1)=3,a(3)=3,a(5)=3a(3:)=5!a(3)以及之后的所有元素赋值为5a(1:3)=b(4:6)!类似于这种的要求左右数组元素个数相同a(:)=b(:,2)!a(1)=

22、b(1,2),a(2)=b(2,2),以此类推6、WHEREwhere形式上类似于if,但只用于设置数组。设有两个同样类型、维数和大小的数组a,b34mwhere(a<3)b=a!a中小于3的元素赋值给b对应位置的元素endwhere再如：where(a(1:3)/=0)c=a!略去了endwhere,因为只跟了一行where可嵌套，也!可类似do循环有署名标签7、FORALL有点像C中的for循环：forall(triplet1,triplet2,triplet3,mask)其中triplet形如i=2:6:2,表示循环，最后一个数字省略则增量为1例如：forall(i=1:5,j=1

23、:5,a(i,j)<10)a(i,j)=1endforall乂如：forall(i=1:5,j=1:5,a(i,j)/=0)a(i,j)=1/a(i,j)forall也可以嵌套使用，好比C中for循环的嵌套。七、函数Fortran中函数分两类：子程序(subroutine)和自定义函数(function)。自定义函数本质上就是数学上的函数，一般要传递自变量给自定义函数，返回函数值。子程序不一定是这样，可以没有返回值。传递参数要注意类型的对应，这跟C是一样的1、子程序目的：把某一段经常使用的有特定功能的程序独立出来，可以方便调用。习惯上一般都把子程序放在主程序结束之后。形式：subrout

24、inename(parameterl,parameter2)!给子程序起一个有意义的名字。可以传递参数，这样可以有返回值。括号也可以空着，代表不传递参数。implicitnoneinteger:parameter1,parameter2!需要定义一下接收参数的类型。!接下来的程序编写跟主程序没有任何差别。mreturn!跟C不同，这里表示子程序执行后回到调用它的地方继续执行下面的程序。不一定!放在最后。可以放在子程序的其他位置，作用相同；子程序中return之后的部分不再执行。endsubroutinename调用：使用call命令直接使用，不需要声明。在调用处写:callsubroutine

25、name(parameter1,parameter2)汪息点：子程序之间也可相互调用。直接调用就是了，像在主程序中调用子程序一样传递参数的原理和C中不同。Fortran里是传址调用(callbyaddress/reference),就是传递时用的参数和子程序中接收时用的参数使用同一个地址，尽管命名可以不同。这样如果子程序的执行改变子程序中接收参数的值，所传递的参数也相应发生变化。子程序各自部定义的变量具有独立性，类似于G各自的行代码也具有独立性。因此各个子程序和主程序中有相同的变量名、行代码号，并不会相互影响。2、自定义函数和子程序的明显不同在于：需要在主程序中声明之后才能使用。调用方式也有差

26、别。另外按照惯例调用函数不去改变自变量的值。如果要改变传递参数的值，习惯上用子程序来做声明方式：real,external:function_name一般自定义函数也是放在主程序之后。形式：functionfunction_name(parameter1,parameter2)implicitnonereal:parameter1,parameter2!声明函数参数类型，这是必需的real:function_name!声明函数返回值类型，这是必需的function_name=-.!返回值的表达式returnend也可以这样直接声明返回值类型，简洁些：realfunctionfunction_n

27、ame(parameter1,parameter2)implicitnonereal:parameterl,parameter2!这个还是必需的function_name=-.!返回值表达式returnend调用：function_name(parameter1,parameter2)不需要call命令。自定义函数可以相互调用。调用时也需要事先声明。总之，调用自定义函数前需要做声明，调用子程序则不需要3、关于函数中的变量(1)注意类型的对应。Fortran中甚至可以传递数值常量，但只有跟函数定义的参数类型对应才会得到想要的结果。如callShowReal(1.0)就必须用1.0而不是1。(2)

28、传递数组参数，也跟C一样是传地址，不过不一定是数组首地址，而可以是数组某个指定元素的地址。比如有数组a(5),调用callfunction(a)则传递a(1)的地址，调用callfunction(a(3)贝U传递a(3)的地址。多维数组作为函数参数，跟C相反的是，最后一维的大小可以不写，其他维大小必须写。这决定于Fortran中数组元素columnmajor的存放方式。在函数中，如果数组是接收用的参数，则在声明时可以用变量赋值它的大小，甚至可以不指定大小。例如：subroutineArray(num,size)implicitnoneinteger:sizeintegernum(size)!可

29、以定义一个数组，其大小是通过传递过来的参数决定的这个很实用returnendsave命令：将函数中的变量值在调用之后保留下来，下次调用此函数时该变量的值就是上次保留的值。只要在定义时加上save就行：integer,save:a=1传递函数(包括自定义函数、库函数、子程序都是可以的)。类似于C中的函数指针。需要在主程序和调用函数的函数中都声明作为参数传递的函数。如real,external:function!自定义函数real,intrinsic:sin!库函数externalsub!子程序(3) 函数使用接口(interface):一段程序模块。以下情况必需:a. 函数返回值为数组指定参数位

30、置来传递参数时所调用的函数参数个数不固定输入指标参数时函数返回值为指针时。具体用法结合例子容易看懂。例子都很长。看书吧。4、全局变量功能就不用说了。原理：根据声明时的相对位置关系而取用，不同与C中根据变量名使用。如果在主程序中定义：integer:a,bcommona,b!就是这样定义全局变量的在子程序或自定义函数中定义：integer:c,dcommonc,d则a和c共用相同存,b和d共用相同存。全局变量太多时会很麻烦。可以把它们人为归类，只需在定义时在common!面加上区问名。如common/groupel/a,common/group2/b。这样使用时就不必把所有全局变量都列出来了，再

31、声明common/groupel/c就可以用a、c全局变量了。可以使用blockdata程序模块。在主程序和函数中不能直接使用前面提到的data命令给全局变量赋初值。可以给它们各自赋初值；如果要使用data命令必须要这样:blockdatanameimplicitnoneintegera,b,creald,ecommonabccommon/group1/d,edataa,b,c,d,e/1,2,3,4.0,5.0/endblockdataname5、ModuleModule不是函数。它用于封装程序模块，一般是把具有相关功能的函数及变量封装在一起。用法很简单，但能提供很多方便，使程序变得简洁，比

32、如使用全局变量不必每次都声明一长申，写在Module里调用就行了。Module一般写在主程序开始之前形式：modulemodule_nameendmodulemodule_name使用：在主程序或函数中使用时，需要在声明之前先写上一行：usemodule_name37mModule中有函数时必须在contains命令之后（即在某一行写上contains然后下面开始写函数，多所有函数都写在这个contains之后）。并且module中定义过的变量在module里的函数中可以直接使用，函数之间也可以直接相互调用，连module中的自定义函数在被调用时也不用先声明的。6>include放在需要的任何地方，插入另外的文件（必须在同一目录下）。如:include'funcion.f90'八、文件1、文本文件Fortran里有两种读取文件的方式，对应于两种文件顺序读取：用于文本文件直接读取：用于二进制文件这里只摘录关于文本文件的读取。一般模式如下。character(len=20):filenamein="in.txt",filenameout="out.txt"!文件名logicalaliveinteger:fileidin=10,fileidout=20!10,20是给文件编