NESC0924 SUBDOSA.

1. NAME OR DESIGNATION OF PROGRAM:  SUBDOSA.

3. DESCRIPTION OF PROBLEM OR FUNCTION

SUBDOSA calculates external gamma and beta doses to individuals from the accidental release of radionuclides to the atmosphere. Doses are calculated as a function  of: quantity released, duration of release, atmospheric conditions during the release, and horizontal distance from the release point.  Doses from both gamma and beta radiation are calculated as a
function of depth in tissue, summed and reported as skin, eye, gonodal, and total body dose. Doses are calculated within each of several release time intervals with separate nuclide inventories and atmospheric dispersion conditions used for each interval. Radioactive decay using a chain decay scheme with branching to account for transitions to and from isomeric states is considered during the release. The program computes normalized air concentrations at ground level as a function of distance from the point of release.

4. METHOD OF SOLUTION

SUBDOSA divides the gamma spectrum into energy  groups. The incremental dose is calculated separately as a function  of photon energy and integrated over the cloud volume to obtain dose rate factors for each energy group. These dose rate factors are coupled with radionuclide release to give individual dose rates. The total dose is calculated as the time integral of short range of beta particles in air, the radionuclide cloud is assumed to be semi- infinite in dimension. With this assumption, the beta dose to surface tissue is calculated as the product of the effective beta energy for the nuclide being considered and the time-integrated air  concentration at the receptor, taking into account body self- shielding. The total beta dose is the sum of dose contributions from each nuclide.

5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

Maxima of -
   500 nuclides
    10 downwind distances
     6 release time intervals
     4 tissue depths
SUBDOSA does not consider the dose resulting from radionuclides deposited in the body and its organs via inhalation.

6. TYPICAL RUNNING TIME

NESC executed the SUBDOSA sample problem in 2 minutes of CP time on a UNIVAC1100/44.

The program included in this package has been run at NEA-DB on a UNIVAC 1100/92 computer.

7. UNUSUAL FEATURES OF THE PROGRAM:

8. RELATED AND AUXILIARY PROGRAMS

Other programs which calculate organ dose are DACRIN (NESC Abstract 923), which calculates organ dose from acute or chronic radionuclide inhalation, AERIN (NESC Abstract 908), which calculates organ and tissue burdens resulting from acute exposure to a radioactive aerosol, and ARRRG,FOOD (NESC Abstract 925) and PABLM (NESC Abstract 926), which calculate radiation dose to humans from radionuclides in the environment.

10. REFERENCES:

- D.L. Strenge, E.C. Watson and J.R. Houston :
  SUBDOSA : A Computer Program for Calculating External Doses from
  Accidental Atmospheric Releases of Radionuclides
  BNWL-B-351 UC-11 (June 1975)
- SUBDOSA Tape Description BNWL-B-351 Errata and Sample Problem
  Output
  NESC Note 81-81 (April 2, 1981)

11. MACHINE REQUIREMENTS

43,000 words of memory are needed for execution of SUBDOSA. BIVAR, an auxiliary program which calculates dose rate factors as a function of energy, requires 13,000 words of  memory.

13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED:  EXEC 8.

EXEC8 39R7 (UNIVAC 1100/92).

14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS

This
program should be compiled with the ASCII FORTRAN (FTN) compiler.

15. NAME AND ESTABLISHMENT OF AUTHOR

             D.L. Strenge*, E.C. Watson, and J.R. Houston
             Ecological Sciences Department
             Pacific Northwest Laboratories
             Battelle
             P.O. Box 999
             Richland, Washington 99352
* Contact

File name	File description	Records
NESC0924_01.001	Information file	47
NESC0924_01.002	JCL and control information	56
NESC0924_01.003	SUBDOSA source	1612
NESC0924_01.004	SUBDOSA sample problem input	26
NESC0924_01.005	Radionuclide decay data library (RMDLIB)	288
NESC0924_01.006	Gamma energy data library (GISLIB)	800
NESC0924_01.007	Beta energy data library (RNDBET)	497
NESC0924_01.008	Dose rate convers. factor library (BIVLIB)	579
NESC0924_01.009	BIVAR source	796
NESC0924_01.010	BIVAR sample problem input	9
NESC0924_01.011	SUBDOSA sample problem output	182
NESC0924_01.012	BIVAR sample problem output	408