IF(FP(I>.GT.FMAXI902,901 

 90a F!4AX«FP(I) 



901 IF(FP( I >,LT,FMIN>903, 12 

 903 FMIN=FP( I ) 



12 SF = SF«-FP{I) 



13 CONTINUE 

 SF-SF/XJ 



CALCULATE SUMS OF SQUARES AND PRODUCTS 

 SFFxO. 

 SFX=0. 

 SFY=0. 

 SFO'O* 

 SF2=0. 



0T»3. 141 5926536/ (25. ♦T » 

 A--OT 



00 15 1=1, J 

 AsA^-OT 

 Xttl«COS (Al 



e«2.«A 



V( t )3C0S(B) 



0( I |zSIN(BI 



C'SINIAI 



Z(II-C«ABS (CI 



F( I »aFP( I l-SF 



SFF«SFF»F( I )*F< I » 



SFX«SFX*F( I )*X( I J 



SFV=SFy+F( I >*Y( I ) 



SFO«SFQ«-F( I )«0( I > 



SFZ«SFZ*F( I )*Z< I ) 

 15 CONTINUE 



SXX».5*XJ 



SVV».5*XJ 



SOO«.S«X J 



iZZ««3^5*XJ 



AX-SFX/SXX 



Ar"SFV/SW 



AOzSFa/SQQ 



AZ«SFZ/SZZ 



VX>AX*SFX 



yVsAV*SFV 



VQ>AO«SFQ 



VZ»AZ«SFZ 



VR=SFF-VX-VY-VZ-VQ 

 : CALCULATE COEFFICIENTS AND PARAMETERS PHI AND K 



PHI»2a.6478 9*ATAN2<AQ,AV) 



IF (PHI.LT.-45. ) 799 9,8999 

 7999 PMI«PHI*-1')0. 

 8999 CONTINUE 



yA»SQRT( AO*AQ*AV*Ay ) 



CLV»2.*YA/FLV 



ANGz( ANGLE'^PI )/lSO. 



CLVA«CLV/COS( ANGI 



CLVU=CLVA/COS( ANG> 



CMV»-AX/FMV 



CDV»-AZ/FOV 



XK>.5-( SF/(CLV*FLV) ) 



147 



