219 



33050 

 33100 

 33150 

 33200 

 33250 

 33300 

 33350 

 33400 

 33450 

 33500 

 33550 

 33600 

 33650 

 33700 

 33750 

 33800 

 33850 

 33900 

 33950 

 34000 

 34050 

 34100 



34 150 

 34200 

 34250 

 34300 

 34350 

 34400 

 34450 

 34500 

 34550 

 34600 

 34650 

 34700 

 34750 

 34800 

 34850 

 34900 

 34950 

 35000 

 35050 



35 100 

 35150 

 35200 

 35250 

 35300 

 35350 

 35400 

 35450 

 3550O 

 35550 

 35600 

 35650 

 35700 

 35750 

 35800 

 35850 

 35900 

 35950 

 36000 



tanX = tan4 - tan3 

 B2d2d1(1) = 2.0/(pi * wvnr(i)**2) 



* sqrt(d1(1)) * sqrt(d2(1)) *(cos( tanu)+s in( tanU) ) 

 B2d2d3(1) = 2.0/(pi * wvnr(i)**2) 



* sqrt(d3(i)) * sqrt(d2(i)) *(cos( tanV )+s 1 n( tanV ) ) 

 B2d4di(i) = 2.0/(pi * wvnr(i)**2) 



* sqrt(d1(1)) * sqrt(d4(i)) * (cos( tanW)+s in( tanW) ) 

 B2d4d3(1) = 2.0/(pi ♦ wvnr(1)**2) 



* sqrt(d3(i)) * sqrt(d4(1)) * ( cos( tanX )+s i n( tanX ) ) 



D TO 4002 

 tany2 = 

 tanz2 = 

 tany4 = 

 tanz4 = 

 B3d2(i) = 



tanyy - tan2 

 tan2 - tanab 

 tanyy - tan4 

 tan4 - tanab 

 -4 .0/(p1 * wvnr( i )^ 



►3) 



sqrt(dyy(1)) *sin(tany2) 



2 * sqrt(ab(i)) * sqrt(d2(i)) 



1 * (-sqrt(d2( i )) 



2 + 0.75 * wvnr( i ) 



3 * s1n(tan22)) 



B3d4(1) = -4.0/(p1 * wvnr(1)**3) 



1 * (-sqrt(d4( 1 ) ) * sqrt(dyy(1)) * sin(tany4) 



2 + 0.75 * wvnr(1)**2 * sqrt(ab(1)) * sqrt(d4(1)) 



3 * s1n(tanz4)) 

 I Continue 



) CONTINUE 



DO 511 1=2, icut 

 C IF ALL GAGES ARE WORKING. USE THIS LOOP: 

 GO TO 7 



A1(i) = (A1d1(1) + A1d3(1))/2.0 

 A2(1) = (A2d1d2(1) + A2d1d4(i) + A2d3d2(i) 

 1 + A2d3d4( 1 ) ) / 4.0 



A3(1) = (A3d1(1) + A3d3(i))/2.0 

 B1(i) = (B1d2(1) + B1d4(1))/2.0 

 B2(i) = (B2d2d1(1) + B2d2d3( i ) + B2d4d1(i) 

 1 + B2d4d3( 1 ))/ 4.0 



B3(i) = (B3d2(1) + B3d4(i)) / 2.0 

 7 Continue 



C IF ONLY TWO DIFFERENTIAL CHANNELS ARE WORKING. USE: 

 C GO TO 3077 



A1( i ) = A1d3( 1 ) 



A2( 1 ) = A2d3d2( 1 ) 



A3( 1 ) = A3d3( 1 ) 



B1( i ) = B1d2( i ) 



B2( i ) = B2d2d3( i ) 



83(1) = B3d2( 1 ) 

 3077 Continue 



Calculate MEAN direction for n=1,2,3 

 THTAKi) = atan2(B1( i ).A1( i )) 

 THTA2(i) = atan2(B2( i ) ,A2( i )) 

 THTA3(i) = atan2(B3( i ).A3( i )) 



Cent inue 



180. /3. 14159 

 180. /3. 14159 

 180. /3. 14159 



