PROGRAM — RAPLOT II 



229 610 CONTINUE 395" 



230 IF (CCR(1«N).LE.0.0) 60 TO 620 396 

 £31 SUM=SUM+CCR(1.N) 397 



232 SUM2=SUM2+CCR(2»N) 398 

 C 399 

 C COMPUTE WEIGHTED MEAN AND STD. OEV. OF ACTIVITY LOCATION 400 

 C 401 



233 EBAR=EBAR+(tCORD(N)-ECORD(l))*CCR(l»N) 402 



234 NeAR=NBAR+(NC0Ru(N)-NC0RD(l))*CCR(l»N) 403 



235 620 CONTINUE 404 



236 NBAR=NC0KD(1)+N6AR/SUM 405 



237 E3AR=EC0R0(l)+EbAR/SUM 1*06 

 23B 5DNRTH=0.0 407 



239 SDEAST=0.0 408 



240 DO 640 N=1.LINE 409 



241 IF (CCR(1»N>> 640.640.630 410 



242 630 SDNRTH=SDNRTH+(NC0R0(N)-NBAR)»»2*CCR(1»N> 411 



243 SDEAST=S0EAST+(EC0R0(N)-EBAR)**2»CCR(1»N) 412 



244 640 CONTINUE 413 



245 SDNRTH=buRT(SDNPTH/SUM> 414 



246 S0EAST=SQRT(SDEA3T/SUM) 415 



247 WRITE (6.950) SUM.SUM2 416 



248 WRITE (6.960) NBAR.EBAR.SDNRTH.SDEAST 417 

 C 418 

 C COMPUTE AND PRINT 95 PC. CONFIDENCE LIMITS OF MEAN RADIATION 419 

 C LOCATION. 420 

 C 421 



249 RTSUy,=SuRT(SUM/BGCR(D) 422 



250 CFIDN=1.96*S0NRTH/RTSUM 423 



251 CFIDE=1.96*SDEAST/RTSUM 424 



252 WRITE (6.970) CFIDN.CFIDE 425 



253 650 DO 660 1=1.2 426 



254 CCR(I»L1NE+1)=BGCR(I> 427 



255 660 CCR(I.LINE+2)=SIGMA(I> 428 



256 NWRD=(LINE+2>*2 429 



257 CALL NTRAN ( ITAPE. 1 »NWRD»CCR.L<4> > 430 



258 NENTRY=3 431 



259 IF (N0PT(3).E0.1) CALL RDPLT2 (NCORD.ECORO.CCR) 432 



260 WRITE (6.980) NMAX.EMAX 433 



261 WRITE (6.990) NMIN.EMIN 434 

 C 435 

 C WRITE OUT THE NUMBER. COORDINATES AND ACTIVITY OF EACH DATA POINT 436 

 C 437 



262 K0UNT=50 438 



263 DO 720 N=1.LINE 439 



264 DO 690 J=1.3 440 



265 IF (FLD(3b-J.l.IERR(N>)> 670.670.680 441 



266 b70 ISiGN«J>=BLANK 442 



267 GO TO 690 443 



268 680 ISIGN(J)=AS 444 



269 690 CONTINUE 445 



270 IF (KOUNT-50) 710.700.700 446 



271 700 WRITE (.6.1000) LEGEND ,4.47 



272 KOUNT=0 448 



273 710 WRITE (6.1010) NMBR(N) .TIME(N) »D< l.N) » ISIGN(l) .D(2. N) » ISIGN(2) . (RA 449 



1D(I.N).I=1.2).N0RTH(N).EAST(N).ISIGN(3).NC0RD(N).EC0RD(N).(CCR(I.N 450 



2). 1=1. 2) 451 



274 720 KOUNT=KOUNT+l 452 



275 WRITE (6.1020) 453 



276 IF (.NOT.JOBENd) GO TO 10 454 



277 730 CONTINUE 455 

 C 456 

 C END-FILE PLOT TAPE 457 

 C 458 



278 END FILE 9 4585 



279 DO 760 1=1.4 459 



280 740 IF (L(I).NE.-l) GO TO 750 460 



281 GO TO 740 461 



282 750 IF (L(I).LT.-l) GO TO 770 462 



263 760 CONTINUE 463 



264 CALL NTRAN (ITAPE.9) 464 



285 770 STOP 465 

 C 466 



286 7B0 FORMAT (3F3.0 .2F7.0. 1X.2F2.0.F3.0.F5.2.3F2.0.4I 1X.F7.0) . 13) 467 



287 790 FORMAT (41 1 .F10.0. 1X.3I2. 1X.F10.0.2I 1X.F10.O) .IX. A6) 468 



288 800 FORMAT (13A6.L2) 469 



289 810 FORMAT ( 10X. 13A6. 10X.L2) 470 



290 820 FORMAT (/5X.UH8EAC0N 1 .F10.0. 1HN.F10.0. 1HE.5X.8HBEAC0N 2.F10.0.1HN 471 



1.F10.0.1HE//) 472 



35 



