

C 





C 





C 





c 





c 



173 





174 



460 



175 





176 





177 



470 



178 





179 



480 



160 





181 



490 



182 





183 



500 





C 





C 





C 



184 





165 



505 



186 





187 





188 





189 



510 





C 





C 





c 





c 





c 



190 





191 





192 





193 





19* 





195 



515 



196 





197 



520 



198 





199 





200 





201 





202 





203 





204 





205 



530 



206 



540 



207 





208 



550 



209 





210 





211 





212 



560 



213 



570 



214 



580 



215 



590 



216 







C 





C 





C 



217 





218 





219 







C 





C 





C 



220 





221 





222 





223 





224 





225 





226 





227 





228 



600 



PROGRAM — RAPLOT II 



DEN0H=SuRT ( (NORTH(N)-NCORD(N-l) )**2+(EAST(N)-EC0RD(N-l) )*»2> 318 



319 



THE CORRECTION FOR THE DISTANCE BETWEEN VESSEL AND THE DETECTOR 320 



IS EQUAL TO 'CABLE' UNLESS THE VESSEL IS LESS THAN 'CABLE' FEET 321 



AWAY FROM THE LAST COMPUTED POSITION OF THE DETECTOR VEHICLE. IN 322 



THIS EVENT. THE NEW COMPUTED DETECTOR POSITION IS THE SAME AS 323 



THE LAST DETECTOR POSITION 324 



325 



IF (DENOM-CABLE) 470.470.460 326 



NCORD(N)=NORTH(N)-CABLE«(NORTH(N>-NCORD(N-l> )/DEN0M 327 



EC0RD(N)=EAST(N)-CABLE*(EAST(N)-EC0RD(N-1> )/DENOM 328 



GO TO 480 329 



NC0RD(N)=NC0RD(N-1) 330 



ECORD(N)=ECORD(N-l) 331 



CONTINUE 332 



CALL NTRAN ( ITAPE.l. LINE. NCORD. LAG. l.LINE.ECORD.L(S) ) 333 



IF (L(4).NE.-1) GO TO 500 334 



GO TO 490 335 



IF (L(41.LT.-1) GO TO 770 336 



337 



CONVERT RADIATION READINGS TO COUNTS PER SECOND 338 



339 



IF (dKG(ll) 650. 650. 505 3395 



CCR(2.1)=RAD(2.1)/SEC 340 



CCR(1.1)=RAD(1.1)/SEC 341 



DO 510 N=2.LINE 342 



DO 510 1=1.2 343 



CCR(I.N)=RAD(I»N)/SEC 344 



345 



COMPUTE THE MEAN AND STANDARD DEVIATION OF THE BACKGROUND COUNT 346 



RATE FROM THE RADIATION OATA THAT LIES WITtUN THE LIMITS OF THE 347 



ESTIMATEO BACKGROUND COUNT RATE SET BY CHAUVENET'S CRITERION. 348 



349 



IF (CCR(1.1).GT.2.*BKG(1>) CCRC1. 1)=BKG(1) 350 



IF (CCR(2'1).GT.2.*BKG(2>> CCR (2. 1 )=BKG<2) 351 



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



CHVR=TINORM(ALPHA»*515> 353 



GO TO 520 354 



CHVR=5. 355 



WRITE (6.930) ALPHA 356 



DO 590 1=1.2 357 



BGCR(I)=0.0 358 



SIGMA(I)=0.0 359 



NCOUNT=0 360 



GATE1=BKG ( I ) -CHVR*SQRT ( BKG ( I ) /SEC ) 361 



GATE2=BKG ( I ) +CHVR*SQRT(BKG 1 1 ) /SEC) 362 



DO 550 N=1.LINE 363 



IF (CCR(I.N)-GATEl) 550.550.530 364 



IF <CCR(I»N)-GATE2) 540.550.550 365 



NC0UNT=NC0UNT+1 366 



BGCR(I)=BGCR(I)+CCR(I.N) 367 



CONTINUE 368 



BGCR ( I ) =BGCR ( I ) /FLOAT t NCOUNT ) 369 



DO 580 N=1.LINE 370 



IF (CCR(I.N)-GATEl) 580.580.560 371 



IF <CCR(I.N)-!GATE2) 570.580.580 372 



SIGMA(I)=SIGMA(I)+(CCR(I.N)-BGCR(I))**2 373 



CONTINUE 374 



SIGMA(I)=SQRT(SIGMA(I)/FL0AT(NC0UNT)) 375 



WRITE (6.940) BKG. BGCR. SIGMA 376 



377 



CALL THE SUBROUTINE FOR PLOTTING UNCORRECTED RADIATION VALUES 378 



379 



NENTRY=2 380 



IF (N0PT(2).EQ.1> CALL RDPLT2 (NCORD. ECORD. CCR) 381 



IF (NOPT(3).NE.l) GO TO 650 382 



383 



CORRECT RADIATION VALUES FOR BACKGROUND COUNT RATE AND TIME-DECAY 384 



385 



SUM=U.O 386 



SUM2=0.U 387 



NBAR=0.0 388 



EBAR=0.0 389 



00 620 N=1.LINE 390 



DO 610 1=1.2 391 



CCR<I.N)=CCR(I»N)-BGCR(I) 392 



IF (CCR(LN) - 3.*SIGMA(I)> 610.610.600 393 



CCR(I.N)=(CCR(I.N)-3.*SIGMA(I))*EXP(DECAY*(DELAY+TIME(N)/3600.)) 394 



1 + 3.*SIGMA(I) 3945 



34 



