Main page >>>

SHIELDOSE 2. Program for radiation dose calculation.

SHIELDOSE is a computer code for space-shielding radiation dose calculations. It determines the absorbed dose as a function of depth in aluminum shielding material of spacecraft, given the electron and proton fluences encountered in orbit. The code makes use of precalculated, mono-energetic depth-dose data for an isotropic, broad-beam fluence of radiation incident on uniform aluminum plane media. These precalculated values are the result of detailed electron and electron-bremsstrahlung Monte Carlo calculations. The present version of SHIELDOSE calculates, for arbitrary proton and electron incident spectra, the dose absorbed in small volumes of the detector materials Al, H2O (tissue-equivalent detector), Si, and SiO2, in the following aluminum shield geometries: (1) in a semi-infinite plane medium, as a function of depth; (2) at the transmission surface of a plane slab, as a function of slab thickness; and (3) at the center of a solid sphere, as a function of sphere radius.

Main page of the program is modelweb.gsfc.nasa.gov

This package compiled using WATCOM FORTRAN/77 Ver.1.4 on WindowsXP SP2 after fixing minor bugs.

After starting program will ask you name of input file. It is customary to name input file with INP extension, for example AP8.INP. The input file had to contain all of required information. Information is the set of records delimited by spaces of carriage return symbols.

First two records are names of two output files. It's customary to give OUT and ARR extension to these files.

Example:
AP8.OUT
AP8.ARR

After that goes the number indicated detector type, for example 1:

1. AL DETECTOR
2. GRAPHITE DETECTOR
3. SI DETECTOR
4. AIR DETECTOR
5. BONE DETECTOR
6. CALCIUM FLUORIDE DETECTOR
7. GALLIUM ARSENIDE DETECTOR
8. LITHIUM FLUORIDE DETECTOR
9. SILICON DIOXIDE DETECTOR
10. TISSUE DETECTOR
11. WATER DETECTOR

Next number indicate the nuclear interaction account, for example 1:

1. NO NUCLEAR ATTENUATION FOR PROTONS IN AL.
2. NUCLEAR ATTENUATION, LOCAL CHARGED-SECONDARY ENERGY DEPOSITION.
3. NUCLEAR ATTENUATION, LOCAL CHARGED-SECONDARY ENERGY DEPOSITION, AND APPROX EXPONENTIAL DISTRIBUTION OF NEUTRON DOSE.

Next number indicate the number of depths of shieldings, for example 17.

After that number shows the unit of sheildings:
1. mils
2. g/cm2
3. mm

Next set of records are the shieldings in the same units as in previous point.

Next pair of numbers is about solar protons, the min and max energy of solar protons spectrum, default values are 0.1 и 10000 MeV.

Next pair of numbers is about trapped protons, the min and max energy of trapped protons spectrum, default values are 0.1 и 10000 MeV.

Next number is the number of spectrum points which divides spectrum for integration. Default value is 1001.

Next pair of numbers is about trapped electrons, the min and max energy of trapped electrons spectrum, default values are 0.05 и 10 MeV and number of specrtum points. Default value is 1001.

Next 72-symbol line can contain arbitrary text information for example the information about falling particles and detector type.

Next number is the number of points in falling spectrum of solar protons, trapped protons and trapped electrons.

After that the unit for conversation of flux. If your flux is in 1/(cm2 s keV), this unit is 1000. If your flux is 1/(cm2 s MeV), this unit is 1. Source files says: "EUNIT IS CONVERSION FACTOR FROM /ENERGY TO /MEV, E.G., EUNIT = 1000 IF FLUX IS /KEV.".

Next number is time of dose accumulation in seconds. Default value is 3.15360E+07 - year, You can indicate 86400 - one day.

Next numbers are spectra. If the spectrum of solar protons contans non-zero points the energies in MeV are shown. After that the fluxes are shown. After that the trapped proton spectrum and trapped electron spectrum in the condition of non-zero points number of spectrum.

This is the end of input file description.

After the program work finish the OUT file will be created.

Download program.

Main page >>>
Hosted by uCoz