1、续行,放在该行末或下行初。以下都是讨论Fortran 90。3、Fortran的一些特点,和C的一些不同其实很多,在下面涉及具体方面时可以看到。这里只是大致提一些。(1)不分大小写(2)每句末尾不必要写分号(3)程序代码命令间的空格没有意义(4)不像C,Fortran不使用 (5)数据类型多出了复数和逻辑判断类型。比如复数类型complex : a !声明复数的方法。复数显然方便了科学计算,满足了工程方面的需求a=(1.0,2.0) ! a=1+i(6)多出了乘幂运算(*)。乘幂除了整数还可以是实数形式。如开方,开立方a=4.0*0.5,a=8.0*(1.0/3.0)。(7)数组有一些整体操作
2、的功能;可以方便的对部分元素进行操作(8)有些情况下可以声明大小待定的数组,很实用的功能4、Fortran的基本程序结构先看一看所谓的Hello Fortran程序。program main !程序开始,main是program的名字,完全自定义write(*,*) Hello !主程序stop !终止程序end programmain !end用于封装代码,表示代码编写完毕。 中的内容可省略,下同。 再看一段实用一些的程序,好有点感性认识。程序用于计算圆柱的表面积,要求输入底面半径和高。其中展示了Fortran的一些特色用法。程序摘自维基。其实是一个叫的网上引用的维基的网页。推荐去看看!能查
3、到不少有意思的东西。program cylinder !给主函数起个名字 Calculate the area of a cylinder. Declare variables and constants. constants=pi variables=radius squared and height implicit none ! Require all variables to be explicitly declared这个一般都是要写上的。下面会进一步说明。 integer : ierr character : yn real : radius, height, area real,
4、 parameter : pi = 3.1415926536 !这是常量的声明方法 interactive_loop: do !do循环,Fortran中的循环可以加标签,如do前面的 interactive_loop就是标签 Prompt the user for radius and height and read them. write (*,*) Enter radius and height.屏幕输出 read (*,*,iostat=ierr) radius,height !键盘输入。isotat的值用于判断输入成功否。 If radius and height could not
5、 be read from input, then cycle through the loop. if (ierr /= 0) then write(*,*) Error, invalid input. cycle interactive_loop !cycle 相当于C里的continue end if Compute area. The * means raise to a power. area = 2 * pi * (radius*2 + radius*height) ! 指数运算比C方便 Write the input variables (radius, height) and
6、output (area) to the screen. write (*,(1x,a7,f6.2,5x,a7,f6.2,5x,a5,f6.2) &表示续行。这里还显示了格式化输出 radius=,radius,height=,height,area=,area yn = yn_loop:内嵌的另一个do循环Perform another calculation? yn read(*,(a1) yn if (yn=y .or. yn=Y) exit yn_loopnN) exit interactive_loop end do yn_loop !结束内嵌do循环 end do interact
7、ive_loop end program cylinder Fortran程序的主要结构就是这样了。一般还会有些module的部分在主函数前,函数在主函数后。三、数据类型及基本输入输出1、数据类型,声明及赋初值(1)integer: 短整型kind=2, 长整型kind=4 integer(kind=2) : a=3 如果声明成integer: a,则默认为长整型。 在声明并同时赋初值时必须要写上;类型名后面有形容词时也必须有;其他情况可略去所谓形容词,可以看一下这个。比如声明常数 real,parameter : pi=3.1415926 。parameter就是形容词。(2)real:单精
8、度kind=4(默认),双精度kind=8 real(kind=8) : a=3.0 还有指数的形式,如1E10为单精度,1D10为双精度(3)complex 单精度和双精度 complex(kind=4) b(4)character character(len=10) c !len为最大长度(5)logical logical*2 : d=.ture. (等价于logical(2):d=.ture.)(6)自定义类型type:类似于C中的structFortran 77中给变量赋初值常用DATA命令,可同时给多个变量赋初值data a,b,string /1, 2.0, fortran/与C
9、不同的是,Fortran中变量不声明也能使用,即有默认类型(跟implicit命令有关)。按照默认的规定,以i,j,k,l,m,n开头的变量被定义为integer,其余为real。取消该设置需在程序声明部分之前加implicit none。彭国伦建议一般都使用该语句。另一点关于声明的不同是Fortran有等价声明:integer a,bequivalence(a,b)使得a,b使用同一块内存。这样可以节省内存;有时可精简代码。如:equivalence(很长名字的变量比如三维数组的某个元素,a),之后使用a来编写程序就简洁多了。2、基本输入输出输入:read(*,*) a !从键盘读入 输出:
10、text在屏幕上输出。Fortran 77用 text。Fortran 90 一般 和都可 print *,只能屏幕输出(*,*)完整写为(unit=*,fmt=*)。其中unit为输入/输出位置,如屏幕,文件等;fmt为格式。如果这两项都写成*,则按默认的方式进行,即上面描述的。print后面的*表示按默认格式输出。四、流程控制1、运算符(1)逻辑运算符= /= = = !Fortran 90用法.EQ. .NE. .GT. .GE. .LT. .LE. !Fortran 77用法(2)涉及相互关系的集合运算符.AND. .OR. .NOT. .EQV. .NEQV. 仅.NOT.连接一个表
11、达式,其余左右两边都要有表达式(可以是logical类型的变量).EQV.:当两边逻辑运算值相同时为真, .NEQV.:当两边逻辑运算值不同时为真2、IF(1) 基本 :if(逻辑判断式) then end if 如果then后面只有一句,可写为if(逻辑判断式) !then和end if可省略(2) 多重判断:else ifelseend if(3) 嵌套: if(逻辑判断式) then else if(逻辑判断式) then else end if !没必要跟C语句多时用了,因为有end if(4) 算术判断:program exampleimplicit nonereal cwrite
12、(*,*) input a numberread (*,*) cif(c) 10,20,30 !10,20和30为行代码,根据c小于/等于/大于0,执行10/20/30行的程序10 write (*,*) Agoto 40 !goto可实现跳到任意前面或后面的行代码处,但用多了破坏程序结构20 write (*,*) Bgoto 4030 write (*,*) C40 stopend3、SELECT CASE类似于C的switch语句select case(变量)case(数值1) ! 比如case(1:5)代表1=变量=5会执行该模块 !case(1,3,5)代表变量等于1或3或5会执行该模块case(数值2) !括号中数值只能是integer,character或logical型常量,不能是real型case defaultend case4、PAUSE, CONTINUEpause暂停程序执行,按enter可继续执行continue貌似没什么用处,可用作封装程序的标志五、循环1、DOdo counter=初值, 终值, 增/减量 !counter的值从初值到终值按增/减量变
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1