      SUBROUTINE PASS                                                   00161290
      IMPLICIT REAL*8(A-H,O-Z)                                          00161300
      REAL*8 NPAR                                                       00161310
      COMMON /MISC/ NBLOCK,NEQB,LL,NF,LB                                00161320
      COMMON /JUNK/ DUK(200),KKK(4),NDYN,NRJUNK(49)                     R0161330
      COMMON /BAND/ KOPT,NRBAND(7)                                      R0161340
      DIMENSION T(3)                                                    00161350
      COMMON /ELPAR/ NPAR(14),NUMNP,MBAND,NELTYP,N1,N2,N3,N4,N5,MTOT,NEQ00161360
     & ,RRELPA(24)                                                      R0161361
      COMMON /OUT/ IDUMM(4),IOSIG,IODISP,NROUT(4)                       R0161370
      COMMON /HEADIN/ TITLE1(8),TITLE2(2),TITLE3(4),RRHEAD(21)          R0161380
      COMMON A(1)                                                       00161390
      COMMON /TRASH/ TIM(400),RRTRAS(90)                                R0161400
      COMMON /PREP/ XD(2),KSKIP,RRPREP(8)                               R0161410
        LL=NF                                                           00161420
        NEQB=NEQ                                                        00161430
        LB=1                                                            00161440
100     FORMAT(8F10.0)                                                  00161450
      WRITE(6, 110) LL                                                  00161460
110   FORMAT (1X ,20X,38HTHERE WILL BE OUTPUT AT THE FOLLOWING ,I5,7H TI00161470
     1MES.//)                                                           00161480
      IF (LL.LE.400) GO TO 130                                          00161490
      WRITE(6, 120)                                                     00161500
120   FORMAT (1H0,20X,34HTOO MANY TIMES HAVE BEEN SELECTED.//)          00161510
      KSKIP=1                                                           00161520
130   READ(5, 100) (TIM(I),I=1,LL)                                      00161530
      TL=-1.0                                                           00161540
      DO 140 I=1,LL                                                     00161550
         IF (TIM(I).LT.0.0) GO TO 150                                   00161560
         IF (TIM(I).LE.TL) GO TO 150                                    00161570
140   TL=TIM(I)                                                         00161580
      GO TO 170                                                         00161590
150   WRITE(6, 160)                                                     00161600
160   FORMAT (/20X,49HERROR - TIMES MUST BE POS. AND IN ASCENDING ORDER/00161610
     1)                                                                 00161620
      KSKIP=1                                                           00161630
170   WRITE(6, 180) (TIM(I),I=1,LL)                                     00161640
180   FORMAT (20X,6G15.6)                                               00161650
        IF(KSKIP.EQ.1) RETURN                                           00161660
      NT14=14                                                           00161670
        ND=1                                                            00161680
        NV=N1+NEQ                                                       00161690
      REWIND NT14                                                       00161700
      NREC=2*((NEQ-1)/5+1)                                              00161710
      NT=2                                                              00161720
      REWIND NT                                                         00161730
      LY=0                                                              00161740
      TH=-1.0                                                           00161750
      TL=-10.0                                                          00161760
      TIME=-10.0                                                        00161770
        LFL=0                                                           00161780
      DO 260 I=1,LL                                                     00161790
         IF (I.EQ.1) GO TO 190                                          00161800
         IF (TIM(I).LE.TIM(LY)) GO TO 260                               00161810
        IF((TIM(I)-TIM(LY)).LT.(TH-TL))  TIM(I)=TH*(2.0+1.0E-8)-TL      00161820
190     CONTINUE                                                        00161830
        IF(LFL.EQ.1) GO TO 200                                          00161840
        TL=TIME                                                         00161850
         CALL RDDIS (TIME,1,NT14,KFL,0,I)                               00161860
         IF (KFL.EQ.1) GO TO 270                                        00161870
        TH=TIME                                                         00161880
200      CALL RDDIS (A(N1),NEQ,NT14,KFL,0,I)                            00161890
         IF (KFL.EQ.1) GO TO 320                                        00161900
         CALL RDDIS (A(NV),NEQ,NT14,KFL,0,NREC)                         00161910
         CALL RDDIS (A(NV),NEQ,NT14,KFL,0,NREC)                         00161920
         IF (KFL.EQ.1) GO TO 320                                        00161930
        TL=TIME                                                         00161940
         CALL RDDIS (TIME,1,NT14,KFL,0,I)                               00161950
         IF (KFL.EQ.1) GO TO 270                                        00161960
         TH=TIME                                                        00161970
         IF (I.GT.1) GO TO 210                                          00161980
         IF (TIM(I).LT.TL) GO TO 230                                    00161990
210      DIFF1=TIM(I)-TL                                                00162000
         DIFF2=TIME-TIM(I)                                              00162010
         IF (TIM(I).GE.TL.AND.TIM(I).LE.TH) GO TO 220                   00162020
         GO TO 200                                                      00162030
220      IF (DIFF1.GT.DIFF2) GO TO 240                                  00162040
230     IF(LY.EQ.0) GO TO 231                                           00162050
        IF(TIM(LY).EQ.TL) GO TO 260                                     00162060
231     CONTINUE                                                        00162070
        LY=LY+1                                                         00162080
         TIM(LY)=TL                                                     00162090
        LFL=1                                                           00162100
         GO TO 250                                                      00162110
240     JFL=0                                                           00162120
        IF(LY.EQ.0) GO TO 241                                           00162130
        IF(TIM(LY).EQ.TH) JFL=1                                         00162140
        IF(JFL.EQ.1) GO TO 242                                          00162150
241     LY=LY+1                                                         00162160
         TIM(LY)=TH                                                     00162170
242     CONTINUE                                                        00162180
        LFL=0                                                           00162190
         CALL RDDIS (A(N1),NEQ,NT14,KFL,0,I)                            00162200
         IF (KFL.EQ.1) GO TO 320                                        00162210
         CALL RDDIS (A(NV),NEQ,NT14,KFL,0,NREC)                         00162220
         CALL RDDIS (A(NV),NEQ,NT14,KFL,0,NREC)                         00162230
         IF (KFL.EQ.1) GO TO 320                                        00162240
        IF(JFL.EQ.1) GO TO 260                                          00162250
250      CALL RDWRT (NT,A(N1),NEQ,13,I)                                 00162260
260   CONTINUE                                                          00162270
      GO TO 280                                                         00162280
270     CONTINUE                                                        00162290
      CALL RDWRT (NT,A(N1),NEQ,13,I)                                    00162300
      LY=LY+1                                                           00162310
      TIM(LY)=TL                                                        00162320
280   REWIND NT                                                         00162330
        IF(LY.EQ.1) GO TO 281                                           00162340
        IF(TIM(LY).EQ.TIM(LY-1)) LY=LY-1                                00162350
281     CONTINUE                                                        00162360
      LL=LY                                                             00162370
        NBLOCK=LL                                                       00162380
      CALL SECOND (T(1))                                                00162390
      N2=N1+NUMNP*3                                                     00162400
      N3=N2+6                                                           00162410
      IF (IODISP.EQ.1) TITHOL=TITLE3(3)                                 00162420
      CALL DPRNTD (A(N1),A(N2),A(N3),NUMNP,LL,NEQ,NT)                   00162430
      IF (IODISP.EQ.1) REWIND 6                                         00162440
      IF (IODISP.EQ.1) TITLE3(3)=TITHOL                                 00162450
      CALL SECOND (T(2))                                                00162460
      NADD=1                                                            00162470
      IF (KOPT.GT.0) NADD=NUMNP                                         00162480
      N2=N1                                                             00162490
      N2A=N2+NEQ                                                        00162500
      N3=N2A+NADD                                                       00162510
      NDYN=0                                                            00162520
      KL=1                                                              00162530
      IF (IOSIG.EQ.1) TITHOL=TITLE3(3)                                  00162540
      CALL DSTRES (A(N3),LL,NEQ,A(KL),A(N2A),NADD)                      00162550
      IF (IOSIG.EQ.1) REWIND 6                                          00162560
      IF (IODISP.EQ.1) WRITE (6,300)                                    00162570
      IF (IOSIG.EQ.1) WRITE (6,310)                                     00162580
      IF (IOSIG.EQ.1) TITLE3(3)=TITHOL                                  00162590
      CALL SECOND (T(3))                                                00162600
      T(1)=T(2)-T(1)                                                    00162610
      T(2)=T(3)-T(2)                                                    00162620
      WRITE (6,290) T(1),T(2)                                           00162630
      RETURN                                                            00162640
290   FORMAT (27H1....TIME LOG (CPU SECONDS)///33H  PRINT DISPLACEMENTS.00162650
     1.......... ,F8.2//33H  COMPUTE STRESSES.............. ,F8.2//)    00162660
300   FORMAT (///20X,34(1H*)/20X,34HDISPLACEMENTS WILL NOT BE PRINTED./200162670
     10X,34(1H*)//)                                                     00162680
310   FORMAT (///20X,29(1H*)/20X,29HSTRESSES WILL NOT BE PRINTED./20X,2900162690
     1(1H*)//)                                                          00162700
320   WRITE(6, 330) NT14                                                00162710
330   FORMAT (///20X,35HAN ERROR WAS FOUND ON RESTART TAPE ,I4///)      00162720
      KSKIP=1                                                           00162730
      RETURN                                                            00162740
      END                                                               00162750
      SUBROUTINE RDDIS (A,N,NT,K,J,NREC)                                00196520
        IMPLICIT REAL*8(A-H,O-Z)                                        00196530
      DIMENSION A(N)                                                    00196540
      K=0                                                               00196550
      IF (J.GT.0) GO TO 40                                              00196560
      READ (NT,10,END=20) A                                             00196570
        GO TO 30                                                        00196580
10    FORMAT (5E16.9)                                                   00196590
20    K=1                                                               00196600
30    RETURN                                                            00196610
40    DO 50 I=1,NREC                                                    00196620
         READ (NT,10,END=20)X                                           00196630
50    CONTINUE                                                          00196640
      RETURN                                                            00196650
      END                                                               00196660
      SUBROUTINE DPRNTD (ID,D,B,NUMNP,LL,NEQ,NT)                        00062880
      IMPLICIT REAL*8(A-H,O-Z)                                          00062890
      REAL*8 ID                                                         00062900
      COMMON /PREP/ XMX,XAD,KSKIP,NDYN,I1,RRPREP(7)                     R0062910
      INTEGER DUMMY                                                     00062920
      COMMON /OUT/ NRES,NSTR,NDIS,DUMMY,IOSIG,IODISP,NROUT(4)           R0062930
      COMMON/ELPAR/XPAR(14),NDUM(8),MTOT,NDUM1(6),NUMEL,NUMEL2,NELPA(41)R0062940
      COMMON A(1)                                                       00062950
      DIMENSION ID(NUMNP,3), B(NEQ), D(6)                               00062960
      IF (NDIS.LT.0) RETURN                                             00062970
      REWIND NT                                                         00062980
      REWIND 8                                                          00062990
      READ (8) ID                                                       00063000
      REWIND 17                                                         00063010
      REWIND 18                                                         00063020
      NDPBLK=(MTOT-6-NUMNP)/6                                           00063030
      NBLK=(NUMNP-1)/NDPBLK+1                                           00063040
      KK=1                                                              00063050
      IF (NDIS.GT.0) WRITE (NDIS,20) KK,LL,NDYN,NDIS,NSTR,NUMNP,NUMEL,NU00063060
     1MEL2                                                              00063070
      DO 30 LD=1,LL                                                     00063080
        READ (NT) B                                                     00063090
      DO 30 N=1,NUMNP                                                   00063100
         DO 10 I=1,6                                                    00063110
        D(I)=0.0                                                        00063120
            CALL UNPKID (ID,NUMNP,W,WX,2,N,I)                           00063130
            NNN=W                                                       00063140
            IF (NNN.LT.1) GO TO 10                                      00063150
            D(I)=B(NNN)                                                 00063160
10       CONTINUE                                                       00063170
20    FORMAT (8I5)                                                      00063180
         WRITE (18) D                                                   00063190
30    CONTINUE                                                          00063200
      K=1                                                               00063210
      N2=K+NUMNP                                                        00063220
      N3=N2+6                                                           00063230
      N4=N3+6*NDPBLK                                                    00063240
      CALL WRDISD (A(K),A(N2),A(N3),NUMNP,LL,NDPBLK,NDIS,NBLK)          00063250
      RETURN                                                            00063260
      END                                                               00063270
      SUBROUTINE WRDISD (NORD,A,B,NUMNP,LL,NDPBLK,NDIS,NBLK)            00327500
      REAL*8 A,B                                                        00327510
      DIMENSION NORD(NUMNP), A(6), B(NDPBLK,6)                          00327520
      COMMON /TRASH/ TIM(400),RRTRAS(580)                               R0327530
      COMMON /OUT/ KDUMMY(9),KROT                                       00327540
      COMMON /BAND/ KOPT,NRBAND(7)                                      R0327550
      NT1=17                                                            00327560
      NT2=18                                                            00327570
        REWIND NT2                                                      00327580
      IF (KOPT.GT.0) READ (NT1) NORD                                    00327590
      DO 70 M=1,LL                                                      00327600
         WRITE(6, 10) TIM(M)                                            00327610
10    FORMAT (1X ,/20X,24HDISPLACEMENTS AT TIME = ,G15.6//)             00327620
         WRITE (6,110)                                                  00327630
         KSHF2=0                                                        00327640
         KSHF=1-NDPBLK                                                  00327650
         KNT=0                                                          00327660
         DO 60 I=1,NBLK                                                 00327670
            KOUNT=0                                                     00327680
            KK=NUMNP+1                                                  00327690
            KSHF=KSHF+NDPBLK                                            00327700
            KSHF2=KSHF2+NDPBLK                                          00327710
            DO 30 KK=1,NUMNP                                            00327720
               READ (NT2) A                                             00327730
               KCH=KK                                                   00327740
               IF (KOPT.GT.0) KCH=NORD(KK)                              00327750
               IF (KCH.LT.KSHF.OR.KCH.GT.KSHF2) GO TO 30                00327760
               KNT=KNT+1                                                00327770
               KOUNT=KOUNT+1                                            00327780
               NSHFT=KCH-KSHF+1                                         00327790
               DO 20 K=1,6                                              00327800
20             B(NSHFT,K)=A(K)                                          00327810
               IF (KOUNT.EQ.NDPBLK.OR.KNT.EQ.NUMNP) GO TO 40            00327820
30          CONTINUE                                                    00327830
40          KNT1=(I-1)*NDPBLK                                           00327840
            DO 50 J=1,NDPBLK                                            00327850
               KNT1=KNT1+1                                              00327860
               IF (KNT1.GT.NUMNP) GO TO 70                              00327870
               WRITE (6,80) KNT1,(B(J,K),K=1,6)                         00327880
               IF (NDIS.GT.0) WRITE (NDIS,90) KNT1,(B(J,K),K=1,3)       00327890
      IF(NDIS.GT.0.AND.KROT.EQ.2) WRITE(NDIS,90) KNT1,(B(J,K),K=4,6)    00327900
        WRITE(33,200)KNT1                                               00327910
        WRITE(33,201)(B(J,K),K=1,3)                                     00327920
200     FORMAT(I5)                                                      00327930
201     FORMAT(3(F20.10))                                               00327940
50          CONTINUE                                                    00327950
60       CONTINUE                                                       00327960
        KBL=0                                                           00327970
        WRITE(33,200)KBL                                                00327980
70    IF (KOPT.GT.0) WRITE (6,100)                                      00327990
80    FORMAT (1H0,I4,1P3E12.3,3E11.2)                                   00328000
90    FORMAT (I10,7E10.4/(8E10.4))                                      00328010
100   FORMAT (//3X,46H*** NOTE *** NODE NUMBERS ARE ORIGINAL NUMBERS//) 00328020
110   FORMAT (40H0.......NODE DISPLACEMENTS AND ROTATIONS//5H NODE,7X,  00328030
     11HX,11X,1HY,11X,1HZ,9X,2HXX,9X,2HYY,9X,2HZZ)                      00328040
      RETURN                                                            00328050
      END                                                               00328060
      SUBROUTINE DSTRES (D,LL,NEQ,SA,NORD,NADD)                         00064960
      IMPLICIT REAL*8(A-H,O-Z)                                          00064970
      REAL*8 NPAR                                                       00064980
      DIMENSION D(NEQ)                                                  00064990
      DIMENSION SA(1)                                                   00065000
      COMMON /ELPAR/ NPAR(14),NUMNP,MBAND,NELTYP,N1,N2,N3,N4,N5,MTOT,MEQ00065010
     1,N2P,N3P,NMRI,NTRI,N1P,NUMEL,NUMEL2,KZ(20,2),NEMN                 00065020
      INTEGER DUMMY                                                     00065030
      COMMON /TRASH/ TIM(400),RRTRAS(90)                                R0065040
      COMMON /OUT/ NRES,NSTR,NDIS,DUMMY,IOSIG,IODISP,NROUT(4)           R0065050
      DIMENSION NORD(NADD)                                              00065060
      COMMON /JUNK/ SIG(200),MM,L,KDU,NTAG,NDYN,NRJUNK(49)              R0065070
      COMMON /GPS/ NEQ4(10),NRGPS(10)                                   R0065080
      IF (NSTR.LT.0) RETURN                                             00065090
      NT1=1                                                             00065100
      NT3=3                                                             00065110
      NT9=17                                                            00065120
      IF (NADD.GT.1) REWIND NT9                                         00065130
      IF (NADD.GT.1) READ (NT9) NORD                                    00065140
      NT=LL                                                             00065150
      LB=1                                                              00065160
      LH=0                                                              00065170
      CALL RDWRT (NT3,SA,1,6,J)                                         00065180
      NELAST=0                                                          00065190
      N=2                                                               00065200
      IF (NSTR.GT.0) WRITE (NSTR,10) N,LL,NDYN,NDIS,NSTR,NUMNP,NUMEL,NUM00065210
     1EL2                                                               00065220
10    FORMAT (8I5)                                                      00065230
        NTU=2                                                           00065240
        REWIND NTU                                                      00065250
      DO 260 II=1,NT                                                    00065260
         WRITE(6, 20) TIM(II)                                           00065270
20    FORMAT (1X ,/20X,19HSTRESSES AT TIME = ,G15.6//)                  00065280
         LT=LH+1                                                        00065290
         LLT=1-LT                                                       00065300
         LH=LT+LB-1                                                     00065310
         IF (LH.GT.LL) LH=LL                                            00065320
         IF (NSTR.GT.0) WRITE (NSTR,10) LT,LH                           00065330
        READ (NTU) D                                                    00065340
         LK=LH-LT+1                                                     00065350
         CALL RDWRT (NT1,SA,1,6,J)                                      00065360
         NPAR(1)=0                                                      00065370
         NUME=NUMEL+NUMEL2                                              00065380
      DO 260 MM=1,NUME                                                  00065390
         CALL RDWRT (NT1,SA,NEMN,0,KOUNT)                               00065400
         Z=SA(KOUNT)                                                    00065410
         IF (Z.NE.7.OR.MM.GT.NUMEL) GO TO 60                            00065420
         NELAST=7                                                       00065430
         WRITE (6,30) MM                                                00065440
30    FORMAT (/20X,7HELEMENT,I5,19H IS A BLANK ELEMENT/)                00065450
         IF (NDYN.EQ.3) CALL RDWRT (NT3,Z,1,1,J)                        00065460
         DO 40 L=LT,LH                                                  00065470
40       IF (NSTR.GT.0) WRITE (NSTR,50) L                               00065480
50    FORMAT (3X,1H1,I2,2X,2H 7,6G10.4)                                 00065490
         GO TO 260                                                      00065500
60       CONTINUE                                                       00065510
         NS1=SA(KOUNT-1)                                                00065520
         ND1=SA(KOUNT-2)                                                00065530
         K=0                                                            00065540
         KSTR=ND1*NS1+ND1                                               00065550
         DO 250 L=LT,LH                                                 00065560
            K=K+1                                                       00065570
            LMIN1=(L-1)*NS1+KSTR                                        00065580
            DO 70 N=1,NS1                                               00065590
               NPN=N+ND1                                                00065600
               SIG(N)=0.0                                               00065610
            DO 70 J=1,ND1                                               00065620
               NELM=NPN+(J-1)*NS1                                       00065630
               JJ=SA(J)                                                 00065640
               IF (JJ.LE.0) GO TO 70                                    00065650
               SIG(N)=SIG(N)+SA(NELM)*D(JJ)                             00065660
70          CONTINUE                                                    00065670
            MTYPE=Z                                                     00065680
            NTAG=1                                                      00065690
            IF (MTYPE.NE.NELAST) NTAG=0                                 00065700
            NELAST=MTYPE                                                00065710
            GO TO (80,90,100,110,120,130,140,150,160,210,170,170,170,23000065720
     1      ), MTYPE                                                    00065730
80          CALL DTRUSS                                                 00065740
            GO TO 240                                                   00065750
90          CONTINUE                                                    00065760
            SIG(200)=NS1                                                00065770
            CALL DBEAM                                                  00065780
            GO TO 240                                                   00065790
100         CALL DPLANE                                                 00065800
            GO TO 240                                                   00065810
110         CALL DAXIS                                                  00065820
            GO TO 240                                                   00065830
120         NELM=KSTR                                                   00065840
            SIG(13)=SA(NELM+4)                                          00065850
            SIG(14)=SA(NELM+6)                                          00065860
            SIG(15)=NS1                                                 00065870
            CALL DTHRED                                                 00065880
            GO TO 240                                                   00065890
130         CALL DSHELL                                                 00065900
            GO TO 240                                                   00065910
140     SIG(150)=DABS(SA(KSTR+1))                                       00065920
        SIG(151)=DABS(SA(KSTR+2))                                       00065930
        CALL DBOUND(NORD,NADD)                                          00065940
            GO TO 240                                                   00065950
150         CALL DPLANE                                                 00065960
            GO TO 240                                                   00065970
160         CONTINUE                                                    00065980
            SIG(200)=NS1                                                00065990
        CALL DELBST                                                     00066000
            GO TO 240                                                   00066010
170         NELM=KSTR-1                                                 00066020
            JJ=(NS1-16)/4                                               00066030
            IF (JJ.LE.0) JJ=1                                           00066040
            KK=NS1                                                      00066050
            DO 180 N=1,JJ                                               00066060
               NELM=NELM+4                                              00066070
               KK=KK+1                                                  00066080
180         SIG(KK)=SA(NELM)                                            00066090
            SIG(150)=NS1                                                00066100
            IF (NS1.LE.4) GO TO 200                                     00066110
            JJ=JJ+1                                                     00066120
            JF=NS1/4                                                    00066130
            DO 190 N=JJ,JF                                              00066140
               NELM=NELM+4                                              00066150
               KK=KK+1                                                  00066160
               SIG(KK)=SA(NELM)                                         00066170
               KK=KK+1                                                  00066180
190         SIG(KK)=SA(NELM+1)                                          00066190
200         SIG(151)=Z                                                  00066200
            CALL ST2DD (NORD,NADD)                                      00066210
            GO TO 240                                                   00066220
210         NELM=KSTR-2                                                 00066230
            JJ=NS1/6                                                    00066240
            KK=NS1                                                      00066250
            DO 220 N=1,JJ                                               00066260
               NELM=NELM+6                                              00066270
               KK=KK+1                                                  00066280
220         SIG(KK)=SA(NELM)                                            00066290
            SIG(150)=NS1                                                00066300
            CALL ST3DD (NORD,NADD)                                      00066310
            GO TO 240                                                   00066320
230         CALL SIXSTD                                                 00066330
240         CONTINUE                                                    00066340
250      CONTINUE                                                       00066350
260   CONTINUE                                                          00066360
      RETURN                                                            00066370
      END                                                               00066380
      SUBROUTINE DTRUSS                                                 00066920
      IMPLICIT REAL*8  (A-H,O-Z)                                        00066930
        REAL*8 NPAR                                                     00066940
      COMMON A(1)                                                       00066950
      COMMON /ELPAR/ NPAR(14),NUMNP,MBAND,NELTYP,N1,N2,N3,N4,N5,MTOT,NEQ00066960
     1,N2P,N3P,NMRI,NTRI,N1P,NUMEL,NUMEL2,KZ(20,2),NEMN                 00066970
      COMMON /JUNK/ SIG(12),EXTRA(188),MM,L,K,NTAG,NDYN,NRJUNK(49)      R0066980
      COMMON /OUT/ NRES,NSTR,NDIS,NROUT(7)                              R0066990
10    IF (NTAG.EQ.0) WRITE (6,30)                                       00067000
      WRITE (6,40) MM,L,SIG(1),SIG(2)                                   00067010
      NTAG=1                                                            00067020
      IF (NSTR.GT.0) WRITE (NSTR,20) L,SIG(1),SIG(2)                    00067030
20    FORMAT (3X,1H2,I2,2X,2H 1,6G10.4)                                 00067040
      RETURN                                                            00067050
30    FORMAT (23H1  TRUSS MEMBER ACTIONS//46H0 MEMBER    LOAD         ST00067060
     1RESS          FORCE/)                                             00067070
40    FORMAT (2I8,F15.0,F15.0)                                          00067080
      END                                                               00067090
      SUBROUTINE DBEAM                                                  00052470
      IMPLICIT REAL*8  (A-H,O-Z)                                        00052480
      COMMON /JUNK/ SIG(26),EXTRA(174),MM,L,K,NTAG,NDYN,NRJUNK(49)      R0052490
      COMMON /OUT/ NRES,NSTR,NDIS,NBMSTR,NROUT(6)                       R0052500
10    IF (NTAG.EQ.0) WRITE (6,30)                                       00052510
      IF (NTAG.EQ.0.AND.NBMSTR.EQ.1) WRITE (6,40)                       00052520
      IF (NTAG.EQ.0) WRITE (6,50)                                       00052530
      NS=EXTRA(174)                                                     00052540
      IF (NS.GT.12) WRITE (6,60) MM,L,(SIG(I1),I1=1,6),(SIG(I2),I2=13,1900052550
     1),(SIG(I3),I3=7,12),(SIG(I4),I4=20,26)                            00052560
      IF (NS.LE.12) WRITE (6,70) MM,L,(SIG(I1),I1=1,12)                 00052570
      NTAG=1                                                            00052580
      IF (NSTR.GT.0) WRITE (NSTR,20) NS,L,(SIG(I),I=1,NS)               00052590
20    FORMAT (I4,I2,2X,2H 2,7G10.4/(8G10.4))                            00052600
      RETURN                                                            00052610
30    FORMAT (/59H1....BEAM FORCES(LBS), MOMENTS(IN-LBS), AND STRESSES(P00052620
     1SI)  //1X,7HELEMENT,2X,4HLOAD,2X,7HSTATION,15X,5HAXIAL,2(6X,5HSHEA00052630
     2R),4X,9HTORSION  ,2(5X,13HB E N D I N G,4X)/2X,6HNUMBER,2X,4HCASE,00052640
     311X,5HFORCE,11X,2HR1,9X,2HR2,9X,2HR3,9X,2HM1,20X,2HM2,20X,2HM3)   00052650
40    FORMAT (25X,6HSTRESS,7X,5HR1/A1,3X,8HAT Q3/B3,4X,8HAT Q2/B2,5X,4H-00052660
     1 - ,1X,2X,9H   AT +C3,2X,9H   AT -C3,2X,9H   AT +C2,2X,9H   AT -C200052670
     2)                                                                 00052680
50    FORMAT (1X)                                                       00052690
60    FORMAT (4X,I4,2X,I4,4X,5HEND-I,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X00052700
     1,F11.0/23X,2X,6HSTRESS,1X,3F11.0,11X,4F11.0/18X,5HEND-J,2X,6HFORCE00052710
     2 ,1X,4F11.0,11X,F11.0,11X,F11.0/23X,2X,6HSTRESS,1X,3F11.0,11X,4F1100052720
     3.0)                                                               00052730
70    FORMAT (4X,I4,2X,I4,4X,5HEND-I,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X00052740
     1,F11.0/18X,5HEND-J,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X,F11.0)     00052750
      END                                                               00052760
      SUBROUTINE DPLANE                                                 00062090
      IMPLICIT REAL*8  (A-H,O-Z)                                        00062100
      COMMON /OUT/ NRES,NSTR,NDIS,NROUT(7)                              R0062110
      COMMON /JUNK/ SIG(12),EXTRA(188),MM,L,K,NTAG,NDYN,NRJUNK(49)      R0062120
10    IF (NTAG.EQ.0) WRITE (6,40)                                       00062130
      CC=(SIG(1)+SIG(2))/2.0                                            00062140
      BB=(SIG(1)-SIG(2))/2.                                             00062150
      CR=DSQRT(BB**2+SIG(4)**2)                                         00062160
      SIG(5)=CC+CR                                                      00062170
      SIG(6)=CC-CR                                                      00062180
      SIG(7)=0.0                                                        00062190
      EF=(SIG(3)-SIG(5))**2+(SIG(5)-SIG(6))**2+(SIG(6)-SIG(3))**2       00062200
      EF=DSQRT(EF/2.)                                                   00062210
      IF ((BB.EQ.0.0).AND.(SIG(4).EQ.0.0)) GO TO 20                     00062220
      SIG(7)=28.648D0*DATAN2(SIG(4),BB)                                 00062230
20    WRITE (6,50) MM,L,(SIG(I),I=1,7),EF                               00062240
      NTAG=1                                                            00062250
      IF (NSTR.GT.0) WRITE (NSTR,30) L,SIG(5),SIG(6),EF,(SIG(I),I=1,4)  00062260
30    FORMAT (3X,1H7,I2,4H12 3,7G10.4)                                  00062270
      RETURN                                                            00062280
40    FORMAT (1X ,37H PLANE STRESS/STRAIN ELEMENT STRESSES//            00062290
     112H0EL.NO. LOAD,7X,8HR-STRESS,7X,                                 00062300
     18HZ-STRESS,7X,8HT-STRESS,6X,9HRZ-STRESS,5X,10HMAX-STRESS,5X,      00062310
     210HMIN-STRESS,3X,5HANGLE,5X,6HSIG-EF/)                            00062320
50    FORMAT (2I6,6F15.0,F8.3,F10.0)                                    00062330
      END                                                               00062340
      SUBROUTINE DAXIS                                                  00052180
      IMPLICIT REAL*8  (A-H,O-Z)                                        00052190
        REAL*8 NPAR                                                     00052200
      COMMON /JUNK/ SIG(200),MM,L,K,NTAG,NDYN,NRJUNK(49)                R0052210
      COMMON /ELPAR/ NPAR(14),NUMNP,MBAND,NELTYP,N1,N2,N3,N4,N5,MTOT,NEQ00052220
     1,N2P,N3P,NMRI,NTRI,N1P,NUMEL,NUMEL2,KZ(20,2),NEMN                 00052230
      COMMON /OUT/ NRES,NSTR,NDIS,NROUT(7)                              R0052240
10    IF (NTAG.EQ.0) WRITE (6,40)                                       00052250
      CC=(SIG(1)+SIG(2))/2.0                                            00052260
      BB=(SIG(1)-SIG(2))/2.                                             00052270
      CR=DSQRT(BB**2+SIG(4)**2)                                         00052280
      SIG(5)=CC+CR                                                      00052290
      SIG(6)=CC-CR                                                      00052300
      SIG(7)=0.0                                                        00052310
      EF=(SIG(3)-SIG(5))**2+(SIG(5)-SIG(6))**2+(SIG(6)-SIG(3))**2       00052320
      EF=DSQRT(EF/2.)                                                   00052330
      IF ((BB.EQ.0.0).AND.(SIG(4).EQ.0.0)) GO TO 20                     00052340
      SIG(7)=28.648D0*DATAN2(SIG(4),BB)                                 00052350
20    WRITE (6,50) MM,L,(SIG(I),I=1,7),EF                               00052360
      IF (NSTR.GT.0) WRITE (NSTR,30) L,SIG(5),SIG(6),EF,(SIG(I),I=1,4)  00052370
30    FORMAT (3X,1H7,I2,4H12 4,7G10.4)                                  00052380
      NTAG=1                                                            00052390
      RETURN                                                            00052400
40    FORMAT (1X ,30H AXISYMMETRIC ELEMENT STRESSES//                   00052410
     112H0EL.NO. LOAD,7X,8HR-STRESS,                                    00052420
     17X,8HZ-STRESS,7X,8HT-STRESS,6X,9HRZ-STRESS,5X,10HMAX-STRESS,5X,   00052430
     210HMIN-STRESS,3X,5HANGLE,5X,6HSIG-EF/)                            00052440
50      FORMAT(2I6,6F15.0,F8.3,F10.0)                                   00052450
      END                                                               00052460
      SUBROUTINE DTHRED                                                 00066390
      IMPLICIT REAL*8(A-H,O-Z)                                          00066400
      REAL*8 NS                                                         00066410
      REAL*8 IS1,IS2,IS3,IS4                                            00066420
      DIMENSION SPR(6)                                                  00066430
      DIMENSION EFS(2)                                                  00066440
      COMMON /JUNK/ SIG(15),EXRA(185),MM,L,K,NTAG,NDYN,NRJUNK(49)       R0066450
      COMMON /OUT/ NRES,NSTR,NDIS,NROUT(7)                              R0066460
      EQUIVALENCE (IS3,SIG(13)), (IS4,SIG(14)), (NS,SIG(15))            00066470
10    CONTINUE                                                          00066480
      IF (NTAG.EQ.0) WRITE (6,50)                                       00066490
      NG=NS                                                             00066500
      IS1=IS3                                                           00066510
      IS2=IS4                                                           00066520
      IF (IS1.LT.1) IS1=IS1*100.+0.001                                  00066530
      IF (IS2.LT.1) IS2=IS2*100.+0.001                                  00066540
      NF1=IS1                                                           00066550
      NF2=IS2                                                           00066560
      IS1=NF1                                                           00066570
      IS2=NF2                                                           00066580
      IF (NF1.EQ.0) NF1=7                                               00066590
      IF (NF2.EQ.0) NF2=7                                               00066600
      NF=10*NF1+NF2                                                     00066610
      IF (NF.EQ.77) NF=71                                               00066620
      EFS(2)=0.0                                                        00066630
      SPR(4)=0.0                                                        00066640
      SPR(5)=0.0                                                        00066650
      CALL DPRIST (IS1,IS2,SIG,SPR,NS)                                  00066660
      II=NG/6                                                           00066670
      DO 30 J=1,II                                                      00066680
         M=6*J-6                                                        00066690
         EFS(J)=(SIG(M+1)-SIG(M+2))**2+(SIG(M+2)-SIG(M+3))**2+(SIG(M+3)-00066700
     1   SIG(M+1))**2                                                   00066710
         DO 20 I=4,6                                                    00066720
            MI=M+I                                                      00066730
20       EFS(J)=EFS(J)+6.*SIG(MI)**2                                    00066740
30    EFS(J)=DSQRT(EFS(J)/2.)                                           00066750
      WRITE (6,60) MM,L,IS1,(SIG(I),I=1,6),(SPR(I),I=1,3),EFS(1)        00066760
      IF (NG.EQ.12) WRITE (6,70) IS2,(SIG(I),I=7,12),(SPR(I),I=4,6),EFS(00066770
     12)                                                                00066780
      II=6+NG                                                           00066790
      IF (NSTR.GT.0) WRITE (NSTR,40) II,L,NF,SPR(1),SPR(2),EFS(1),SPR(4)00066800
     1,SPR(5),EFS(2),(SIG(I),I=1,NG)                                    00066810
40    FORMAT (I4,I2,I2,2H 5,7G10.4/(8G10.4))                            00066820
      RETURN                                                            00066830
50    FORMAT (1X ,35H .....8-NODE SOLID ELEMENT STRESSES//              00066840
     113H ELEM. LOAD F,                                                 00066850
     162HACE    SIG-XX    SIG-YY    SIG-ZZ    SIG-XY    SIG-YZ   SIG-ZX,00066860
     238H   SIG-MAX   SIG-MIN   S2 OR    SIG-EF/16H   NO.  NO.  NO.,84X,00066870
     35HANGLE/)                                                         00066880
60    FORMAT (I5,I4,F5.0,1X,10F10.0)                                    00066890
70    FORMAT (9X,F5.0,1X,10F10.0)                                       00066900
      END                                                               00066910
      SUBROUTINE DPRIST (IS1,IS2,SIG,SPR,NS)                            00062350
      IMPLICIT REAL*8  (A-H,O-Z)                                        00062360
      REAL*8 IS1,IS2,NS                                                 00062370
      DIMENSION SIG(12), SPR(6), IS(2), SG(6)                           00062380
      IS(1)=IS1                                                         00062390
      IS(2)=IS2                                                         00062400
      NNS=1                                                             00062410
      IF (NS.EQ.12.) NNS=2                                              00062420
      DO 60 N=1,NNS                                                     00062430
         IN=3*N-3                                                       00062440
         II=IN*2                                                        00062450
         IF (IS(N).EQ.0) GO TO 10                                       00062460
         CC=(SIG(II+1)+SIG(II+2))/2.                                    00062470
         BB=(SIG(II+1)-SIG(II+2))/2.                                    00062480
         CR=DSQRT(BB**2+SIG(II+4)**2)                                   00062490
         SPR(IN+1)=CC+CR                                                00062500
         SPR(IN+2)=CC-CR                                                00062510
         SPR(IN+3)=0.                                                   00062520
         IF (BB.NE.0.) SPR(IN+3)=28.648*DATAN2(SIG(II+4),BB)            00062530
         GO TO 60                                                       00062540
10       CC=(SIG(II+1)+SIG(II+2)+SIG(II+3))/3.                          00062550
         DO 20 I=1,3                                                    00062560
            SG(I)=SIG(II+I)-CC                                          00062570
20       SG(I+3)=SIG(II+I+3)                                            00062580
         C2=(SG(1)**2+SG(2)**2+SG(3)**2)*.5+SG(4)**2+SG(5)**2+SG(6)**2  00062590
         C3=SG(1)*(SG(2)*SG(3)-SG(5)*SG(5))+SG(4)*(SG(5)*SG(6)-SG(4)*SG(00062600
     1   3))+SG(6)*(SG(4)*SG(5)-SG(2)*SG(6))                            00062610
         IF (C2.EQ.0.0) C2=1.0E-08                                      00062620
         T=DSQRT(C2/1.5)                                                00062630
         A=C3*1.414214/T**3                                             00062640
         AXQ=A+1.0                                                      00062650
         IF (AXQ.GT.0.0.AND.AXQ.LT.2.0) A=DACOS(A)/3.0                  00062660
         IF (AXQ.GE.2.0) A=0.0                                          00062670
         IF (AXQ.LE.0.0) A=1.0471976                                    00062680
         T=T*1.414214                                                   00062690
         SPR(IN+1)=T*DCOS(A)                                            00062700
         SPR(IN+2)=T*DCOS(A+2.0944)                                     00062710
         SPR(IN+3)=T*DCOS(A-2.0944)                                     00062720
         DO 30 I=2,3                                                    00062730
            IF (SPR(IN+1).GT.SPR(IN+I)) GO TO 30                        00062740
            C3=SPR(IN+1)                                                00062750
            SPR(IN+1)=SPR(IN+I)                                         00062760
            SPR(IN+I)=C3                                                00062770
30       CONTINUE                                                       00062780
         IF (SPR(IN+2).LE.SPR(IN+3)) GO TO 40                           00062790
         C3=SPR(IN+2)                                                   00062800
         SPR(IN+2)=SPR(IN+3)                                            00062810
         SPR(IN+3)=C3                                                   00062820
40       DO 50 I=1,3                                                    00062830
50       SPR(IN+I)=SPR(IN+I)+CC                                         00062840
60    CONTINUE                                                          00062850
      RETURN                                                            00062860
      END                                                               00062870
      SUBROUTINE DBOUND(NORD,NADD)                                      00052770
      IMPLICIT REAL*8  (A-H,O-Z)                                        00052780
      COMMON /JUNK/ SIG(12),EXRA(188),MM,L,K,NTAG,NDYN,NRJUNK(49)       R0052790
        DIMENSION NORD(NADD)                                            00052800
10    IF (NTAG.EQ.0) WRITE (6,20)                                       00052810
        NN=EXRA(138)*10000.+.001                                        00052820
        IF(EXRA(138).GT.1.0) NN=EXRA(139)*10000.+.001                   00052830
        IF(NADD.GT.1.AND.NN.LE.NADD.AND.NN.GT.0) NN=NORD(NN)            00052840
        IF(NN.LE.0) NN=0                                                00052850
        WRITE(6,30) MM,NN,L,SIG(1),SIG(2)                               00052860
      RETURN                                                            00052870
20      FORMAT(1X ,39H   CONSTRAINT FORCE - BOUNDARY ELEMENTS //        00052880
     $77H   EL. NUMBER     NODE N     LOAD CASE          FORCE(LBS.)   M00052890
     $OMENT(IN-LBS.) /)                                                 00052900
30      FORMAT(1X,I10,I12,I11,4X,2F15.0)                                00052910
      END                                                               00052920
      SUBROUTINE DELBST                                                 00055930
      IMPLICIT REAL*8(A-H,O-Z)                                          00055940
      COMMON /JUNK/ SIG(39),EXR(161),MM,L,K,NTAG,RRJUNK(25)             R0055950
      COMMON /OUT/ N,NSTR,NDIS,NBMSTR,NROUT(6)                          R0055960
      IF (NTAG.EQ.0) WRITE (6,20)                                       00055970
      IF (NTAG.EQ.0.AND.NBMSTR.EQ.1) WRITE (6,30)                       00055980
      IF (NTAG.EQ.0) WRITE (6,10)                                       00055990
10    FORMAT (1X)                                                       00056000
20    FORMAT (85H1E L B O W - F O R C E S(LBS), M O M E N T S(IN-LBS), A00056010
     1 N D  S T R E S S E S(PSI)    //1X,7HELEMENT,2X,4HLOAD,2X,7HSTATIO00056020
     2N,15X,5HAXIAL,2(6X,5HSHEAR),4X,9HTORSION  ,2(5X,13HB E N D I N G,400056030
     3X)/2X,6HNUMBER,2X,4HCASE,11X,5HFORCE,11X,2HRX,9X,2HRY,9X,2HRZ,9X,200056040
     4HMX,20X,2HMY,20X,2HMZ)                                            00056050
30    FORMAT (25X,6HSTRESS,7X,5HRX/A1,3X,8HAT Q3/B3,4X,8HAT Q2/B2,5X,4H-00056060
     1 - ,1X,2X,9H   AT +C3,2X,9H   AT -C3,2X,9H   AT +C2,2X,9H   AT -C200056070
     2)                                                                 00056080
      NS=EXR(161)                                                       00056090
      IF (NS.GT.18) WRITE (6,50) MM,L,(SIG(I1),I1=1,6),(SIG(I2),I2=19,2500056100
     1),(SIG(I3),I3=7,12),(SIG(I4),I4=26,32),(SIG(I5),I5=13,18),(SIG(I6)00056110
     2,I6=33,39)                                                        00056120
      IF (NS.LE.18) WRITE (6,60) MM,L,(SIG(I1),I1=1,18)                 00056130
      NTAG=1                                                            00056140
      IF (NSTR.GT.0) WRITE (NSTR,40) NS,L,(SIG(I),I=1,NS)               00056150
40    FORMAT (I4,I2,2X,2H 9,7G10.4/(8G10.4))                            00056160
50    FORMAT (4X,I4,2X,I4,4X,5HEND-I,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X00056170
     1,F11.0/23X,2X,6HSTRESS,1X,3F11.0,11X,4F11.0/17X,6HCENTER,2X,6HFORC00056180
     2E ,1X,4F11.0,11X,F11.0,11X,F11.0/23X,2X,6HSTRESS,1X,3F11.0,11X,4F100056190
     31.0/18X,5HEND-J,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X,F11.0/23X,2X,600056200
     4HSTRESS,1X,3F11.0,11X,4F11.0)                                     00056210
60    FORMAT (4X,I4,2X,I4,4X,5HEND-I,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X00056220
     1,F11.0/17X,6HCENTER,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X,F11.0/18X,00056230
     25HEND-J,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X,F11.0)                00056240
      RETURN                                                            00056250
      END                                                               00056260
      SUBROUTINE ST2DD (NORD,NADD)                                      00262790
      IMPLICIT REAL*8(A-H,O-Z)                                          00262800
      COMMON /JUNK/ SIG(200),MM,L,K,NTAG,RRJUNK(25)                     R0262810
      COMMON /BAND/ KOPT,NRBAND(7)                                      R0262820
      COMMON /OUT/ NRES,NSTR,NROUT(8)                                   R0262830
      COMMON /PREP/ XMX,XAD,NRPREP(17)                                  R0262840
      DIMENSION NORD(NADD)                                              00262850
      COMMON /QTSARG/ SIGEX(100),X(2),RRQTSA(898)                       R0262860
      DIMENSION C(4), D(4), Q(3)                                        00262870
      DATA Q/4HNODE,4HC.G.,4HI.P./                                      00262880
      JJ=SIG(150)/4                                                     00262890
      MTYP=SIG(151)                                                     00262900
      KL=0                                                              00262910
      IF (NTAG.GT.0) GO TO 70                                           00262920
      IF (MTYP.EQ.11.AND.NTAG.EQ.0) WRITE (6,10)                        00262930
      IF (MTYP.EQ.12.AND.NTAG.EQ.0) WRITE (6,20)                        00262940
      IF (MTYP.EQ.13.AND.NTAG.EQ.0) WRITE (6,30)                        00262950
10    FORMAT (1H0,32H .... AXISYMMETRIC ELEMENTS ..../)                 00262960
20    FORMAT (1H0,32H .... PLANE STRESS ELEMENTS ..../)                 00262970
30    FORMAT (1H0,32H .... PLANE STRAIN ELEMENTS ..../)                 00262980
      IF (NTAG.EQ.0.AND.MTYP.NE.11) WRITE (6,40)                        00262990
      IF (NTAG.EQ.0.AND.MTYP.EQ.11) WRITE (6,50)                        00263000
40    FORMAT (12H0EL.NO. LOAD,12X,                                      00263010
     1                          90HY-STRESS   Z-STRESS  YZ-STRESS   T-ST00263020
     1RESS MAX-STRESS MIN-STRESS  ANGLE   SIG-EF   LOCATION/)           00263030
50    FORMAT (12H0EL.NO. LOAD,12X,                                      00263040
     1                          90HR-STRESS   Z-STRESS  RZ-STRESS   T-ST00263050
     1RESS MAX-STRESS MIN-STRESS  ANGLE   SIG-EF   LOCATION/)           00263060
      IF (NTAG.EQ.0.AND.KOPT.GT.1.AND.NADD.EQ.1) WRITE (6,60)           00263070
60    FORMAT (5X,39HNOTE - NODES ARE NEW RENUMBERED NODES  /)           00263080
70    KK=SIG(150)                                                       00263090
      JL=1000                                                           00263100
      IF (JJ.GT.1) JL=(KK-16)/4                                         00263110
      DO 140 I=1,JJ                                                     00263120
         KK=KK+1                                                        00263130
         FACE=SIG(KK)*10000.+.001                                       00263140
         NF=FACE                                                        00263150
         LL=NF                                                          00263160
         A=Q(2)                                                         00263170
         IF (I.EQ.1) GO TO 80                                           00263180
         IF (NADD.GT.1.AND.I.LE.JL) NF=NORD(LL)                         00263190
         A=Q(1)                                                         00263200
         IF (I.EQ.1) A=Q(2)                                             00263210
         IF (I.LE.JL) GO TO 80                                          00263220
         A=Q(3)                                                         00263230
         NF=I-JL                                                        00263240
         LL=SIG(KK)                                                     00263250
         IF (SIG(KK).LT.0.0) LL=LL-1                                    00263260
         X(1)=(SIG(KK)-LL-XAD)*XMX                                      00263270
         KK=KK+1                                                        00263280
         LL=SIG(KK)                                                     00263290
         IF (SIG(KK).LT.0.0) LL=LL-1                                    00263300
         X(2)=(SIG(KK)-LL-XAD)*XMX                                      00263310
80       CONTINUE                                                       00263320
         LL=4*(I-1)                                                     00263330
         DO 90 J=1,4                                                    00263340
            LL=LL+1                                                     00263350
90       D(J)=SIG(LL)                                                   00263360
         CALL MXMND (D,C(1),C(2),C(3),C(4))                             00263370
         DO 100 K=1,4                                                   00263380
            KL=KL+1                                                     00263390
100      SIGEX(KL)=C(K)                                                 00263400
         IF (I.EQ.1) WRITE (6,120) MM,L,D,C,A                           00263410
         IF (I.GT.1.AND.I.LE.JL) WRITE (6,130) D,C,A,NF                 00263420
         IF (I.GT.JL) SIGEX(KL-1)=X(1)                                  00263430
         IF (I.GT.JL) SIGEX(KL)=X(2)                                    00263440
         IF (I.GT.JL) GO TO 140                                         00263450
         IF (I.GT.JL.AND.L.EQ.1) WRITE (6,110) X(1),X(2),D,C,A,NF       00263460
         IF (I.GT.JL.AND.L.GT.1) WRITE (6,130) D,C,A,NF                 00263470
110   FORMAT (3H Y=,F8.2,3H Z=,F8.2,F10.0,5F11.0,F7.2,F9.0,2X,A4,I5)    00263480
120   FORMAT (2I6,9X,6F11.0,F7.2,F9.0,2X,A4,I5)                         00263490
130   FORMAT (21X,6F11.0,F7.2,F9.0,2X,A4,I5)                            00263500
140   CONTINUE                                                          00263510
      KK=SIG(150)                                                       00263520
      LL=JJ*4                                                           00263530
      NS=KK+LL                                                          00263540
      IF (NSTR.GT.0) WRITE (NSTR,150) NS,L,MTYP,(SIGEX(I),I=1,LL),(SIG(I00263550
     1),I=1,KK)                                                         00263560
150   FORMAT (I4,I2,2X,I2,7G10.4/(8G10.4))                              00263570
      RETURN                                                            00263580
      END                                                               00263590
      SUBROUTINE MXMND (STRESS,P1,P2,AG,EF)                             00152110
      IMPLICIT REAL*8 (A-H,O-Z)                                         00152120
      DIMENSION STRESS(1)                                               00152130
      CC=(STRESS(1)+STRESS(2))*0.5E0                                    00152140
      BB=(STRESS(1)-STRESS(2))*0.5E0                                    00152150
      CR=DSQRT(BB**2+STRESS(3)**2)                                      00152160
      P1=CC+CR                                                          00152170
      P2=CC-CR                                                          00152180
      AG=45.0E0                                                         00152190
      IF (DABS(BB).LT.1.0D-8) GO TO 10                                  00152200
      AG=28.648D0*DATAN2(STRESS(3),BB)                                  00152210
10    EF=(P1-P2)**2+(P2-STRESS(4))**2+(STRESS(4)-P1)**2                 00152220
      EF=DSQRT(EF/2.0)                                                  00152230
      RETURN                                                            00152240
      END                                                               00152250
      SUBROUTINE ST3DD (NORD,NADD)                                      00264960
      IMPLICIT REAL*8(A-H,O-Z)                                          00264970
      COMMON /JUNK/ SIG(200),MM,L,K,NTAG,RRJUNK(25)                     R0264980
      COMMON /BAND/ KOPT,NRBAND(7)                                      R0264990
      COMMON /OUT/ NR,NSTR,NROUT(8)                                     R0265000
      COMMON /QTSARG/ SIGEX(90),RRQTSA(910)                             R0265010
      DIMENSION C(6), D(12), Q(2)                                       00265020
      DIMENSION NORD(NADD)                                              00265030
      DATA Q/4HNODE,4HC.G./                                             00265040
      JJ=SIG(150)/6.                                                    00265050
      IF (NTAG.GT.0) GO TO 40                                           00265060
      IF (NTAG.EQ.0) WRITE (6,10)                                       00265070
10    FORMAT (1H0,39H .... 3-D SOLID ELEMENTS 8-20NODE .... /)          00265080
      IF (NTAG.EQ.0.AND.KOPT.GT.1.AND.NADD.EQ.1) WRITE (6,20)           00265090
20    FORMAT (5X,35H NODE NUMBERS ARE RENUMBERED NODES /)               00265100
      IF (NTAG.EQ.0) WRITE (6,30)                                       00265110
30    FORMAT (118H ELEM. LOAD LOC.    SIG-XX    SIG-YY    SIG-ZZ    SIG-00265120
     1XY    SIG-YZ    SIG-ZX    SIG-MAX    SIG-MIN    SIG-MID   SIG-EF/)00265130
40    KK=SIG(150)                                                       00265140
      A=0.                                                              00265150
      B=0.                                                              00265160
      E=6.                                                              00265170
      KL=0                                                              00265180
      DO 100 I=1,JJ                                                     00265190
         KK=KK+1                                                        00265200
         FACE=SIG(KK)*10000.+.001                                       00265210
         NF=FACE                                                        00265220
         LL=NF                                                          00265230
         IF (NADD.GT.1.AND.I.GT.1) NF=NORD(LL)                          00265240
         LL=6*(I-1)                                                     00265250
         DO 50 J=1,6                                                    00265260
            LL=LL+1                                                     00265270
50       D(J)=SIG(LL)                                                   00265280
         CALL DPRIST (A,B,D,C,E)                                        00265290
         EFS=(C(1)-C(2))**2+(C(2)-C(3))**2+(C(3)-C(1))**2               00265300
         EFS=DSQRT(EFS/2.)                                              00265310
         DO 60 K=1,3                                                    00265320
            KL=KL+1                                                     00265330
60       SIGEX(KL)=C(K)                                                 00265340
         KL=KL+1                                                        00265350
         SIGEX(KL)=EFS                                                  00265360
         IF (I.GT.1) GO TO 80                                           00265370
         WRITE (6,70) MM,L,Q(2),(D(J),J=1,6),(C(J),J=1,3),EFS           00265380
70    FORMAT (I5,I4,1X,A4,1X,10F10.0)                                   00265390
         GO TO 100                                                      00265400
80       WRITE (6,90) Q(1),NF,(D(J),J=1,6),(C(J),J=1,3),EFS             00265410
90    FORMAT (5X,A4,I5,1X,10F10.0)                                      00265420
100   CONTINUE                                                          00265430
      KK=SIG(150)                                                       00265440
      LL=JJ*4                                                           00265450
      NS=KK+LL                                                          00265460
      IF (NSTR.GT.0) WRITE (NSTR,110) NS,L,(SIGEX(I),I=1,LL),(SIG(I),I=100265470
     1,KK)                                                              00265480
110   FORMAT (I4,I2,2X,2H10,7G10.4/(8G10.4))                            00265490
      RETURN                                                            00265500
      END                                                               00265510
      SUBROUTINE SIXSTD                                                 00238470
        IMPLICIT REAL*8(A-H,O-Z)                                        00238480
      COMMON /OUT/ NRES,NSTR,NROUT(8)                                   R0238490
      COMMON /JUNK/ SIG(6),EXTRA(194),MM,L,K,NTAG,RRJUNK(25)            R0238500
      IF (NTAG.EQ.0) WRITE (6,20)                                       00238510
      WRITE (6,10) MM,L,SIG                                             00238520
      NTAG=1                                                            00238530
      RETURN                                                            00238540
10    FORMAT (2I6,6F15.1)                                               00238550
20    FORMAT (//10X,48H SIX BY SIX STIFFNESS ELEMENT FORCES AND MOMENTS/00238560
     120X,18H(LOCAL DIRECTIONS)/                                        00238570
     112H0EL.NO. LOAD,5X,10HFORCE - XX,5X,10HFORCE - YY,                00238580
     25X,10HFORCE - ZZ,4X,11HMOMENT - XX,4X,11HMOMENT - YY,4X,          00238590
     311HMOMENT - ZZ)                                                   00238600
      END                                                               00238610
