home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 9 / 09.iso / e / e003 / 2.ddi / SAP6PC1.FOR < prev   
Encoding:
Text File  |  1980-01-04  |  32.4 KB  |  422 lines

  1.       PROGRAM SAP6P1                                                    R0001101
  2.       IMPLICIT REAL*8(A-H,O-Z)                                          00001100
  3.       LOGICAL ELPRT,ELPCH,GENPRT,GENPCH,DEFPCH,GEOST                    00001110
  4.       COMMON/BMDATA/IAISC,ILDTYP,ILDMUL,IBMSEC,NLWAVE,ILOCAL            00001120
  5.       COMMON/EQUILB/NEQIL,NX43                                          00001130
  6.       COMMON/CTL/ELPRT,ELPCH,GENPRT,GENPCH                              00001140
  7.       COMMON / JUNK / DUK(227)                                          00001150
  8.       REAL*8  NPAR                                                      00001160
  9.       COMMON /QTSARG/ QQQ(1000)                                         00001170
  10.       COMMON/DYN3/ NEIG,NAD,ANORM,NVV,NFO                               00001180
  11.       COMMON /DYN4/KDYN,NRESS,NCRD,NCWT,NRESS1                          00001190
  12.       COMMON /DYN5/FRSHFT,FRINIT,FREND,MODEFR,NOSS                      00001200
  13.        COMMON /TAPES/NSTIF,NRED,NL,NR,NT,NMASS                          00001210
  14.       COMMON /EXTRA/MODEX,NT8,N10SV,NT10,KEQB,NY,T(10)                  00001220
  15.       COMMON/GEOSTF/GEOST,NELGEO                                        00001230
  16.       COMMON/MASS/LMASS                                                 00001240
  17.       COMMON/MATL/MATLCO                                                R0001241
  18.       COMMON/SIZE/NDMX,MXDF,NSMX,NTERM,NADND                            00001250
  19.       COMMON/SLVE/NSLAVE                                                00001260
  20.       COMMON /ELPAR/ NPAR(14),NUMNP,MBAND,NELTYP,N1,N2,N3,N4,N5,MTOT,NEQ00001270
  21.      $ ,N2P,N3P,NMRI,NTRI,N1P,NUMEL,NUMEL2,KZ(20,2),NEMN                00001280
  22.       COMMON / MISC / NBLOCK,NEQB,LL,NFREQ,LB                           00001290
  23.       COMMON/AMB/ GRAV,REFT,JROT                                        00001300
  24.       COMMON /PREP/ XMX,XAD,KSKIP,NDYN,I1,I99,POS,PRTCOD                00001310
  25.      $              ,POSSAV,PRTOFF,PRTON,PRTDUM,IDIRC                   00001320
  26.       COMMON /DYN/ IFIL1(11),IFIL2                                      00001330
  27.       COMMON/ELARRY/NELAR(4,20)                                         00001340
  28.       COMMON /ELTEMP/ SET1(103)                                         00001350
  29.      $       /OUT/KSET2(6),KELRST,MAXDF,IFIL3(2)                        00001360
  30.      $       /SQZ/    SET3,LIST,LISTC,LISTB,LISTA                       00001370
  31.      $       /TRASH/  SET4(490)                                         00001380
  32.      $       /GPS/    SET5(10)                                          00001390
  33.      $       /CG/     SET6(4),RFIL1(2)                                  00001400
  34.      $       /TAPES/ SET7(6)                                            00001410
  35.      $       /DYN2/KSET8(3),NFVC,SET8(12)                               00001420
  36.       COMMON /WORDS/ NWDS(30,2)                                         00001430
  37.       COMMON /BAND/  NRNM(3),IRSK,IFIL4(4)                              00001440
  38.       COMMON /HEADIN/TITLE1(20),TITLE2(5),TITLE3(10)                    00001450
  39.       COMMON /SUPEL/ NSELEM,NEQL,NODESE                                 00001460
  40.      $,KSET(3)                                                          00001470
  41.       COMMON/FORCE/ NLC,NELD                                            00001480
  42.         COMMON/DAPG/NQ1,NQX,DUMY(210)                                   00001490
  43.       COMMON/VAXPLT/IVPLT,XVB(6),XVA(4),MINX,MINY,LEN                   00001500
  44.       COMMON /ICM/ICOMP,MMRI,MTRI,M1P,M2P,M3P                           00001510
  45.       COMMON/PLOTH/IPLT,IPLWRT                                          00001520
  46.       COMMON/COMMT/NTYP,NUMET,NNRRC,NNRRC1                              R0001521
  47.       COMMON/COMMT1/NDSSS,KDSSS,NTY,NSLDM,NBLANK,MDYN,NE2B,KRK1,MCB,MLT R0001522
  48.       DIMENSION KZN(20),ZD(31)                                          00001530
  49.       DIMENSION NEXPDT(2),NOWDTE(4)                                     00001540
  50.       COMMON A(1)                                                       00001550
  51.       DATA KZN/2,7,1,2,7,7,10,7,21,1177,6,51,10,21,9,7,8,5,0,0/         00001560
  52.       DATA NEXPDT/78,222/                                               00001570
  53.       DATA CONKEY/8HSUBSTRUC/                                           R0001571
  54.       WRITE (*,990)                                                     R0001611
  55.   990 FORMAT (5X,'*********** PROGRAM SAP6P1 STARTING *********')       R0001612
  56. CC    CALL OPEN1                                                        00001580
  57. CCR    CALL TIME(NOWDTE)                                                00001590
  58. CCR   CALL TODATE                                                       00001600
  59.       READ (5,2001) COKEY1                                              R0001601
  60.       WRITE (6,2001) COKEY1                                             R0001602
  61.  2001 FORMAT (A8)                                                       R0001603
  62.  2002 FORMAT (5X,'***** SUBSTRUCT STARTING *****'/)                     R0001604
  63.       IF (CONKEY .EQ. COKEY1) CALL COMMRW (1)                           R0002241
  64.       IF (CONKEY .NE. COKEY1) REWIND 5                                  R0002242
  65.       IF (CONKEY .EQ. COKEY1) WRITE (6,2002)                            R0002244
  66.       IF (CONKEY .EQ. COKEY1) GO TO 100                                 R0002243
  67.       CALL SIZER1                                                       00001610
  68. CCR    CALL ERRSET(208,0,-1,1,1)                                        00001620
  69. CCR   CALL CPUINT                                                       00001630
  70.       NX43=1700000                                                      00001640
  71. CCR   CALL FDATE(NOWDTE)                                                00001650
  72. CCR   REWIND 31                                                         00001660
  73.       NOWDTE(1) = 7                                                     R0001661
  74.       NOWDTE(2) = 345                                                   R0001662
  75.       WRITE(31,988)NOWDTE(1),NOWDTE(2)                                  00001670
  76.   988 FORMAT(I2,I3)                                                     R0001680
  77.       REWIND 31                                                         00001690
  78.        READ(31,989)N,M                                                  00001700
  79.   989 FORMAT(I2,I3)                                                     R0001710
  80.       ITX1=(N*10000+M*3+1)*2                                            00001720
  81.       IF(NX43.EQ.1700000) NX43=ITX1                                     00001730
  82.       REWIND 31                                                         00001740
  83.       LIST=0                                                            00001750
  84.       LISTC=0                                                           00001760
  85.       LISTB=0                                                           00001770
  86.       LISTA=0                                                           00001780
  87.       NAD=0                                                             00001790
  88.       KSKIP=0                                                           00001800
  89.       XAD=0.5E+00                                                       00001810
  90.                 I1=10000                                                00001820
  91.   100 DO 110 I=1,7                                                      00001830
  92.   110 T(I)=0.0                                                          00001840
  93.       T(1)=-1.0                                                         00001850
  94.       MBAND =1                                                          00001860
  95.       KEQB=99999                                                        00001870
  96.         KZN(1)=2                                                        00001880
  97.         KZN(2)=7                                                        00001890
  98.         KZN(3)=1                                                        00001900
  99.         KZN(4)=2                                                        00001910
  100.         KZN(5)=7                                                        00001920
  101.         KZN(6)=7                                                        00001930
  102.         KZN(7)=10                                                       00001940
  103.         KZN(8)=7                                                        00001950
  104.         KZN(9)=21                                                       00001960
  105.         KZN(10)=1177                                                    00001970
  106.         KZN(11)=6                                                       00001980
  107.         KZN(12)=51                                                      00001990
  108.         KZN(13)=10                                                      00002000
  109.         KZN(14)=21                                                      00002010
  110.         KZN(15)=9                                                       00002020
  111.         KZN(16)=7                                                       00002030
  112.         KZN(17)=8                                                       00002040
  113.         KZN(18)=5                                                       00002050
  114.         KZN(19)=0                                                       00002060
  115.         KZN(20)=0                                                       00002070
  116. CCR   CALL QDBGOF(4HOARC)                                               00002080
  117.       MKZ=20                                                            00002090
  118.         NEXPDT(1)=78                                                    00002100
  119.         NEXPDT(2)=222                                                   00002110
  120.       NCRD=0                                                            00002120
  121.       NCWT=0                                                            00002130
  122.        NSTIF=4                                                          00002140
  123.        NMASS=9                                                          00002150
  124.        NL=62                                                            R0002160
  125.        NR=3                                                             00002170
  126.       NT=10                                                             00002180
  127.       NRED=15                                                           00002190
  128.       M1P=0                                                             00002200
  129.       M2P=0                                                             00002210
  130.         NOSS=0                                                          00002220
  131.       NTYP = 0                                                          R0002221
  132.       CALL SECOND (T(1))                                                00002230
  133.       CALL MODEL(NTYP)                                                  00002240
  134.       CALL COMMRW (0)                                                   R0002241
  135.       WRITE (*,991) NTYP
  136.   991 FORMAT (5X,'*********** SAP6P1 FINISHED ** NTYP ***',I5/)         R0002244
  137.       STOP                                                              R0002245
  138.       END                                                               R0002246
  139.       SUBROUTINE ERROR(I)                                               00086230
  140.       IMPLICIT REAL*8(A-H,O-Z)                                          00086240
  141.       REAL*8  X                                                         00086250
  142.       COMMON /EXTRA/ MODEX,NREXTR(25)                                   R0086260
  143.       COMMON /PREP/ X(2),KSKIP,RRPREP(8)                                R0086270
  144.       KSKIP=1                                                           00086280
  145.       MODEX=1                                                           00086290
  146.       WRITE(6,100)I                                                     00086300
  147.   100 FORMAT (1H0//1X,30HALLOCATED STORAGE EXCEEDED BY   ,I7,6H WORDS)  00086310
  148.       WRITE(6,110)                                                      00086320
  149.   110 FORMAT(/1X, 29HNO EXECUTION WILL BE ALLOWED./)                    00086330
  150.       RETURN                                                            00086340
  151.       END                                                               00086350
  152.       SUBROUTINE CLOSE
  153.       RETURN
  154.       END
  155.       SUBROUTINE EXIT
  156.       WRITE (6,101)
  157.   101 FORMAT (5X,'********  SAP6 PROGRAM STOP  ********')
  158.       STOP
  159.       END
  160.       SUBROUTINE FILES(NN)
  161.       RETURN
  162.       END
  163.       BLOCKDATA                                                         00007790
  164.       IMPLICIT REAL*8(A-H,O-Z)                                          00007800
  165.         COMMON/HEADIN/TITLE1(20),TITLE2(5),TITLE3(10)                   00007810
  166.       COMMON/ELARRY/NELAR(4,20)                                         00007820
  167.       COMMON /GASS/ XK(4,4),WGT(4,4),IPERM(3)                           00007830
  168.       COMMON/GASS2/A5(7,2),W5(7)                                        00007840
  169.       COMMON /PREP/XMX,XAD,KSKIP,NDYN,I1,I99,POS,PRTCOD                 00007850
  170.      1,POSSAV,PRTOFF,PRTON,PRTDUM,IDIRC                                 00007860
  171.       DATA XK /     0.D0,     0.D0,               0.D0,            0.D0,00007870
  172.      $ -.5773502691896D0, .5773502691896D0,     0.D0,              0.D0,00007880
  173.      $ -.7745966692415D0, .0000000000000D0, .7745966692415D0,      0.D0,00007890
  174.      $ -.8611363115941D0,-.3399810435849D0, .3399810435849D0,           00007900
  175.      $.8611363115941D0/                                                 00007910
  176.       DATA WGT / 2.000D0,     0.D0,          0.D0,           0.D0,      00007920
  177.      $ 1.0000000000000D0,1.0000000000000D0,  0.D0,           0.D0,      00007930
  178.      $  .5555555555556D0, .8888888888889D0, .5555555555556D0,0.D0,      00007940
  179.      $  .3478548451375D0, .6521451548625D0, .6521451548625D0,           00007950
  180.      $  .3478548451375D0/                                               00007960
  181.       DATA IPERM / 2,3,1 /                                              00007970
  182.       DATA  A5(1,1)/-0.333333333333D0/,A5(2,1)/-0.88056825640D0/        00007980
  183.       DATA  A5(3,1)/-0.05971587178D0/,A5(4,1)/-0.05971587178D0/         00007990
  184.       DATA  A5(5,1)/ 0.59485397070D0/, A5(6,1)/-0.79742698530D0/        00008000
  185.       DATA  A5(7,1)/-0.79742698530D0/, A5(1,2)/-0.333333333333D0/       00008010
  186.       DATA  A5(2,2)/-0.05971587178D0/, A5(3,2)/-0.88076825640D0/        00008020
  187.       DATA  A5(4,2)/-0.05971587178D0/ ,A5(5,2)/-0.79742698530D0/        00008030
  188.       DATA  A5(6,2)/ 0.59485397070D0/ ,A5(7,2)/-0.79742698530D0/        00008040
  189.       DATA    W5(1)/ 0.225        D0/,  W5(2)/ 0.13239415   D0/         00008050
  190.       DATA    W5(3)/ 0.13239415   D0/,  W5(4)/ 0.13239415   D0/         00008060
  191.       DATA    W5(5)/ 0.12593918   D0/,  W5(6)/ 0.12593918   D0/         00008070
  192.       DATA    W5(7)/ 0.12593918   D0/                                   00008080
  193.       DATA NELAR /                                                      00008090
  194.      $   2,   2,   6,   2,                                              00008100
  195.      $   3,   2,  12,  28,                                              00008110
  196.      $   4,   4,  12,   8,                                              00008120
  197.      $   4,   4,   8,   4,                                              00008130
  198.      $   8,   8,  33,  54,                                              00008140
  199.      $   4,   4,  42,  24,                                              00008150
  200.      $   1,   1,   1,   1,                                              00008160
  201.      $   4,   4,   8,   4,                                              00008170
  202.      $   3,   2,  12,  39,                                              00008180
  203.      $  20,  20,  60,  54,                                              00008190
  204.      $   8,   8,  16,  52,                                              00008200
  205.      $   8,   8,  16,  52,                                              00008210
  206.      $   8,   8,  16,  52,                                              00008220
  207.      $   4,   1,   6,   6,                                              00008230
  208.      $   8,   8,  48,   6,                                              00008240
  209.      $   20*0/                                                          00008250
  210.       DATA TITLE2/4H    ,4HSAP6,4H    ,4HVER.,4H 2.0/                   00008260
  211.       DATA TITLE3(3)/4H LPI/,TITLE3(4)/4HAUTO/,TITLE3(5)/54./           00008330
  212.       DATA POS/3H   /,PRTCOD/3H   /                                     00008340
  213.       DATA POSSAV/3H   /,PRTOFF/3HOFF/,PRTON/3HON-/,PRTDUM/3HDUM/       00008350
  214.       DATA IDIRC/0/                                                     00008360
  215.       END                                                               00008370
  216.       SUBROUTINE SECOND(T)                                              00234270
  217.       IMPLICIT REAL*8 (A-H,O-Z)                                         00234280
  218.       CALL GETTIM(NA,NB,NC,ND)                                          R0234290
  219.       AA = NA * 100.0                                                   R0234291
  220.       CC = NC                                                           R0234292
  221.       CC = CC / 100.0                                                   R0234292
  222.       T  = AA + NB + CC                                                 R0234293
  223.       RETURN                                                            00234300
  224.       END                                                                       
  225.       SUBROUTINE INLIST(NCARDS)                                         00116150
  226.       IMPLICIT REAL*8(A-H,O-Z)                                          00116160
  227.       DIMENSION DATA(20)                                                00116170
  228.       INPUT=5                                                           00116180
  229.       IF(NCARDS.EQ.0) GO TO 20                                          00116190
  230.       DO 10 I=1,NCARDS                                                  00116200
  231.    10 BACKSPACE INPUT                                                   00116210
  232.    20 WRITE(6,100)                                                      00116220
  233.       KARD=0                                                            00116230
  234.       KARDNO=0                                                          00116240
  235.    30 READ (5,110,END=60) DATA                                          00116250
  236.    40 KARDNO=KARDNO+1                                                   00116260
  237.       KARD=KARD+1                                                       00116270
  238.       IF(KARDNO.LE.43) GO TO 50                                         00116280
  239.       WRITE(6,120)                                                      00116290
  240.       WRITE(6,100)                                                      00116300
  241.       KARDNO=1                                                          00116310
  242.    50 WRITE(6,130) KARD,DATA                                            00116320
  243.       GO TO 30                                                          00116330
  244.    60 WRITE(6,120)                                                      00116340
  245.       WRITE(6,140)                                                      00116350
  246.       KARD=KARD-NCARDS+1                                                00116360
  247.       DO 70 I=1,KARD                                                    00116370
  248.    70 BACKSPACE INPUT                                                   00116380
  249.   100 FORMAT(1X ,10X,45HF O L L O W I N G    I S    A    C A R D    I,  00116390
  250.      162HM A G E    L I S T I N G    O F    T H E    I N P U T    D A T,00116400
  251.      22H A//68X,26HC O L U M N    N U M B E R/12X,4HCARD,34X,7H1111111, 00116410
  252.      362H11122222222223333333333444444444455555555556666666666777777777,00116420
  253.      42H78/11X,6HNUMBER,24X,40H1234567890123456789012345678901234567890,00116430
  254.      540H1234567890123456789012345678901234567890/)                     00116440
  255.   110 FORMAT(20A4)                                                      00116450
  256.   120 FORMAT(/12X,4HCARD,34X,39H111111111122222222223333333333444444444,00116460
  257.      132H45555555555666666666677777777778/11X,6HNUMBER,24X,9H123456789, 00116470
  258.      261H0123456789012345678901234567890123456789012345678901234567890, 00116480
  259.      310H1234567890/68X,26HC O L U M N    N U M B E R)                  00116490
  260.   130 FORMAT(I15,26X,20A4)                                              00116500
  261.   140 FORMAT(1H0,10X,34(1H'),37H E N D    O F    I N P U T    L I S T,  00116510
  262.      16H I N G,34(1H'))                                                 00116520
  263.       RETURN                                                            00116530
  264.       END                                                               00116540
  265.       SUBROUTINE HDPRNT(HED,LL,NF,NDYN,KDYN,XMX)                        00111270
  266.       IMPLICIT REAL*8(A-H,O-Z)                                          00111280
  267.       REAL*8  HED,XMX                                                   00111290
  268.       COMMON /SUPEL/ NSELEM,MEQL,NODESE,MATNO,NEADD,NSEL                00111300
  269.         COMMON /PREP/XDM(2),KSKIP,RRPREP(8)                             R0111310
  270.       DIMENSION HED(20)                                                 00111320
  271.       COMMON/DAYTIM/NOWDAY(5),NOWDTE(2)                                 00111330
  272.       WRITE(6,90)                                                       00111340
  273.       KDYN=NDYN                                                         00111350
  274.       NDYN=IABS(NDYN)                                                   00111360
  275.       IF(LL.EQ.0) LL=1                                                  00111370
  276.       IF(NDYN.GT.3.AND.NDYN.LT.8) LL=1                                  00111380
  277.       CALL DATE(NOWDAY)                                                 00111390
  278.       WRITE(6,420)(HED(I),I=1,15),NOWDAY                                00111400
  279.       WRITE(6,220)                                                      00111410
  280.       WRITE(6,260)                                                      00111420
  281.       WRITE(6,270) NOWDAY,NOWDTE,HED                                    00111430
  282.       NDYNMN=0                                                          00111440
  283.       NDYNMX=12                                                         00111450
  284.         IF(NDYN.EQ.10)LL=1                                              00111460
  285.       IF(NDYN.GE.NDYNMN.AND.NDYN.LE.NDYNMX) GO TO 10                    00111470
  286.       WRITE(6,140)                                                      00111480
  287.         KSKIP=1                                                         00111490
  288.       RETURN                                                            00111500
  289.    10 CONTINUE                                                          00111510
  290.       IF(NDYN.GT.0) GO TO 20                                            00111520
  291.       WRITE(6,150)LL                                                    00111530
  292.       GO TO 410                                                         00111540
  293. 20      IF(NDYN.GT.6) GO TO 30                                          00111550
  294.         IF(NF.GT.0)   GO TO 30                                          00111560
  295.         WRITE(6,25)                                                     00111570
  296. 25      FORMAT(///50X,34HNO. OF FREQUENCIES MUST BE NONZERO///)         00111580
  297.         KSKIP=1                                                         00111590
  298.         RETURN                                                          00111600
  299. 30      CONTINUE                                                        00111610
  300.       GO TO (310,320,330,310,320,330,340,350,360,365,1360               00111620
  301.      $,366),NDYN                                                        00111630
  302.   310 WRITE(6,160)                                                      00111640
  303.       GO TO 370                                                         00111650
  304.   320 WRITE(6,170)                                                      00111660
  305.       GO TO 370                                                         00111670
  306.   330 WRITE(6,180)                                                      00111680
  307.       GO TO 370                                                         00111690
  308.   340 WRITE(6,110)                                                      00111700
  309.       GO TO 370                                                         00111710
  310.   350 WRITE(6,230)NF,LL                                                 00111720
  311.       GO TO 370                                                         00111730
  312.   360 WRITE(6,250)LL,NF                                                 00111740
  313.         GO TO 370                                                       00111750
  314.   365 WRITE (6,2600) NF                                                 00111760
  315.       GO TO 370                                                         00111770
  316. 366     WRITE(6,181)                                                    00111780
  317.         GO TO 370                                                       00111790
  318.  1360 CONTINUE                                                          00111800
  319.       IF(NF.EQ.0) NF=1                                                  00111810
  320.       WRITE(6,280)LL,NF                                                 00111820
  321.   370 GO TO (380,380,380,390,390,390,410,400,400,410,410                00111830
  322.      $,390),NDYN                                                        00111840
  323.   380 WRITE(6,190)LL,NF                                                 00111850
  324.       WRITE(6,120)                                                      00111860
  325.       GO TO 410                                                         00111870
  326.   390 WRITE(6,200)NF                                                    00111880
  327.       WRITE(6,130)                                                      00111890
  328.       GO TO 410                                                         00111900
  329.   400 MATNO=IABS(NF)                                                    00111910
  330.       NSEL=0                                                            00111920
  331.       IF(NF.GE.0) GO TO 405                                             00111930
  332.       READ (5,401) NSEL                                                 00111940
  333.   401 FORMAT(I5)                                                        00111950
  334.       WRITE(6,402)MATNO,NSEL                                            00111960
  335.   402 FORMAT(/20X,6HMATRIX,I4,30H WAS ATTACHED TO SUPER ELEMENT,I4/)    00111970
  336.   405 CONTINUE                                                          00111980
  337.       NF=0                                                              00111990
  338.   410 CONTINUE                                                          00112000
  339.   420 FORMAT(/  5H AT**,15A4,7X,'SAP6/AT (87.4) BEIJING ',7X,5A4)    
  340.       RETURN                                                            00112020
  341.  90   FORMAT(1X,5X)                                                     00112030
  342.   110 FORMAT(1H+,51X,41HDYNAMIC ANALYSIS USING DIRECT INTEGRATION)      00112040
  343.   120 FORMAT(52X,51HMODAL EXTRACTION WILL USE THE RAYLIEGH-RITZ METHOD.)00112050
  344.   130 FORMAT(52X,51HMODAL EXTRACTION WILL USE A DETERMINANT SEARCH OR A,00112060
  345.      $29H   SUBSPACE ITERATION METHOD.)                                 00112070
  346.   140 FORMAT(1H+,51X,46HAN ANALYSIS WAS REQUESTED THAT IS NOT POSSIBLE) 00112080
  347.   150 FORMAT(1H+,51X,20HSTATIC ANALYSIS WITH,I3,18H LOAD CONDITION(S))  00112090
  348.   160 FORMAT(1H+,51X,41HDYNAMIC ANALYSIS FOR FREQUENCIES AND MODE,      00112100
  349.      $13H SHAPES   FOR)                                                 00112110
  350.   170 FORMAT(1H+,51X,47HDYNAMIC ANALYSIS INCLUDING TIME-DEPENDENT LOADS,00112120
  351.      $16H IT WILL REQUIRE)                                              00112130
  352.   180 FORMAT(1H+,51X,46HDYNAMIC ANALYSIS INCLUDING A RESPONSE SPECTRUM, 00112140
  353.      $26H ANALYSIS, IT WILL REQUIRE)                                    00112150
  354. 181     FORMAT(1H+,51X,26HDYNAMIC ANALYSIS INCLUDING/52X,9H A STEADY,   00112160
  355.      &  53H STATE SINUSOIDAL RESPONSE ANALYSIS.  IT WILL REQUIRE)       00112170
  356.   190 FORMAT(52X,I2,12H LOAD(S) AND,I3,12H FREQUENCIES)                 00112180
  357.   200 FORMAT(52X,I2,12H FREQUENCIES)                                    00112190
  358.   220 FORMAT(18X,10(1HS),5X,10(1HA),4X,11(1HP),5X,10(1H6),6X,'**',4X,
  359.      $10(1HA),4X,12(1HT),/,
  360.      217X,12(1HS),3X,12(1HA),3X,12(1HP),3X,12(1H6),5X,'**',3X,12(1HA)
  361.      $,3X,12(1HT),/,
  362.      317X,2HSS,8X,2HSS,3X,2HAA,8X,2HAA,3X,2HPP,8X,2HPP,3X,2H66,8X,
  363.      X2H66,4X,'**',4X,2HAA,8X,2HAA,8X,2HTT,/,
  364.      417X,2HSS,10X,3X,2HAA,8X,2HAA,3X,2HPP,8X,2HPP,3X,2H66,13X,'*',6X,
  365.      $2HAA,8X,2HAA,8X,2HTT,/,
  366.      517X,3HSSS,12X,2HAA,8X,2HAA,3X,2HPP,8X,2HPP,3X,2H66,20X,2HAA,8X,
  367.      $2HAA,8X,2HTT,/,
  368.      X18X, 9(1HS),5X,12(1HA),3X,12(1HP),3X,11(1H6),11X,12(1HA),8X,
  369.      72HTT,/,19X,9(1HS),4X,12(1HA),3X,11(1HP),4X,12(1H6),10X,12(1HA),
  370.      88X,2HTT,/,26X,3HSSS,3X,2HAA,8X,2HAA,3X,2HPP,13X,2H66,8X,2H66
  371.      $,10X,2HAA,8X,2HAA,8X,2HTT,/,27X
  372.      9,2HSS,3X,2HAA,8X,2HAA,3X,2HPP,13X,2H66,8X,2H66,10X,2HAA,8X,2HAA,
  373.      $8X,2HTT,/,17X,2HSS,8X,
  374.      A2HSS,3X,2HAA,8X,2HAA,3X,2HPP,13X,2H66,8X,2H66,10X,2HAA,8X,2HAA,8X,
  375.      B2HTT,/,17X,12(1HS),3X,2HAA,8X,2HAA,3X,2HPP,13X,12(1H6),10X,2HAA,8X
  376.      C,2HAA,8X,2HTT,/,18X,10(1HS),4X,2HAA,8X,2HAA,3X,2HPP,14X,10(1H6),
  377.      $11X,2HAA,8X,2HAA,8X,2HTT,////)
  378.   230 FORMAT(1H+,51X,37HFORMATION OF SUPER ELEMENT MATRIX NO.,I4,       00112330
  379.      $5H WITH,I3,7H LOADS.                                              00112340
  380.      $ /52X,55HTHE MATRICES ARE WRITTEN ON TAPE20 AND SHOULD BE STORED, 00112350
  381.      $15H FOR LATER USE.)                                               00112360
  382.   250 FORMAT(1H+,51X,11HRECOVERY OF,I3,29H SETS OF DISPLACEMENTS OR MOD,00112370
  383.      $26HE SHAPES FOR SUPER ELEMENT/52X,6HMATRIX,I4,11H. THE SUPER,     00112380
  384.      $37H ELEMENT MATRICES ARE TO BE ON TAPE20/52X,                     00112390
  385.      $42HAND THE DISPLACEMENTS ARE TO BE ON TAPE27.)                    00112400
  386.   260 FORMAT(46X,41(1H*),/,46X,1H*,39X,1H*,/,46X,1H*,39X,1H*,/,46X,1H*, 00112410
  387.      17X,24HAT VERSION: APRIL , 1987 ,8X,1H*,/,46X,1H*,39X,1H*/46X,     00112420
  388.      21H*,4X,'INSTITUTE OF COMPUTER TECHNOLOGY',3X,1H*/46X,1H*,39X,1H*,/
  389.      346X,1H*,4X,'STATE MACHINE INDUSTRY COMMISION',3X,1H*/46X,1H*,39X,1
  390.      4H*/46X,1H*,10X,'DEPARTMENT OF CAD',12X,1H*/46X,1H*,39X,1H*/
  391.      546X,1H*,4X,'ADDR: 10 CHEDAOGO BEIJIBG',10X,1H*/46X,1H*,39X,1H*/
  392.      646X,1H*,4X,'P.O.BOX 2418-26 BEIJING',12X,1H*/46X,1H*
  393.      7,39X,1H*/46X,1H*,8X,23HTEL. NO. 890581-632    ,8X,1H*/46X,1H*,39X
  394.      8,1H*/46X,1H*,39X,1H*/46X,41(1H*))                                 00112490
  395.   270  FORMAT(//,33X,18HDATE             - ,5A4,5X,2A4,6H HOURS,//      00112500
  396.      133X,18HTITLE OF RUN     - ,1X,20A4//                              00112510
  397.      233X,18HTYPE OF ANALYSIS -)                                        00112520
  398.   280 FORMAT(1H+,51X,17HBUCKLING ANALYSIS,/,52X,                        00112530
  399.      119HSTATIC ANALYSIS FOR,I3,32H LOAD CONDITIONS IS FIRST SOLVED,/,  00112540
  400.      252X,26HTHEN BUCKLING ANALYSIS FOR,I3,25H MODE SHAPES IS PERFORMED)00112550
  401.  2600 FORMAT(1H+,51X,12HRESULTS FOR ,I5,30H TIME STEPS WILL BE RECOVERED00112560
  402.      $.)                                                                00112570
  403.       END                                                               00112580
  404.       SUBROUTINE ELSZ(N)                                                00084970
  405.       IMPLICIT REAL*8(A-H,O-Z)                                          00084980
  406.       REAL*8  X                                                         00084990
  407.       COMMON /PREP/ X(2),KSKIP,RRPREP(8)                                R0085000
  408.       COMMON/SIZE/NDMX,MXDF,NSMX,NTERM,NADND                            00085010
  409.       COMMON/ELARRY/NELAR(4,20)                                         00085020
  410.       MAXEL=15                                                          00085030
  411.       IF(N.GT.0.AND.N.LE.MAXEL) GO TO 120                               00085040
  412.   100 WRITE(6,110)N                                                     00085050
  413.   110 FORMAT (20X, 12HELEMENT TYPE,I10, 15H DOES NOT EXIST/)            00085060
  414.       KSKIP=1                                                           00085070
  415.       RETURN                                                            00085080
  416.   120 IF(NELAR(1,N).EQ.0   ) GO TO 100                                  00085090
  417.       IF(NELAR(1,N).GT.NDMX) NDMX=NELAR(1,N)                            00085100
  418.       IF(NELAR(3,N).GT.MXDF) MXDF=NELAR(3,N)                            00085110
  419.       IF(NELAR(4,N).GT.NSMX) NSMX=NELAR(4,N)                            00085120
  420.       RETURN                                                            00085130
  421.       END                                                               00085140
  422.