218 - 



30050 



30100 



30150 



30200 



30250 



30300 



30350 



30400 



304 50 



30500 



30550 



3060O 



30650 



30700 



30750 



30800 



30850 



30900 



30950 



31000 



31050 



31 100 



31 150 



3 1200 



31250 



3 1300 



31350 



3 1400 



31450 



31500 



31550 



31600 



31650 



31700 



31750 



3 1800 



31850 



3 1900 



31950 



32000 



32050 



32100 



32150 



32200 



32250 



32300 



32350 



32400 



32450 



32500 



32550 



32600 



32650 



32700 



32750 



32800 



32850 



32900 



32950 



33000 



4001 



wvnr(1) = wvnr(1)/12.0 

 tanab = atan2(ab1( 1 ) . abs( i ) ) 

 IF(DP1(I).NE.0.0) 



TAN1 = atan2(D1I(I) ,DP1(I )) 

 1f (dpK 1 ) .eq.0.0) tani =0.0 

 1f(dp2(1 ).ne.O.O) 



tan2 = atan2(d21( i ) ,dp2( 1 )) 

 If (d21( 1 ) .eq.O.O) tan2 = . OO 

 if (dp3( 1 ) .ne.O.O) 



tan3 = atan2(d31 ( 1 ) .dp3( 1 )) 

 if (dp3( 1 ) .eq.O.O) tan3 =0.0 

 if (dp4( 1 ) .ne.O.O) 



tan4 = atan2(d4i ( 1 ) ,dp4( i ) ) 

 if (d41( i ) .eq.0.00) tan4 =0.00 

 if (dpxx( 1 ) .ne.O.O) 



tanxx = atan2(dxxi ( i ) .dpxx( i ) ) 

 if (dpxx( i ) .eq.O.O) tanxx = O.O 

 if (dpyy( 1 ) .ne.O.O) 



tanyy = atan2(dyy i ( 1 ) ,dpyy( i ) ) 

 if (dyyi ( i ) .eq.O.O) tanyy =0.0 



AO(i) = ab(i) / (2.0 * pi) 



tanM = tan1 - tanAB 

 tanN = tan3 - tanAB 

 A1d1(i) = -1.0/(pi * wvnr(i)) ♦ sqrt(ab(i)) * 



sqrt(d1(i)) • sin(tanM) 

 A1d3(i) = -1.0/(pi * wvnr(1)) ♦ sqrt(ab(i)) • 

 sqrt(d3(i)) * sin(tanN) 



A2d1d2{1) = -(d2(i) - d1( i ))/(wvnr(1 )**2*pi) 



A2d1d4(i) = -(d4(i) - d1(i))/(wvnr(i)**2*pi) 



A2d3d2(i) = -(d2(i) - d3( i ) )/(wvnr( i )**2*pi ) 



A2d3d4(i) = -(d4(l) - d3(i ))/(wvnr( i )**2*pi ) 



GO TO 4001 



tanO = tanXX - tani 

 tanP = tanXX - tan3 

 tanO = tani - tanAB 

 tanR = tanS - tanAB 

 A3d1(i) = 4.0/(pi * wvnr(i)**3) 



* (-sqrt(d1( i )) ♦ sqrt(dxx(i)) ♦ sin(tanO) 



+ 0.75 * wvnr(i)**2 * sqrt(ab(i)) * sqrt(d1(i)) 



* sin(tanO)) 



A3d3(i) = 4.0/(pi * wvnr(i)**3) 



* (-sqrt(d3( i )) • sqrt(dxx(i)) * sin(tanP) 



+ 0.75 * wvnr(i)**2 * sqrt(ab(i)) * sqrt(d3(i)) 



* sin(tanR)) 

 Cont inue 



tanS = tan2 - tanAB 

 tanT = tan4 - tanAB 

 B1d2(i) = -1.0/(pi * wvnr(i)) 



* sqrt(d2(i)) ♦ sln(tanS) 

 B1d4(1) = -1.0/(pi * wvnr(i)) 



* sqrt(d4(i)) * sin(tanT) 



tanU = tan2 - tani 

 tanV = tan2 - tan3 

 tanW = tan4 - tani 



* sqrt(ab( 1 )) 



• sqrt(ab( 1 ) ) 



