
       FGR13PAK: Dose and Risk Coefficient File Package for FGR-13
                     K.F. Eckerman and R.W. Leggett
                         Dosimetry Research Group
                      Oak Ridge National Laboratory


INTRODUCTION

FGR13PAK is a collection of data files and supportive software written to 
facilitate use of the dose and risk coefficients from Federal Guidance
Report 13 in radiological assessment codes.  The package includes an 
illustrative FORTRAN module designed to access the data files and return 
the relevant information to the user's assessment code.

DISKETTE CONTENTS

All dosimetric data in FGR13PAK are in SI units.  The SI unit of absorbed
dose is the Gy (1 rad = 0.01 rad); the unit of activity is the Bq (1 Ci is
equivalent to 3.7E+10 Bq).  When the effective dose is shown it is indicated
to be in units of Sv (1 Sv = 100 rem). The following data files, all direct 
access files, are included in FGR13PAK:

1. FGR13CD.NDX (record length 173)
   The records of this file contain the record numbers (or keys) for the 
   radionuclides in all other files.  In addition, this file contains 
   information on the half-lives, decay mode, branching fractions, etc. of 
   the radionuclide. The file is used by FGR13PAK to assemble the information 
   on the decay chain of the parent.

2. FGR13ING.GDB (record length 354)
   The records of this file contain the committed absorbed dose coefficients 
   for 23 organs and the effective dose for intakes by ingestion. Coefficients
   are presented for the six ages considered in the ICRP Publication 56
   series. That is, the infant, 1-, 5-, 10-, 15-year-old, and the adult.  The
   low and high LET components of the absorbed dose are included in the file 
   and are displayed by FGR13PAK along with the effective dose.  The file 
   CHEMDATA.DAT contains information to identify multiple chemical forms of 
   the element.

3. FGR13INH.HDB (record length 363)
   The records of this file contain the committed absorbed dose coefficients
   for 23 organs and the effective dose for intakes by inhalation.  
   Coefficients are presented for the six ages considered in the ICRP 
   Publication 56 series.  That is, the infant, 1-, 5-, 10-, 15-year-old, and 
   the adult.  The low and high LET components of the absorbed dose are 
   included in the file.  The present file contains coefficients for only an 
   AMAD of 1.0E-06 m; the final version will contain data for additional 
   particle sizes. The file CHEMDATA.DAT contains information to identify 
   multiple chemical forms of the element.

4. F12TIII1.EXT (record length 250)
   The records of this file contain the dose coefficient for exposure to the
   radionuclide in the air; so called submersion.  The units are absorbed dose
   per time integrated air concentration (Gy per Bq s m^-3).  These data are 
   from Federal Guidance Report 12 (see Table III.1).

5. F12TIII3.EXT (record length 250)
   The records of this file contain the dose coefficient for exposure to the
   radionuclide on the ground surface.  The units are absorbed dose per time
   integrated surface concentration (Gy per Bq s m^-2).  These data are from
   Federal Guidance Report 12 (see Table III.3).

6. F12TIII7.EXT (record length 250)
   The records of this file contain the dose coefficient for exposure to the
   radionuclide distributed to an infinite depth in the soil.  The units are 
   absorbed dose per time integrated concentration (Gy per Bq s m^-3).  These
   data are from Federal Guidance Report 12 (see Table III.7).

   The above dose files contain data for the following organs.

            Adrenals     Bld Wall     B Surface    Brain 
            Breasts      Esophagus    St Wall      SI Wall
            ULI Wall     LLI Wall     Kidneys      Liver
            Lungs        Muscle       Ovaries      Pancreas
            R Marrow     Skin         Spleen       Testes
            Thymus       Thyroid      Uterus

   In addition the effective dose coefficient, denoted as e (based on ICRP 
   Publication 60) is also included in the files.

7. FGR13EXT.RBS (record length 118)
   The records of this file contain the risk (mortality and morbidity) 
   coefficients for external exposure pathways.  Data are presented for 
   exposure to radionuclides in the air (submersion), on the ground surface, 
   and within the volume on soil. Coefficients are presented in this file, and
   the other risk files, for the following cancer sites:

          Esophagus      Stomach       Colon       Liver
          Lung           Bone          Skin        Breast
          Ovary          Bladder       Kidneys     Thyroid
          Leukemia       Residual      Total


8. FGR13ING.RBS (record length 123)
   The records of this file contain the risk coefficients for ingestion
   intakes of the radionuclides in drinking water and in the diet.

9. FGR13INH.RBS (record length 124)
   The records of this file contain the risk coefficients for inhalation
   intakes of the radionuclides.

10. CHEMDATA.DAT (sequential file)
   File contains the identification information for multiple chemical
   forms of the elements ingested or inhaled.

The CD ROM includes the following software:

1. FGR13PAK.FOR
   This file is the FORTRAN source code for the access module.

2. PAKPARM.CMN
   This file contains information on the dimensioning of arrays in FGR13PAK.
   It should be included in the user's code with the include file FGR13PAK.CMN.

3. FGR13PAK.CMN
   This include file contains the common blocks which will contain the 
   information extracted from the data files per the user's request.  This
   file should be included in the user's code with PAKPARM.CMN.

4. IOLIST.CMN
   This include file contains the unit designations for the IO used by 
   FGR13PAK.  This is for information purposes.

5. BATCH.CMN
   This include file defines a common block with the logical variable DBATCH.
   If DBATCH is .true. then FGR13PAK runs silently as might be expected if 
   these routines are included in an assessment code.

6. FGR13PAK.INI
   This data file informs FGR13PAK of the location of all data files and the 
   specification of the RECL parameter in the FORTRAN open statement.  The 
   file is intended to enable the user to place the data files in an 
   appropriate directory structure where the data files are isolated from the 
   directory within which the software is being run.  It is necessary that 
   this file be present in the directory within which the software is being 
   run.  See the notes in FGR13PAK.INI regarding compiler-dependent RECL.

7. READEM.FOR
   This is the source code for a module that accesses FGR13PAK and returns to
   the screen the decay chain and dose coefficients for the members based on
   the users input.

8. READEM.EXE 
   A executable of the above with FGR13PAK not in the silent mode.


USE OF FGR13PAK

The FGR13PAK routines are accessed by the call

            call getcof(nuke, ipath)

where the variables have the following meaning:

  NUKE
     This is a character*7 variable representing the name of the radionuclide.
     The naming convention follows normal practice; i.e, Co-60, Cs-137,
     Tc-99m.  FGR13PAK includes a character*7 function NUKEOK that checks for
     proper notation in the nuclide name and for its existence within the
     data files.  See READEM.FOR for examples.
  
  IPATH
    This logical array (size 3) serves as a flag to which dose coefficients
    are of interest to the user.  This should be set to true to return 
    coefficients from all exposure routes (1 - inhalation, 2- ingestion, and
    3 - external)

The above call to GETCOF will return the information in the common blocks
defined in the included file FGR13PAK.CMN.  The following is a brief
description to the common blocks.

   common /dfacts/chemfinh(mspec,10),chemfing(mspec,10),namage(mage),
  :       cancer(mcan),organ(morg), dfinh(minh,mspec,mage,morg,mlet),
  :       LETh(mspec,mlet), dfing(ming,mspec,mage,morg,mlet),
  :       LETg(mspec,mlet), f1inh(minh,mspec,mage),
  :       f1ing(ming,mspec,mage), dfext(mspec, mext, morg),
  :       rinh(minh,mspec,mcan,2), ring(ming,mspec,mcan,2,2),
  :       rext(mspec,mcan,mext,2), NLET(mspec),Type(minh),
  :       iflag(mspec, mfact), ntypes(mspec), nfings(mspec),
  :       nfinh(mspec), nfing(mspec), nint, next

      chemfinh   character*25 array of chemical form names for inhalation
      chemfing   character*25 array of chemical form names for ingestion
      nameage    character*9 array of names for the mage ages
      cancer     character*9 array of the names of the mcan cancer sites
      organ      character*9 array of the morg organ names
      dfinh      inhalation absorbed dose coefficient array by type, chain
                 member, age, organ, and let
      LETh       character array of LET designation (L or H) for ispec for
                 each LET type for inhalation
      dfinh      ingestion absorbed dose coefficient array by chemical form, 
                 chain member, age, organ, and let
      LETg       character array of LET designation (L or H) for ispec for
                 each LET type for ingestion
      f1inh      real array of f_1 values for inhalation by type, chain
                 member, and age
      f1ing      real array of f_1 values for ingestion by chemical form,
                 chain member and age
      dfext      external absorbed dose coefficient array by chain member,
                 external pathway, and organ
      rinh       inhalation risk coefficient by type, chain member, cancer
                 site, and type of risk (mortality and morbidity)
      ring       ingestion risk coefficient by chemical form, chain member,
                 cancer site, pathway, and type of risk (mortality and
                 morbidity)
      rext       external risk coefficient by chain member, cancer site,
                 exposure pathway, and risk (mortality and morbidity)
      NLET       array of the number of LET types for each chain member
      type       character*1 array of the Type designation for each chain
                 member
      iflag      logical array indicating the presence of dose and risk
                 coefficients for each chain member and exposure route
      ntypes     array containing the number of inhalation types for each
                 chain member
      nfings     array containing the number of ingestion chemical forms
                 of each chain member
      iflag      logical flags for dose factor by chain member and pathway
      nfinh      array of number of special chemical forms for inhalation
      nfing      array of number of special chemical forms for ingestion
      nint       length of chain for internal factors
      next       length of chain for external factors
                                                                             
Information on the decay chain is contained in the following common block.

     common /radat/ thalf(mspec), iu(mspec), nucnam(mspec),                  
    :               branch(mspec, mspec), lmr(mspec), ibr(mspec,mspec),      
    :               nbr(mspec), nspec                                        
                                                                             
       thalf      A character (length 8) array of the half-life of each chain
                  member.
       iu         A character (length 2) array of the units of the half-life
                  of each chain member.
       nucnam     A character (length 7) array of the names of the chain
                  members.
       branch     A two dimension real array of branching fraction representing
                  the fraction of the decays of the ith chain member forming
                  the jth chain member.  An upper triangular matrix.
       lmr        A double precision array of decay constants (/d) for the
                  chain members.
       nbr        An integer array containing the number of precursors of each 
                  chain member.
       ibr        An integer array containing indices of the precursor
                  nuclides for each chain member. (ibr(i, ispec), i = 1,
                  nbr(ispec)) is the set of indices of radionuclides which decay 
                  directly to chain member ispec. 
       nspec      The length of the decay chain.

FGR13PAK truncates decay chains for both internal and external exposures.  For
internal exposure, the chain is truncated at member j if the sum of the energy 
of each radiation type (alpha, electron, and photon) emitted during a 100 year 
period by the first j chain members is greater than 99% of the cumulative 
values for the entire chain.  The truncation for external radiation 
considers only electrons and photon.

EXECUTABLE CODE

The READEM code can be used to examine the package.
