home *** CD-ROM | disk | FTP | other *** search
/ HAM Radio 1 / HamRadio.cdr / misc / pc_ham3 / falout.bas < prev    next >
BASIC Source File  |  1986-12-14  |  2KB  |  62 lines
  1. 20 REM
  2. 30 REM FALOUT - A NUCLEAR FALLOUT CALC PROGRAM
  3. 40 REM (C) JOHN R. FANCHI, JULY 1985
  4. 50 REM MICROSOFT BASIC ON CP/M-80 OPERATING SYSTEM
  5. 60 REM CHR$(26) CLEARS SCREEN
  6. 70 PRINT CHR$(26)
  7. 80 PRINT " FALOUT - ESTIMATING THE DISTRIBUTION OF NUCLEAR FALLOUT"
  8. 90 PRINT
  9. 100 PRINT "(C) J.R. FANCHI, JULY 1985
  10. 110 PRINT: PRINT: PRINT 
  11. 120 PRINT "ESTIMATING RADIATION DOSAGE"
  12. 130 PRINT
  13. 140 INPUT "ENTER THE YIELD OF THE NUCLEAR BLAST IN MEGATONS: ", Y
  14. 150 PRINT
  15. 160 INPUT "ENTER YOUR DISTANCE FROM THE BLAST IN MILES: ", D
  16. 170 PRINT
  17. 180 REMS=250*1000*Y/(16*3.1416*D*D)
  18. 190 PRINT "ESTIMATED RADIATION DOSAGE IN REMS: ",REMS
  19. 200 PRINT
  20. 210 PRINT: PRINT: PRINT
  21. 220 PRINT "ESTIMATING FALLOUT DISTRIBUTION"
  22. 230 PRINT 
  23. 240 PRINT "THE DISTRIBUTION OF FALLOUT IS TREATED AS A RANDOM WALK PROCESS"
  24. 250 PRINT "SIMILAR TO BROWNIAN MOTION."
  25. 260 PRINT
  26. 270 INPUT "ENTER THE LINE-OF-SIGHT WIND SPEED (MPH): ",VPAR
  27. 280 VPAR=VPAR*24!
  28. 290 PRINT
  29. 300 INPUT "ENTER THE LINE-OF-SIGHT DISPERSION (SQ MILES/D): ",DPAR
  30. 310 PRINT
  31. 320 INPUT "ENTER THE TRANSVERSE WIND SPEED (MPH): ",VTRANS
  32. 330 VTRANS=VTRANS*24!
  33. 340 PRINT
  34. 350 INPUT "ENTER THE TRANSVERSE DISPERSION (SQ MILES/D): ",DTRANS
  35. 360 PRINT
  36. 370 REM
  37. 380 PRINT "RADIATION DOSAGE AT YOUR LOCATION AS A FUNCTION OF TIME"
  38. 390 PRINT "DAYS","REMS","NORM","HRS"
  39. 400 REM
  40. 410 ALPHA=D*D/(4*DPAR)
  41. 420 BETA=-2*D*VPAR/(4*DPAR)
  42. 430 GAMMA=VPAR*VPAR/(4*DPAR) + VTRANS*VTRANS/(4*DTRANS)
  43. 440 TFAC=SQR(1+4*ALPHA*GAMMA)
  44. 450 TYMMAX=(1+TFAC)/(2*GAMMA)
  45. 460 FACMAX=((D-VPAR*TYMMAX)^2)/(4*DPAR*TYMMAX)
  46. 470 FACMAX=FACMAX+VTRANS*VTRANS*TYMMAX/(4*DTRANS)
  47. 480 FOR I=1 TO 20 STEP 1
  48. 490 TYM=TYMMAX*(1-.02*(10-I))
  49. 500 FACPAR=0! : FACPER=0!
  50. 510 DISPAR=D-VPAR*TYM
  51. 520 DISPER=VTRANS*TYM
  52. 530 FACPAR=DISPAR*DISPAR/(4!*DPAR*TYM)
  53. 540 IF DTRANS<>0! THEN FACPER=DISPER*DISPER/(4!*DTRANS*TYM)
  54. 550 FACTOR=FACPAR+FACPER
  55. 560 COEF=TYMMAX/TYM
  56. 570 EFFREM=REMS*COEF*EXP(FACMAX-FACTOR)
  57. 580 NORM=EFFREM/REMS
  58. 590 TYMHRS=TYM*24
  59. 600 PRINT TYM,EFFREM,NORM,TYMHRS
  60. 610 NEXT I
  61. 620 END
  62.