44 DAVE=DAVE+X(I)-X(J) 

 0AVE=0.3*(ABS(DAVE/AN)) 

 CALL S0RTY(AY,AX,AZ,Y,X,Z,N,1,1,DAVE) 

 IF(JDIR.LT.l) 60 TO 92 

 DO 93 I=1,N 

 Z(I)=AZ(I) 

 X(I)=AY(I) 



93 Y(I)=AX(I) 

 60 TO 94 



92 DO 95 1=1, N 



Z(I)=AZ(I) 



X(I)=AX(I) 

 95 Y(I)=AY(I) 



94 CONTINUE 



C RESET 0RI6IN TO FIRST MAPPED POINT 



BLAT=BLAT+Y(1) 



BL0N6=BL0N6-X(1) 



DX=X(1) 



0Y=Y(1) 



DO 71 I=1,N 



X(I)=X(I)-DX 

 71 Y(I)=Y(I)-DY 

 C NOW INTERPOLATE ALONG TRACK AT DESIRED GRID SPACIN6 

 C COMPUTE APPROX. LENGTH OF LONGITUDE FOR THIS TRACK 

 C USE APPROX. MIDLATITUDE AS BASIS AND TABLE 6(B0WDITCH) 



NN=N/2 



AL=(Y(NN)+BLAT)*RAD 



A=(111415. 13*C0S (AL )-94. 55*C0S ( 3.*AL )+. 012*C0S ( 5.*AL ) )/60. 

 C UNITS OF A ARE METERS/MINUTE OF LONGITUDE 

 C COMPUTE APPROX. NAUTICAL MILES/MINUTE OF LONGITUDE 



19 B=A*(1. 0/1852.0) 

 C CONVERT X COORDINATE OF MAPPED POSITION TO MILES AND COMPUTE DISTANCE 

 C DOWN TRACK. 



00 25 I=1,N 



25 DIST( I )=SQRT( (X( I )*B)**2+Y( I )**2) 

 IF(ABS(X(N)).LT. 0.00001) X(N)=0. 00001 

 THETA=ATAN2(Y(N ) , (X(N )*B) ) 



CALL GINT(GRID,N,0.0,DIST(N),ICT,DIST,Z,AX,AZ) 

 C STORE NEW POSIT OF INTERPOLATED PT.IN AX (L0N6,AY(LAT) DECIMAL DE6. 

 C NEW INTERPOLATED VALUE OF Z IS STORED IN AZ 



DO 26 I=1,ICT 



AY(I)=(AX(I)*SIN(THETA)+BLAT)/60.0 



26 AX(I)=(BL0NG-(AX(I)*C0S(THETA)/B))/60.0 

 RETURN 



END 



161 



