FOETRAN 三角形节点程序

C 主程序 
C 
       DIMENSION C(500000),IA(50000),EK(36) 
  

       CHARACTER*12 IN,OUT 
C  IN和OUT为输入文件和输出文件的文件名 
       WRITE(*,*)' ' 
       WRITE(*,*)' PLEASE INPUT THE INPUT-FILE NAME (A<12)' 
       WRITE(*,*)' ' 
       READ(*,5) IN 
C  输入输入文件的文件名 
       WRITE(*,*)' ' 
       WRITE(*,*)' PLEASE INPUT THE OUTPUT-FILE NAME (A<12)' 
       WRITE(*,*)' ' 
       READ(*,5) OUT 
C  输入输出文件的文件名 
  5     FORMAT(A12) 
  8     OPEN(5,FILE=IN, STATUS='OLD') 
       OPEN(6,FILE=OUT,STATUS='UNKNOWN') 
C  打开对应的输入和输出文件 
 10    READ(5,*) NG 
       IF(NG.EQ.0) STOP 
  write(6,15)
 15    FORMAT(//,5X,
     & '*******************  PLANE PROBLEM  *******************',/)    
C  输入结构结点数;如果结点数为0则停止运行 
       READ(5,*)NE,MC,NX,NB,ND,EO,VO,T 
 WRITE(6,20)NG,NE,MC,NX,NB,ND,EO,VO,T
 20    FORMAT(//,10X,'NG=',I4,5X,'NE=',I4,5X,'MC=',I4,5X,'NX=',I4,/,

     &        10X,'NB=',I4,5X,'ND=',I4,/,10X,'EO=',E10.4,5X,'VO=',F7.2,
     &         5X,'T=',F7.2)
C  按顺序输入结构单元数,问题类型参数,载荷组数,给定位移个数 
C   结构刚度阵的半带宽,弹性模量,泊松比和结构厚度 
       READ(5,*)NWA,NWE,NWK,NWP,NWD 
WRITE(6,25)NWA,NWE,NWK,NWP,NWD
  25   FORMAT(/10X,'NWA=',I2,3X,'NWE=',I2,3X,'NWK=',I2,3X,
     &        'NWP=',I2,3X,'NWD=',I2)
C  按顺序输入各输出控制参数 
       NT=2*NG 
C  确定总刚度矩阵阶数NT 

C 
C  计算变界数组的下限 
C 
       M1=3*NE+1 
       M2=M1+2*NB 
       N1=2*NG+1 
       N2=N1+NB 
       N3=7*NE+N2 
       N4=N3+NT*ND 
       N5=N4+NT 
C 得到各变界数组在一维大数组中的起始元素编号 

C 
C 检验实型数组C的下限 
C 
       NEND=N5 
       IF(NEND.LE.500000) GOTO 35 
       WRITE(*,*)'*** EXCEED THE LIMIT OF ARRAY C(IN THE MIDDLE)!! ***'