1 C THIS IS THE MAIN PART OF NLAYNM 



2 IMPLICIT DOUBLE PRECISION (A-H,0-Z) 



3 DOUBLE PRECISION LAMBDA , LAMBD I 



4 INTEGER COL 



5 REAL R ATTEN, T RE, RX 



6 DIMENSION L0SPCH(5,25) 



7 COMMON /HAN/ H2R , H2 I , H1 R , HI I ,H2PR ,H2PI , HI PR , HI PI , R 



8 COMMON/INPUT/ Z(12), N, OMEGA. V, VI, C0N(12), GSQ(12), 



9 1 CAY(12), LAMBDA, LAMBDI, G(12) 



10 2,RHO(12), GI(12), G S0I(12), CAYI(12) 



11 COMMON /FXPO/ EXSUM, CNTR, RATI0(25) 



12 COMMON/DETMNT/ A(25,4), 0(25,4) 



13 COMMON/PARTS/ ZT { t 2 ) , ZT I( 1 2 ) , ZB ( 1 2 ) . ZBI ( 1 2 ) 



14 COMMON/REFL/ AF(12,200), AG(12,200), BF(12,200), BG(12,200), 



15 2 EIGEN(350), EIGENI(350), B (25,4), BI(25,4), CB(12), CBI(12), 



16 3 CAYS0(12), CAYS0t(12), NN 



17 DIMENSION D(350) , DI(350), F(100), FI(IOO), HZER02(350), 



18 * DA(350), SRES(350), GAMMAI(12), BLPK(12), 



19 4 HZER2I(350), DPK(12), GCU(12), GCUI(12), CI(12), 



20 3 PHASE V(350), PHASI V(350), UU(2000), UUI(2000) 



21 COMMON /LIMIT/ TLIM, EXPONT , SLIM 



22 DIMENSION LOSS(101 ) 



23 DIMENSION C(12), DEPTH(52), DBLOSS(350), COL(120), 



24 1CONTR(10), EF(2), FMAG(350), FANG(100), 



25 2GAMMA( 12),JSMBL(10) , JCaUNT(5) , JC0U(5) , LEVEL (41 ) ,PLEV(5) ,RLOSS( 100} 



26 3 , RLOS( 101 ) ,RECVRS(51 ) ,TEST(3) , ING( 11 ) 



27 EQUIVALENCE ( FF , E F ( 1) ),( DEPTH ( 1 ), SOURCE ) , (DEPTH(2) ,RECVRS(1)) , 



28 1 (RL0S(2) ,RL0SS(1 )) 



29 COMMON /AHZERO/ HZ EROR ,HZERQI 



30 DATA ( CONTR(I), 1=1,4) / 1 1 0. DO , 95 . DO , 80 . DO , -1 000 . DO/ , 



31 1 (J SMBL(I), 1=1,3) /1H1 ,1H*,1H8/, 



32 *(ING(I), I = 1 , 10)/1H0,1H1 , 1H2, 1H3,1H4,1H5, 1H6,1H7,1H8, 1H9/, 



33 2 (TEST(I), 1=1,3) / . 2D0 , 1 . DO , 5 . DO/ 



34 TLIM = 25. 



35 SLIM = -8.54988 



36 C READ IN PARAMETERS 



37 1 READ 11, K1 , K2, K3, K4, K5, K6 , K7 , K8, K9 



38 C KEYSS 1-DEPT FN PRINTOUT, 2-LOSS PRINTOUT, 3-REFLECTION COEFF PRINTOUT 



39 C 4-CHANGE CONTOURS, 5-CONTINUE MODES 



40 11 FORMAT (1014) 



41 i-HiiNir I J, M, K2, K3, K4, K5, K6 , K7 , K8 , K9 



42 13 FORMAT (6H1KEYS , 1014) 



43 IF (K1 . LT. 7) GO TO 5 



44 READ 11 , M 



45 READ 434, (SRES(I), I = 1,M) 



46 434 FORMAT (5D16.7) 



47 5 IF (K8 .NE. 1 ) GO TO 8 



48 READ 20, T LIM 



49 SLIM = -DCBRT(TLIM**2) 



50 PRINT 30, TLIM, SLIM 



51 8 EXPONT = DEXP((TLIM + TLIM) / 3.) 

 52 

 53 

 54 

 55 

 56 



KB = K6 



+ 



1 











INSUR = 

















IF (K2 , 



. LT 





10) 



GO 



TO 



16 



K2 = K2 



- 



10 











INSUR = 



1 













44 



