Trailing-Edge
-
PDP-10 Archives
-
-
There are no other files named in the archive.
00100 DIMENSION NP(10),XT(12),YT(12),PT(12),CHAR(10),
00200 . LT(10),IX(4),IY(10),YS(2000),IDLIN(60)
00300 DIMENSION VARN(4,10),IVARN(4),VARNM(10)
00400 EQUIVALENCE(VARN(1,1),VARNM(1)),(IVARN(1),IVARNM)
00500 DIMENSION LNLBL(4,10),IWHX(10),IPLTN(4),SPX(2),SPY(2)
00600 COMMON/XOUT/XZ,XQ,PXP(2000)
00700 COMMON/YOUT/YZ,YQ,PYP(2000)
00800 COMMON/GRID/XGR,YGR
00900 INTEGER XT,YT,PT,CHAR,PNUM
01000 LOGICAL TS
01100 DIMENSION XBUF(15),XS(800)
01200 DATA CHAR/10*' '/
01250 DIMENSION IOLD(2),INEW(2)
01300 DATA IOLD/' ',' '/
01400 INTEGER VARNM,VARN
01500 DATA SX,SY/.5,.5/
01600 DIMENSION XD(1)
01700 EQUIVALENCE(HMIN,IHMIN),(HMAX,IHMAX)
01800 C
01900 C OPEN PLOT FILE
02000 CALL HEREIS(XD)
02100 ASSIGN 994 TO KRET
02200 TYPE 100
02300 100 FORMAT(' NAME OF PLOTTER FILE :'$)
02400 ACCEPT 200,IPLT
02500 200 FORMAT(A5,2I)
02600 CALL PLOTS(II,IPLT)
02700 C
02800 C SET SOME CONSTANTS
02900 PNUM=0
03000 2 KS=0
03100 KKS=1
03200 YMIN=9999999999.
03300 YMAX=-999999999.
03400 XMIN=99999999999.
03500 XMAX=-999999999999.
03600 NSINES=0
03700 TS=.FALSE.
03800 5 PNUM=PNUM+1
03900 6 K=KS
04000 KK=KKS
04100 NLINES=NSINES
04200 TYPE 101,PNUM
04300 101 FORMAT(' PLOT ',I1/)
04400 C
04500 C DETERMINE DATA SOURCE
04600 7 TYPE 102
04700 102 FORMAT(' DATA SOURCE - P(LANETS) OR O(THER) :'$)
04800 ACCEPT 202,IPLT
04900 202 FORMAT(A1)
05000 C
05100 C IPT=1 IF PLANETS, IPT=2 IF FILEN, IPT=3 IF TTY
05200 IPTN=1
05300 IF (IPLT .NE. 'P') GO TO 20
05400 10 TYPE 103
05500 103 FORMAT(' PLANETS DATA BANK NAME :'$)
05600 ACCEPT 200,II
05700 CALL SYMCK1(II,$30,$6,$7)
05800 REREAD 203,IPLTN
05900 203 FORMAT(4A4)
06000 11 TYPE 1030
06100 1030 FORMAT(' PPN :'$)
06200 ACCEPT 212,IPR,IPG
06300 IPP=IPR*262144+IPG
06500 GO TO 22
06600 20 TYPE 104
06700 IPTN=2
06800 104 FORMAT(' TTY OR FILENAME :'$)
06900 ACCEPT 200,II
07000 CALL SYMCK1(II,$30,$6,$7)
07100 REREAD 204,INEW
07150 204 FORMAT(2A5)
07200 IF (INEW(1) .EQ. 'TTY') IPTN=IPTN+1
07350 IF (IPTN .EQ. 2) GO TO 11
07375 22 IPT=IPTN
07400 C
07500 C DETERMINE TIME SERIES OR REGULAR
07600 30 TYPE 108
07700 108 FORMAT(3X,'TIME SERIES PLOT (Y OR N):'$)
07800 ACCEPT 200,IANS
07900 CALL SYMCK1(IANS,$55,$6,$30)
08000 TS=IANS .EQ. 'Y'
08050 55 K=K+1
08100 IF (TS) GO TO 31
08200 GO TO (32,61,61),IPT
08300 C
08400 C GET PLOT PERIOD
08500 31 TYPE 105
08600 105 FORMAT(' PLOT PERIOD :'$)
08700 ACCEPT 200,II
08800 CALL SYMCK1(II,$32,$6,$10)
08900 REREAD 205,IFPL,ILPL
09000 205 FORMAT(3I)
09100 C
09200 C IF NOT PLANETS GET FREQUENCY
09300 IF (IPT .EQ. 1) GO TO 32
09400 TYPE 1050
09500 1050 FORMAT(' DATA FREQUENCY (1,4,6, OR 12) :'$)
09600 ACCEPT 200,II
09700 CALL SYMCK1(II,$32,$6,$10)
09800 REREAD 205,NFREQ
09900 NPRS=KINDEX(ILPL,IFPL,NFREQ)
09910 TYPE 1051,NPRS
09920 1051 FORMAT(' THE PLOT PERIOD REPRES. ',I5,' OBS. '/)
10000 C
10100 C GET # OF DEP VARS THIS X.
10200 32 TYPE 106
10300 106 FORMAT(' NUMBER OF DEPENDENT VARS :'$)
10400 ACCEPT 200,II
10500 CALL SYMCK1(II,$33,$6,$31)
10600 REREAD 205,NL
10700 33 KFIN=NL+K-1
10800 IP=IFPL
10850 IFF=ILPL
10900 NLINES=NLINES+NL
11000 IF (IPT .GT. 1) GO TO 620
11100 C
11200 C PLANETS - CHECK IF DB IS THERE
11300 CALL PLBANK(IPLTN,IPP,NFREQ,IERR)
11400 IF (IERR .GT. 1) GO TO 5000
11500 34 DO 51 I=K,KFIN
11600 IWHX(I)=KK
11800 J=I-K+1
11900 351 TYPE 107,J
12000 107 FORMAT(3X,'NAME OF DEP. VAR ',I1,':'$)
12100 ACCEPT 200,II
12200 CALL SYMCK1(II,$51,$6,$35)
12300 REREAD 203,(VARN(J,I),J=1,4)
12400 C
12500 C READ IN PLANETS VARS.
12600 ISW=1
12700 CALL ALLOC(INDEX,200,KRET)
12800 IY(I)=INDEX
12900 CALL PLREAD(VARN(1,I),IP,IFF,XD(IY(I)),IERR,ISW)
13000 IF (IERR .LT. -1) GO TO 993
13001 IF (IERR .NE. -1) GO TO 352
13002 TYPE 9933,(VARN(J,I),J=1,4)
13004 GO TO 351
13100 352 NP(I)=KINDEX(IFF,IP,NFREQ)
13200 IWHX(I)=1
13300 51 CONTINUE
13400 GO TO 52
13500 35 NLINES=NLINES-NL
13600 GO TO 32
13700 C
13800 52 IF (TS) GO TO 75
13900 C
14000 TYPE 109
14100 109 FORMAT(' NAME OF INDEP. VAR :'$)
14200 ACCEPT 200,II
14300 CALL SYMCK1(II,$75,$6,$52)
14400 REREAD 203,(IVARN(J),J=1,4)
14500 IP=IFPL
14600 IFF=ILPL
14700 ISW=1
14800 CALL ALLOC(INDEX,200,KRET)
14900 IX(1)=INDEX
15000 CALL PLREAD(IVARN,IP,IFF,XD(IX(1)),IERR,ISW)
15100 IF (IERR .LT. 0) GO TO 993
15200 GO TO 75
15300 C
15400 C FILE NAME OR TTY
15600 61 TYPE 110
15700 110 FORMAT(3X,'# OF X,Y PAIRS EACH LINE :'$)
15800 ACCEPT 200,II
15900 CALL SYMCK1(II,$62,$6,$63)
16000 REREAD 205,NPRS
16100 GO TO 62
16200 63 K=K-1
16300 GO TO 20
16400 62 TYPE 111
16500 111 FORMAT(3X,'# OF LINES THIS INDEP VAR :'$)
16600 ACCEPT 200,II
16700 CALL SYMCK1(II,$64,$6,$61)
16800 REREAD 205,NL
16900 NLINES=NLINES+NL
17000 620 IF (NLINES .GT. 10) GO TO 992
17100 KFIN=K+NL-1
17200 DO 621 I=K,KFIN
17300 IWHX(I)=KK
17400 621 NP(I)=NPRS
17500 64 GO TO (990,65,70),IPT
17600 C
17700 C FILE NAME
17800 65 TYPE 112
17900 112 FORMAT(3X,'DEP VAR. #''S. (SEP. BY ,) :'$)
18000 CALL SETCOR(' ',IDLIN(1),IDLIN(60))
18100 ACCEPT 206,(IDLIN(J),J=1,60)
18200 206 FORMAT(60A1)
18300 DO 650 I=1,60
18400 650 CALL SYMCK1(IDLIN(I),$652,$6,$651)
18500 REREAD 212,(VARNM(J),J=K,KFIN)
18600 GO TO 652
18700 651 NLINES=NLINES-NL
18800 GO TO 62
18900 212 FORMAT(10I)
19000 C
19100 C FIND MAX VAR WANT
19200 652 MV=0
19300 DO 66 J=K,KFIN
19400 66 MV=MAX0(MV,VARNM(J))
19500 IF (TS) GO TO 68
19600 TYPE 113
19700 113 FORMAT(3X,'INDEP VAR. # :'$)
19800 ACCEPT 200,II
19900 CALL SYMCK1(II,$67,$6,$65)
20000 REREAD 205,IVARNM
20100 67 IF (IVARNM .GT. MV) MV=IVARNM
20200 C
20300 C READ IN FROM FILE
20400 68 IF(IOLD(1).NE.INEW(1).OR.IOLD(2).NE.INEW(2))
20450 . CALL IFIL(22,INEW(1),INEW(2),IPP)
20500 IOLD(1)=INEW(1)
20550 IOLD(2)=INEW(2)
20600 DO 700 I=K,KFIN
20700 CALL ALLOC(INDEX,NPRS,KRET)
20800 700 IY(I)=INDEX
20900 CALL ALLOC(INDEX,NPRS,KRET)
21000 IX(KK)=INDEX
21100 DO 690 I=1,NPRS
21300 READ(22) (XBUF(J),J=1,MV)
21350 IF (TS) GO TO 681
21400 XD(IX(KK)+I-1)=XBUF(IVARNM)
21500 681 DO 69 J=K,KFIN
21600 69 XD(IY(J)+I-1)=XBUF(VARNM(J))
21700 690 CONTINUE
21800 GO TO 75
21900 C
22000 C FROM TTY
22100 70 IF (TS) GO TO 72
22200 TYPE 114
22300 114 FORMAT(3X,'TYPE X VALUES (10 PER LINE SEP. BY ,)'/' *'$)
22400 ACCEPT 200,II
22500 CALL SYMCK1(II,$72,$6,$651)
22600 CALL ALLOC(INDEX,NPRS,KRET)
22700 IX(KK)=INDEX
22800 MMMAX=9
22900 MLEFT=0
23000 IF (NPRS .GT. 10) MLEFT=1
23100 REREAD 214,(XD(IX(KK)+I),I=0,MMMAX)
23200 IF (MLEFT .GT. 0) ACCEPT 214,(XD(IX(KK)+I),I=10,NPRS-1)
23300 214 FORMAT(10F)
23400 72 TYPE 115
23500 115 FORMAT(3X,'TYPE Y VALUES (10 PER LINE SEP. BY ,)'/)
23600 DO 71 I=K,KFIN
23700 L=I-K+1
23800 TYPE 116,L
23900 116 FORMAT(' *VAR ',I1,' :'$)
24000 ACCEPT 220,II
24100 CALL SYMCK1(II,$71,$6,$72)
24200 CALL ALLOC(INDEX,NP(I),KRET)
24300 IY(I)=INDEX
24400 MMMAX=9
24500 MLEFT=0
24600 IF (NP(I) .GT. 10) MLEFT=1
24700 REREAD 214,(XD(IY(I)+J),J=0,MMMAX)
24800 IF (MLEFT .GT. 0) ACCEPT 214,(XD(IY(I)+J),J=10,NP(I)-1)
24900 71 CONTINUE
25000 C
25100 C FIND MIN AND MAX OF X AND Y VALUES
25200 75 DO 73 I=K,KFIN
25300 DO 73 J=1,NP(I)
25400 YMIN=AMIN1(YMIN,XD(IY(I)+J-1))
25500 YMAX=AMAX1(YMAX,XD(IY(I)+J-1))
25600 73 CONTINUE
25700 IF (TS) GO TO 80
25800 DO 74 I=1,NPRS
25900 XMIN=AMIN1(XMIN,XD(IX(KK)+I-1))
26000 XMAX=AMAX1(XMAX,XD(IX(KK)+I-1))
26100 74 CONTINUE
26200 C
26300 C GET LINE TYPE FOR EACH LINE, SEP BY COMMAS
26400 80 TYPE 117
26500 117 FORMAT(3X,'LINE TYPE FOR EACH (SEP. BY ,):'$)
26600 ACCEPT 200,II
26700 CALL SYMCK1(II,$82,$6,$81)
26800 REREAD 212,(LT(J),J=K,KFIN)
26900 GO TO 82
27000 81 GO TO (52,65,70),IPT
27100 C
27200 C GET CHAR FOR EACH LINE, SEP BY COMMAS - NO BLANKS
27300 82 TYPE 118
27400 ACCEPT 200,II
27500 CALL SYMCK1(II,$84,$6,$80)
27600 118 FORMAT(3X,'CHAR FOR EACH LINE (SEP. BY ,) :'$)
27700 REREAD 218,(CHAR(J),J=K,KFIN)
27800 218 FORMAT(10(A1,1X))
27900 C
28000 C GET TITLES FOR EACH DEP VARIABLE
28100 84 TYPE 123
28200 123 FORMAT(3X,'UP TO 20 CHAR ID FOR EACH DEP VAR.'/)
28300 DO 91 I=K,KFIN
28400 L=I-K+1
28500 GO TO (913,92,914),IPT
28600 913 TYPE 124,(VARN(J,I),J=1,4)
28700 124 FORMAT(1X,4A4,' :'$)
28800 GO TO 93
28900 92 TYPE 125,VARNM(I)
29000 125 FORMAT(3X,'VAR ',I1,' :'$)
29100 GO TO 93
29200 914 TYPE 125,L
29300 93 ACCEPT 220,II
29400 CALL SYMCK(II,LNLBL(1,I),4,$91,$6)
29500 REREAD 220,(LNLBL(J,I),J=1,4)
29600 91 CONTINUE
29700 IF (TS) GO TO 90
29800 C
29900 C IF NOT FROM PLANETS, MAY BE MORE LINES THIS PLOT - ANOTHER INDEP VAR.
30000 TYPE 119
30100 119 FORMAT(3X,'MORE LINES THIS PLOT (Y OR N) :'$)
30200 ACCEPT 200,IANS
30300 CALL SYMCK1(IANS,$990,$6,$82)
30400 IF (IANS .NE. 'Y') GO TO 90
30500 K=NLINES
30600 KK=KK+1
30700 IF (KK .GT. 4) GO TO 991
30800 IF (K .GT. 10) GO TO 992
30900 GO TO 20
31000 C
31100 C GET X AXIS TITLE
31200 90 TYPE 120
31300 120 FORMAT(3X,'X-AXIS TITLE (<= 60 CHAR)'/' *'$)
31400 ACCEPT 220,II
31500 220 FORMAT(12A5)
31600 CALL SYMCK(II,XT,12,$912,$6)
31700 REREAD 220,XT
31800 C
31900 C GET Y AXIS TITLE
32000 912 TYPE 121
32100 121 FORMAT(3X,'Y-AXIS TITLE (AS ABOVE)'/' *'$)
32200 ACCEPT 220,II
32300 CALL SYMCK(II,YT,12,$910,$6)
32400 REREAD 220,YT
32500 C
32600 C GET PLOT TITLE
32700 910 TYPE 122
32800 122 FORMAT(3X,'PLOT TITLE (AS ABOVE) '/' *'$)
32900 ACCEPT 220,II
33000 CALL SYMCK(II,PT,12,$915,$6)
33100 REREAD 220,PT
33200 C
33300 C GET LENGTH OF X AXIS
33400 915 IF (.NOT. TS)TYPE 126,XMIN,XMAX
33500 126 FORMAT(3X,'HORIZ. RANGE IS : ',2G12.3)
33600 1262 FORMAT(3X,'X-AXIS LENGTH ( > 1 IN.) :'$)
33700 IF (TS) TYPE 1261,IP,IFF
33800 1261 FORMAT(3X,'HORIZ. RANGE IS : ',2I8)
33900 TYPE 1262
34000 ACCEPT 200,II
34100 CALL SYMCK1(II,$97,$6,$90)
34200 IF (TS) REREAD 2141,XIN,IHMIN,IHMAX
34300 IF (.NOT. TS) REREAD 214,XIN,HMIN,HMAX
34400 2141 FORMAT(F,2I)
34500 C
34600 C FIND Y AXIS LENGTH
34700 97 TYPE 127,YMIN,YMAX
34800 127 FORMAT(3X,'VERT. RANGE IS : ',2G12.3/3X,'Y-AXIS LENGTH '
34900 . ,'(2-10 IN.) :'$)
35000 ACCEPT 200,II
35100 CALL SYMCK1(II,$98,$6,$915)
35200 REREAD 214,YIN,VMIN,VMAX
35300 C
35400 C DETERMINE IS LOG SCALE
35500 98 TYPE 128
35600 128 FORMAT(3X,'LOG SCALE - X,Y, OR B(OTH) :'$)
35700 ACCEPT 202,LOGS
35800 IF (LOGS .EQ. '$') GO TO 6
35900 LX=0
36000 LY=0
36100 IF (LOGS .EQ. 'X' .OR. LOGS .EQ. 'B') LX=1
36200 IF (LOGS .EQ. 'Y' .OR. LOGS .EQ. 'B') LY=1
36300 C
36400 C DETERMINE IF GRID LINES
36500 TYPE 129
36600 129 FORMAT(3X,'GRID LINES - X,Y, OR B(OTH) :'$)
36700 ACCEPT 202,IGRID
36800 IF (IGRID .EQ. '$') GO TO 6
36900 XGR=0.
37000 YGR=0.
37100 IF (IGRID .EQ. 'X' .OR. IGRID .EQ. 'B') XGR=1.
37200 IF (IGRID .EQ. 'Y' .OR. IGRID .EQ. 'B') YGR=1.
37300 IF (TS ) GO TO 95
37400 C
37500 C DETERMINE IF ZERO EXCLUDED FROM SCALING
37600 TYPE 130
37700 130 FORMAT(3X,'ZERO EXCL. FROM SCALE - X,Y, OR B(OTH) :'$)
37800 94 ACCEPT 202,IZR
37900 IF (IZR .EQ. '$') GO TO 6
38000 ZX=1.
38100 ZY=1.
38200 IF (IZR .EQ. 'X' .OR. IZR .EQ. 'B') ZX=0.
38300 IF (IZR .EQ. 'Y' .OR. IZR .EQ. 'B') ZY=0.
38400 GO TO 96
38500 95 TYPE 131
38600 131 FORMAT(3X,'ZERO EXCL. FROM Y-AXIS (Y OR N) :'$)
38700 GO TO 94
38800 96 CONTINUE
38900 C
39000 C DONE INITIALIZING, PREPARE POINTS FOR PLOTTING
39100 KP=1
39200 DO 300 I=1,NLINES
39300 L=IWHX(I)
39400 DO 301 J=1,NP(I)
39500 IF (.NOT. TS) XS(KP)=XD(IX(L)+J-1)
39600 YS(KP)=XD(IY(I)+J-1)
39700 KP=KP+1
39800 301 CONTINUE
39900 300 CONTINUE
40000 XS(KP)=XMIN
40100 YS(KP)=YMIN
40200 XS(KP+1)=XMAX
40300 YS(KP+1)=YMAX
40400 KP=KP+2
40500 C
40600 C SEE IF USER SPECIFIED MAX AND MIN
40700 IF (TS) GO TO 302
40800 IF (HMIN .EQ. 0.0 .AND. HMAX .EQ. 0.0) GO TO 302
40900 IF (HMIN .LT. XMIN) XS(KP-2)=HMIN
41000 IF (HMAX .GT. XMAX) XS(KP-1)=HMAX
41100 302 IF (VMIN .EQ. 0. .AND. VMAX .EQ. 0) GO TO 303
41200 IF (VMIN .LT. YMIN) YS(KP-2)=VMIN
41300 IF (VMAX .GT. YMAX) YS(KP-1)=VMAX
41400 303 N=KP-1
41500 IF (TS ) GO TO 350
41600 C
41700 C REMEMBER: LAST 2 POINTS ARENT TO BE PLOTTED
41800 C ONLY SCALED
41900 C
42000 C NOT A TIME SERIES
42100 CALL XYAXIS(LX,LY,XS,YS,N,XT,YT,XIN,YIN,SX,SY,ZX,ZY)
42200 IST=0
42300 IST1=0
42400 C
42500 C ADJUST LINE TYPES
42600 304 DO 310 I=1,NLINES
42700 IST=IST+1
42800 IST1=IST1+1
42900 IF (LT(I) .GT. 1) LT(I)=LT(I)+2
43000 IF (LT(I) .EQ. 1) LT(I)=2
43100 IF (LT(I) .LE. 0) GO TO 305
43200 CALL LINE(PXP(IST1),PYP(IST),NP(I),1,LT(I))
43300 305 IF (CHAR(I) .EQ. ' ') GO TO 307
43400 CALL POINTS(PXP(IST1),PYP(IST),NP(I),1,CHAR(I),.09)
43500 307 IST=IST+NP(I)-1
43600 IST1=IST
43700 IF (TS) IST1=0
43800 310 CONTINUE
43900 GO TO 390
44000 C
44100 C TIME SERIES PLOTS
44200 350 CALL TAXIS(SX,SY,XIN,IP,IFF,NFREQ,XDEL)
44300 XQ=SX+XIN+.75
44400 CALL YAXIS(LY,YS,N,YT,YIN,SX,SY,ZY)
44500 PXP(1)=SX+.75
44600 DO 360 I=2,NP(1)
44700 360 PXP(I)=PXP(I-1)+XDEL
44800 CALL PLACET(XT,BX,SZ,NCHARS,IBL,.09,XQ,SX+.75)
44900 CALL SYMBOL(BX,SY+.18,.09,XT(IBL),0.,NCHARS)
45000 C
45100 C DRAW ZERO LINE
45200 IF (SY+.75 .EQ. YZ) GO TO 304
45300 CALL PLOT(SX+.75,YZ,3)
45400 CALL PLOT(XQ,YZ,2)
45500 GO TO 304
45600 C
45700 C ALL DONE WITH PLOTS - DO TITLE
45800 390 CALL PLACET(PT,BX,SZ,NCHARS,IBL,.09,XQ,SX+.75)
45900 IF (YQ+.55 .GT. 10.7) YQ=YQ-.3
46000 CALL SYMBOL(BX,YQ+.55,.1,PT(IBL),0.,NCHARS)
46100 C
46200 C PLOT INDIVIDUAL LINE LABEL
46300 BY=NLINES*.25+SY+1.
46400 CALL PLOT(XQ+1,BY,3)
46500 BX=XQ+1.
46600 DO 395 I=1,NLINES
46700 SPX(1)=BX
46800 SPX(2)=BX+1.5
46900 SPY(1)=BY
47000 SPY(2)=BY
47100 IF (LT(I) .LE. 0) GO TO 393
47200 CALL LINE(SPX,SPY,2,1,LT(I))
47300 393 IF (CHAR(I) .EQ. ' ') GO TO 394
47400 CALL SYMBOL(SPX(2)+.3,BY,.12,CHAR(I),0.,1)
47500 394 CALL SYMBOL(SPX(2)+.5,BY,.12,LNLBL(1,I),0.,20)
47600 BY=BY-.25
47700 395 CONTINUE
47800 CALL PLOT(SX+XIN+6.,0.,-3)
47900 TYPE 132
48000 132 FORMAT(3X,'PLOT COMPLETED - ANOTHER PLOT (Y OR N) :'$)
48100 ACCEPT 232,IANS
48200 232 FORMAT(A1)
48300 IF (IANS .NE. 'Y') GO TO 400
48400 GO TO 2
48500 400 CALL CPLOT
48600 410 CALL EXIT
48700 C
48800 C ERROR MESSAGES
48900 990 TYPE 9901
49000 9901 FORMAT(' IMPOSSIBLE ERROR '/)
49100 GO TO 410
49200 991 TYPE 9911
49300 9911 FORMAT(' ONLY 4 INDEP VBLS. ALLOWED PER PLOT '/)
49400 GO TO 400
49500 992 TYPE 9921
49600 9921 FORMAT(' ONLY 10 DEP. VBLS. ALLOWED PER PLOT '/)
49700 GO TO 400
50000 9933 FORMAT(' ',4A4,' NOT IN DATA BANK '/)
50200 993 TYPE 9931,IERR
50300 9931 FORMAT(' PLANETS DATA BANK PROBLEM. FROM PLREAD IERR = ',
50400 . I4/)
50500 GO TO 400
50600 994 TYPE 9941
50700 9941 FORMAT(' NOT ENUF CORE TO ALLOCATE '/)
50800 GO TO 400
50900 5000 PAUSE 'PLANETS DATA BANK NOT THERE'
51000 GO TO 10
51100 END
51200 C
51300 C::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
51400 C
51500 SUBROUTINE SYMCK(II,IA,L,$,$)
51600 DIMENSION IA(1)
51700 IF (II .EQ. '>>') RETURN 4
51800 IF (II .EQ. '$') RETURN 5
51900 CALL SETCOR(' ',IA(1),IA(L))
52000 RETURN
52100 END
52200 C
52300 C::::::::::::::::::::::::::::::::::::::::::::::::::::::::
52400 C
52500 SUBROUTINE SYMCK1(II,$,$,$)
52600 IF (II .EQ. '>>') RETURN 2
52700 IF (II .EQ. '$') RETURN 3
52800 IF (II .EQ. '<<') RETURN 4
52900 RETURN
53000 END