C ***************************************************************** 00000001 C PRIM PRIMARY STRUCTURE INPUT 00000002 C ***************************************************************** 00000003 C READ(12,.) - FROM CARDS, TAPE OR DISC 00000008 C ================================================================= 00000009 SUBROUTINE PRIM(NL1) 00000010 COMMON MTN,MLP,N1,N,N2/INOU/KPROU1,AM(40),EVSTR(22),JA(507),B(500)00000011 1/CT/NZON,NN,NNN/VECT/V(2580)/ENER/LINST(5,150,4)/ENER1/V2(1800) 00000012 2/DT/X(160),VW(20),NV/NDA/HYDR(20),AM1(40),NV1,NSTRT,NINT,NN0,NN1, 00000013 3KEND,NBEG0,NKUS,NCL1,NCL,ND,NN11,NV11,NKC1,NKN1,NNAA,NSTRT1 00000014 DIMENSION PLIN(80),ANAME(6),SECST(4,500),JSEC(3,4), 00000015 *LINE(6),TEXT(33),LC(3,150),V10(500),V11(500),V12(500), 00000016 *V13(500),V14(500) 00000017 EQUIVALENCE (V(1),V10(1),NW),(V(2),ANAME(1)),(V(501),SECST(1)), 00000018 *(V(8),JSEC(1)),(V(20),LC(1)),(V11(1),V(501)),(V12(1),V(1001)) 00000019 *,(V13(1),V(1501)),(V14(1),V(2001)),(PLIN(1),V(2501)) 00000020 *,(V(273),NAC),(V(274),KON) 00000021 DATA TEXT/4H1 ,4H2 ,4H3 ,4H4 ,4H5 ,4H6 ,4H7 ,4H8 ,00000022 14H9 ,4H- ,4H? ,4HH ,4HS ,4H& ,4HX ,4HB ,4HQ ,4HN00000023 2 ,4HY ,4HA ,4HP ,4HM ,4H| ,4H ,4H. ,4H' ,4H: 00000024 3,4H* ,4H+ ,4H= ,4H( ,4H) ,4H/ /,TX/' '/,XX/'='/, 00000025 4 K45/5/,K73/4/,K22/4/,K23/3/,I21/7/,I25/20/,K20/0/,YX/'*'/ 00000026 5 FORMAT(80A1) 00000045 READ (12,25) ANAME,NDOM,NPR,NAN,NV 00000046 25 FORMAT(6A4,36X,4I5) 00000047 READ (12,203) N1,NN,NNN,NPR1,NDOM1 00000048 203 FORMAT(I3,I4,I2,1X,I5,I5) 00000051 C 00000052 NZAP=0 00000055 444 NZAP1=NZAP/22+1 00000056 NZAP=NZAP+1 00000057 K=NZAP-NZAP1*22+22 00000058 K2=K*80 00000059 K1=K2-79 00000060 READ (12,5) (V2(I),I=K1,K2) 00000061 IF(V2(K1).EQ.XX.OR.K2.GE.1760) CALL RITE(NZAP1,1800,V2) 00000062 IF(V2(K1).NE.XX) GOTO 444 00000063 C +++ 00000064 PRINT 200,ANAME,NDOM,NPR,NAN,NV,N1,NN,NNN,NPR1,NDOM1 00000065 200 FORMAT(3X,1H|,80(1H=),37H|=INPUT & INPUT DATA & CONTROL=======/ 00000066 14H |,16(5H----+),1H|,36(1H-)/4H |,6A4,36X,4I5,27H| ANAME, NDOM00000067 2, NPR, NAN, NV/4H |,80X,1H|,36(1H-)/4H |,I3,I4,I2,I6,I5,60X, 00000068 329H| N1, NN, NNN, NPR1, NDOM1 /4H |,80X,1H|,36(1H-)) 00000069 IF(NPR1.GT.0) NPR=NPR1 00000070 IF(NPR.LT.1) NPR=1 00000071 IF(NDOM1.GT.0) NDOM=NDOM1 00000072 IF(NDOM.LT.1) NDOM=1 00000073 IF(NAN.NE.0) NAN=ISIGN(1,NAN) 00000074 IF(NN.EQ.0.AND.NNN.EQ.0) NNN=NAN 00000075 C CONTROL NV,N1,NN 00000076 IF(NV.LT.N1) PRINT 26,NV,N1 00000077 26 FORMAT(/' ||| CHAIN LENGTH NV=',I5,' WAS LESS THAN FRAGMENT L00000078 1ENGTH N1=',I5,'; NV=N1 IS TAKEN |'/) 00000079 IF(NV.LT.N1) NV=N1 00000080 IF(NN.LT.0.OR.NN.GE.NV) NN=0 00000081 N=NV-NN 00000082 IF(N1.LE.0) N1=N 00000083 IF(N.LT.N1) PRINT 27,N1,NV,NN,N 00000084 27 FORMAT(/' ||| FRAGMENT LENGTH N1=',I5,' WAS GREATER THAN ONE OF00000085 1 REMAINING CHAIN: N=NV-NN=',I5,'-',I5,'=',I5,'; N1=N IS TAKEN |'/)00000086 IF(N.LT.N1) N1=N 00000087 NN11=N1 00000088 ND=50 00000089 ND2=NL1/5 00000090 IF(ND.GT.ND2) ND=ND2 00000091 NKUS=(N1-2*ND-1)/(NL1-2*ND)+1 00000092 K=(N1+2*ND*(NKUS-1)-1)/NKUS+1 00000093 K2=(N1+2*ND*NKUS-1)/NKUS+1 00000094 IF(K2.LE.NL1.AND.NKUS.GT.2) K=K2 00000095 NCL1=K 00000096 IF(K2.LE.NL1.AND.NKUS.GT.2) K=K-ND 00000097 NKC1=N1-K+2*ND-(NKUS-2)*(NCL1-2*ND) 00000098 NKN1=K 00000099 NV1=NN+N1 00000100 NBEG0=NN 00000101 NNAA=NN 00000102 NV11=NV 00000103 CALL REW(K,1800,V2) 00000104 GOTO 450 00000105 C 00000106 ENTRY PRIMM(NL1) 00000109 NN=NNAA 00000110 NSTRT=NSTRT1 00000111 NNN=1 00000112 450 NV=NV11 00000113 NCL=NCL+1 00000114 N1=NKN1 00000115 IF(NCL.GT.1) N1=NCL1 00000116 IF(NCL.GT.1.AND.NCL.EQ.NKUS) N1=NKC1 00000117 C NEXT FRAGMENT 00000118 NNAA=NN+N1-2*ND 00000119 IF(NSTRT.GT.ND) N1=N1-10 00000120 IF(NSTRT.GT.ND) NN=NN+10 00000121 IF(NSTRT.GT.ND) NSTRT=NSTRT-10 00000122 NINT=0 00000123 IF(NCL.LT.NKUS) NINT=ND 00000124 C 00000125 IF(NCL.EQ.NKUS) GOTO 29 00000126 NSTRT1=ND 00000127 K1=NNAA+ND 00000128 K1=K1-(K1/20)*20 00000129 C 00000130 IF(K1.LT.10) N1=N1-K1 00000131 IF(K1.LT.10) NSTRT1=ND-K1 00000132 C 00000133 IF(K1.GE.10) NINT=NINT-20+K1 00000134 IF(K1.GE.10) NSTRT1=ND+20-K1 00000135 29 KZAP=1760 00000136 NZAP=0 00000137 C 00000138 NZON=0 00000142 N2=NN+1 00000143 MLP=-NN-NN11 00000144 MTN=NV 00000145 M=-MLP 00000146 DO 1 K=1,20 00000147 1 VW(K)=0. 00000148 N=0 00000149 201 FORMAT(4H |,80A1,10H| SEQUENCE) 00000150 C ++ 00000151 4 NZAP2=NZAP/22+1 00000152 IF(KZAP.GE.1680) CALL RID1(NZAP1,1800,V2) 00000153 NZAP=NZAP+1 00000154 K=NZAP-NZAP2*22+22 00000155 KZAP=(K-1)*80 00000156 DO 446 I=1,80 00000157 J=I+KZAP 00000158 446 PLIN(I)=V2(J) 00000159 C ++ 00000160 IF(NCL.EQ.1) PRINT 201,PLIN 00000161 C SEQUENCE INPUT BLANKS ARE IGNORED 00000162 DO 8 K=1,80 00000163 Q=PLIN(K) 00000164 IF(Q.EQ.TX.OR.N.GE.NV) GOTO 8 00000165 N=N+1 00000166 II=21 00000167 DO 7 J=1,40 00000168 IF(Q.NE.AM1(J)) GOTO 7 00000169 II=J 00000170 GOTO 6 00000171 7 CONTINUE 00000172 PRINT 19,N,Q 00000173 6 I1=II 00000174 IF(II.GT.20) I1=II-20 00000175 IF(NNN.GE.-1.AND.NCL.EQ.1) VW(I1)=VW(I1)+1./FLOAT(NV) 00000176 I=N-NN 00000177 IF(I.LT.1.OR.I.GT.N1) GOTO 8 00000178 B(I)=Q 00000179 JA(I)=II 00000180 IF(NNN.LT.-1.AND.NCL.EQ.1) VW(I1)=VW(I1)+1./FLOAT(NN11) 00000181 8 CONTINUE 00000182 IF(N.LT.NV) GOTO 4 00000183 IF(NCL.GT.1) GOTO 451 00000184 NV=INT(FLOAT(NV*NPR)/FLOAT(NDOM)) 00000185 IF(NNN.LT.-1) NV=NN11 00000186 IF(NNN.NE.0) NNN=ISIGN(1,NNN) 00000187 19 FORMAT(1X,119('-') /' *||* WRONG RESIDUE',I5,' NAME: ',A1 /) 00000188 451 DO 48 K=1,4 00000189 DO 46 K1=1,3 00000190 46 JSEC(K1,K)=0 00000191 DO 47 K1=1,5 00000192 DO 47 K2=1,100 00000193 47 LINST(K1,K2,K)=0 00000194 DO 48 K1=1,N1 00000195 48 SECST(K,K1)=TX 00000196 DO 49 K=1,3 00000197 DO 49 K1=1,100 00000198 49 LC(K,K1)=0 00000199 IF(K73.GE.4) GOTO 500 00000200 K2=K73+1 00000201 DO 499 K=K2,4 00000202 JSEC(1,K)=1 00000203 499 JSEC(2,K)=1 00000204 500 CONTINUE 00000205 IF(NCL.EQ.1) PRINT 242 00000206 2 IF(PLIN(1).EQ.XX) GOTO 382 00000212 IF(PLIN(1).EQ.YX) GOTO 383 00000213 C ++ 00000214 NZAP2=NZAP/22+1 00000215 IF(KZAP.GE.1680) CALL RID1(NZAP1,1800,V2) 00000216 NZAP=NZAP+1 00000217 K=NZAP-NZAP2*22+22 00000218 KZAP=(K-1)*80 00000219 DO 447 I=1,80 00000220 J=I+KZAP 00000221 447 PLIN(I)=V2(J) 00000222 C ++ 00000223 IF(NCL.EQ.1) PRINT 211,PLIN 00000224 GOTO 2 00000225 211 FORMAT(4H |,80A1,24H| ADDITIONAL INFORMATION) 00000226 242 FORMAT(4H |,16(5H----+),1H|,36(1H-)) 00000227 383 IF(NCL.EQ.1) PRINT 242 00000228 205 FORMAT(4H |,80A1,36H| SEC.STR: L&G,AUTHOR,H-BONDS,KABSCH) 00000229 206 FORMAT(4H |,18A1,62X,36H| CONNECTION: '1'-SS,'2'-COORD.,ETS.) 00000230 K5=0 00000231 K6=0 00000232 C INPUT OF X-RAY SECONDARY STRUCTURE (FOR COMPARISON) 00000233 C * BEFORE THESE CARDS: ONE CARD WITH * IN THE 1-ST POSITION 00000234 C = AFTER ALL CARDS: ONE CARD WITH = IN THE 1-ST POSITION 00000235 C SECONDARY STRUCTURE CARDS: 00000236 C * HELIX : |H.NNNN-EEEE........... .. .. | 00000237 C ------------------ --,,,-- --,,,-- --,,,-- 00000238 C * TURN |T.NNNN-EEEE........... .. .. | 00000239 C ------------------ --,,,-- --,,,-- --,,,-- 00000240 C * STRAND: |S.NNNN-EEEE*K1*K2*K3.. .. .. | 00000241 C ------------------ --,,,-- --,,,-- --,,,-- 00000242 C * CONTACT: |C.NNNN-EEEE..K.................................| 00000243 C | | | | | | | | | | | | | 00000244 C COLUMN 1 3 7 11 14 17 20 23 40 43 60 63 80 00000245 C 00000246 C 'NNNN'-START,'EEEE'-FINISH, 00000247 C CONTACTS OF STRANDS: '*'-> + PARALLEL; '*'-> - ANTIPARALLEL 00000248 C WITH: 'GIVEN STRAND - K1', 'GIVEN STRAND -K2', ... 00000249 C (E.G. K1=1: WITH PRECEEDING IN THE LIST) 00000250 C K - TYPE OF CONTACT : 1 - SS-BOND, 2 - COORD. BOND, ETC. 00000251 C YOU CAN USE UP TO 4 SETS OF STRUCTURE (E.G. 'LEVITT'S','AUTHOR'S',00000252 C ETC.) THE COLUMNS 23-40, 43-60, 63-80 FOR S, H, T : 00000253 C ARE USED LIKE COLUMNS 3-20) 00000254 C ++ 00000255 30 NZAP2=NZAP/22+1 00000256 IF(KZAP.GE.1680) CALL RID1(NZAP1,1800,V2) 00000257 NZAP=NZAP+1 00000258 K=NZAP-NZAP2*22+22 00000259 KZAP=(K-1)*80 00000260 DO 448 I=1,80 00000261 J=I+KZAP 00000262 448 PLIN(I)=V2(J) 00000263 C ++ 00000264 Q1=PLIN(1) 00000265 K4=0 00000266 IF(Q1.EQ.AM1(18)) K4=1 00000267 IF(Q1.EQ.TEXT(12))K4=2 00000268 IF(Q1.EQ.TEXT(13))K4=K45 00000269 K5=K4 00000270 IF(K5.EQ.1) K5=2 00000271 IF(Q1.EQ.AM1(2)) K5=3 00000272 IF(K5.EQ.0.AND.Q1.EQ.XX) GOTO 37 00000273 IF(K5.EQ.0) GOTO 371 00000274 IF(NCL.GT.1) GOTO 452 00000275 IF(Q1.EQ.AM1(2)) PRINT 206,(PLIN(K),K=1,18) 00000276 IF(Q1.EQ.AM1(18).OR.Q1.EQ.TEXT(12).OR.Q1.EQ.TEXT(13)) PRINT 205, 00000277 *PLIN 00000278 C DECODING: SECONDARY STRUCTURE 00000279 452 I1=1+K4/K45 00000280 I2=I21 00000281 K7=K73 00000282 IF(K5.EQ.3) K7=1 00000283 DO 31 K=1,K7 00000284 K21=K22 00000285 I3=I2 00000286 IF(PLIN(I2).EQ.TEXT(11).AND.K4.GE.2) JSEC(I1,K)=1 00000287 DO 32 K1=1,K5 00000288 M=0 00000289 J1=1 00000290 IF(K1.GT.2) K21=K23 00000291 C 00000292 DO 33 K2=1,K21 00000293 J=0 00000294 L1=I3-K2 00000295 P=PLIN(L1) 00000296 DO 34 K3=1,9 00000297 Q=TEXT(K3) 00000298 34 IF(Q.EQ.P) J=K3 00000299 M=M+J1*J 00000300 33 J1=J1*10 00000301 IF(P.EQ.TEXT(10)) M=-M 00000302 IF(K1.EQ.1) I3=I3+1+K22 00000303 IF(K1.GT.1) I3=I3+K20+K23 00000304 32 LINE(K1)=M 00000305 I2=I2+I25 00000306 K1=MIN0(IABS(LINE(1)),IABS(LINE(2)))-NN 00000307 K2=MAX0(IABS(LINE(1)),IABS(LINE(2)))-NN 00000308 IF(K2.LT.1.OR.K1.GT.N1) GOTO 31 00000309 C 00000310 IF(K2.GT.N1.OR.K1.LT. 1.OR.K5.NE.3) GOTO 50 00000311 K6=K6+1 00000312 IF(K6.GT.150) PRINT 23 00000313 IF(K6.GT.150) GOTO 31 00000314 LC(1,K6)=K1 00000315 LC(2,K6)=K2 00000316 LC(3,K6)=LINE(3) 00000317 50 IF(K5.EQ.3) GOTO 31 00000318 IF(K1.LT.1) K1=1 00000319 IF(K2.GT.N1) K2=N1 00000320 IF(I1.EQ.2) GOTO 35 00000321 P=TEXT(14) 00000322 IF(K4.EQ.1) P=AM1(18) 00000323 C 00000324 IF(K1.EQ.1.OR.K4.EQ.1) GOTO 361 00000325 K0=K1-1 00000326 IF(SECST(K,K0).EQ.P) SECST(K,K0)=TEXT(33) 00000327 361 DO 36 J=K1,K2 00000328 36 SECST(K,J)=P 00000329 GOTO 31 00000330 C 00000331 35 M4=JSEC(3,K)+1 00000332 IF(M4.GT.150) PRINT 23 00000333 IF(M4.GT.150) GOTO 31 00000334 JSEC(3,K)=M4 00000335 LINST(1,M4,K)=K1 00000336 LINST(2,M4,K)=K2 00000337 DO 310 K3=3,K45 00000338 J=LINE(K3) 00000339 IF(J.EQ.0) GOTO 310 00000340 K1=ISIGN(1,J)+4 00000341 LINST(K1,M4,K)=1 00000342 LINST(4,M4,K)=LINST(4,M4,K)+1 00000343 C 00000344 K2=M4-IABS(J) 00000345 IF(K2.LE.0) GOTO 310 00000346 LINST(K1,K2,K)=1 00000347 LINST(4,K2,K)=LINST(4,K2,K)+1 00000348 310 CONTINUE 00000349 31 CONTINUE 00000350 GOTO 30 00000351 371 IF(NCL.EQ.1) PRINT 22,PLIN 00000352 C 00000353 37 IF(NCL.EQ.1) PRINT 211,PLIN 00000354 DO 38 K=1,K73 00000355 K3=JSEC(3,K) 00000356 IF(K3.LT.1) GOTO 38 00000357 DO 381 K4=1,K3 00000358 K1=LINST(1,K4,K) 00000359 K2=LINST(2,K4,K) 00000360 J=LINST(4,K4,K) 00000361 IF(J.GT.2) J=2 00000362 J1=15 00000363 IF(J.GT.0) J1=4*(J-1)+2*LINST(5,K4,K)+LINST(3,K4,K)+J1 00000364 IF(K1.EQ.1) GOTO 362 00000365 K0=K1-1 00000366 DO 363 J=15,22 00000367 IF(SECST(K,K0).EQ.TEXT(J)) SECST(K,K0)=TEXT(29) 00000368 363 CONTINUE 00000369 362 DO 381 J=K1,K2 00000370 381 SECST(K,J)=TEXT(J1) 00000371 38 CONTINUE 00000372 382 IF(NCL.EQ.1) PRINT 242 00000373 NW=N1 00000374 NAC=NBEG0+1 00000375 KON=NBEG0+NN11 00000376 IF(NCL.EQ.1) CALL RITE(41,500,V10) 00000377 CALL RITE(42,500,V11) 00000378 CALL RITE(43,500,V12) 00000379 CALL RITE(44,500,V13) 00000380 CALL RITE(45,500,V14) 00000381 IF(NCL.GT.1) PRINT 453,NCL 00000382 23 FORMAT(/' ||| NO AREA FOR THIS INFORMATION |||') 00000383 453 FORMAT(/1X,55('='),' PART',I3,1X,55('=')/) 00000384 22 FORMAT(4H |,80A1,22H| ??? DATA END ???/) 00000385 RETURN 00000386 END 00000387 C ********************* 00000001 C * SUBROUTINE DATIN * CONTROL INFORMATION INPUT 00000002 C ********************* 00000003 SUBROUTINE DATIN(MINB,MAXB,MAXBP,MUP,KPRU1) 00000004 COMMON /ENER2/KNSP,KSIG,KTES,KPRIN,KPRU2/NDA/MARKA(67),NKUS 00000005 1/PART/IP1,IP2,GN(2),JY(6),AY(6)/ENER/ESIG(9),SIG(9),SIG1(4),SIGBA 00000006 2,SIGP(5),FI,DSHI,DSBI,DSBP,TSHI,TSHE,TSP, TSI,TELP,TSBI,TSBE 00000007 3,TSPS,TSPB,TSBP,TSBPE,TSPBP,PH,YMY,E,T,TSPSP,SH(20),SB(20),SBP(20)00000008 4,ZNEL(20),IZ(40),YJ(1650)/DT/XXX(140),VW(40),NVV 00000009 * /SURF1/SETPOT(4),NIT,ITUR1,MXB1,IS1,IPR1,MU100000010 5,NIT2,ITUR2,MXB2,IS2,IPR2,MU2,PROB1,QRAT,DSTAB,PA1,PB1,PSTR1,HSF1,00000011 6HEL1,HTW1,BSF1,BEL1,BTW1,PROB2,QRAT2,DSTAB2,PA2,PB2,PSTR2,HSF2, 00000012 7HEL2,HTW2,BSF2,BEL2,BTW2,HIN1,SIN1,TIN1,EU1,HIN2,SIN2,TIN2,EU2 00000013 DIMENSION VAD1(6,4),SZG(4),SUG(2),SIGP1(3), V2(1800),IKN(5) 00000014 *,IGN(12),EXG(8),VADA(24) 00000015 EQUIVALENCE (VAD1(1),VADA(1),PROB1) 00000016 *,(IGN(1),NIT),(EXG(1),HIN1),(SIG(2),SZG(1)),(SIG(8),SUG(1)) 00000017 *,(IKN(1),KNSP) 00000018 *,(SIGP(3),SIGP1(1)),(SIGAL,SIG(1)),(SIGTN,SIG(7)), 00000019 *(V2(1),ESIG(1)), (SIGLP,SIG(6)),(SIGBP,SIGP(2)) 00000020 C=== SEE 'MANUAL FOR USER', MAINPROGRAM 00000021 323 FORMAT(4F10.2) 00000022 READ 323,PH,YMY,E,T 00000023 PRINT 304,PH,YMY,E,T 00000024 304 FORMAT(4H |,4F10.2,40X,15H| PH, YMY, E, T,19X,3H* 1 00000025 */4H |,80X,1H|,36(1H-)) 00000026 READ 299,IP1,IP2,NIT,MUP,MINB,MAXB,MAXBP,QRAT,DSTAB 00000027 299 FORMAT(7I3,2F5.1) 00000028 PRINT 300,IP1,IP2,NIT,MUP,MINB,MAXB,MAXBP,QRAT,DSTAB 00000029 300 FORMAT(4H |,7I3,2F5.1,49X,37H| IP1,IP2,NIT,MUP,MINB,MAXB,MAXBP,*00000030 1 2/4H |,80X,12H| QRAT,DSTAB,22X,1H*/4H |,80X,1H|,36(1H-)) 00000031 301 FORMAT(6I1,4X,3(2X,F8.7),2(2X,F3.2)) 00000032 READ 301,KTES,KSIG,KNSP,KPRIN,KPRU1,KPRU2 00000033 *,SIGAL,SIGTN,SIGLP,SIGBA,SIGBP 00000034 PRINT 302,KTES,KSIG,KNSP,KPRIN,KPRU1,KPRU2 00000035 *,SIGAL,SIGTN,SIGLP,SIGBA,SIGBP 00000036 302 FORMAT(4H |,6I1,4X,3F10.7,2F5.2,30X,37H| KTES,KSIG,KNSP,KPRIN,KP00000037 *RU1,KPRU2* 3 /4H |,80X,37H| SIGAL,SIGTN,SIGLP,SIGBA,SIGBP * 00000038 1/4H |,80X,1H|,36(1H-)) 00000039 IF(KTES.EQ.0) GOTO 320 00000040 306 FORMAT(20F4.2) 00000041 READ 306,SH,SB,SBP 00000042 307 FORMAT(20F4.1) 00000043 READ 307,ZNEL 00000044 PRINT 308,SH,SB,SBP,ZNEL 00000045 308 FORMAT(4H |,20F4.2,14H| SH(L),L=1,20,18X,5H* 3A1/ 00000046 1 4H |,20F4.2,14H| SB(L),L=1,20,18X,5H* 3A2/ 00000047 2 4H |,20F4.2,15H| SBP(L),L=1,20,17X,5H* 3A3/ 00000048 3 4H |,20F4.1,16H| ZNEL(L),L=1,20,16X,5H* 3A4/ 00000049 4 4H |,80X,1H|,36(1H-)) 00000050 320 IF(KSIG.EQ.0) GOTO 321 00000051 309 FORMAT(4(2X,F8.7),2(2X,F3.2)) 00000052 1309 FORMAT(3(2X,F8.7),2X,F3.2) 00000053 2309 FORMAT(2X,F8.7,1X,F4.2,7X,F8.7,1X,F4.2) 00000054 READ 309,SZG,SUG 00000055 READ 1309,SIG1 00000056 READ 2309,SIGP(1),SIGP1 00000057 PRINT 310,SZG,SUG,SIG1,SIGP(1),SIGP1 00000058 310 FORMAT(4H |,4F10.7,2F5.2,30X,18H| SIG(I),I=2-5,8-9,14X,5H* 3B1/ 00000059 1 4H |,3F10.7, F5.2,45X,15H| SIG1(I),I=1,4,17X,5H* 3B2/ 00000060 2 4H |,F10.7,F5.2,5X,F10.7,F5.2,45X,20H| SIGP(I),I=1,3,4,5 00000061 3,12X,5H* 3B3/ 4H |,80X,1H|,36(1H-)) 00000062 321 IF(KNSP.EQ.0) GOTO 322 00000063 303 FORMAT(10F6.2) 00000064 READ 303,DSHI,TSHI,TSHE,TSP,TSI,TELP 00000065 READ 303,DSBI,TSBI,TSBE,TSPS,TSPB 00000066 READ 303,DSBP,TSBP,TSBPE,TSPSP,TSPBP 00000067 PRINT 305,DSHI,TSHI,TSHE,TSP,TSI,TELP,DSBI,TSBI,TSBE,TSPS,TSPB, 00000068 1DSBP,TSBP,TSBPE,TSPSP,TSPBP 00000069 305 FORMAT(4H |,6F6.2, 44X,37H| DSHI,TSHI,TSHE,TSP,TSI,TELP * 3C1/00000070 1 4H |,5F6.2, 50X,37H| DSBI,TSBI,TSBE,TSPS,TSPB * 3C2/00000071 2 4H |,5F6.2, 50X,37H| DSBP,TSBP,TSBPE,TSPSP,TSPBP * 3C3/00000072 3 4H |,80X,1H|,36(1H-)) 00000073 322 CONTINUE 00000074 312 FORMAT(2I3,4F6.1) 00000075 READ 312,ITUR1,MXB1,PROB1,PA1,PB1,PSTR1 00000076 PRINT 313,ITUR1,MXB1,PROB1,PA1,PB1,PSTR1 00000077 313 FORMAT(4H |,2I3,4F6.1,50X,32H| ITUR1,MXB1,PROB1,PA1,PB1,PSTR1, 00000078 *5H * 4 /4H |,80X,1H|36(1H-)) 00000079 DO 319 I=1,6 00000080 IF(VAD1(I,1).EQ.0.) VAD1(I,1)=VAD1(I,2) 00000081 319 CONTINUE 00000082 314 FORMAT(2I3,2X,9F6.3) 00000083 READ 314,IS1,IPR1,HIN1,SIN1,TIN1,HSF1,HEL1,HTW1,BSF1,BEL1,BTW1 00000084 PRINT 315,IS1,IPR1,HIN1,SIN1,TIN1,HSF1,HEL1,HTW1,BSF1,BEL1,BTW1 00000085 315 FORMAT(4H |,2I3,2X,9F6.3, 18X,26H| IS1,IPR1,HIN1,SIN1,TIN1,00000086 *,8X,3H* 5/4H |,80X,37H| HSF1,HEL1,HTW1,BSF1,BEL1,BTW1 * 00000087 */4H |,80X,1H|36(1H-)) 00000088 READ 312,ITUR2,MXB2,PROB2,PA2,PB2,PSTR2 00000089 PRINT 317,ITUR2,MXB2,PROB2,PA2,PB2,PSTR2 00000090 317 FORMAT(4H |,2I3,4F6.1,50X,32H| ITUR2,MXB2,PROB2,PA2,PB2,PSTR2, 00000091 *5H * 6/4H |,80X,1H|36(1H-)) 00000092 DO 326 I=1,6 00000093 IF(VAD1(I,3).EQ.0.) VAD1(I,3)=VAD1(I,4) 00000094 326 CONTINUE 00000095 READ 314,IS2,IPR2,HIN2,SIN2,TIN2,HSF2,HEL2,HTW2,BSF2,BEL2,BTW2 00000096 PRINT 318,IS2,IPR2,HIN2,SIN2,TIN2,HSF2,HEL2,HTW2,BSF2,BEL2,BTW2 00000097 318 FORMAT(4H |,2I3,2X,9F6.3, 18X,26H| IS2,IPR2,HIN2,SIN2,TIN2,00000098 *,8X,3H* 7/4H |,80X,37H| HSF2,HEL2,HTW2,BSF2,BEL2,BTW2 * 00000099 */4H |,80X,1H|36(1H-)) 00000100 2000 FORMAT(2I2,2I3,I4,I6,2X,2F8.6,2F8.4,2F8.2) 00000101 READ 2000,JY,AY 00000102 PRINT 2001,JY,AY 00000103 2001 FORMAT(4H |2I2,2I3,I4,I6,2X,2F8.6,2F8.4,2F8.2,10X,37H| (JY(I),I=00000104 *1,6), (AY(I),I=1,6) * 8) 00000105 C===== 00000106 PRINT 324 00000107 324 FORMAT(3X,1H|,16(5H----+),1H|,36(1H-)/ 00000108 * 3X,1H|,80(1H=),37H|=INPUT & INPUT DATA & CONTROL=======/) 00000109 NIT=IABS(NIT+1) 00000110 QRAT2=QRAT 00000111 DSTAB2=DSTAB 00000112 IP1=IABS(IP1+1) 00000113 IF(IP2.LE.0) IP2=0 00000114 IF(IP1.EQ.0.AND.IP2.EQ.0) IP2=1 00000115 IF(IP1.EQ.0.AND.IP2.EQ.0.AND.NIT.EQ.0) NIT=1 00000116 NIT2=NIT 00000117 IF(PSTR1.LT.0.) PSTR1=0. 00000118 IF(PSTR2.LT.-100.) PSTR2=-100. 00000119 IF(NKUS.EQ.1) GOTO 1009 00000120 IZ(1)=MINB 00000121 IZ(2)=MAXB 00000122 IZ(3)=MAXBP 00000123 IZ(4)=MUP 00000124 IZ(5)=KPRU1 00000125 DO 1000 K=1,5 00000126 I=K+5 00000127 1000 IZ(I)=IKN(K) 00000128 IZ(11)=IP1 00000129 IZ(12)=IP2 00000130 DO 1001 K=1,6 00000131 I=K+12 00000132 1001 IZ(I)=JY(K) 00000133 DO 1002 K=1,12 00000134 I=K+18 00000135 1002 IZ(I)=IGN(K) 00000136 IZ(31)=NVV 00000137 DO 1003 K=1,6 00000138 I=K 00000139 1003 YJ(I)=AY(K) 00000140 DO 1004 K=1,4 00000141 I=K+6 00000142 1004 YJ(I)=SETPOT(K) 00000143 DO 1005 K=1,24 00000144 I=K+10 00000145 1005 YJ(I)=VADA(K) 00000146 DO 1006 K=1,8 00000147 I=K+34 00000148 1006 YJ(I)=EXG(K) 00000149 DO 1007 K=1,40 00000150 I=K+42 00000151 1007 YJ(I)=VW(K) 00000152 K=NKUS 00000153 CALL RITE(K,1800,V2) 00000154 1009 RETURN 00000155 ENTRY DATIN1(MINB,MAXB,MAXBP,MUP,KPRU1) 00000156 K=NKUS 00000157 CALL RID1(K,1800,V2) 00000158 MINB=IZ(1) 00000159 MAXB=IZ(2) 00000160 MAXBP=IZ(3) 00000161 MUP=IZ(4) 00000162 KPRU1=IZ(5) 00000163 DO 1010 K=1,5 00000164 I=K+5 00000165 1010 IKN(K)=IZ(I) 00000166 IP1=IZ(11) 00000167 IP2=IZ(12) 00000168 DO 1011 K=1,6 00000169 I=K+12 00000170 1011 JY(K)=IZ(I) 00000171 DO 1012 K=1,12 00000172 I=K+18 00000173 1012 IGN(K)=IZ(I) 00000174 NVV=IZ(31) 00000175 DO 1013 K=1,6 00000176 I=K 00000177 1013 AY(K)=YJ(I) 00000178 DO 1014 K=1,4 00000179 I=K+6 00000180 1014 SETPOT(K)=YJ(I) 00000181 DO 1015 K=1,24 00000182 I=K+10 00000183 1015 VADA(K)=YJ(I) 00000184 DO 1016 K=1,8 00000185 I=K+34 00000186 1016 EXG(K)=YJ(I) 00000187 DO 1017 K=1,40 00000188 I=K+42 00000189 1017 VW(K)=YJ(I) 00000190 RETURN 00000191 END 00000192