home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 9 / 09.iso / e / e003 / 10.ddi / FRFREQT.FOR next >
Encoding:
Text File  |  1987-06-23  |  57.7 KB  |  721 lines
  1.       SUBROUTINE FRFREQ                                                 00095590
  2.         IMPLICIT REAL*8 (A-H,O-Z)                                       00095600
  3.         REAL DIR(3)                                                     00095610
  4.         COMMON A(1)                                                     00095620
  5.         COMMON /MISC/NBLOCK,NEQB,LL,NF,LB                               R0095630
  6.         COMMON /ELPAR/ APAR(14),NUMNP,MBAND,NLTP,N1,MMM(4),MTOT,NEQ     00095640
  7.      & ,RRELPA(24)                                                      R0095641
  8.         DATA FIRST /8HSINE-RES/                                         00095650
  9.         DATA DIR/1HX,1HY,1HZ/                                           00095660
  10.         READ(5,999)TIT                                                  00095670
  11.         READ(5,1000)NGM,IDIR,NFN,IPRINT,XREF,YREF,ZREF                  00095680
  12.         NGMX=NGM                                                        00095690
  13.         WRITE(6,1010)                                                   00095700
  14.         IF(NGM.NE.2)WRITE(6,2000)DIR(IDIR),NFN                          00095710
  15.         IF(NGM.EQ.2)WRITE(6,2001)DIR(IDIR),XREF,YREF,ZREF,NFN           00095720
  16.         IF(IPRINT.EQ.1)WRITE(6,2010)                                    00095730
  17.         NFN=NFN+NF                                                      00095740
  18.         N2=N1+NUMNP*3                                                   00095750
  19.         N3=N2+NEQB                                                      00095760
  20.         CALL FEMID(A(N1),A(N2),NUMNP,NEQB)                              00095770
  21.         N2=N1+NF                                                        00095780
  22.         N3=N2+NF                                                        00095790
  23.         N4=N3+NEQB*NF                                                   00095800
  24.         N5=N4+NEQB                                                      00095810
  25.         N6=N5+NEQB                                                      00095820
  26.         CALL FRMPF(A(N1),A(N2),A(N3),A(N4),A(N5),IDIR                   00095830
  27.      &  ,NEQB,NF,NBLOCK,NUMNP)                                          00095840
  28.         N4=N3+NF                                                        00095850
  29.         N5=N4+NFN*NF                                                    00095860
  30.         N6=N5+NFN                                                       00095870
  31.         N7=N6+NFN                                                       00095880
  32.         CALL FRESP(A(N1),A(N2),A(N3),A(N4),A(N5),A(N6)                  00095890
  33.      &  ,NF,NFN,IPRINT)                                                 00095900
  34.         NSB=NEQB*NBLOCK                                                 00095910
  35.         N2=N1+8*NF                                                      00095920
  36.         N3=N2+NF*NFN                                                    00095930
  37.         N4=N3+NFN                                                       00095940
  38.         N5=N4+NUMNP                                                     00095950
  39.         IF(N5.GT.MTOT)CALL ERROR(N5-MTOT)                               00095960
  40.         CALL FRDSPR(A(N1),A(N1),A(N2),A(N2),A(N3)                       00095970
  41.      &  ,NEQB,NF,NFN,NUMNP,NBLOCK,NSB,A(N4))                            00095980
  42.         N2=N1+1                                                         00095990
  43.         N3=N2+8*NF                                                      00096000
  44.         N4=N3+NSB*NF                                                    00096010
  45.         N5=N3+NF*NFN                                                    00096020
  46.         N4=MAX0(N4,N5)                                                  00096030
  47.         N5=N4+NFN                                                       00096040
  48.         IF(N5.GT.MTOT)CALL ERROR(N5-MTOT)                               00096050
  49.         CALL FRSTRS(A(N1),A(N2),A(N3),A(N3),A(N4),NF                    00096060
  50.      &  ,NSB,NFN,NEQB,NBLOCK,A(1))                                      00096070
  51.         RETURN                                                          00096080
  52. 999     FORMAT(A8)                                                      00096090
  53. 1000    FORMAT(I4,I1,2I5,3F10.0)                                        00096100
  54. 1010    FORMAT(39H1S I N U S O I D A L    R E S P O N S E,              00096110
  55.      &  19H    A N A L Y S I S//)                                       00096120
  56. 2000    FORMAT(15H BASE MOTION IN,A1,10H DIRECTION/                     00096130
  57.      &  31H NUMBER OF OUTPUT FREQUENCIES =,I5)                          00096140
  58. 2001    FORMAT(15H NOT IN SERVICE)                                      00096150
  59. 2010    FORMAT(46H ADDITIONAL PRINTING OUTPUT HAS BEEN REQUESTED)       00096160
  60.         END                                                             00096170
  61.       SUBROUTINE FEMID(ID,MASS,NUMNP,NEQB)                              00086470
  62.       IMPLICIT REAL*8(A-H,O-Z)                                          00086480
  63.       REAL*8  MASS                                                      00086490
  64.       REAL*8  ID                                                        00086500
  65.       COMMON/PREP/XMX,XAD,KSKIP,NDY ,I1,RRPREP(7)                       R0086510
  66.       DIMENSION ID(NUMNP,3),MASS(NEQB)                                  00086520
  67.       IWRITE=0                                                          00086530
  68.       REWIND 3                                                          00086540
  69.       REWIND 8                                                          00086550
  70.       READ (8) ID                                                       00086560
  71.       DO 100 L=1,NEQB                                                   00086570
  72.   100 MASS(L)=0.0D0                                                     00086580
  73.       NT=1                                                              00086590
  74.       DO 140 N=1,NUMNP                                                  00086600
  75.       DO 130 I=1,6                                                      00086610
  76.       NEQBS=NEQB*(NT-1)                                                 00086620
  77.       NEQBE=NEQB*NT                                                     00086630
  78.       CALL UNPKID(ID,NUMNP,W,WX,2,N,I)                                  00086640
  79.       NNN=W                                                             00086650
  80.       IF(NNN.LE.0) GO TO 130                                            00086660
  81.       IF(NNN.GT.NEQBS.AND.NNN.LE.NEQBE) GO TO 110                       00086670
  82.       IF(NNN.LE.NEQBS) GO TO 130                                        00086680
  83.       NT=NT+1                                                           00086690
  84.       DO 105 M=1,NEQB                                                   00086700
  85.   105 MASS(M)=0.0D0                                                     00086710
  86.   110 IF(I.GT.3) GO TO 120                                              00086720
  87.       L=NNN-(NT-1)*NEQB                                                 00086730
  88.       MASS(L)=I                                                         00086740
  89.   120 IF(NNN.EQ.NEQBE) WRITE(3) MASS                                    00086750
  90.       IF(NNN.EQ.NEQBE) IWRITE=IWRITE+1                                  00086760
  91.   130 CONTINUE                                                          00086770
  92.   140 CONTINUE                                                          00086780
  93.       IF(IWRITE.LT.NT) WRITE(3) MASS                                    00086790
  94.       RETURN                                                            00086800
  95.       END                                                               00086810
  96.       SUBROUTINE FRMPF(GAM,FR,F,XM,MASS,IDIR,NEQB                       00096180
  97.      &  ,NF,NBLOCK,NUMNP)                                               00096190
  98.         IMPLICIT REAL*8 (A-H,O-Z)                                       00096200
  99.         REAL*8 MASS                                                     00096210
  100.         DIMENSION GAM(1),FR(NF),F(NEQB,NF),XM(NEQB),MASS(NEQB)          00096220
  101.          L10RC=NEQB*NF*4                                                00096230
  102.         DO 100 I=1,NF                                                   00096240
  103. 100     GAM(I)=0.                                                       00096250
  104.         REWIND 9                                                        00096260
  105.         REWIND 3                                                        00096270
  106.         DO 130 N=1,NBLOCK                                               00096280
  107.         BACKSPACE 10                                                    00096290
  108.         READ (10) F                                                     00096300
  109.         BACKSPACE 10                                                    00096310
  110.         READ(3)MASS                                                     00096320
  111.         READ(9)XM                                                       00096330
  112.         DO 120 I=1,NEQB                                                 00096340
  113.         J=MASS(I)                                                       00096350
  114.         IF(J.NE.IDIR)GO TO 120                                          00096360
  115.         IF(J.LE.0)GO TO 120                                             00096370
  116.         DO 110 L=1,NF                                                   00096380
  117. 110     GAM(L)=GAM(L)+F(I,L)*XM(I)                                      00096390
  118. 120     CONTINUE                                                        00096400
  119. 130     CONTINUE                                                        00096410
  120.         BACKSPACE 10                                                    00096420
  121.         READ (10) FR                                                    00096430
  122.         WRITE (6,2000)                                                  00096440
  123.         TPI=6.2831853072D0                                              00096450
  124.         DO 230 I=1,NF                                                   00096460
  125.         FR(I)=FR(I)/TPI                                                 00096470
  126. 230     WRITE(6,2010)FR(I),GAM(I)                                       00096480
  127.         RETURN                                                          00096490
  128. 2000    FORMAT(11H0 FREQ, HZ.,5X,26HMODAL PARTICIPATION FACTOR)         00096500
  129. 2003    FORMAT(36H  BASE TRANSLATION, DIRECTION CODE =,I2)              00096510
  130. 2010    FORMAT(1X,1PE10.3,5X,E10.3)                                     00096520
  131.         END                                                             00096530
  132.       SUBROUTINE FRESP(GAM,FR,PHI,Q,FF,DAMP,NF,NFN,IPRINT)              00093300
  133.         IMPLICIT REAL*8(A-H,O-Z)                                        00093310
  134.         REAL HED                                                        00093320
  135.         COMMON /AMB/GRAV,REFT,JROT                                      R0093330
  136.         COMMON /PREP/XMX,XAD,MODEX,RRPREP(8)                            R0093340
  137.         DIMENSION GAM(NF),FR(NF),PHI(NF),Q(NF,NFN),HED(15),G(4)         00093350
  138.      &  ,F(8),FF(NFN),DAMP(NFN)                                         00093360
  139.         DO 10 I=1,NFN                                                   00093370
  140. 10      FF(I)=0.D0                                                      00093380
  141.         IF(NFN.EQ.NF)GO TO 300                                          00093390
  142.         READ(5,1000)ITYPE,HED                                           00093400
  143.         WRITE(6,2000)HED                                                00093410
  144.         GO TO (100,200),ITYPE                                           00093420
  145. 100     I=0                                                             00093430
  146.         J=0                                                             00093440
  147.         NF1=NFN-NF                                                      00093450
  148. 110     I=I+8                                                           00093460
  149.         READ(5,1010)F                                                   00093470
  150.         DO 120 N=1,8                                                    00093480
  151.         J=J+1                                                           00093490
  152.         IF(J.GT.NF1)GO TO 130                                           00093500
  153.         FF(J)=F(N)                                                      00093510
  154. 120     CONTINUE                                                        00093520
  155.         GO TO 110                                                       00093530
  156. 130     IF(MODEX.EQ.1)GO TO 250                                         00093540
  157.         DO 140 N=1,NF                                                   00093550
  158.         FF(J)=FR(N)                                                     00093560
  159.         J=J+1                                                           00093570
  160. 140     CONTINUE                                                        00093580
  161.         GO TO 250                                                       00093590
  162. 200     NF1=NFN-NF                                                      00093600
  163.         READ(5,1010)F1,F2                                               00093610
  164.         D=DLOG(F2/F1)/(NF1-1)                                           00093620
  165.         DO 210 I=1,NF1                                                  00093630
  166.         FF(I)=F1*DEXP(DBLE  (I-1)*D)                                    R0093640
  167. 210     CONTINUE                                                        00093650
  168.         J=NF1                                                           00093660
  169.         IF(MODEX.EQ.1)GO TO 250                                         00093670
  170.         DO 220 I=1,NF                                                   00093680
  171.         J=J+1                                                           00093690
  172.         FF(J)=FR(I)                                                     00093700
  173. 220     CONTINUE                                                        00093710
  174. 250     N2=NFN-1                                                        00093720
  175.         DO 270 I=1,N2                                                   00093730
  176.         I1=I+1                                                          00093740
  177.         DO 260 J=I1,NFN                                                 00093750
  178.         IF(FF(J).GE.FF(I))GO TO 260                                     00093760
  179.         SWAP=FF(J)                                                      00093770
  180.         FF(J)=FF(I)                                                     00093780
  181.         FF(I)=SWAP                                                      00093790
  182. 260     CONTINUE                                                        00093800
  183. 270     CONTINUE                                                        00093810
  184.         GO TO 400                                                       00093820
  185. 300     IF(MODEX.EQ.1)GO TO 400                                         00093830
  186.         DO 310 I=1,NFN                                                  00093840
  187. 310     FF(I)=FR(I)                                                     00093850
  188. 400     WRITE(6,2020)(I,FF(I),I=1,NFN)                                  00093860
  189.         IF(MODEX.EQ.1)WRITE(6,2180)                                     00093870
  190.         NRD=0                                                           00093880
  191. 405     READ(5,1000)NFMX,HED                                            00093890
  192.         WRITE(6,2030)NFMX,HED                                           00093900
  193.         IF(NFMX.LT.2)WRITE(6,2050)                                      00093910
  194.         READ(5,1010)(F(L),G(L),L=1,4)                                   00093920
  195.         WRITE(6,2040)(F(L),G(L),L=1,4)                                  00093930
  196.         IF(MODEX.EQ.0)GO TO 408                                         00093940
  197.         NRD=NRD+4                                                       00093950
  198.         IF(NRD.LT.NFMX)GO TO 405                                        00093960
  199.         GO TO 500                                                       00093970
  200. 408     FOLD=F(1)                                                       00093980
  201.         GOLD=G(1)                                                       00093990
  202.         I=1                                                             00094000
  203.         J=0                                                             00094010
  204.         K=4                                                             00094020
  205.         NFMX3=NFMX*3                                                    00094030
  206. 410     J=J+1                                                           00094040
  207.         GO TO 430                                                       00094050
  208. 420     I=I+1                                                           00094060
  209.         FOLD=F(I)                                                       00094070
  210.         GOLD=G(I)                                                       00094080
  211. 430     IF(FF(J).GE.FOLD)GO TO 440                                      00094090
  212.         WRITE(6,2050)                                                   00094100
  213.         STOP                                                            00094110
  214. 440     IF(I.LT.4)GO TO 460                                             00094120
  215.         K=K+4                                                           00094130
  216.         IF(K.LE.NFMX3)GO TO 450                                         00094140
  217.         WRITE(6,2050)                                                   00094150
  218.         STOP                                                            00094160
  219. 450     READ(5,1010)(F(L),G(L),L=1,4)                                   00094170
  220.         WRITE(6,2040)(F(L),G(L),L=1,4)                                  00094180
  221.         I=0                                                             00094190
  222. 460     IF(FF(J).GT.F(I+1))GO TO 420                                    00094200
  223.         DF=DLOG(F(I+1))-DLOG(FOLD)                                      00094210
  224.         DG=DLOG(G(I+1))-DLOG(GOLD)                                      00094220
  225.         IF(DF)470,420,480                                               00094230
  226. 470     WRITE(6,2060)                                                   00094240
  227.         STOP                                                            00094250
  228. 480     SLOPE=DG/DF                                                     00094260
  229.         GO=GRAV*(SLOPE*(DLOG(FF(J))-DLOG(FOLD))+GOLD)                   00094270
  230.         DO 485 L=1,NF                                                   00094280
  231. 485     Q(L,J)=GO                                                       00094290
  232.         IF(J.LT.NFN)GO TO 410                                           00094300
  233.         IF(IPRINT.EQ.1)                                                 00094310
  234.      &  WRITE(6,2070)(L,FF(L),Q(1,L),L=1,NFN)                           00094320
  235. 490     IF(K.GE.NFMX)GO TO 500                                          00094330
  236.         READ(5,1010)F                                                   00094340
  237.         K=K+4                                                           00094350
  238.         GO TO 490                                                       00094360
  239. 500     NRD=0                                                           00094370
  240. 505     READ(5,1000)NFMX,HED                                            00094380
  241.         WRITE(6,2080)NFMX,HED                                           00094390
  242.         IF(NFMX.LT.2)WRITE (6,2090)                                     00094400
  243.         READ(5,1010)(F(L),G(L),L=1,4)                                   00094410
  244.         WRITE(6,2040)(F(L),G(L),L=1,4)                                  00094420
  245.         IF(MODEX.EQ.0)GO TO 508                                         00094430
  246.         NRD=NRD+4                                                       00094440
  247.         IF(NRD.LT.NFMX)GO TO 505                                        00094450
  248.         GO TO 600                                                       00094460
  249. 508     FOLD=F(1)                                                       00094470
  250.         DOLD=G(1)                                                       00094480
  251.         I=1                                                             00094490
  252.         J=0                                                             00094500
  253.         K=4                                                             00094510
  254.         NFMX3=NFMX+3                                                    00094520
  255. 510     J=J+1                                                           00094530
  256.         GO TO 530                                                       00094540
  257. 520     I=I+1                                                           00094550
  258.         FOLD=F(I)                                                       00094560
  259.         DOLD=G(I)                                                       00094570
  260. 530     IF(FF(J).GE.FOLD)GO TO 540                                      00094580
  261.         WRITE(6,2090)                                                   00094590
  262.         STOP                                                            00094600
  263. 540     IF(I.LT.4)GO TO 560                                             00094610
  264.         K=K+4                                                           00094620
  265.         IF(K.LE.NFMX3)GO TO 550                                         00094630
  266.         WRITE(6,2090)                                                   00094640
  267.         STOP                                                            00094650
  268. 550     READ(5,1010)(F(L),G(L),L=1,4)                                   00094660
  269.         WRITE(6,2040)(F(L),G(L),L=1,4)                                  00094670
  270.         I=0                                                             00094680
  271. 560     IF(FF(J).GT.F(I+1))GO TO 520                                    00094690
  272.         DF=F(I+1)-FOLD                                                  00094700
  273.         DD=G(I+1)-DOLD                                                  00094710
  274.         IF(DF)570,520,580                                               00094720
  275. 570     WRITE(6,2060)                                                   00094730
  276.         STOP                                                            00094740
  277. 580     SLOPE=DD/DF                                                     00094750
  278.         DAMP(J)=(FF(J)-FOLD)*SLOPE+DOLD                                 00094760
  279.         IF(J.LT.NFN)GO TO 510                                           00094770
  280.         IF(IPRINT.EQ.1)                                                 00094780
  281.      &  WRITE(6,2110)(L,FF(L),DAMP(L),L=1,NFN)                          00094790
  282. 590     IF(K.GE.NFMX)GO TO 600                                          00094800
  283.         READ(5,1010)F                                                   00094810
  284.         K=K+4                                                           00094820
  285.         GO TO 590                                                       00094830
  286. 600     IF(MODEX.EQ.1)RETURN                                            00094840
  287.         IF(IPRINT.EQ.1)WRITE(6,2160)                                    00094850
  288.         DO 630 N=1,NFN                                                  00094860
  289.         F0=FF(N)                                                        00094870
  290.         XS=0.D0                                                         00094880
  291.         YS=0.D0                                                         00094890
  292.         DO 610 I=1,NF                                                   00094900
  293.         F1=FR(I)                                                        00094910
  294.       RNUM=2.0*DAMP(I)*F1*F0                                            00094920
  295.       RDEN=F1*F1-F0*F0                                                  00094930
  296.       IF(DABS(RNUM).LT.1.0E-10.AND.DABS(RDEN).LT.1.0E-10) GO TO 605     00094940
  297.         PHI(I)=DATAN2(2.*DAMP(I)*F1*F0,F1*F1-F0*F0)                     00094950
  298.         DEN=DSQRT((F1*F1-F0*F0)**2+(2.*DAMP(I)*F1*F0)**2)               00094960
  299.         Q(I,N)=Q(I,N)*GAM(I)/(DEN*39.4784176D0)                         00094970
  300.       GO TO 606                                                         00094980
  301.   605 WRITE(6,2190)                                                     00094990
  302.       Q(I,N)=1.0E10                                                     00095000
  303.       PHI(I)=0.0D0                                                      00095010
  304.   606 CONTINUE                                                          00095020
  305.         XS=XS+Q(I,N)*DCOS(PHI(I))                                       00095030
  306.         YS=YS+Q(I,N)*DSIN(PHI(I))                                       00095040
  307. 610     CONTINUE                                                        00095050
  308.         IF(IPRINT.EQ.1)                                                 00095060
  309.      &  WRITE(6,2150)F0,(I,PHI(I),Q(I,N),I=1,NF)                        00095070
  310.         PHI0=DATAN2(YS,XS)                                              00095080
  311.         IF(IPRINT.EQ.1)WRITE(6,2170)PHI0                                00095090
  312.         DO 620 I=1,NF                                                   00095100
  313. 620     Q(I,N)=Q(I,N)*DCOS(PHI(I)-PHI0)                                 00095110
  314. 630     CONTINUE                                                        00095120
  315.         IF(IPRINT.EQ.0)GO TO 650                                        00095130
  316.         WRITE(6,2120)                                                   00095140
  317.         DO 640 N=1,NFN                                                  00095150
  318. 640     WRITE(6,2130)FF(N),(Q(I,N),I=1,NF)                              00095160
  319. 650     REWIND 4                                                        00095170
  320.         WRITE(4)Q,FF                                                    00095180
  321. 1000    FORMAT(I5,5X,15A4)                                              00095190
  322. 1010    FORMAT(8F10.0)                                                  00095200
  323. 1020    FORMAT(2F10.0,I5)                                               00095210
  324. 2000    FORMAT(1X ,15A4)                                                00095220
  325. 2010    FORMAT(32H0FREQUENCY OUTPUT REQUESTED FROM,E10.4,7H HZ. TO      00095230
  326.      &  ,E10.4,4H HZ./5H WITH,I5,15H LOG INCREMENTS)                    00095240
  327. 2020    FORMAT(33H0REQUESTED FREQUENCIES FOR OUTPUT/                    00095250
  328.      &18H NO.    FREQ.(HZ.)/(I5,1PE12.4))                               R0095260
  329. 2030    FORMAT(34H0TABLE OF LOG(G) VS LOG(FREQ) WITH,I5,7H POINTS/      00095270
  330.      &  1X,15A4/22H FREQ, HZ.      PEAK G)                              00095280
  331. 2040    FORMAT((1PE10.3,3X,E10.3))                                      00095290
  332. 2050    FORMAT(47H0***ERROR*** NOT ENOUGH DATA POINTS FOR G LEVEL,      00095300
  333.      &  6H TABLE)                                                       00095310
  334. 2060    FORMAT(40H0***ERROR*** FREQUENCY DATA OUT OF ORDER)             00095320
  335. 2070    FORMAT(45H0FINAL ACCELERATION APPLIED AT EACH FREQUENCY/        00095330
  336.      &  34H  NO.     FREQ(HZ)    ACCELERATION/(I5,1P2E13.3))            00095340
  337. 2080    FORMAT(30H0TABLE OF DAMPING VS FREQ WITH,I5,7H POINTS/          00095350
  338.      &  1X,15A4/26H FREQ(HZ)    DAMPING RATIO)                          00095360
  339. 2090    FORMAT(47H0***ERROR*** NOT ENOUGH DATA POINTS FOR DAMPING,      00095370
  340.      &  6H TABLE)                                                       00095380
  341. 2110    FORMAT(39H0FINAL DAMPING RATIO FOR EACH FREQUENCY/              00095390
  342.      &  35H  NO.     FREQ(HZ)    DAMPING RATIO/(I5,1PE13.3,0PF13.5))    00095400
  343. 2120    FORMAT(43H0FINAL AMPLITUDE (WITH RESPECT TO REFERENCE,          00095410
  344.      &  13H PHASE ANGLE),37H OF EACH MODE FOR EACH DESIRED OUTPUT,      00095420
  345.      &  10H FREQUENCY//8H APPLIED,12X,27HAMPLITUDES FOR MODE NUMBERS/   00095430
  346.      &  6H  FREQ,12X,1H1,12X,1H2,12X,1H3,12X,1H4,12X,6H5.....)          00095440
  347. 2130    FORMAT(1P10E13.4/(13X,9E13.4))                                  00095450
  348. 2140    FORMAT(36H0***ERROR*** IMPROPER OPTION, ITYPE=,I12)             00095460
  349. 2150    FORMAT(19H0APPLIED FREQUENCY=,E10.3,3HHZ./                      00095470
  350.      &  (25H   MODE ANGLE   AMPLITUDE)/5(13H    NO  (RAD),12X)/         00095480
  351.      &  (1X,5(I5,F7.3,E13.4)))                                          00095490
  352. 2160    FORMAT(43H1TABLE OF RESPONSE (PHASE AND AMPLITUDE) OF,          00095500
  353.      &  10H EACH MODE,26H AT EACH APPLIED FREQUENCY)                    00095510
  354. 2170    FORMAT(23H0REFERENCE PHASE ANGLE=,F6.3,8H RADIANS)              00095520
  355. 2180    FORMAT(44H0DURING THE DATA CHECK,THE STRUCTURE NATURAL,         00095530
  356.      &  29H FREQUENCIES ARE SET TO ZERO.)                               00095540
  357.  2190 FORMAT(5X,39HRESONANCE IS ACHIEVED, DAMPING IS ZERO.,             00095550
  358.      12X,       48HAMPLITUDE IS SET TO AN ARBITRARY VALUE OF 1.0E10)    00095560
  359.         RETURN                                                          00095570
  360.         END                                                             00095580
  361.       SUBROUTINE FRDSPR(ID,F,FI,X,FF,NEQB,NF,NDS,NUMNP,NBLOCK,NSB,ISIR) 00092390
  362.       IMPLICIT REAL*8(A-H,O-Z)                                          00092400
  363.       REAL*8  ID                                                        00092410
  364.       INTEGER*2 ISIR                                                    00092420
  365.       DIMENSION ID(NUMNP,3),F(8,NF),FI(NSB ,NF),X(NF,NDS),FF(NDS),NUM(1)R0092430
  366.       COMMON/JUNK/JJ, NP,IC(6),D(8),L,  II,MSB,NS,NE,N,DDT,TIME,        00092440
  367.      $                M,J,K,MM,KD(3,8),DD,XUM,IEQ,NRD                   00092450
  368.      $                ,DM(8),TM(8),NRJUNK(354)                          R0092460
  369.       COMMON / DYN / NT,NOT,DAMP,DT,RRDYN(3)                            R0092470
  370.       COMMON/PREP/XMX,XAD,KSKIP,NDYN,I1,RRPREP(7)                       R0092480
  371.       COMMON /BAND/ KOPT,NRBAND(7)                                      R0092490
  372.       COMMON /SSIT/ NV                                                  00092500
  373.       DIMENSION ISIR(NUMNP)                                             00092510
  374.       LBKS=4*NEQB*NF                                                    00092520
  375.       LRDS=4*NEQB*NV                                                    00092530
  376.       DT=DT*NOT                                                         00092540
  377.       READ (5,230) KKK,ISP                                              00092550
  378.       IF(KKK.EQ.0) RETURN                                               00092560
  379.       REWIND 9                                                          00092570
  380.       REWIND 8                                                          00092580
  381.       READ (8) ID                                                       00092590
  382.       NT1=17                                                            00092600
  383.       IF(KOPT.GT.0) REWIND NT1                                          00092610
  384.       IF(KSKIP.EQ.1) GO TO 5                                            00092620
  385.       IF(KOPT.GT.0) READ (NT1)                                          00092630
  386.       IF(KOPT.GT.0) READ (NT1) ISIR                                     00092640
  387.     5 CONTINUE                                                          00092650
  388.       L=0                                                               00092660
  389.       NUM(1) = 0                                                        R0092670
  390.       WRITE (6,220)                                                     00092680
  391.   100 READ (5,230) NP,IC                                                00092690
  392.       IF(NP.GT.0) WRITE(6,240)NP,IC                                     00092700
  393.       IF(NP.GT.0) GO TO 120                                             00092710
  394.       IF(L.EQ.0) GO TO 160                                              00092720
  395.       IF (KSKIP.EQ.1) GO TO 110                                         00092730
  396.       WRITE (9) KD,L                                                    00092740
  397.   110 CONTINUE                                                          00092750
  398.       NUM(1) = NUM(1) + 1                                               R0092760
  399.       GO TO 160                                                         00092770
  400.   120 DO 150 I=1,6                                                      00092780
  401.       II=IC(I)                                                          00092790
  402.       IF(II.EQ.0) GO TO 100                                             00092800
  403.   130 L=L+1                                                             00092810
  404.       KD(1,L)=NP                                                        00092820
  405.       N=NP                                                              00092830
  406.       IF(KOPT.GT.0)   NP=ISIR(N)                                        00092840
  407.       KD(2,L)=II                                                        00092850
  408.       CALL UNPKID (   ID,NUMNP,W      ,WX      ,2,NP,II)                00092860
  409.       NP=N                                                              00092870
  410.       NNN=W                                                             00092880
  411.       KD(3,L)=NNN                                                       00092890
  412.       IF(NNN.LE.0) L=L-1                                                00092900
  413.       IF(L.LT.8) GO TO 150                                              00092910
  414.       IF (KSKIP.EQ.1) GO TO 140                                         00092920
  415.       WRITE (9) KD,L                                                    00092930
  416.   140 CONTINUE                                                          00092940
  417.       NUM(1) = NUM(1) + 1                                               R0092950
  418.       L=0                                                               00092960
  419.   150 CONTINUE                                                          00092970
  420.       GO TO 100                                                         00092980
  421.   160 IF(NUM(1).EQ.0) RETURN                                            R0092990
  422.       WRITE (6,250) KKK,ISP                                             00093000
  423.       IF (KSKIP.EQ.1)RETURN                                             00093010
  424.       REWIND 3                                                          00093020
  425.       REWIND 9                                                          00093030
  426.       REWIND 10                                                         00093040
  427.       READ (10)                                                         00093050
  428.       NE=NSB                                                            00093060
  429.       NS=NE+1-NEQB                                                      00093070
  430.       DO 170 I=1,NBLOCK                                                 00093080
  431.       READ (10)((FI(J,K),J=NS,NE),K=1,NF)                               00093090
  432.       NS=NS-NEQB                                                        00093100
  433.   170 NE=NE-NEQB                                                        00093110
  434.       NUMRR = NUM(1)                                                            
  435.       DO 190 N=1,NUMRR                                                  R0093120
  436.       READ (9) KD,L                                                     00093130
  437.       DO 180 I=1,L                                                      00093140
  438.       II=KD(3,I)                                                        00093150
  439.       DO 180 J=1,NF                                                     00093160
  440.   180 F(I,J)=FI(II,J)                                                   00093170
  441.   190 WRITE (3) L,KD,F                                                  00093180
  442.   200 CONTINUE                                                          00093190
  443.       CALL FRDSPL (X,F,FF,NF,NDS,NUM,1,KKK,2,ISP)                       00093200
  444.   210 RETURN                                                            00093210
  445.   220 FORMAT (35H1DISPLACEMENT COMPONENTS FOR WHICH ,                   00093220
  446.      $   19H OUTPUT IS REQUIRED  //                                     00093230
  447.      $ 31H  NODE  DISPLACEMENT COMPONENTS  /)                           00093240
  448.   230 FORMAT (7I5)                                                      00093250
  449.   240 FORMAT (I5,4X,6I3)                                                00093260
  450.   250 FORMAT (/16H OUTPUT TYPE....,I1/                                  00093270
  451.      $         16H PLOT SPACING...,I1)                                  00093280
  452.       END                                                               00093290
  453.       SUBROUTINE FRDSPL (X,F,FF,NF,NDS,NUM,NN,KKK,ISD,ISP)              00091450
  454.       IMPLICIT REAL*8(A-H,O-Z)                                          00091460
  455.       REAL*8  NPAR                                                      00091470
  456.       DIMENSION FF(NDS),X(NF,NDS),F(8,NF),NUM(NN)                       00091480
  457.       COMMON / JUNK /TM(8),DM(8),D(8),KD(3,8),RRJUNK(191)               R0091490
  458.       COMMON / DYN / NT,NOT,DAMP,DT,RRDYN(3)                            R0091500
  459.       COMMON / ELPAR / NPAR(14),RRELPA(29)                              R0091510
  460.       REWIND 3                                                          00091520
  461.       REWIND 4                                                          00091530
  462.       READ (4) X,FF                                                     00091540
  463.       DO 320 N=1,NN                                                     00091550
  464.       REWIND 2                                                          00091560
  465.       REWIND 9                                                          00091570
  466.       MM=NUM(N)                                                         00091580
  467.   100 IF(MM.EQ.0) GO TO 320                                             00091590
  468.       DO 270 M=1,MM                                                     00091600
  469.       READ (3) L,KD,F                                                   00091610
  470.       GO TO (110,160,140),KKK                                           00091620
  471.   110 IF(ISD.EQ.1) GO TO 120                                            00091630
  472.       WRITE (6,330) M                                                   00091640
  473.       GO TO 130                                                         00091650
  474.   120 WRITE (6,390) M                                                   00091660
  475.   130 WRITE (6,400) (KD(1,I),KD(2,I),I=1,L)                             00091670
  476.       GO TO 160                                                         00091680
  477.   140 IF(M.GT.1) GO TO 160                                              00091690
  478.       IF(ISD.EQ.1) GO TO 150                                            00091700
  479.       WRITE (6,340)                                                     00091710
  480.       WRITE (6,450)                                                     00091720
  481.       GO TO 160                                                         00091730
  482.   150 WRITE (6,410)                                                     00091740
  483.       WRITE (6,430)                                                     00091750
  484.   160 DO 170 I=1,L                                                      00091760
  485.       TM(I)=0.                                                          00091770
  486.   170 DM(I)=0.                                                          00091780
  487.       DO 230 K=1,NDS                                                    00091790
  488.       TIME=FF(K)                                                        00091800
  489.       DO 200 I=1,L                                                      00091810
  490.       DD=0.                                                             00091820
  491.       DO 180 J=1,NF                                                     00091830
  492.   180 DD = DD + F(I,J)*X(J,K)                                           00091840
  493.       AD= DABS(DD)                                                      00091850
  494.       IF(AD-DM(I)) 200,200,190                                          00091860
  495.   190 DM(I)=AD                                                          00091870
  496.       TM(I)=TIME                                                        00091880
  497.   200 D(I)=DD                                                           00091890
  498.       GO TO (210,220,230),KKK                                           00091900
  499.   210 WRITE (6,350) TIME,(D(I),I=1,L)                                   00091910
  500.       GO TO 230                                                         00091920
  501.   220 WRITE (9) D                                                       00091930
  502.   230 CONTINUE                                                          00091940
  503.       GO TO (240,250,260),KKK                                           00091950
  504.   240 WRITE (6,360) (DM(I),I=1,L)                                       00091960
  505.       WRITE (6,370) (TM(I),I=1,L)                                       00091970
  506.       GO TO 270                                                         00091980
  507.   250 WRITE (2) KD,DM,TM,L                                              00091990
  508.       GO TO 270                                                         00092000
  509.   260 WRITE (6,380) (KD(1,I),KD(2,I),DM(I),TM(I),I=1,L)                 00092010
  510.   270 CONTINUE                                                          00092020
  511.       IF(KKK.NE.2) GO TO 320                                            00092030
  512.       REWIND 2                                                          00092040
  513.       REWIND 9                                                          00092050
  514.       DO 310 M=1,MM                                                     00092060
  515.       GO TO (280,290),ISD                                               00092070
  516.   280 WRITE (6,420) M                                                   00092080
  517.       WRITE (6,430)                                                     00092090
  518.       GO TO 300                                                         00092100
  519.   290 WRITE (6,440) M                                                   00092110
  520.       WRITE (6,450)                                                     00092120
  521.   300 CALL FPLOTD (2,9,NDS,ISP)                                         00092130
  522.   310 CONTINUE                                                          00092140
  523.   320 CONTINUE                                                          00092150
  524.       RETURN                                                            00092160
  525.   330 FORMAT (1X ,      43H NUMBER OF SELECTED DISPLACEMENT COMPONENTS, 00092170
  526.      $ 5H..... ,I3//10X,40HNODE NUMBERS AND DISPLACEMENT COMPONENTS )   00092180
  527. 340     FORMAT(1X ,35H MAXIMUM DISPLACEMENT VALUES FROM A,              00092190
  528.      &42H STEADY STATE SINUSOIDAL RESPONSE ANALYSIS//)                  00092200
  529.   350 FORMAT (F8.3,2X,1P8E12.3)                                         00092210
  530.   360 FORMAT ( /24H MAXIMUM ABSOLUTE VALUES //                          00092220
  531.      $ 10H MAXIMUM   ,1P8E12.3)                                         00092230
  532.   370 FORMAT (10H FREQ      ,1P8E12.3)                                  00092240
  533.   380 FORMAT (I6,I13,1PE18.3,E12.3,5X,2HNA)                             00092250
  534.   390 FORMAT (1X ,       37H NUMBER OF SELECTED STRESS COMPONENTS,      00092260
  535.      $ 5H..... , I3//10X,37H ELEMENT AND STRESS COMPONENT NUMBERS )     00092270
  536.   400 FORMAT (8H    FREQ , 2X, 8(I8,1H-,I2,1X))                         00092280
  537.   410 FORMAT (                                                          00092290
  538.      $1X ,58H MAXIMUM STRESS VALUES FROM SINUSOIDAL RESPONSE ANALYSIS//)00092300
  539.   420 FORMAT (                                                          00092310
  540.      $ 39H1NORMALISED PLOT OF STRESS HISTORIES...,I3/)                  00092320
  541.   430 FORMAT(58H  ELEMENT      STRESS       MAXIMUM      FREQ AT     PLO00092330
  542.      $T / 58H  NUMBER     COMPONENT       VALUE       MAXIMUM    SYMBOL)00092340
  543.   440 FORMAT (46H1NORMALISED PLOT OF DISPLACEMENT HISTORIES....,I3/)    00092350
  544.   450 FORMAT(58H   NODE     DISPLACEMENT    MAXIMUM      FREQ AT     PLO00092360
  545.      $T / 58H  NUMBER     COMPONENT       VALUE       MAXIMUM    SYMBOL)00092370
  546.       END                                                               00092380
  547.       SUBROUTINE FPLOTD (IT,JT,NDS,ISP)                                 00090550
  548.       IMPLICIT REAL*8 (A-H,O-Z)                                         00090560
  549.       COMMON /QTSARG/ PP(101),KD(3,8),XM(8),TM(8),IP(8),X(8),RRQTSA(859)R0090570
  550.       COMMON /DYN/    NT,NOT,DAMP,DT,RRDYN(3)                           R0090580
  551.       DIMENSION SM(8)                                                   00090590
  552.       DATA SM /1H1,1H2,1H3,1H4,1H5,1H6,1H7,1H8 /                        00090600
  553.       DATA BL /1H /,V /1HX/,AST /1H*/                                   00090610
  554.       LINE=53                                                           00090620
  555.       READ (IT) KD,XM,TM,L                                              00090630
  556.       WRITE (6,260) (KD(1,I),KD(2,I),XM(I),TM(I),I,I=1,L)               00090640
  557.       DO 100 K=1,L                                                      00090650
  558.       TT = XM(K)                                                        00090660
  559.       IF(  DABS(TT).GT.1.0D-8 ) XM(K) = 50.0D0/ TT                      00090670
  560.   100 CONTINUE                                                          00090680
  561.       TT=0.E0                                                           00090690
  562.       WRITE (6,220)                                                     00090700
  563.   110 FORMAT(/4X,47HNOTE - AN * INDICATES MORE THAN ONE PLOT SYMBOL)    00090710
  564.       WRITE (6,230)                                                     00090720
  565.       WRITE (6,240) TT,(V,I=1,101),TT                                   00090730
  566.       LLCT=9                                                            00090740
  567.       K=1                                                               00090750
  568.       DO 120 I=2,100                                                    00090760
  569.   120 PP(I)=BL                                                          00090770
  570.       DO 210 N=1,NDS                                                    00090780
  571.       READ (JT) X                                                       00090790
  572.       PP(1)=V                                                           00090800
  573.       PP(51)=V                                                          00090810
  574.       PP(101)=V                                                         00090820
  575.   130 II=ISP                                                            00090830
  576.   140 IF(II.LE.0) GO TO 150                                             00090840
  577.       WRITE (6,250) PP                                                  00090850
  578.       LLCT=LLCT+1                                                       00090860
  579.       IF(LLCT.EQ.LINE) WRITE(6,250) (V,I=1,101)                         00090870
  580.       IF(LLCT.EQ.LINE) WRITE(6,230)                                     00090880
  581.       IF(LLCT.EQ.LINE) WRITE(6,110)                                     00090890
  582.       IF(LLCT.EQ.LINE) WRITE(6,220)                                     00090900
  583.       IF(LLCT.EQ.LINE) WRITE(6,230)                                     00090910
  584.       IF(LLCT.EQ.LINE) WRITE(6,250) (V,I=1,101)                         00090920
  585.       IF(LLCT.EQ.LINE) LLCT=9                                           00090930
  586.       II=II-1                                                           00090940
  587.       GO TO 140                                                         00090950
  588.   150 TT=TT+DT                                                          00090960
  589.       DO 170 I=1,L                                                      00090970
  590.       XX=XM(I)*X(I)                                                     00090980
  591.       M=XX                                                              00090990
  592.       M=M+51                                                            00091000
  593.       IP(I)=M                                                           00091010
  594.       IF(PP(M).EQ.V .OR. PP(M).EQ.BL) GO TO 160                         00091020
  595.       PP(M) = AST                                                       00091030
  596.       GO TO 170                                                         00091040
  597.   160 PP(M) = SM(I)                                                     00091050
  598.   170 CONTINUE                                                          00091060
  599.       IF(K.LT.10)  GO TO 180                                            00091070
  600.       K=1                                                               00091080
  601.       WRITE (6,240) TT,PP,TT                                            00091090
  602.       LLCT=LLCT+1                                                       00091100
  603.       IF(LLCT.EQ.LINE) WRITE(6,250) (V,I=1,101)                         00091110
  604.       IF(LLCT.EQ.LINE) WRITE(6,230)                                     00091120
  605.       IF(LLCT.EQ.LINE) WRITE(6,110)                                     00091130
  606.       IF(LLCT.EQ.LINE) WRITE(6,220)                                     00091140
  607.       IF(LLCT.EQ.LINE) WRITE(6,230)                                     00091150
  608.       IF(LLCT.EQ.LINE) WRITE(6,250) (V,I=1,101)                         00091160
  609.       IF(LLCT.EQ.LINE) LLCT=9                                           00091170
  610.       GO TO 190                                                         00091180
  611.   180 WRITE (6,250) PP                                                  00091190
  612.       LLCT=LLCT+1                                                       00091200
  613.       IF(LLCT.EQ.LINE) WRITE(6,250) (V,I=1,101)                         00091210
  614.       IF(LLCT.EQ.LINE) WRITE(6,230)                                     00091220
  615.       IF(LLCT.EQ.LINE) WRITE(6,110)                                     00091230
  616.       IF(LLCT.EQ.LINE) WRITE(6,220)                                     00091240
  617.       IF(LLCT.EQ.LINE) WRITE(6,230)                                     00091250
  618.       IF(LLCT.EQ.LINE) WRITE(6,250) (V,I=1,101)                         00091260
  619.       IF(LLCT.EQ.LINE) LLCT=9                                           00091270
  620.       K=K+1                                                             00091280
  621.   190 DO 200 I=1,L                                                      00091290
  622.       M=IP(I)                                                           00091300
  623.   200 PP(M)=BL                                                          00091310
  624.   210 CONTINUE                                                          00091320
  625.       TT=TT+DT                                                          00091330
  626.       WRITE (6,240) TT,(V,I=1,101),TT                                   00091340
  627.       WRITE (6,230)                                                     00091350
  628.       WRITE(6,110)                                                      00091360
  629.       RETURN                                                            00091370
  630.   220 FORMAT (1X ,57X,15HO R D I N A T E )                              00091380
  631.   230 FORMAT ( / 1H ,3X,7HT I M E,2X,4H-1.0,21X,4H-0.5,22X,3H0.0,22X,   00091390
  632.      $        3H0.5,22X,3H1.0,4X,7HT I M E, 1X)                         00091400
  633.   240 FORMAT (1H ,F10.4,4X,101A1,F12.4)                                 00091410
  634.   250 FORMAT (1H ,14X,101A1)                                            00091420
  635.   260 FORMAT (I8,12X,I3,1P2E14.4,3X,I6)                                 00091430
  636.       END                                                               00091440
  637.       SUBROUTINE FRSTRS (T  ,SF,FI,X,FF,NF,NSB,NDS,NEQB,NBLOCK,SA)      00096540
  638.       IMPLICIT REAL*8(A-H,O-Z)                                          00096550
  639.       REAL*8  NPAR                                                      00096560
  640.       DIMENSION NUM(1),SF(8,NF),FI(NSB,NF),X(NF,NDS),FF(NDS)            R0096570
  641.       DIMENSION SA(1)                                                   00096580
  642.       COMMON /ELPAR/ NPAR(14),NUMNP,MBAND,NELTYP,N1,N2,N3,N4,N5,MTOT,NEQ00096590
  643.      $ ,N2P,N3P,NMRI,NTRI,N1P,NUMEL,NUMEL2,KZ(20,2),NEMN                00096600
  644.       COMMON/PREP/XMX,XAD,KSKIP,NDYN,I1,RRPREP(7)                       R0096610
  645.       COMMON / JUNK / N,NEL,IS(12),M,I,L,KS(3,8),II,K,J,SS,JJ           00096620
  646.      $                ,NUME,NE,NRJUNK(405)                              R0096630
  647.       NT1=1                                                             00096640
  648.       READ (5,190) KKK,ISP                                              00096650
  649.       IF(KKK.EQ.0) RETURN                                               00096660
  650.       N=1                                                               00096670
  651.       IF (KSKIP.EQ.1) GO TO 105                                         00096680
  652.       REWIND 10                                                         00096690
  653.       READ (10)                                                         00096700
  654.       NE=NSB                                                            00096710
  655.       NS1=NE+1-NEQB                                                     00096720
  656.       DO 100 I=1,NBLOCK                                                 00096730
  657.       READ (10)((FI(J,K),J=NS1,NE),K=1,NF)                              00096740
  658.       NS1=NS1-NEQB                                                      00096750
  659.   100 NE=NE-NEQB                                                        00096760
  660.   105 CONTINUE                                                          00096770
  661.       CALL RDWRT(NT1,SA,1,6,JK)                                         00096780
  662.       REWIND 3                                                          00096790
  663.       WRITE (6,210)                                                     00096800
  664.       WRITE (6,220)                                                     00096810
  665.       READ (5,190) NEL,IS                                               00096820
  666.       WRITE (6,200) NEL,IS                                              00096830
  667.       NUME=NUMEL+NUMEL2                                                 00096840
  668.       L=0                                                               00096850
  669.       NUM(N)=0                                                          00096860
  670.       DO 170 M=1,NUME                                                   00096870
  671.       IF (KSKIP.EQ.1) GO TO 110                                         00096880
  672.       CALL RDWRT(NT1,SA,NEMN,0,KOUNT)                                   00096890
  673.   110 CONTINUE                                                          00096900
  674.       NS1=SA(KOUNT-1)                                                   00096910
  675.       ND1=SA(KOUNT-2)                                                   00096920
  676.       IF(NEL.NE.M) GO TO 170                                            00096930
  677.       IF(KSKIP.EQ.1) GO TO 155                                          00096940
  678.       DO 150 I=1,NS1                                                    00096950
  679.       II=IS(I)                                                          00096960
  680.       IF(II.EQ.0) GO TO 160                                             00096970
  681.       NPN=ND1+II                                                        00096980
  682.       L=L+1                                                             00096990
  683.       KS(1,L)=NEL                                                       00097000
  684.       KS(2,L)=II                                                        00097010
  685.       DO 140 K=1,NF                                                     00097020
  686.       SS=0.                                                             00097030
  687.       DO 130 J=1,ND1                                                    00097040
  688.       NELM=NPN+(J-1)*NS1                                                00097050
  689.       JJ=SA(J)                                                          00097060
  690.       IF(JJ) 130,130,120                                                00097070
  691.   120 SS = SS + SA(NELM)*FI(JJ,K)                                       00097080
  692.   130 CONTINUE                                                          00097090
  693.   140 SF(L,K)=SS                                                        00097100
  694.       IF(L.LT.8) GO TO 150                                              00097110
  695.       WRITE (3) L,KS,SF                                                 00097120
  696.       L=0                                                               00097130
  697.       NUM(N)=NUM(N) + 1                                                 00097140
  698.   150 CONTINUE                                                          00097150
  699.   155 CONTINUE                                                          00097160
  700.   160 READ  (5,190) NEL,IS                                              00097170
  701.       IF(NEL.GT.0) WRITE(6,200)NEL,IS                                   00097180
  702.   170 CONTINUE                                                          00097190
  703.       IF(L.EQ.0) GO TO 180                                              00097200
  704.       IF (KSKIP.EQ.1) GO TO 180                                         00097210
  705.       WRITE (3) L,KS,SF                                                 00097220
  706.       NUM(N)=NUM(N) + 1                                                 00097230
  707.   180 CONTINUE                                                          00097240
  708.       WRITE (6,230) KKK,ISP                                             00097250
  709.       IF (KSKIP.EQ.1)RETURN                                             00097260
  710.       NELTYP=1                                                          00097270
  711.       CALL FRDSPL (X,SF,FF,NF,NDS,NUM,NELTYP,KKK,1,ISP)                 00097280
  712.       RETURN                                                            00097290
  713.   190 FORMAT (13I5)                                                     00097300
  714.   200 FORMAT (I6,4X,12I3)                                               00097310
  715. 210     FORMAT(27H1STRESS COMPONENTS REQUIRED)                          00097320
  716.   220 FORMAT (/16H ELEMENT TYPE .. ,F4.2//                              00097330
  717.      $ 38H ELEMENT     DESIRED STRESS COMPONENTS )                      00097340
  718.   230 FORMAT (/16H OUTPUT TYPE....,I1/                                  00097350
  719.      $         16H PLOT SPACING...,I1)                                  00097360
  720.       END                                                               00097370
  721.