实战fortran77基础语法2

由于大量的接触fortran77，先学习fortran77吧

1. 1.fortran不区分大小写，fortran77一般全部大写，感觉没必要，不好读，还是用大小写混合着写吧。后缀 .for 为固定格式的代码，比如fortran77的代码语法就是。另一种以 .f90 为后缀，是自由格式。

      Program ex1
      print *,'Hello world!'
      end

• 　最好在程序开头写上implicit none语句，作用是不按照默认的数据类型。要是默认的则I和N之间所有字母开头的变量都代表这个变量是整数类型。

         Program ex1
         implicit none
         real x,y,z,aver
         x=90.5
         y=87.0
         z=60.5
         aver=(x+y+z)/
         print *,'the three number are:',x,y,z
         print *,'average=',aver
         end

  但是注意所有变量要在所有执行语句之前声明。

• 调用子函数：

在子函数中声明与子函数同名的变量，记得要返回的值赋值给这个与子函数同名的变量。

 c     Main program,to calculate the area of the circle
       Program ex2
       implicit none
       real x,y,z,areaofcircle
       x=90.5
       y=87.0
       z=60.5
       print *,'r=',x,'area of circle=',areaofcircle(x)
       print *,'r=',y,'area of circle=',areaofcircle(y)
       print *,'r=',z,'area of circle=',areaofcircle(z)
       end
 
       function  areaofcircle(r)
       implicit none
       real r ,areaofcircle
       areaofcircle=3.1415926*r**
       return
       end

从上面可以看出，子函数更像数学里一个公式，其实areaofcircle可以用一个更简单的f代替，就更好理解啦。

• 子例行函数

 c     Main program,to calculate the area of the circle
       Program ex2
       implicit none
       real x,y,z,xarea,yarea,zarea
       x=90.5
       y=87.0
       z=60.5
 c     必须有call
       call areaofcircle(x,xarea)
       call areaofcircle(y,yarea)
       call areaofcircle(z,zarea)
       print *,'r=',x,'area of circle=',xarea
       print *,'r=',y,'area of circle=',yarea
       print *,'r=',z,'area of circle=',zarea
       end
 c     通过传址方式回传数值，利用了area的地址
       subroutine  areaofcircle(r,area)
       implicit none
 c     含有implicit none语句，那么本子程序所有变量都要显示声明
       real r,area
       area=3.1415926*r**
       end

2. 算法

• if：

       program readGrade
       implicit none
       integer num,grade,i,j
       i=
       j=
       print *,"请输入",i,"个人的学号和成绩："
     read *,num,grade
       if(grade .GT. ) then
           print *,num,grade
       end if
       j=j+
       if (j .LE. i) goto
       end

尽量不要用goto这种跳转，多了逻辑会很乱。可以用do 代替。

由于visual studio 2013安装的inter 的fortran编译器兼容性比较好，在.for中写fortran90的语法也不报错，所以，大于号也可以直接写 >,比如：      if(grade > 80) then

判断是否是闰年：

 C     判断是否是闰年，这种c开头的注释是fortran77语法
       Program isYear
       implicit none
       integer year!year是整数类型； 以 ! 开头的语句或!之后的部分是fortran90语法的注释
 !下面是执行语句
       year=
       !下面是默认的读入方式
       print *,'请输入年份，为整数：'
       read *,year
       if ( MOD(year,) .EQ. ) then
             if (MOD(year,) .NE. )  then
                print *,year,'年是闰年,1'
             else
                print *,year,'年不是闰年,1'
             endif
       else
             if (MOD(year,) .EQ. ) then
                   print *,year,"年是闰年,2"
             else
                   print *,year,'年不是闰年,2'
             endif
       endif
 
       end

• 计算1-1/2+1/3-...-1/100：

       program computeSum
       implicit none
       real sum
       integer i
       sum=0.0
       do  i=,
             if (mod(i,) .EQ. ) then
                   sum=sum-1.0/i
             else
                   sum=sum+1.0/i
     end if
       print *,sum
 
       end

       program computeSum
       implicit none
       real sum
       integer i
       sum=0.0
       do , i=, ,!10后面的逗号要不要都行
             if (mod(i,) .EQ. ) then
                   sum=sum-1.0/i
             else
                   sum=sum+1.0/i
           end if!fortran77中不让这样，但也可以执行，兼容了
       print *,sum
             sum=
       do , i=, ,!10后面的逗号要不要都行
                   if (mod(i,) .EQ. ) then
                         sum=sum-1.0/i
                   else
                         sum=sum+1.0/i
                   end if
     continue!最好这样
             print *,sum
 
       end

优化

--------------

       program mymath
       !double 比 real 精度要高
       double precision a,b!fortran77语法
       a=11111.1
       b=1111.11
       print *,a*b
       print *,sin(1.0)
       print *,asin(1.0)
       print *,log(3.0)
       print *,int(8.6)
       print *,mod(,)
       print *,sign(8.5,-2.0)!将后一个数的符号给前一个数
       print *,max(1.0,3.0,9.0,11.0,22.0,100.0)!求一系列数的最大值，但类型必须相同，都是实数
       print *,max(,**,,,,)!类型必须相同，都是实数,**10代表2的10次方
       print *,min(,,,,,)
       !real 只能表示7位数
       print *,0.001+123456.0-123455.0
       print *,123456.0-123455.0+0.001
       print *,1.0/3.0*3.0
       print *,,
       print *,sqrt(4.0),sqrt(9.0)!必须为实型数据
       print *,sqrt(a),sqrt(b)
       write(*,*) ''!整数占11列，两个整数占23列，说明中间有个空格
 
       end

--------------------------

       program arrayStudy
       double precision  a(:)!一维数组，包含10个整形数组
       double precision b(:,:)!二维数组，3行3列
       !dimension b(:,:)!和上面声明效果一样，不过是多写了一行
       !double precision b
       DATA C,D,E,F/-1.0,2.0,3.0,4.0/
       DATA G/2.0/,H/3.0/
 
       integer i
       do  ,i=,
             a(i)=i**
             !print *,a(i)
     continue
       print *,a()
       end

　

       Program ex1
       print *,'Hello world!'
       end

fortran77.for