1 



2 





3 





4 





5 





6 





7 





8 





9 





10 





11 





12 



4 



13 



C-- 



14 



c** 



15 





16 





17 



10 



18 





19 





20 





21 





22 





23 





24 





25 



30 



26 



c** 



27 





28 





29 



C*t 



30 



c** 



31 





32 





33 



C** 1 



34 





35 





36 



C** 1 



37 





38 





39 



C** 1 



40 





41 





42 



C** 1 



43 





44 





45 



128 



46 





47 





48 



c** 



49 





50 





51 



1 19 



52 



21 



53 





54 





55 



c** : 



56 



107 



PROGRAM MFOLLO 



IMPLICIT DOUBLE PRECISION (A-H.Q-Z) 



COMMON /INPUT/ Z(1 0) , N .OMEGA , V , VI ,GCU( 10) ,GSO( 10) ,CAY( 10) , LAMBDA, L 

 1 AMBDI ,G(10),RHO(10),GI(10) ,GSOI( 10) ,CAYI(10) 

 COMMON /DETMNT/ A ( 2 1 , 4 ) ,0 ( 21 , 4 ) 



REAL INCA, INCB, INCC, INCD, INCE, LAMBDA, LAMBDI 

 DIMENSION T(4), PV(4), W(8), WI(8), CB(10), CBI(IO), 0(10), 



1 CAY S0(10), GAMMA(IO), DPK(IO), GCUI(IO), CI(10), OR ( 1 0) , PVI (4) 



2 , CAYSOI(IO) , SR(4) , SI(4) 

 CHNG = 1 . / 8192. 

 CHNGI = 0. 

 CONTINUE 



;;0 - TOTAL. STEP LIMIT, K1, K2 PRINT KEYS, K3 = i r\ttpS SAivit kkumlc 

 FOR NEXT RUN. 



READ 10, KO, K1 , K2, K6, K3, TLIM, BLIM, RATIO, EX 



PRINT 10, KO, K1, K2, K6, K3, TLIM, BLIM, RATIO, EX 



FORMAT (514, 4E10. 1 ) 



IF (TLIM .EO. 0.) TLIM = 1 . E-5 



IF (BLIM .EQ. 0.) BLIM = 1 . E-2 



IF (EX .EO. 0. ) EX = 28. 



RLIM = 10.+*EX 



IF (RATIO .EO. 0. ) RATIO = 2. 



IF (KO .EO. 0) KO = 300 



IF (K3 .NE. 0) GO TO 128 



READ 1240, N,FREQ ,ATTEN 

 STOP IF N = 0. THIS IS THE ONLY PROGRAMED STOP. 



IF (N.EO.O) GO TO 1200 



PRINT 1250, N,FREO 

 PARAMETERS READ IN BELOW ARE THOSE AT THE TOP OF EACH LAYER. 

 READ IN VELOCITIES. 



READ 1260, (C( I) , 1 = 1 ,N) 



PRINT 1280, (C( I J , 1 = 1 ,N) 

 READ IN DEPTHS. 



READ 1260, (Z( I ) , 1=1 ,N) 



PRINT 1 280, ( Z( I ) , 1=1 ,N) 

 READ IN GRADIENTS 



READ 1260, (GA;«1MA( I ) , 1 = 1 ,N) 



PRINT 1280, (GAMMA( I ) , 1 = 1 ,N) 

 READ IN ATTENUATION FACTOR IN LOSS PER KILOMETER. 



READ 1260, (DPK( I ) , 1 = 1 ,N) 



PRINT 1280, (DPK( I ) , 1 = 1 ,N) 

 READ IN DENSITIES (BLANK INPUT IMPLIES SEA WATER DENSITY). 



READ 1260, (RHO( I ) , 1 = 1 ,N) 



PRINT 1280, ( RHO( I ) , 1 = 1 ,N) 



CONTINUE 



NUMBER = 1 



dX = 

 NX = VARIABLE, NY = LAYER NUMBER, NZ = CONTINUITY 



READ 119, NX, NY, NZ, PK, VALL.DP, V, VI, STEP, STEPI 



PRINT 21, NX, NY, NZ, PK, VALL.DP, V, VI, STEP, STEPI 



FORMAT (312, 4X, 7010.2) 



FORMAT (10H VARIABLE ,12, 1 0H LAYER NO , 1 2 , 1 2H CONTIWIITtv 

 * 12, / 7G15.5) 



PK = PK - DP 

 START NEW CYCLE BY INCREMENTING VARIABLE. 



PK = PK + DP 



94 



