gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS


/tabNuclearMass [ Modules ]

[ Top ] [ Modules ]

NAME

module tabNuclearMass

PURPOSE

  • Provide database structure for nuclides (mass, binding energy, ground state spin, name)
  • intitialize nuclear decay properties (unused/untested in this context)
  • provide access to single data points without loading the entire 200kB data file to RAM
  • Read table "path_To_Input/nndc.dat", which contains the nuclide properties: "Ame2003 atomic mass evaluation (II)"
  • G.Audi, A.H. Wapsta, C. Thibautl, NPA729(2003)337-676, original file: mass.mas03
  • Based on Theo Gaitanos 'testRun/auswerteTools/ClusterCode/code/InitGEM.f90'

AUTHORS

   Theo Gaitanos, Ivan Lappo-Danilevski


tabNuclearMass/Nuclide [ Types ]

[ Top ] [ tabNuclearMass ] [ Types ]

NAME

type Nuclide

SOURCE

  type Nuclide
     integer                            :: ZahlN        !neutron number
     integer                            :: ZahlZ        !proton number
     real                               :: MassEx       !mass excess
     real                               :: Mass         !mass in GeV
     real                               :: Bind         !binding energy
     real                               :: spin         !groundstate spin
     character(7)                       :: Name         !name of element
     integer, allocatable, dimension(:) :: Ejectile     !index-position of ejectile
     integer, allocatable, dimension(:) :: Daughter     !index-position of daughter
     logical, allocatable, dimension(:) :: ChannelStatus !status of decay channel
  End Type Nuclide

NOTES

Constructor for nuclide properties


tabNuclearMass/NucMax [ Global module-variables ]

[ Top ] [ tabNuclearMass ] [ Global module-variables ]

NAME

integer, parameter :: NucMax

PURPOSE

Maximum number of nuclides considered in evaporation model, number of entries in 'nndc.dat'

NOTES

DO NOT CHANGE THIS NUMBER!!!


tabNuclearMass/Elements [ Global module-variables ]

[ Top ] [ tabNuclearMass ] [ Global module-variables ]

NAME

type(Nuclide), Allocatable,dimension (:),save :: Elements

PURPOSE

Database for elements properties

USAGE

  • To get e.g. the mass excess of the nuclide IN type "Elements(IN)%MassEx". The ID's are choosen according to NNDC-Table.
  • The array indizes of "Elements" run from NucMax=1n up to NucMax=208Pb (IN=1,...,NucMax).


tabNuclearMass/NucPosition [ Global module-variables ]

[ Top ] [ tabNuclearMass ] [ Global module-variables ]

NAME

Character, Allocatable,dimension(:,:),save :: NucPosition

PURPOSE

(N,Z)-Matrix providing with the position (in the nndc.dat-file) of a stable (according N/Z-ratio) nuclide. Negative Value of a matrix element indicates not existing nulcides.


tabNuclearMass/read_element [ Subroutines ]

[ Top ] [ tabNuclearMass ] [ Subroutines ]

NAME

subroutine read_element

PURPOSE

INPUTS

  • integer, intent(in) :: Ni -- Neutron number of requested element
  • integer, intent(in) :: Zi -- Proton number

OUTPUT

  • type(Nuclide),intent(out) :: element -- information about the nucleus

NOTES

returns an element with %name='missing', when element is not found in table


tabNuclearMass/read_mass [ Functions ]

[ Top ] [ tabNuclearMass ] [ Functions ]

NAME

function read_mass

INPUTS

OUTPUT

  • real :: read_mass

PURPOSE

reads the tabulated mass (in GeV) for a given N,Z


tabNuclearMass/init_Mass_Table [ Subroutines ]

[ Top ] [ tabNuclearMass ] [ Subroutines ]

NAME

subroutine init_Mass_Table

PURPOSE

NOTES

  • Initialization of (1:MaxN,1:MaxZ)-Matrix "NucPosition" is done here.
  • Variables "MaxN" and "MaxZ" give the max. possible value of neutron- and proton-content of a cluster. Its status (stable/unstable) according N/Z-ratio is initialized here.
  • NucPosition is needed in checking the correct N/Z-ratio of produced clusters, before starting the evaporation procedure.


tabNuclearMass/SetDecays [ Subroutines ]

[ Top ] [ tabNuclearMass ] [ Subroutines ]

NAME

subroutine SetDecays

PURPOSE

  • Defines valid decay channels for each nuclide
  • Important for Monte-Carlo decision of decaying cluster
  • The list of elements is taken from the mass.mas03-table (cf. init_Mass_Table)

NOTES

Tested only with these lines of code in GiBUU

     call Init_Mass_Table
     do i = 1,NucMax
       write (*,*) Elements(i)%name
     end do
     call SetDecays
     call End_Mass_Table
     write (*,*) 'allocation/deallocation successfull'
     STOP

Thus allocation and deallocation works. Fully tested in other context by Theo Gaitanos