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



169 IF (DENOM-CABLE) 420.420.410 315 



170 410 NC0RD(N)=N0RTH(N>-CABLE*(N0RTH<N>-NC0R0(N-1)>/DEN0M 316 



171 EC0RD(N>=EAST(N>-CABLE*(EAST(N)-EC0RD<N-1)>/DEN0M 317 



172 GO TO 430 318 



173 420 NC0RD(N)=NC0RD(N-1> 319 



174 EC0RD<N)=EC0RD(N-1> 320 



175 430 CONTINUE 321 

 C 322 

 C CONVERT RADIATION READINGS TO COUNTS PER SECOND 323 

 C 324 



176 CCR(1)=KAD(1)/SEC 325 



177 DO 450 N=2.LINE 326 



178 IF (NMBR(N).NE.NMBR(N-1)+1) GO TO U««0 327 



179 CCR(N)=RAD(N)/(TIME(N)-TIME(N-1>) 328 



180 GO TO 450 329 



181 WO CCR(N)=RAD(N)/SEC 330 



182 450 CONTINUE 331 

 C 332 

 C COMPUTE THE MEAN AND STANDARD DEVIATION OF THE BACKGROUND COUNT 333 

 C RATE FROM THE RADIATION DATA THAT LIES WITHIN THE LIMITS OF THE 334 

 C ESTIMATED BACKGROUND COUNT RATE SET BY CHAUVENET'S CRITERION. 335 

 c 336 



183 IF (CCRll).GT.2.*BKG) CCR(1>=BKG 337 



184 ALPHA=1.-1./FL0AT(2«LINE) 338 



185 CALL TINORM (CHVRr ALPHA. S455) 339 



186 GO TO 4b0 340 



187 455 CHVR=5. 341 



188 WRITE (6(770) ALPHA 342 



189 460 8GCR=0.0 343 



190 SIGMA=0.0 344 



191 NCOUNT=0 3 '' 5 



192 GATE1=BKG-CHVR*SQRT(BKG/SEC) 346 



193 GATE2=BKG+CHVR*S0RT(BKG/5EC) 347 



194 DO 490 N=1.LINE 3I » 8 



195 IF (CCR(N)-GATEl) 490.490.470 3"»9 



196 470 IF (CCR(N)-GATE2) 480.490.490 350 



197 480 NC0UNT=NC0UNT+1 351 



198 13GCR=BGCR + CCR(N) 352 



199 490 CONTINUE 353 



200 BGCR=BGCR/FLOAT(NCOUNT) 354 



201 DO 520 N=1»LINE 355 



202 IF (CCR(N)-GATEl) 520>520.500 356 



203 500 IF (CCR(N)-GATE2) 510.520.520 357 



204 510 SIGMA=SIGMA+(CCR(N)-BGCR>**2 358 



205 520 CONTINUE 359 



206 SIGMA=SGRT(SIGMA/FLOAT(NCOUNT>) 360 



207 WRITE (6>780) BKG.BGCR. SIGMA 361 

 C 362 

 C CALL THE SUBROUTINE FOR PLOTTING UNCORRECTED RADIATION VALUES 363 

 C 364 



208 NENTRY=2 365 



209 IF (NOPT(2).EQ.l) CALL RADPLT (NCORD.ECORD.CCR) 366 



210 IF (NOPT(3).NE.l) GO TO 570 367 

 C 368 

 C CORRECT RADIATION VALUES FOR BACKGROUND COUNT RATE AND TIME-DECAY 369 

 C 370 



211 SUM=0.0 371 



212 NBAR=0.0 372 



213 ERAR=0.0 373 



214 SDNRTH=0.0 374 



215 SDEAST=0.0 375 



216 DO 540 N=1.LINE 376 



217 CCRU)=CCR(N)-BGCR 377 



218 IF(CCR(N)-3.*SIGMA) 540.540.530 378 



219 530 CCR(N)=(CCR(N)-3.*SIGMA)*EXP(DECAY*(DELAY+TIME(N)/3600.)> 379 



1 + 3.»SIGMA 3795 



220 SUM=SUM+CCR(N) 380 



221 NRAR=NBAR+(NCORD(N)-NCORD(l))*CCR(N) 381 



222 EBAR=EBAR+(ECORD(N)-ECORD(l)>*CCR(N) 382 



223 540 CONTINUE 383 



224 WRITE (6.790) SUM 3B4 

 C 385 

 C COMPUTE WEIGHTED MEAN AND STD. DEV. OF ACTIVITY LOCA 1 * .ON 386 

 C 387 



225 NBAR=NCORD(l)+NfaAR/SUM 388 



226 EBAR=ECORD(l)+EBAR/SUM 389 



227 DO 560 N=1.LINE 390 



228 IF (CCR(N)) 560.560.550 391 



52 



