      PROGRAM SAP6P1                                                    R0001101
      IMPLICIT REAL*8(A-H,O-Z)                                          00001100
      LOGICAL ELPRT,ELPCH,GENPRT,GENPCH,DEFPCH,GEOST                    00001110
      COMMON/BMDATA/IAISC,ILDTYP,ILDMUL,IBMSEC,NLWAVE,ILOCAL            00001120
      COMMON/EQUILB/NEQIL,NX43                                          00001130
      COMMON/CTL/ELPRT,ELPCH,GENPRT,GENPCH                              00001140
      COMMON / JUNK / DUK(227)                                          00001150
      REAL*8  NPAR                                                      00001160
      COMMON /QTSARG/ QQQ(1000)                                         00001170
      COMMON/DYN3/ NEIG,NAD,ANORM,NVV,NFO                               00001180
      COMMON /DYN4/KDYN,NRESS,NCRD,NCWT,NRESS1                          00001190
      COMMON /DYN5/FRSHFT,FRINIT,FREND,MODEFR,NOSS                      00001200
       COMMON /TAPES/NSTIF,NRED,NL,NR,NT,NMASS                          00001210
      COMMON /EXTRA/MODEX,NT8,N10SV,NT10,KEQB,NY,T(10)                  00001220
      COMMON/GEOSTF/GEOST,NELGEO                                        00001230
      COMMON/MASS/LMASS                                                 00001240
      COMMON/MATL/MATLCO                                                R0001241
      COMMON/SIZE/NDMX,MXDF,NSMX,NTERM,NADND                            00001250
      COMMON/SLVE/NSLAVE                                                00001260
      COMMON /ELPAR/ NPAR(14),NUMNP,MBAND,NELTYP,N1,N2,N3,N4,N5,MTOT,NEQ00001270
     $ ,N2P,N3P,NMRI,NTRI,N1P,NUMEL,NUMEL2,KZ(20,2),NEMN                00001280
      COMMON / MISC / NBLOCK,NEQB,LL,NFREQ,LB                           00001290
      COMMON/AMB/ GRAV,REFT,JROT                                        00001300
      COMMON /PREP/ XMX,XAD,KSKIP,NDYN,I1,I99,POS,PRTCOD                00001310
     $              ,POSSAV,PRTOFF,PRTON,PRTDUM,IDIRC                   00001320
      COMMON /DYN/ IFIL1(11),IFIL2                                      00001330
      COMMON/ELARRY/NELAR(4,20)                                         00001340
      COMMON /ELTEMP/ SET1(103)                                         00001350
     $       /OUT/KSET2(6),KELRST,MAXDF,IFIL3(2)                        00001360
     $       /SQZ/    SET3,LIST,LISTC,LISTB,LISTA                       00001370
     $       /TRASH/  SET4(490)                                         00001380
     $       /GPS/    SET5(10)                                          00001390
     $       /CG/     SET6(4),RFIL1(2)                                  00001400
     $       /TAPES/ SET7(6)                                            00001410
     $       /DYN2/KSET8(3),NFVC,SET8(12)                               00001420
      COMMON /WORDS/ NWDS(30,2)                                         00001430
      COMMON /BAND/  NRNM(3),IRSK,IFIL4(4)                              00001440
      COMMON /HEADIN/TITLE1(20),TITLE2(5),TITLE3(10)                    00001450
      COMMON /SUPEL/ NSELEM,NEQL,NODESE                                 00001460
     $,KSET(3)                                                          00001470
      COMMON/FORCE/ NLC,NELD                                            00001480
        COMMON/DAPG/NQ1,NQX,DUMY(210)                                   00001490
      COMMON/VAXPLT/IVPLT,XVB(6),XVA(4),MINX,MINY,LEN                   00001500
      COMMON /ICM/ICOMP,MMRI,MTRI,M1P,M2P,M3P                           00001510
      COMMON/PLOTH/IPLT,IPLWRT                                          00001520
      COMMON/COMMT/NTYP,NUMET,NNRRC,NNRRC1                              R0001521
      COMMON/COMMT1/NDSSS,KDSSS,NTY,NSLDM,NBLANK,MDYN,NE2B,KRK1,MCB,MLT R0001522
      DIMENSION KZN(20),ZD(31)                                          00001530
      DIMENSION NEXPDT(2),NOWDTE(4)                                     00001540
      COMMON A(1)                                                       00001550
      DATA KZN/2,7,1,2,7,7,10,7,21,1177,6,51,10,21,9,7,8,5,0,0/         00001560
      DATA NEXPDT/78,222/                                               00001570
      DATA CONKEY/8HSUBSTRUC/                                           R0001571
      WRITE (*,990)                                                     R0001611
  990 FORMAT (5X,'*********** PROGRAM SAP6P1 STARTING *********')       R0001612
CC    CALL OPEN1                                                        00001580
CCR    CALL TIME(NOWDTE)                                                00001590
CCR   CALL TODATE                                                       00001600
      READ (5,2001) COKEY1                                              R0001601
      WRITE (6,2001) COKEY1                                             R0001602
 2001 FORMAT (A8)                                                       R0001603
 2002 FORMAT (5X,'***** SUBSTRUCT STARTING *****'/)                     R0001604
      IF (CONKEY .EQ. COKEY1) CALL COMMRW (1)                           R0002241
      IF (CONKEY .NE. COKEY1) REWIND 5                                  R0002242
      IF (CONKEY .EQ. COKEY1) WRITE (6,2002)                            R0002244
      IF (CONKEY .EQ. COKEY1) GO TO 100                                 R0002243
      CALL SIZER1                                                       00001610
CCR    CALL ERRSET(208,0,-1,1,1)                                        00001620
CCR   CALL CPUINT                                                       00001630
      NX43=1700000                                                      00001640
CCR   CALL FDATE(NOWDTE)                                                00001650
CCR   REWIND 31                                                         00001660
      NOWDTE(1) = 7                                                     R0001661
      NOWDTE(2) = 345                                                   R0001662
      WRITE(31,988)NOWDTE(1),NOWDTE(2)                                  00001670
  988 FORMAT(I2,I3)                                                     R0001680
      REWIND 31                                                         00001690
       READ(31,989)N,M                                                  00001700
  989 FORMAT(I2,I3)                                                     R0001710
      ITX1=(N*10000+M*3+1)*2                                            00001720
      IF(NX43.EQ.1700000) NX43=ITX1                                     00001730
      REWIND 31                                                         00001740
      LIST=0                                                            00001750
      LISTC=0                                                           00001760
      LISTB=0                                                           00001770
      LISTA=0                                                           00001780
      NAD=0                                                             00001790
      KSKIP=0                                                           00001800
      XAD=0.5E+00                                                       00001810
                I1=10000                                                00001820
  100 DO 110 I=1,7                                                      00001830
  110 T(I)=0.0                                                          00001840
      T(1)=-1.0                                                         00001850
      MBAND =1                                                          00001860
      KEQB=99999                                                        00001870
        KZN(1)=2                                                        00001880
        KZN(2)=7                                                        00001890
        KZN(3)=1                                                        00001900
        KZN(4)=2                                                        00001910
        KZN(5)=7                                                        00001920
        KZN(6)=7                                                        00001930
        KZN(7)=10                                                       00001940
        KZN(8)=7                                                        00001950
        KZN(9)=21                                                       00001960
        KZN(10)=1177                                                    00001970
        KZN(11)=6                                                       00001980
        KZN(12)=51                                                      00001990
        KZN(13)=10                                                      00002000
        KZN(14)=21                                                      00002010
        KZN(15)=9                                                       00002020
        KZN(16)=7                                                       00002030
        KZN(17)=8                                                       00002040
        KZN(18)=5                                                       00002050
        KZN(19)=0                                                       00002060
        KZN(20)=0                                                       00002070
CCR   CALL QDBGOF(4HOARC)                                               00002080
      MKZ=20                                                            00002090
        NEXPDT(1)=78                                                    00002100
        NEXPDT(2)=222                                                   00002110
      NCRD=0                                                            00002120
      NCWT=0                                                            00002130
       NSTIF=4                                                          00002140
       NMASS=9                                                          00002150
       NL=62                                                            R0002160
       NR=3                                                             00002170
      NT=10                                                             00002180
      NRED=15                                                           00002190
      M1P=0                                                             00002200
      M2P=0                                                             00002210
        NOSS=0                                                          00002220
      NTYP = 0                                                          R0002221
      CALL SECOND (T(1))                                                00002230
      CALL MODEL(NTYP)                                                  00002240
      CALL COMMRW (0)                                                   R0002241
      WRITE (*,991) NTYP
  991 FORMAT (5X,'*********** SAP6P1 FINISHED ** NTYP ***',I5/)         R0002244
      STOP                                                              R0002245
      END                                                               R0002246
      SUBROUTINE ERROR(I)                                               00086230
      IMPLICIT REAL*8(A-H,O-Z)                                          00086240
      REAL*8  X                                                         00086250
      COMMON /EXTRA/ MODEX,NREXTR(25)                                   R0086260
      COMMON /PREP/ X(2),KSKIP,RRPREP(8)                                R0086270
      KSKIP=1                                                           00086280
      MODEX=1                                                           00086290
      WRITE(6,100)I                                                     00086300
  100 FORMAT (1H0//1X,30HALLOCATED STORAGE EXCEEDED BY   ,I7,6H WORDS)  00086310
      WRITE(6,110)                                                      00086320
  110 FORMAT(/1X, 29HNO EXECUTION WILL BE ALLOWED./)                    00086330
      RETURN                                                            00086340
      END                                                               00086350
      SUBROUTINE CLOSE
      RETURN
      END
      SUBROUTINE EXIT
      WRITE (6,101)
  101 FORMAT (5X,'********  SAP6 PROGRAM STOP  ********')
      STOP
      END
      SUBROUTINE FILES(NN)
      RETURN
      END
      BLOCKDATA                                                         00007790
      IMPLICIT REAL*8(A-H,O-Z)                                          00007800
        COMMON/HEADIN/TITLE1(20),TITLE2(5),TITLE3(10)                   00007810
      COMMON/ELARRY/NELAR(4,20)                                         00007820
      COMMON /GASS/ XK(4,4),WGT(4,4),IPERM(3)                           00007830
      COMMON/GASS2/A5(7,2),W5(7)                                        00007840
      COMMON /PREP/XMX,XAD,KSKIP,NDYN,I1,I99,POS,PRTCOD                 00007850
     1,POSSAV,PRTOFF,PRTON,PRTDUM,IDIRC                                 00007860
      DATA XK /     0.D0,     0.D0,               0.D0,            0.D0,00007870
     $ -.5773502691896D0, .5773502691896D0,     0.D0,              0.D0,00007880
     $ -.7745966692415D0, .0000000000000D0, .7745966692415D0,      0.D0,00007890
     $ -.8611363115941D0,-.3399810435849D0, .3399810435849D0,           00007900
     $.8611363115941D0/                                                 00007910
      DATA WGT / 2.000D0,     0.D0,          0.D0,           0.D0,      00007920
     $ 1.0000000000000D0,1.0000000000000D0,  0.D0,           0.D0,      00007930
     $  .5555555555556D0, .8888888888889D0, .5555555555556D0,0.D0,      00007940
     $  .3478548451375D0, .6521451548625D0, .6521451548625D0,           00007950
     $  .3478548451375D0/                                               00007960
      DATA IPERM / 2,3,1 /                                              00007970
      DATA  A5(1,1)/-0.333333333333D0/,A5(2,1)/-0.88056825640D0/        00007980
      DATA  A5(3,1)/-0.05971587178D0/,A5(4,1)/-0.05971587178D0/         00007990
      DATA  A5(5,1)/ 0.59485397070D0/, A5(6,1)/-0.79742698530D0/        00008000
      DATA  A5(7,1)/-0.79742698530D0/, A5(1,2)/-0.333333333333D0/       00008010
      DATA  A5(2,2)/-0.05971587178D0/, A5(3,2)/-0.88076825640D0/        00008020
      DATA  A5(4,2)/-0.05971587178D0/ ,A5(5,2)/-0.79742698530D0/        00008030
      DATA  A5(6,2)/ 0.59485397070D0/ ,A5(7,2)/-0.79742698530D0/        00008040
      DATA    W5(1)/ 0.225        D0/,  W5(2)/ 0.13239415   D0/         00008050
      DATA    W5(3)/ 0.13239415   D0/,  W5(4)/ 0.13239415   D0/         00008060
      DATA    W5(5)/ 0.12593918   D0/,  W5(6)/ 0.12593918   D0/         00008070
      DATA    W5(7)/ 0.12593918   D0/                                   00008080
      DATA NELAR /                                                      00008090
     $   2,   2,   6,   2,                                              00008100
     $   3,   2,  12,  28,                                              00008110
     $   4,   4,  12,   8,                                              00008120
     $   4,   4,   8,   4,                                              00008130
     $   8,   8,  33,  54,                                              00008140
     $   4,   4,  42,  24,                                              00008150
     $   1,   1,   1,   1,                                              00008160
     $   4,   4,   8,   4,                                              00008170
     $   3,   2,  12,  39,                                              00008180
     $  20,  20,  60,  54,                                              00008190
     $   8,   8,  16,  52,                                              00008200
     $   8,   8,  16,  52,                                              00008210
     $   8,   8,  16,  52,                                              00008220
     $   4,   1,   6,   6,                                              00008230
     $   8,   8,  48,   6,                                              00008240
     $   20*0/                                                          00008250
      DATA TITLE2/4H    ,4HSAP6,4H    ,4HVER.,4H 2.0/                   00008260
      DATA TITLE3(3)/4H LPI/,TITLE3(4)/4HAUTO/,TITLE3(5)/54./           00008330
      DATA POS/3H   /,PRTCOD/3H   /                                     00008340
      DATA POSSAV/3H   /,PRTOFF/3HOFF/,PRTON/3HON-/,PRTDUM/3HDUM/       00008350
      DATA IDIRC/0/                                                     00008360
      END                                                               00008370
      SUBROUTINE SECOND(T)                                              00234270
      IMPLICIT REAL*8 (A-H,O-Z)                                         00234280
      CALL GETTIM(NA,NB,NC,ND)                                          R0234290
      AA = NA * 100.0                                                   R0234291
      CC = NC                                                           R0234292
      CC = CC / 100.0                                                   R0234292
      T  = AA + NB + CC                                                 R0234293
      RETURN                                                            00234300
      END                                                                       
      SUBROUTINE INLIST(NCARDS)                                         00116150
      IMPLICIT REAL*8(A-H,O-Z)                                          00116160
      DIMENSION DATA(20)                                                00116170
      INPUT=5                                                           00116180
      IF(NCARDS.EQ.0) GO TO 20                                          00116190
      DO 10 I=1,NCARDS                                                  00116200
   10 BACKSPACE INPUT                                                   00116210
   20 WRITE(6,100)                                                      00116220
      KARD=0                                                            00116230
      KARDNO=0                                                          00116240
   30 READ (5,110,END=60) DATA                                          00116250
   40 KARDNO=KARDNO+1                                                   00116260
      KARD=KARD+1                                                       00116270
      IF(KARDNO.LE.43) GO TO 50                                         00116280
      WRITE(6,120)                                                      00116290
      WRITE(6,100)                                                      00116300
      KARDNO=1                                                          00116310
   50 WRITE(6,130) KARD,DATA                                            00116320
      GO TO 30                                                          00116330
   60 WRITE(6,120)                                                      00116340
      WRITE(6,140)                                                      00116350
      KARD=KARD-NCARDS+1                                                00116360
      DO 70 I=1,KARD                                                    00116370
   70 BACKSPACE INPUT                                                   00116380
  100 FORMAT(1X ,10X,45HF O L L O W I N G    I S    A    C A R D    I,  00116390
     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
     22H A//68X,26HC O L U M N    N U M B E R/12X,4HCARD,34X,7H1111111, 00116410
     362H11122222222223333333333444444444455555555556666666666777777777,00116420
     42H78/11X,6HNUMBER,24X,40H1234567890123456789012345678901234567890,00116430
     540H1234567890123456789012345678901234567890/)                     00116440
  110 FORMAT(20A4)                                                      00116450
  120 FORMAT(/12X,4HCARD,34X,39H111111111122222222223333333333444444444,00116460
     132H45555555555666666666677777777778/11X,6HNUMBER,24X,9H123456789, 00116470
     261H0123456789012345678901234567890123456789012345678901234567890, 00116480
     310H1234567890/68X,26HC O L U M N    N U M B E R)                  00116490
  130 FORMAT(I15,26X,20A4)                                              00116500
  140 FORMAT(1H0,10X,34(1H'),37H E N D    O F    I N P U T    L I S T,  00116510
     16H I N G,34(1H'))                                                 00116520
      RETURN                                                            00116530
      END                                                               00116540
      SUBROUTINE HDPRNT(HED,LL,NF,NDYN,KDYN,XMX)                        00111270
      IMPLICIT REAL*8(A-H,O-Z)                                          00111280
      REAL*8  HED,XMX                                                   00111290
      COMMON /SUPEL/ NSELEM,MEQL,NODESE,MATNO,NEADD,NSEL                00111300
        COMMON /PREP/XDM(2),KSKIP,RRPREP(8)                             R0111310
      DIMENSION HED(20)                                                 00111320
      COMMON/DAYTIM/NOWDAY(5),NOWDTE(2)                                 00111330
      WRITE(6,90)                                                       00111340
      KDYN=NDYN                                                         00111350
      NDYN=IABS(NDYN)                                                   00111360
      IF(LL.EQ.0) LL=1                                                  00111370
      IF(NDYN.GT.3.AND.NDYN.LT.8) LL=1                                  00111380
      CALL DATE(NOWDAY)                                                 00111390
      WRITE(6,420)(HED(I),I=1,15),NOWDAY                                00111400
      WRITE(6,220)                                                      00111410
      WRITE(6,260)                                                      00111420
      WRITE(6,270) NOWDAY,NOWDTE,HED                                    00111430
      NDYNMN=0                                                          00111440
      NDYNMX=12                                                         00111450
        IF(NDYN.EQ.10)LL=1                                              00111460
      IF(NDYN.GE.NDYNMN.AND.NDYN.LE.NDYNMX) GO TO 10                    00111470
      WRITE(6,140)                                                      00111480
        KSKIP=1                                                         00111490
      RETURN                                                            00111500
   10 CONTINUE                                                          00111510
      IF(NDYN.GT.0) GO TO 20                                            00111520
      WRITE(6,150)LL                                                    00111530
      GO TO 410                                                         00111540
20      IF(NDYN.GT.6) GO TO 30                                          00111550
        IF(NF.GT.0)   GO TO 30                                          00111560
        WRITE(6,25)                                                     00111570
25      FORMAT(///50X,34HNO. OF FREQUENCIES MUST BE NONZERO///)         00111580
        KSKIP=1                                                         00111590
        RETURN                                                          00111600
30      CONTINUE                                                        00111610
      GO TO (310,320,330,310,320,330,340,350,360,365,1360               00111620
     $,366),NDYN                                                        00111630
  310 WRITE(6,160)                                                      00111640
      GO TO 370                                                         00111650
  320 WRITE(6,170)                                                      00111660
      GO TO 370                                                         00111670
  330 WRITE(6,180)                                                      00111680
      GO TO 370                                                         00111690
  340 WRITE(6,110)                                                      00111700
      GO TO 370                                                         00111710
  350 WRITE(6,230)NF,LL                                                 00111720
      GO TO 370                                                         00111730
  360 WRITE(6,250)LL,NF                                                 00111740
        GO TO 370                                                       00111750
  365 WRITE (6,2600) NF                                                 00111760
      GO TO 370                                                         00111770
366     WRITE(6,181)                                                    00111780
        GO TO 370                                                       00111790
 1360 CONTINUE                                                          00111800
      IF(NF.EQ.0) NF=1                                                  00111810
      WRITE(6,280)LL,NF                                                 00111820
  370 GO TO (380,380,380,390,390,390,410,400,400,410,410                00111830
     $,390),NDYN                                                        00111840
  380 WRITE(6,190)LL,NF                                                 00111850
      WRITE(6,120)                                                      00111860
      GO TO 410                                                         00111870
  390 WRITE(6,200)NF                                                    00111880
      WRITE(6,130)                                                      00111890
      GO TO 410                                                         00111900
  400 MATNO=IABS(NF)                                                    00111910
      NSEL=0                                                            00111920
      IF(NF.GE.0) GO TO 405                                             00111930
      READ (5,401) NSEL                                                 00111940
  401 FORMAT(I5)                                                        00111950
      WRITE(6,402)MATNO,NSEL                                            00111960
  402 FORMAT(/20X,6HMATRIX,I4,30H WAS ATTACHED TO SUPER ELEMENT,I4/)    00111970
  405 CONTINUE                                                          00111980
      NF=0                                                              00111990
  410 CONTINUE                                                          00112000
  420 FORMAT(/  5H AT**,15A4,7X,'SAP6/AT (87.4) BEIJING ',7X,5A4)    
      RETURN                                                            00112020
 90   FORMAT(1X,5X)                                                     00112030
  110 FORMAT(1H+,51X,41HDYNAMIC ANALYSIS USING DIRECT INTEGRATION)      00112040
  120 FORMAT(52X,51HMODAL EXTRACTION WILL USE THE RAYLIEGH-RITZ METHOD.)00112050
  130 FORMAT(52X,51HMODAL EXTRACTION WILL USE A DETERMINANT SEARCH OR A,00112060
     $29H   SUBSPACE ITERATION METHOD.)                                 00112070
  140 FORMAT(1H+,51X,46HAN ANALYSIS WAS REQUESTED THAT IS NOT POSSIBLE) 00112080
  150 FORMAT(1H+,51X,20HSTATIC ANALYSIS WITH,I3,18H LOAD CONDITION(S))  00112090
  160 FORMAT(1H+,51X,41HDYNAMIC ANALYSIS FOR FREQUENCIES AND MODE,      00112100
     $13H SHAPES   FOR)                                                 00112110
  170 FORMAT(1H+,51X,47HDYNAMIC ANALYSIS INCLUDING TIME-DEPENDENT LOADS,00112120
     $16H IT WILL REQUIRE)                                              00112130
  180 FORMAT(1H+,51X,46HDYNAMIC ANALYSIS INCLUDING A RESPONSE SPECTRUM, 00112140
     $26H ANALYSIS, IT WILL REQUIRE)                                    00112150
181     FORMAT(1H+,51X,26HDYNAMIC ANALYSIS INCLUDING/52X,9H A STEADY,   00112160
     &  53H STATE SINUSOIDAL RESPONSE ANALYSIS.  IT WILL REQUIRE)       00112170
  190 FORMAT(52X,I2,12H LOAD(S) AND,I3,12H FREQUENCIES)                 00112180
  200 FORMAT(52X,I2,12H FREQUENCIES)                                    00112190
  220 FORMAT(18X,10(1HS),5X,10(1HA),4X,11(1HP),5X,10(1H6),6X,'**',4X,
     $10(1HA),4X,12(1HT),/,
     217X,12(1HS),3X,12(1HA),3X,12(1HP),3X,12(1H6),5X,'**',3X,12(1HA)
     $,3X,12(1HT),/,
     317X,2HSS,8X,2HSS,3X,2HAA,8X,2HAA,3X,2HPP,8X,2HPP,3X,2H66,8X,
     X2H66,4X,'**',4X,2HAA,8X,2HAA,8X,2HTT,/,
     417X,2HSS,10X,3X,2HAA,8X,2HAA,3X,2HPP,8X,2HPP,3X,2H66,13X,'*',6X,
     $2HAA,8X,2HAA,8X,2HTT,/,
     517X,3HSSS,12X,2HAA,8X,2HAA,3X,2HPP,8X,2HPP,3X,2H66,20X,2HAA,8X,
     $2HAA,8X,2HTT,/,
     X18X, 9(1HS),5X,12(1HA),3X,12(1HP),3X,11(1H6),11X,12(1HA),8X,
     72HTT,/,19X,9(1HS),4X,12(1HA),3X,11(1HP),4X,12(1H6),10X,12(1HA),
     88X,2HTT,/,26X,3HSSS,3X,2HAA,8X,2HAA,3X,2HPP,13X,2H66,8X,2H66
     $,10X,2HAA,8X,2HAA,8X,2HTT,/,27X
     9,2HSS,3X,2HAA,8X,2HAA,3X,2HPP,13X,2H66,8X,2H66,10X,2HAA,8X,2HAA,
     $8X,2HTT,/,17X,2HSS,8X,
     A2HSS,3X,2HAA,8X,2HAA,3X,2HPP,13X,2H66,8X,2H66,10X,2HAA,8X,2HAA,8X,
     B2HTT,/,17X,12(1HS),3X,2HAA,8X,2HAA,3X,2HPP,13X,12(1H6),10X,2HAA,8X
     C,2HAA,8X,2HTT,/,18X,10(1HS),4X,2HAA,8X,2HAA,3X,2HPP,14X,10(1H6),
     $11X,2HAA,8X,2HAA,8X,2HTT,////)
  230 FORMAT(1H+,51X,37HFORMATION OF SUPER ELEMENT MATRIX NO.,I4,       00112330
     $5H WITH,I3,7H LOADS.                                              00112340
     $ /52X,55HTHE MATRICES ARE WRITTEN ON TAPE20 AND SHOULD BE STORED, 00112350
     $15H FOR LATER USE.)                                               00112360
  250 FORMAT(1H+,51X,11HRECOVERY OF,I3,29H SETS OF DISPLACEMENTS OR MOD,00112370
     $26HE SHAPES FOR SUPER ELEMENT/52X,6HMATRIX,I4,11H. THE SUPER,     00112380
     $37H ELEMENT MATRICES ARE TO BE ON TAPE20/52X,                     00112390
     $42HAND THE DISPLACEMENTS ARE TO BE ON TAPE27.)                    00112400
  260 FORMAT(46X,41(1H*),/,46X,1H*,39X,1H*,/,46X,1H*,39X,1H*,/,46X,1H*, 00112410
     17X,24HAT VERSION: APRIL , 1987 ,8X,1H*,/,46X,1H*,39X,1H*/46X,     00112420
     21H*,4X,'INSTITUTE OF COMPUTER TECHNOLOGY',3X,1H*/46X,1H*,39X,1H*,/
     346X,1H*,4X,'STATE MACHINE INDUSTRY COMMISION',3X,1H*/46X,1H*,39X,1
     4H*/46X,1H*,10X,'DEPARTMENT OF CAD',12X,1H*/46X,1H*,39X,1H*/
     546X,1H*,4X,'ADDR: 10 CHEDAOGO BEIJIBG',10X,1H*/46X,1H*,39X,1H*/
     646X,1H*,4X,'P.O.BOX 2418-26 BEIJING',12X,1H*/46X,1H*
     7,39X,1H*/46X,1H*,8X,23HTEL. NO. 890581-632    ,8X,1H*/46X,1H*,39X
     8,1H*/46X,1H*,39X,1H*/46X,41(1H*))                                 00112490
  270  FORMAT(//,33X,18HDATE             - ,5A4,5X,2A4,6H HOURS,//      00112500
     133X,18HTITLE OF RUN     - ,1X,20A4//                              00112510
     233X,18HTYPE OF ANALYSIS -)                                        00112520
  280 FORMAT(1H+,51X,17HBUCKLING ANALYSIS,/,52X,                        00112530
     119HSTATIC ANALYSIS FOR,I3,32H LOAD CONDITIONS IS FIRST SOLVED,/,  00112540
     252X,26HTHEN BUCKLING ANALYSIS FOR,I3,25H MODE SHAPES IS PERFORMED)00112550
 2600 FORMAT(1H+,51X,12HRESULTS FOR ,I5,30H TIME STEPS WILL BE RECOVERED00112560
     $.)                                                                00112570
      END                                                               00112580
      SUBROUTINE ELSZ(N)                                                00084970
      IMPLICIT REAL*8(A-H,O-Z)                                          00084980
      REAL*8  X                                                         00084990
      COMMON /PREP/ X(2),KSKIP,RRPREP(8)                                R0085000
      COMMON/SIZE/NDMX,MXDF,NSMX,NTERM,NADND                            00085010
      COMMON/ELARRY/NELAR(4,20)                                         00085020
      MAXEL=15                                                          00085030
      IF(N.GT.0.AND.N.LE.MAXEL) GO TO 120                               00085040
  100 WRITE(6,110)N                                                     00085050
  110 FORMAT (20X, 12HELEMENT TYPE,I10, 15H DOES NOT EXIST/)            00085060
      KSKIP=1                                                           00085070
      RETURN                                                            00085080
  120 IF(NELAR(1,N).EQ.0   ) GO TO 100                                  00085090
      IF(NELAR(1,N).GT.NDMX) NDMX=NELAR(1,N)                            00085100
      IF(NELAR(3,N).GT.MXDF) MXDF=NELAR(3,N)                            00085110
      IF(NELAR(4,N).GT.NSMX) NSMX=NELAR(4,N)                            00085120
      RETURN                                                            00085130
      END                                                               00085140
