lSlNA/SINA/,6X,eMVARI*tMCF I 

 PPtNT JO 1 I SFF , Vr, VQ, VX, VZ , VR 



301 FCBM*T (6F ie,6/ ) 

 PRINT 310 



31 F0RMAT(26X,2HAY,13X,2HAQ,l3X,2HAXfl3X,2MAZ, 13X,2HrA ) 

 PRINT 102 , AY, AO, AX , AZ I YA 



302 FOfiMAT( 15X, 5Flf., 6/ ) 

 PRINT 308 



308 F0PMAT(4 1X, 3MFLV,12Xt 3HFMV, I2X,3HFDV) 

 PRINT 30 Jf FLV,FMV,FOV 



303 FOWVIAT ( 30X, JFie.e//// ) ■ 

 PRINT lOft 



104 FOfiMATC lOX, 4H CLV,10X,5h CLVA,10X,5H CLVU,11X,5H CMV ,10X,5H COV , 

 1 1 1 X,2H K , I2X,4H PHI ) 

 PRINT 305,CLV,CLVA,CLVLl,CMV|COV,XK,PHl 

 305 F0RNIAT(5F15.3,F15,4,F15.2////) 



PRINT 309 

 30 9 FORMAT (38X, 3HFAVG(L8) , 7X,8HFMAX<LB ) i 7Xt BHFMI N( L9 1 ) 



PRINT 304 ,SF ,FMAX,FMIN 

 3C4 FORMAT ( 30X,F 1£. 6 ,2F15.S///'/ ) 



PUNCH 387,CL,DI A,ANGLE,T,H,XL,U,CLV,CLVA,CLVU,PHI ,XtC,CMV,CDV 

 387 F0RMATCF4,3 |F5. 3 ,F4. i F5. 2 , F5. 3, F6. 2 ,F6,4 , F5. 2 i 2F6 . 2 t F7 ,2i F 7,4 , F6. 

 12,F8.2) 

 PLOT ORIGINAL DATA AND RESULTS FCR COMPARISON 



PRINT 26 

 26 F0RMAT{4X,7H FP(LB),3X,7H FV(L8),2XiaH RES(L9)/) 

 DO 31 L=l,101 

 G(L)=H 



31 CONTINUE 



00 25 I=li40 



FV( I ) = AX*X( I )*AQ*Q( I )*AY*Y( I )»AZ*Z( I > i-SF 

 ReS( I )=FP( I )-FV( I ) 

 25 CONTINUE 



e»ABS(FP( 1 1 ) 

 DO 59 1=2,40 

 CrA9S(FV( I ) ) 

 A«ABS(FP( I ) ) 

 IF(C-B)57,57,56 



56 e«c 



57 IF( A-9)59,59,58 



58 B=A 



59 CONTINUE 

 BB=49./8 

 IF (63-80 0. ) 61 ,43,4 3 



61 IF(B8-400, 162,44,44 



62 IF(BB-200. >63, 45,45 



63 IF(8B-100. ) 64,46 ,46 



64 BB=50. 

 GO TO 47 



43 BB=800. 

 GO TO 4 7 



44 BB=400. 

 GO TO 47 



AS BB=200. 

 GO TO 47 



46 88=100. 



47 CONTINUE 



00 32 1=1,40 

 G(51 ) = 1HI 

 J = S1.*^B*FP( I ) 

 K=51. ♦BB«FV( I ) 

 L = 5l.*r3B*RES( I > 

 G( J )=1H* 

 G<K)=1H* 

 G(L1=1H. 



PRINT 100,FP( I ) tFV< II,.RES^< I),G 

 100 FORMAT( lHZ,3F10.5,101AI/» 

 G( J ) = 1H 

 G(K) =IH 

 6(L)=1H 



32 CONTINUE 

 GO TO a 



10 CONTINUE 

 END 



SUBRCUTJNE XAVFL ( Tf D, XL ) 

 B«32.2*T*T/6.283185 

 TP0=6.283185«O 



OEEp"-ATER Initial estimate for wavelength 



2 XL*B 



SHALLQW^KATER IN I TI Al. ESTIMATE FOR WAVELEMSTM 



3 XL=T*SQHT(0*32.2) 



4 XLX=XL 

 XL=8»TANH( TPD/XLX) 

 IFIABS (XLX-XL •-•005)5t4»4 



5 RETURN 

 END 



143 



