实战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)/3
        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**2
      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**2
      end

2. 算法

• if：

      program readGrade
      implicit none
      integer num,grade,i,j
      i=3  
      j=1
      print *,"请输入",i,"个人的学号和成绩："  
10    read *,num,grade
      if(grade .GT. 80) then
          print *,num,grade
      end if
      j=j+1
      if (j .LE. i) goto 10     
      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=0
      !下面是默认的读入方式
      print *,'请输入年份，为整数：'
      read *,year
      if ( MOD(year,4) .EQ. 0) then
            if (MOD(year,100) .NE. 0)  then
               print *,year,'年是闰年,1'
            else
               print *,year,'年不是闰年,1'     
            endif
      else
            if (MOD(year,400) .EQ. 0) 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 10 i=1, 100
            if (mod(i,2) .EQ. 0) then
                  sum=sum-1.0/i
            else
                  sum=sum+1.0/i
 10   end if
      print *,sum

      end

      program computeSum
      implicit none
      real sum
      integer i
      sum=0.0
      do 10, i=1, 100,1!10后面的逗号要不要都行
            if (mod(i,2) .EQ. 0) then
                  sum=sum-1.0/i
            else
                  sum=sum+1.0/i
 10         end if!fortran77中不让这样，但也可以执行，兼容了
      print *,sum
            sum=0
      do 20, i=1, 100,1!10后面的逗号要不要都行
                  if (mod(i,2) .EQ. 0) then
                        sum=sum-1.0/i
                  else
                        sum=sum+1.0/i
                  end if
20    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(8,3)
      print *,sign(8.5,-2.0)!将后一个数的符号给前一个数
      print *,max(1.0,3.0,9.0,11.0,22.0,100.0)!求一系列数的最大值，但类型必须相同，都是实数
      print *,max(1,2**10,9,11,22,100)!类型必须相同，都是实数,2**10代表2的10次方
      print *,min(1,3,9,11,22,100)
      !real 只能表示7位数
      print *,0.001+123456.0-123455.0
      print *,123456.0-123455.0+0.001
      print *,1.0/3.0*3.0
      print *,1,2
      print *,sqrt(4.0),sqrt(9.0)!必须为实型数据
      print *,sqrt(a),sqrt(b)
      write(*,*) '123456789012345678901234567890'!整数占11列，两个整数占23列，说明中间有个空格

      end

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

 

      program arrayStudy
      double precision  a(1:100000000)!一维数组，包含10个整形数组
      double precision b(1:3,4:6)!二维数组，3行3列
      !dimension b(1:3,4:6)!和上面声明效果一样，不过是多写了一行
      !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 10 ,i=1,100000000
            a(i)=i**2
            !print *,a(i)
   10 continue
      print *,a(100000000)
      end

　

      Program ex1
      print *,'Hello world!'
      end