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



CONTINUE 237 



IF (LAG.EQ.O) GO TO 330 238 



ASSIGN 330 TO KEY 239 



240 



ELIMINATE DATA SETS FOR WHICH AN 2<«1 



IMAGIN ART FIX WAS OBTAINED 242 



243 



LAG=0 >44 



DO 320 N=1»LINE 2<»5 



IF (NORTH(N)) 290t290i300 246 



LAG=LAG+1 2*7 



GO TO 320 2KB 



NN=N-LAG 249 



NMBR(NN)=NMBR(N) 250 



TIME(NN)=TIME(N> 251 



DO 310 1 = 1.2 252 



D(I>NN)=D(ItN) 253 



RAD(NN)=RAD(N) 254 



EAST(NN)=EAST(N> 255 



NORTH(NN)=NORTH(N> 256 



CONTINUE 257 



LINE=LINE"LAG 258 



GO TO KEY. (330.400) 259 



260 



CHECK NORTH AND EAST COORDINATES FOR 261 



EXTREME VALUES BY CHAUVENET'S CRITERION 262 



263 



NF1AR=-1. 264 



CALL STUEV (NORTH. LINE.NBAR.SDNRTH) 265 



EBAR=-1. 266 



CALL STDEV (EAST. LINE. EBAR.SDEAST) 267 



ALPHA=1.-1./FL0AT(2*LINE) 268 



CALL TINORM (CHVR. ALPHA. *335) 269 



GO TO 340 270 



CHVR=5. 271 



WRITE (6.770) ALPHA 272 



GATE1=E8AR-CHVR*SDEAST 273 



GATE2=EdAK+CHVR*SDEAST 274 



GATEN1=NBAR-CHVR»SDNRTH 275 



GATE.I2=NBAR + CHVR»SDNRTH 276 



LAG=0 277 



DO 390 n=l.LINE 278 



IF (EAST(N)-GATEl) 380.350.350 279 



IF (EAST(N)-GATE2) 360.360.380 280 



IF (NORTH(N)-GATENl) 380.370.370 281 



IF (N0RTH(N)-GATEN2) 390.390.380 282 



NORTH(N)=-l. 283 



LAG=LAG+1 28<» 



CONTINUE 285 



ASSIGN 400 TO KEY 286 



C 287 



C ELIMINATE ANY DATA SETS THAT HAVE AN EXTREME 288 



C VALUES OF THE NORTH OR EAST COORDINATES 289 



C 290 



159 IF <LAG.GT.O) GO TO 280 291 

 C 292 

 C CALL SUBROUTINE FOR PLOTTING THE TRACK OF THE SURVEY VESSEL. 293 

 C 294 



160 400 NENTRY=1 295 



161 IF (NOPT(l).EO.l) CALL TRACK (NORTH. EAST. NMBR) 296 

 C 297 

 C APPLY A CORRECTION TO ALLOW FOR THE DISTANCE THE DETECTOR IS TOWED 298 

 C ASTERN OF THE SURVEY SHIP. 299 

 C 300 



162 DN0=NORTH(l)-(NORTH(2)-NORTH(D) 301. 



163 DE0=EAST(1)-(EAST(2)-EAST(1>) 302 



164 OENOM=SGRT((NORTH(l)-DN0)**2+(EAST(l)-DE0)**2> 303 



165 NCORD(l)=N0RTH(l)-CABLE* (NORTH! l)-DN0)/DENOM 304 



166 ECORD(l)=EAST(l)-CABLE*(EAST(l)-DE0)/DENOM 305 



167 DO 430 N=2.LINE 306 



168 DENOM=SQRT ( ( NORTH ( N) -NCORD ( N-l ) ) **2+ (EAST ( N ) -ECORD (N-l ) ) **2 ) 307 

 C 308 

 C THE CORRECTION FOR THE DISTANCE BETWEEN VESSEL AND THE DETECTOR 309 

 C IS EQUAL TO 'CABLE' UNLESS THE VESSEL IS LESS THAN 'CABLE' FEET 310 

 C AWAY FROM THE LAST COMPUTED POSITION OF THE DETECTOR VEHICLE. IN 311 

 C THIS EVENT, THE NEW COMPUTED DETECTOR POSITION IS THE SAME AS 312 

 C THE LAST DETECTOR POSITION 313 

 C 314 



51 



117 



270 



118 





119 







C 





C 





C 





C 



120 



280 



121 





122 





123 



290 



124 





125 



300 



126 





127 





128 





129 



310 



130 





131 





132 





133 



320 



134 





135 







C 





c 





c 





c 



136 



330 



137 





138 





139 





140 





141 





142 





143 



335 



144 





145 



340 



146 





147 





148 





149 





150 





151 





152 



350 



153 



360 



154 



370 



155 



380 



156 





157 



390 



158 





