64 



100 



CONTINUE 



65 





N=U-1 



66 





LINE=N-1 



67 



C 

 C 

 C 



IF (ISKIP.LE.O) GO TO 130 





SKIP LEADING CARD IMAGES ' 



68 



NSTART=ISKIP-t-l 



69 





00 120 N=NSTART»LINE 



70 





NN=N-ISKIP 



71 





NK.BR<NN)=NMbR(N> 



72 





TIME(NN)=TU'E<N> 



73 





DO 110 1=1.2 



74 



110 



D(I»NN)=D(I»N> 



75 





RAD(NN)=RAD(N> 



76 



120 



FATH(NN)=FATH(N> 



77 





LINE=LINE"ISKIP 



78 



130 



MST0P=LINE-1 



PROGRAM — RAPLOT III MODIFIED FOR IBM-7094 WITH SC-4060 CRT 



15H 

 1S9 

 160 

 161 

 162 

 :0NTAIN BAD OATA. 163 



164 

 165 

 166 

 167 

 168 

 169 

 170 

 171 

 172 

 173 

 174 

 175 

 C 176 



C CHECK DISTANCES TO RADAR BEACONS FOR ERRORS. IF DISTANCE/TIME 177 

 C FOR SUCCESSIVE BEACON RANGES INDICATE A SHIP SPEED .GT. 6 KNOTS 178 

 C (3.09 METERS/SEC) .RANGE IS IN ERROR. 179 



C 180 



79 DO 240 1=1.2 181 



80 DO 170 M=l.MSTOP 182 



81 IF (0(1. M)) 170.170»1<»0 183 



82 1*0 NSTART=M+1 184 



83 DO 160 N=NSTART.LINE 185 



84 IF (At3S(D(I»N>-D(I»M>>-(TIME<N>-TIME(M>)*3. 08865) 170.170.150 186 



85 150 D(I.N>=-1. 187 



86 160 CONTINUE 18ft 

 67 170 CONTINUE 189 

 88 DO 230 M=l.MSTOP 190 

 69 IF (0(1. Ml) 180.180.230 191 



90 180 NSTART=M 192 

 C 193 

 C CORRECT ERRONEOUS BEACON RANGES BY LINEAR INTERPOLATION (ON TIME) 19<» 

 C BETWEEN NON-ERRONEOUS RANGES. 195 

 C 196 



91 DO 200 N=NSTART.LINE 197 



92 IF (0(1. N)) 200.200.190 19S 



93 190 NSTOP=N 199 



94 GO TO 210 200 



95 200 CONTINUt ^01 



96 210 0TIME=TIME(NST0P)-TIME(NSTART-1) 202 



97 0D1=D(I.NST0P)-D(I.NSTART-1> 203 



98 N=NSTART 204 



99 220 0(I>N)=D(I.NSTART-1)+DD1*(TIME(N)-TIME(NSTART-1))/0TIME 205 



100 N=N+1 206 



101 IF (N-NSTOP) 220.230.230 207 



102 230 CONTINUE 208 



103 2"»0 CONTINUE 209 

 C 210 

 C COMPUTE POSITION OF SHIP FROM DISTANCES FROM THE TWO BEACONS 211 

 C 212 



104 LAG=0 213 



105 DO 270 N=1.LINE 214 

 C 215 

 C MAKE CONSTANT CORRECTION FOR CUBIC AUTOTAPE INTERROGATOR 216 

 C AND CONVERT TO FEET 217 

 C 218 



106 DFT1=(D(1.N)+C0RR1>*3. 28083 219 



107 DFT2=(D(2»N)+C0RR2>*3. 28083 220- 



108 DX1=(SQUSTB+DFT1*DFT1-DFT2»DFT2)/(2.*0ISTB) 221 



109 DY1=DFT1*DFT1-DX1»DX1 222 

 C 223 

 C CHECK FOR IMAGINARY ROOT* 224 

 C 225 



110 IF (UY1) 250.250.260 226 



111 250 N0RTH(N)=-1. 227 



112 LAG=LAG+1 228 



113 GO TO 270 229 



114 260 DY1=-SQRT(DY1) 230 

 C 231 

 C ROTATE COORDINATES AND TRANSLATE INTO CALIFORNIA LAMBERT COORDINAT 232 

 C SYSTEM 233 

 C 234 



115 EAST(N)=DX1*C0SINE-DY1*SINE+BEAC1E 235 



116 N0RTH(N)=DX1*SINE+DY1*C0SINE+BEAC1N 236 



50 



