home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 9 / 09.iso / e / e003 / 14.ddi / MODEL.FOR < prev    next >
Encoding:
Text File  |  1987-04-18  |  33.6 KB  |  428 lines

  1. CC*********************************************************************IMOD00001
  2. CC THE SAP USERS GROUP MAKES NO WARRANTY THAT THE MODEL PROGRAM IS     IMOD00002
  3. CC FREE FROM ERRORS OR DEFECTS AND ASSUMES NO LIABILITY FOR THE        IMOD00003
  4. CC PROGRAM. THE SAP USERS GROUP DISCLAIMS ANY EXPRESS WARRANTY  OR     IMOD00004
  5. CC FITNESS FOR ANY INTENDED USE OR PURPOSE. THE MODEL PROGRAM IS       IMOD00005
  6. CC CONSTANTLY BEING DEVELOPED, MODIFIED AND CHECKED AND ANY KNOWN      IMOD00006
  7. CC ERRORS SHOULD BE REPORTED TO THE SAP USERS GROUP.                   IMOD00007
  8. CC*********************************************************************IMOD00008
  9. CC.....................................................................IMOD00009
  10. CC***************************************************************I      MOD00010
  11. CC                                                               I      MOD00011
  12. CC                                                               I      MOD00012
  13. CC      +++++++++++++++++++++++++++++++++++++++++++++++++        I      MOD00013
  14. CC      +    MAIN PROGRAM ** MODEL ** V0-3 JAN-1982     +        I      MOD00014
  15. CC      +                                               +        I      MOD00015
  16. CC      +++++++++++++++++++++++++++++++++++++++++++++++++        I      MOD00016
  17. CC                                                               I      MOD00017
  18. CC      GENERATE MESH FOR ANY 2 OR 3 DIMENSIONAL SHAPES          I      MOD00018
  19. CC      CREATE THE INPUT FOR THE FOLLOWING PROGRAMS :            I      MOD00019
  20. CC      SAP6 , SAP7 , POST , ?                                   I      MOD00020
  21. CC                                                               I      MOD00021
  22. CC...............................................................I      MOD00022
  23. CC      THIS PROGRAM CAN BE OBTAINED FROM S A P USER GROUP       I      MOD00023
  24. CC      UNIVERSITY OF SOUTHERN CALIFORNIA,UNIVERSITY PARK        I      MOD00024
  25. CC      LOS ANGELES,90007   TEL 213-743-5508                     I      MOD00025
  26. CC                                                               I      MOD00026
  27. CC      COPY RIGHT 1982 SAP USER GROUP    ALL RIGHTS RESERVED    I      MOD00027
  28. CC                                                               I      MOD00028
  29. CC...............................................................I      MOD00029
  30. CC      ACCEPTABL COMMANDS ;                                     I      MOD00030
  31. CC      EN= END                                                  I      MOD00031
  32. CC      SU= SURFACE,INCS,INCT,PERS,PERT                          I      MOD00032
  33. CC      VO= VOLUME,INCS,INCT,INCR,PERS,PERT,PERR                 I      MOD00033
  34. CC      NO= NODE,N,X,Y,Z                                         I      MOD00034
  35. CC      FI= FIXED,N1,N2,INCR,DX,DY,DZ,NM,ND,NR,RX,RY,RZ          I      MOD00035
  36. CC      NW= NWRITE  ;WRITE NODE COORDINATES                             MOD00036
  37. CC      EW= EWRITE  ;WRITE ELEMENT CONNECTION                           MOD00037
  38. CC      EL= ELEMENTS,M,ITYP,IPS,BETA.THICK,MTYP,IEL    ...FOR 2D I      MOD00038
  39. CC      EL= ELEMENTS,M,ITYP,IPS,MTYP,IST,IELD,IELX,NDLS  .FOR 3D I      MOD00039
  40. CC      RE= REPEAT,NTIM,IEL1,IEL2,INC,IADD                       I      MOD00040
  41. CC      RO= ROTATE,NTIM,N1,N2,INC,IADD,DEGREE,AXIS               I      MOD00041
  42. CC      TR= TRANSLATE,NTIM,N1,N2,INC,IADD,DX,DY,DZ               I      MOD00042
  43. CC      SY= SYMET,N1,N2,INC,IADD,X,Y,Z                           I      MOD00043
  44. CC      SA= SAP6                                                 I      MOD00044
  45. CC      SC= SCALE,NTIM,N1,N2,INC,IADD,SX,SY,SZ                   I      MOD00045
  46. CC      ME= MERGE,ZONE                                           I      MOD00046
  47. CC      RN= RNUMBER                                              I      MOD00047
  48. CC      CC= COMMENT CARD                                         I      MOD00048
  49. CC      FIL= FILL,NTIM,N1,N2,NFILL,NSTRT,INC,IADD                I      MOD00049
  50. CC      DN= DNODE,N1,N2,INC                                      I      MOD00050
  51. CC      DI= DIGITISE,I                                           I      MOD00051
  52. CC      LO= LOAD ,LCAE,KEY                                       I      MOD00052
  53. CC      QS= QSURFACE,INCS,INCT,PERS,PERT,KEY                     I      MOD00053
  54. CC          N1,N2,......N8                                       I      MOD00054
  55. CC      QV= QVOLUME,INCS,INCT,INCR,PERS,PERT,PERR,KEY            I      MOD00055
  56. CC          N1,N2,......N20                                      I      MOD00056
  57. CC      TS= TSURFACE,N                                           I      MOD00057
  58. CC          N1,N2,.....N6                                        I      MOD00058
  59. CC      CG= CG,IECHO                                             I      MOD00059
  60. CC      RC= RCOOR,N                                              I      MOD00060
  61. CC                                                               I      MOD00061
  62. CC                                                               I      MOD00062
  63. CC      TAPE REQUIRE BY THE PROGRAM                              I      MOD00063
  64. CC      IN=READ DATA INPUT TO MODEL PROGRAM                      I      MOD00064
  65. CC      ITER=READ FROM TERMINAL FOR INTERACTIVE                  I      MOD00065
  66. CC      IS6= OUTPUT FOR SAP6                                     I      MOD00066
  67. CC      IO=WRITE RESULT AND ERROR MESSAGES                       I      MOD00067
  68. CC      IP=WRITE PLOT FILE FOR POST                              I      MOD00068
  69. CC      INP=WRITE INPUT FOR SAP7 (UNFORMATTED)                   I      MOD00069
  70. CC      INPF = WRITE INPUT FOR SAP7 FORMATTED                    I      MOD00070
  71. CC      FOR015.DAT WRITE BOUNDARY CONDITION FOR POST PROGRAM     I      MOD00071
  72. CC      FOR050.DAT IT IS SCRATCH FILE                            I      MOD00072
  73. CC      FOR051.DAT    "     "     "                              I      MOD00073
  74. CC      FOR052.DAT    "     "     "                              I      MOD00074
  75. CC      FOR058.DAT    "     "     "                              I      MOD00075
  76. CC      FOR017.DAT    "     "     "                              I      MOD00076
  77. CC      FOR022.DAT    "     "     "                              I      MOD00077
  78. CC      FOR060.DAT  READ THE MATERIAL PROPERTIES                 I      MOD00078
  79. CC      FOR026.DAT  READ THE DIGITISE NODES                      I      MOD00079
  80. CC      FOR027.DAT  READ NODES AND COORDINATES (FREE FORMAT)     I      MOD00080
  81. CC      FOR057.DAT  NODEL POINT LOADS                            I      MOD00081
  82. CC---------------------------------------------------------------I      MOD00082
  83. CC      ORIGINAL PROGRAMMER :                                    I      MOD00083
  84. CC                                                               I      MOD00084
  85. CC      DR.FIROOZ GHASSEMI        TEL 213-743 8939               I      MOD00085
  86. CC      ------------------        -----------------              I      MOD00086
  87. CC                                                               I      MOD00087
  88. CC***************************************************************I      MOD00088
  89. CC    PROGRAM NMODEL                                                    MOD##088
  90.       INTEGER AGET,AGETW                                                MOD00089
  91.       LOGICAL IGET,RGET                                                 MOD00090
  92.       COMMON  A(50000) 
  93.       COMMON/NONLIN/INON,ITHIS                                          MOD00092
  94.       COMMON/ECHOT/IECHOT                                               MOD00093
  95.       COMMON/SKE/CANG,SANG                                              MOD00094
  96.       COMMON/TERM1/DUMY26(3)                                            MOD00095
  97.       COMMON/AXES/DUMY27(2)                                             MOD00096
  98.       COMMON/BEAM3/IGEN3                                                MOD00097
  99.       COMMON/A/DUM1(5)                                                  MOD00098
  100.       COMMON/BAND/DUM2(4)                                               MOD00099
  101.       COMMON/ELPAR/DUM3(31)                                             MOD00100
  102.       COMMON/FRECNT/DUM4(97)                                            MOD00101
  103.       COMMON/JACHEK/DUM5(2)                                             MOD00102
  104.       COMMON/NONLIN/DUM6(3)                                             MOD00103
  105.       COMMON/OUTS/DUM7(8)                                               MOD00104
  106.       COMMON/PREP/DUM8(5)                                               MOD00105
  107.       COMMON/SAP/DUM9(3)                                                MOD00106
  108.       COMMON/SIZE/DUM10(5)                                              MOD00107
  109.       COMMON/SUPEL/DUM11(5)                                             MOD00108
  110.       COMMON/TRASH/DUM12(472)                                           MOD00109
  111.       COMMON/SEQUEN/ISEQU,ISTART,INCRE                                  MOD00110
  112.       COMMON/TLOAD1/IDUM13(4),DUM00,MTOT1                               MOD00111
  113.       COMMON/ELARRY/NELAR(4,20)                                         MOD00112
  114.       COMMON/ETITLE/TITLE(20)                                           MOD00113
  115.       COMMON/UNIT/IN,IO,IP,INP,INPF,ITER,IS6,I30,I57,I58,I60            MOD00114
  116.      1,I22,I26,I27                                                      MOD00115
  117.       COMMON/TOTAL/MAXELM,MAXNOD,MAXNDM,MAXGRO                          MOD00116
  118.       COMMON/SAP6/ISAP6                                                 MOD00117
  119.       DIMENSION ICOM(70),ICC(20),COOR(8,3),INDE(8)                      MOD00118
  120.       DATA ICOM(1)/1HE/,ICOM(2)/1HN/                                    MOD00119
  121.       DATA ICOM(3)/1HS/,ICOM(4)/1HU/                                    MOD00120
  122.       DATA ICOM(5)/1HV/,ICOM(6)/1HO/                                    MOD00121
  123.       DATA ICOM(7)/1HN/,ICOM(8)/1HO/                                    MOD00122
  124.       DATA ICOM(9)/1HF/,ICOM(10)/1HI/                                   MOD00123
  125.       DATA ICOM(11)/1HG/,ICOM(12)/1HE/                                  MOD00124
  126.       DATA ICOM(13)/1HE/,ICOM(14)/1HL/                                  MOD00125
  127.       DATA ICOM(15)/1HR/,ICOM(16)/1HE/                                  MOD00126
  128.       DATA ICOM(17)/1HR/,ICOM(18)/1HO/                                  MOD00127
  129.       DATA ICOM(19)/1HT/,ICOM(20)/1HR/                                  MOD00128
  130.       DATA ICOM(21)/1HM/,ICOM(22)/1HE/                                  MOD00129
  131.       DATA ICOM(23)/1HS/,ICOM(24)/1HC/                                  MOD00130
  132.       DATA ICOM(25)/1HC/,ICOM(26)/1HC/                                  MOD00131
  133.       DATA ICOM(27)/1HN/,ICOM(28)/1HW/                                  MOD00132
  134.       DATA ICOM(29)/1HE/,ICOM(30)/1HW/                                  MOD00133
  135.       DATA ICOM(31)/1HR/,ICOM(32)/1HN/                                  MOD00134
  136.       DATA ICOM(33)/1HD/,ICOM(34)/1HN/                                  MOD00135
  137.       DATA ICOM(35)/1HD/,ICOM(36)/1HI/                                  MOD00136
  138.       DATA ICOM(37)/1HS/,ICOM(38)/1HY/                                  MOD00137
  139.       DATA ICOM(39)/1HS/,ICOM(40)/1HA/                                  MOD00138
  140.       DATA ICOM(41)/1HL/,ICOM(42)/1HO/                                  MOD00139
  141.       DATA ICOM(43)/1HN/,ICOM(44)/1HE/                                  MOD00140
  142.       DATA ICOM(45)/1HT/,ICOM(46)/1HS/                                  MOD00141
  143.       DATA ICOM(47)/1HR/,ICOM(48)/1HC/                                  MOD00142
  144.       DATA ICOM(49)/1HC/,ICOM(50)/1HG/                                  MOD00143
  145.       DATA ICOM(51)/1HQ/,ICOM(52)/1HS/                                  MOD00144
  146.       DATA ICOM(53)/1HQ/,ICOM(54)/1HV/                                  MOD00145
  147.       DATA ICOM(55)/1HM/,ICOM(56)/1HA/                                  MOD00146
  148.       DATA ICOM(57)/1HJ/,ICOM(58)/1HA/                                  MOD00147
  149.       DATA ICOM(59)/1HB/,ICOM(60)/1HE/                                  MOD00148
  150.       DATA ICOM(61)/1HP/,ICOM(62)/1HU/                                  MOD00149
  151.       DATA ICOM(63)/1HS/,ICOM(64)/1HE/                                  MOD00150
  152.       DATA ICOM(65)/1HE/,ICOM(66)/1HC/                                  MOD00151
  153.       DATA BLANK/1H /,LLL/1HL/,HELP/4HHELP/                             MOD00152
  154.       CALL OPENS                                                        MOD00153
  155.       MAXELM=0                                                          MOD00154
  156.       ISEQU=0                                                           MOD00155
  157.       IECHOT=0                                                          MOD00156
  158.       ISTART=0                                                          MOD00157
  159.       INCRE=0                                                           MOD00158
  160.       MTOT1=0                                                           MOD00159
  161.       IGEN3=0                                                           MOD00160
  162.       IZERO=0                                                           MOD00161
  163.       MAXNOD=0                                                          MOD00162
  164.       MAXGRO=0                                                          MOD00163
  165.       INON=0                                                            MOD00164
  166.       ITHIS=0                                                           MOD00165
  167.       ICGDON=0                                                          MOD00166
  168.       IJADON=0                                                          MOD00167
  169.       MAXNDM=8                                                          MOD00168
  170.       REWIND 50                                                         MOD00169
  171.       REWIND 52                                                         MOD00170
  172.       REWIND 51                                                         MOD00171
  173.       CALL STIME                                                        MOD00172
  174.       CALL GETTIM(NH1,NM1,NS1,NUM1)
  175.       DO 66 I=1,80                                                      MOD00174
  176. 66    DUM4(I)=BLANK
  177.       IN=5                                                              MOD00176
  178.       I22=22                                                            MOD00177
  179.       I26=26                                                            MOD00178
  180.       I27=27                                                            MOD00179
  181.       ITER=5                                                            MOD00180
  182.       IS6=19                                                            MOD00181
  183.       IO=6                                                              MOD00182
  184.       IP=14                                                             MOD00183
  185.       INP=21                                                            MOD00184
  186.       I60=60                                                            MOD00185
  187.       INPF=20                                                           MOD00186
  188.       I30=30                                                            MOD00187
  189.       I57=57                                                            MOD00188
  190.       I58=58                                                            MOD00189
  191.       REWIND I57                                                        MOD00190
  192.       MTOTAL=50000 
  193.       ISAP6=0                                                           MOD00192
  194.       IFILE=0                                                           MOD00193
  195.       READ(IN,90)TITLE                                                  MOD00194
  196. 90    FORMAT(20A4)                                                      MOD00195
  197. CC    IF(TITLE(1).NE.HELP)GO TO 80                                      MOD00196
  198. CC    CALL HELP1(ICOM)                                                  MOD00197
  199. CC    READ(IN,90)TITLE                                                  MOD00198
  200. 80    CALL GETNL(GET001)                                                MOD00199
  201.       IF(IGET(LIMIT))CALL PRTERR(0)                                     MOD00200
  202.       IF(IGET(NUMEL))CALL PRTERR(0)                                     MOD00201
  203.       IF(LIMIT.LE.0)CALL PRTERR(1)                                      MOD00202
  204.       IF(NUMEL.LE.0)CALL PRTERR(1)                                      MOD00203
  205.       N1=1                                                              MOD00204
  206.       N2=N1+LIMIT                                                       MOD00205
  207.       N3=N2+LIMIT                                                       MOD00206
  208.       N4=N3+LIMIT                                                       MOD00207
  209.       N5=N4+LIMIT                                                       MOD00208
  210.       N6=N5+6*LIMIT                                                     MOD00209
  211.       N7=N6+NUMEL                                                       MOD00210
  212.       N8=N7+NUMEL                                                       MOD00211
  213.       N9=N8+NUMEL                                                       MOD00212
  214.       N10=N9+NUMEL                                                      MOD00213
  215.       N11=N10+NUMEL                                                     MOD00214
  216.       N12=N11+20*NUMEL                                                  MOD00215
  217.       N13=N12+LIMIT                                                     MOD00216
  218.       N14=N13+LIMIT                                                     MOD00217
  219.       N15=N14+NUMEL                                                     MOD00218
  220.       N16=N15+NUMEL                                                     MOD00219
  221.       N17=N16+NUMEL                                                     MOD00220
  222.       N18=N17+LIMIT                                                     MOD00221
  223.       N19=N18+NUMEL                                                     MOD00222
  224.       N20=N19+NUMEL                                                     MOD00223
  225.       N21=N20+NUMEL                                                     MOD00224
  226.       N22=N21+NUMEL                                                     MOD00225
  227.       N23=N22+LIMIT*2                                                   MOD00226
  228.       N24=N23+LIMIT                                                     MOD00227
  229.       IF(N24.GT.MTOTAL)CALL PRTERR(8)                                   MOD00228
  230.       N24X=N24-1                                                        MOD00229
  231.       DO 98 I=1,N24X                                                    MOD00230
  232. 98    A(I)=0.                                                           MOD00231
  233. 100   CALL GETNL(GET001)                                                MOD00232
  234.       I=AGETW(AG001)                                                    MOD00233
  235.       K=AGET(2)                                                         MOD00234
  236.       KEY=1                                                             MOD00235
  237.       DO 110 J=1,65,2                                                   MOD00236
  238.       L=J+1                                                             MOD00237
  239.       IF(ICOM(J).EQ.I.AND.ICOM(L).EQ.K)GO TO 120                        MOD00238
  240.       KEY=KEY+1                                                         MOD00239
  241. 110   CONTINUE                                                          MOD00240
  242.       WRITE(IO,111)I,K                                                  MOD00241
  243. 111   FORMAT('  ++ F ++ CONTROL KEY = ',A1,A1,'  IS ILLEGAL')           MOD00242
  244.       STOP                                                              MOD00243
  245. 120   GO TO (1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,1,17,18,19,20          MOD00244
  246.      1,21,22,23,24,1,26,27,28,1,30,31,32,33),KEY                        MOD00245
  247. 1     IF(ICGDON.EQ.1)GO TO 1016                                         MOD00246
  248.       IF(IJADON.EQ.1)GO TO 1016                                         MOD00247
  249.       IF(IFILE.EQ.0)GO TO 214                                           MOD00248
  250. 223   WRITE(IP,233)MAXELM                                               MOD00249
  251. 233   FORMAT(5X,I5)                                                     MOD00250
  252.       REWIND 50                                                         MOD00251
  253.       DO 234 MM=1,MAXELM                                                MOD00252
  254.       READ(50)MT1,MT2,MH,DENS,AREA                                      MOD00253
  255.       WRITE(IP,222)MT1,MT2,DENS,AREA                                    MOD00254
  256. 222   FORMAT(5X,I5,6X,I5,2E12.5)                                        MOD00255
  257.       READ(50)M1,(ICC(KK),KK=1,MH)                                      MOD00256
  258.       WRITE(IP,228)MM,(ICC(KK),KK=1,MH)                                 MOD00257
  259. 228   FORMAT(5X,I5,20(I5))                                              MOD00258
  260. 234   CONTINUE                                                          MOD00259
  261.       WRITE(IP,222)IZERO,IZERO                                          MOD00260
  262.       IF(MAXELM.NE.NUMEL)WRITE(IO,236)MAXELM                            MOD00261
  263.       IF(MAXELM.EQ.NUMEL)WRITE(IO,300)MAXELM                            MOD00262
  264. 300   FORMAT(' . . .    TOTAL ELEMENTS GENERATED ARE = ',I5)            MOD00263
  265. 236   FORMAT(' ++ W ++  TOTAL ELEMENTS GENERATED = ',I5)                MOD00264
  266. 235   FORMAT(A4)                                                        MOD00265
  267.       IF(ISAP6.NE.1)GO TO 2010                                          MOD00266
  268.       WRITE(IO,1020)                                                    MOD00267
  269.       CALL FCOPY1                                                       MOD00268
  270. 2010  IF(KEY.NE.16)GO TO 1000                                           MOD00269
  271.       IF(INON.EQ.1)GO TO 1000                                           MOD00270
  272. 2030  WRITE(IO,250)                                                     MOD00271
  273. 250   FORMAT(1X,/,' . . .    START RENUMBERING THE MESH   . . .')       MOD00272
  274.       CALL RENUM1                                                       MOD00273
  275. 1000  IF(ISAP6.EQ.1)GO TO 1015                                          MOD00274
  276. 1020  FORMAT(' . . . .  CREATING INPUT FOR SAP6')                       MOD00275
  277. 1015  IF(INON.EQ.0)GO TO 1016                                           MOD00276
  278.       WRITE(IO,1017)                                                    MOD00277
  279. 1017  FORMAT(' * * *    START RENUMBERING THE MESH FOR NON-LINEAR       MOD00278
  280.      1ELEMENTS * * *')                                                  MOD00279
  281. CC    CALL NODERE(LIMIT,A(N17))                                         MOD00280
  282. 1016  IF(KEY.EQ.25)GO TO 25                                             MOD00281
  283.       IF(KEY.EQ.29)GO TO 29                                             MOD00282
  284.       IF(ISAP6.EQ.1)GO TO 1027                                          MOD00283
  285. 1010  REWIND 52                                                         MOD00284
  286.       WRITE(I30,1018)TITLE                                              MOD00285
  287. 1018  FORMAT(20A4)                                                      MOD00286
  288.       WRITE(I30,1019)MAXNOD,MAXELM,MAXGRO,MTOT1                         MOD00287
  289. 1019  FORMAT(4I10)                                                      MOD00288
  290.       DO 1021 MB=1,MAXGRO                                               MOD00289
  291.       READ(52)NNEE,NNTYPE,NONTYP                                        MOD00290
  292.       WRITE(I30,1022)NNEE,NNTYPE,NONTYP                                 MOD00291
  293. 1022  FORMAT(2I10,I2)                                                   MOD00292
  294. 1021  CONTINUE                                                          MOD00293
  295. 1027  CONTINUE                                                          MOD00294
  296.       CALL GETTIM(NH2,NM2,NS2,NUS2)
  297.       NH2=NH2-NH1
  298.       NM2=NM2-NM1
  299.       NS2=NS2-NS1
  300.       WRITE(IO,1001)NH2,NM2,MS2
  301. 1001   FORMAT(//,' . . . .  ',I2,':',I2,':',I2,' TIME WERE USED TO GENER
  302.      1ATE THIS MESH ',' . . . .')                                       
  303.       CALL CFILES
  304.       STOP                                                              MOD00301
  305. 14    IFILE=1                                                           MOD00302
  306. 214   CALL PRCOOR(LIMIT,A(N1),A(N2),A(N3),A(N4),A(N5),A(N23))           MOD00303
  307.       IF(IFILE.EQ.0)GO TO 215                                           MOD00304
  308.       GO TO 100                                                         MOD00305
  309. 15    IFILE=1                                                           MOD00306
  310. 215   CALL PRELEM(LIMIT,NUMEL,A(N6),A(N7),A(N8),A(N9),A(N10),A(N11)     MOD00307
  311.      1,A(N14),A(N15),A(N16),A(N18),A(N19),A(N20)                        MOD00308
  312.      &,A(N21),A(N22),A(N23))                                            MOD00309
  313.       IF(IFILE.EQ.0)GO TO 223                                           MOD00310
  314.       GO TO 100                                                         MOD00311
  315. 2     CALL PRE2(2,LIMIT,NUMEL,A(N1),A(N2),A(N3),A(N4),A(N5))            MOD00312
  316.       GO TO 100                                                         MOD00313
  317. 3     CALL PRE2(3,LIMIT,NUMEL,A(N1),A(N2),A(N3),A(N4),A(N5))            MOD00314
  318.       GO TO 100                                                         MOD00315
  319. 4     CALL NODE1(LIMIT,A(N1),A(N2),A(N3),A(N4))                         MOD00316
  320.       GO TO 100                                                         MOD00317
  321. 5     MM=AGET(3)                                                        MOD00318
  322.       IF(MM.EQ.LLL)GO TO 205                                            MOD00319
  323.       CALL READBN(LIMIT,A(N5))                                          MOD00320
  324.       GO TO 100                                                         MOD00321
  325. 205   CALL FILL1(LIMIT,A(N1),A(N2),A(N3),A(N4))                         MOD00322
  326.       GO TO 100                                                         MOD00323
  327. 6     STOP                                                              MOD00324
  328. 7     CALL ELEM1(LIMIT,NUMEL,A(N6),A(N7),A(N8),A(N9),A(N10),A(N11)      MOD00325
  329.      1,A(N14),A(N15),A(N16),A(N17),A(N18),A(N19)                        MOD00326
  330.      &,A(N20),A(N21))                                                   MOD00327
  331.       GO TO 100                                                         MOD00328
  332. 8     CALL REPEA1(LIMIT,NUMEL,A(N6),A(N7),A(N8),A(N9),A(N10),A(N11)     MOD00329
  333.      1,A(N14),A(N15),A(N16),A(N17),A(N18)                               MOD00330
  334.      &,A(N19),A(N20),A(N21))                                            MOD00331
  335.       GO TO 100                                                         MOD00332
  336. 9     CALL ROTAT1(LIMIT,A(N1),A(N2),A(N3),A(N4))                        MOD00333
  337.       GO TO 100                                                         MOD00334
  338. 10    CALL TRANS1(LIMIT,A(N1),A(N2),A(N3),A(N4))                        MOD00335
  339.       GO TO 100                                                         MOD00336
  340. 11    CALL MERGE1(LIMIT,NUMEL,A(N1),A(N2),A(N3),A(N4),A(N5),A(N11)      MOD00337
  341.      1,A(N12),A(N13),A(N17))                                            MOD00338
  342.       GO TO 100                                                         MOD00339
  343. 12    CALL SCALE1(LIMIT,A(N1),A(N2),A(N3),A(N4))                        MOD00340
  344. 13    GO TO 100                                                         MOD00341
  345. 17    CALL DNODE1(LIMIT,NUMEL,A(N1),A(N2),A(N3),A(N4),A(N5),A(N11))     MOD00342
  346.       GO TO 100                                                         MOD00343
  347. 18    IF(IGET(IREST))CALL PRTERR(0)                                     MOD00344
  348.       CALL RECOOR(LIMIT,NUMEL,A(N1),A(N2),A(N3),A(N4),A(N11),A(N12)     MOD00345
  349.      1,A(N13),A(N6),A(N7),A(N8),A(N9),A(N10),A(N14))                    MOD00346
  350.       IF(IREST.EQ.1)GO TO 100                                           MOD00347
  351. CC    CALL DIGIT1(LIMIT,NUMEL,A(N1),A(N2),A(N3),A(N4),A(N12)            MOD00348
  352. CC   1,A(N13),A(N6),A(N7),A(N8),A(N9),A(N10),A(N11),A(N14))             MOD00349
  353. CC    GO TO 100                                                         MOD00350
  354. 19    CALL SYMET1(LIMIT,A(N1),A(N2),A(N3),A(N4))                        MOD00351
  355.       GO TO 100                                                         MOD00352
  356. 20    ISAP6=1                                                           MOD00353
  357.       GO TO 100                                                         MOD00354
  358. 21    MM=N17-1                                                          MOD00355
  359.       DO 279 I=1,MM                                                     MOD00356
  360.       A(I)=0.0                                                          MOD00357
  361. 279   CONTINUE                                                          MOD00358
  362.       CALL LOAD1(LIMIT,NUMEL)                                           MOD00359
  363.       GO TO 100                                                         MOD00360
  364. 22    INON=1                                                            MOD00361
  365.       ITHIS=1                                                           MOD00362
  366.       GO TO 7                                                           MOD00363
  367. 23    CALL TRIAN1(LIMIT,A(N1),A(N2),A(N3),A(N4))                        MOD00364
  368.       GO TO 100                                                         MOD00365
  369. 24    CALL RECO1(LIMIT,A(N1),A(N2),A(N3),A(N4))                         MOD00366
  370.       GO TO 100                                                         MOD00367
  371. 25    IF(IGET(IECHO))CALL PRTERR(0)                                     MOD00368
  372.       CALL CG1(IECHO)                                                   MOD00369
  373.       ICGDON=1                                                          MOD00370
  374.       GO TO 100                                                         MOD00371
  375. 26    CALL PRE2(4,LIMIT,NUMEL,A(N1),A(N2),A(N3),A(N4),A(N5))            MOD00372
  376.       GO TO 100                                                         MOD00373
  377. 27    CALL PRE2(5,LIMIT,NUMEL,A(N1),A(N2),A(N3),A(N4),A(N5))            MOD00374
  378.       GO TO 100                                                         MOD00375
  379. 28    CALL MATRI1                                                       MOD00376
  380.       GO TO 100                                                         MOD00377
  381. 29    CALL JACOB1                                                       MOD00378
  382.       IJADON=1                                                          MOD00379
  383.       GO TO 100                                                         MOD00380
  384. 30    IF(IGET(IGEN3))CALL PRTERR(0)                                     MOD00381
  385.       GO TO 100                                                         MOD00382
  386. 31    CALL PUSH1(LIMIT,A(N2),A(N22))                                    MOD00383
  387.       GO TO 100                                                         MOD00384
  388. 32    IF(IGET(ISTART))CALL PRTERR(0)                                    MOD00385
  389.       IF(IGET(INCRE))CALL PRTERR(0)                                     MOD00386
  390.       IF(ISTART.LE.0)ISTART=1                                           MOD00387
  391.       IF(INCRE.EQ.0)INCRE=1                                             MOD00388
  392.       ISEQU=1                                                           MOD00389
  393.       GO TO 100                                                         MOD00390
  394. 33    IECHOT=1                                                          MOD00391
  395.       GO TO 100                                                         MOD00392
  396.       END                                                               MOD00393
  397.       BLOCK DATA                                                        MOD00394
  398.       INTEGER SECT,POINT,BEGIN,GETSEC,EPOINT,BLANK,AGET,AGETW           MOD00395
  399.       LOGICAL EOL,EOS,EOF,ERROR,INDEX,GETWRD,IGET,RGET,ERR1             MOD00396
  400.       COMMON/ELARRY/NELAR(4,20)                                         MOD00397
  401.       COMMON  /FRECNT/ LINE(80),SECT,EOL,EOS,EOF,                       MOD00398
  402.      1ERROR,BLANK,ICOMMA,POINT,BEGIN,LENGTH,                            MOD00399
  403.      2EPOINT,NSECT,LINENM,MAXSTR,KUNIT,JUNIT,IUNIT                      MOD00400
  404.       COMMON/FRECNM/MULTIP                                              MOD00401
  405.       DATA MULTIP/0/                                                    MOD00402
  406.       DATA BLANK/1H /,ICOMMA/1H,/,ERROR/.FALSE./,                       MOD00403
  407.      2POINT/1/,BEGIN/1/,LENGTH/0/,EPOINT/1/,                            MOD00404
  408.      3LINE/80*1H /,MAXSTR/1/,LINENM/1/                                  MOD00405
  409.       DATA NELAR/2,2,6,2,                                               MOD00406
  410.      13,2,12,26,                                                        MOD00407
  411.      14,4,12,8,                                                         MOD00408
  412.      14,4,8,4,                                                          MOD00409
  413.      18,8,33,54,                                                        MOD00410
  414.      14,4,42,12,                                                        MOD00411
  415.      11,1,1,1,                                                          MOD00412
  416.      14,4,8,4,                                                          MOD00413
  417.      13,2,12,39,                                                        MOD00414
  418.      120,20,60,54,                                                      MOD00415
  419.      18,8,16,52,                                                        MOD00416
  420.      18,8,16,52,                                                        MOD00417
  421.      18,8,16,52,                                                        MOD00418
  422.      14,1,6,6,                                                          MOD00419
  423.      18,8,48,48,                                                        MOD00420
  424.      19,9,27,52,                                                        MOD00421
  425.      120,20,40,40,                                                      MOD00422
  426.      116,16,32,52,8*0/                                                  MOD00423
  427.       END                                                               MOD00424
  428.