TABLE OF CONTENTS
- 1. /densityStatic
- 1.1. densityStatic/staticDensityInit
- 1.2. densityStatic/TabulateZero
- 1.3. densityStatic/TabulateSphere
- 1.4. densityStatic/TabulateDensityWoodsSaxon
- 1.5. densityStatic/TabulateDensityHarmOsc
- 1.6. densityStatic/TabulateDensityLenske
- 1.7. densityStatic/TabulateDensityLuis
- 1.8. densityStatic/TabulateDensityExRTF
- 1.9. densityStatic/SearchMaxVals
- 1.10. densityStatic/staticDensity
- 1.11. densityStatic/densityLuis
- 1.12. densityStatic/TabulateDensityBirger
- 1.13. densityStatic/TabulateDensityBirgerWelke
- 1.14. densityStatic/ReAdjust
/densityStatic [ Modules ]
NAME
module densityStatic
PURPOSE
Collect routines for STATIC density calculations.
densityStatic/staticDensityInit [ Subroutines ]
[ Top ] [ densityStatic ] [ Subroutines ]
NAME
subroutine staticDensityInit (nuc)
PURPOSE
decide, which density parametrisation is used. Then tabulate this and also set the extreme values for the MC decision.
INPUTS
OUTPUT
- type(tNucleus) :: nuc
densityStatic/TabulateZero [ Subroutines ]
[ Top ] [ densityStatic ] [ Subroutines ]
NAME
subroutine TabulateZero(nuc)
PURPOSE
Set the density table to zero
INPUTS
- type(tNucleus) :: nuc
OUTPUT
- type(tNucleus) :: nuc
densityStatic/TabulateSphere [ Subroutines ]
[ Top ] [ densityStatic ] [ Subroutines ]
NAME
subroutine TabulateFermiGas(nuc)
PURPOSE
Tabulate a sphere with constant density.
Tabulates the static density to make it available faster for later use
INPUTS
- type(tNucleus) :: nuc
OUTPUT
- type(tNucleus) :: nuc
densityStatic/TabulateDensityWoodsSaxon [ Subroutines ]
[ Top ] [ densityStatic ] [ Subroutines ]
NAME
subroutine TabulateDensityWoodsSaxon (nuc)
PURPOSE
Tabulate the Woods-Saxon distribution. Tabulates the static density to make it available faster for later use.
INPUTS
- type(tNucleus) :: nuc
OUTPUT
- type(tNucleus) :: nuc
densityStatic/TabulateDensityHarmOsc [ Subroutines ]
[ Top ] [ densityStatic ] [ Subroutines ]
NAME
subroutine TabulateDensityHarmOsc(nuc)
PURPOSE
Tabulate the density distribution according harmonic oscillator shell modell.
Tabulates the static density to make it available faster for later use
INPUTS
- type(tNucleus) :: nuc
OUTPUT
- type(tNucleus) :: nuc
NOTES
The parameter are taken from the FRITIOF package.
densityStatic/TabulateDensityLenske [ Subroutines ]
[ Top ] [ densityStatic ] [ Subroutines ]
NAME
subroutine TabulateDensityLenske (nuc)
PURPOSE
Tabulate the density distribution according to Woods-Saxon distribution but with refined charge radii for proton and neutron according to H. Lenske. Tabulates the static density to make it available faster for later use.
INPUTS
- type(tNucleus) :: nuc
OUTPUT
- type(tNucleus) :: nuc
NOTES
Everything in fm.
densityStatic/TabulateDensityLuis [ Subroutines ]
[ Top ] [ densityStatic ] [ Subroutines ]
NAME
subroutine TabulateDensityLuis(nuc)
PURPOSE
Tabulate the density distribution of matter (p and n) and the density of centers ( p and n number densities ) following J. Nieves et al., Pionic atoms... NPA554
Tabulates the static density to make it available faster for later use
INPUTS
- type(tNucleus) :: nuc
OUTPUT
- type(tNucleus) :: nuc
NOTES
everything in fm
densityStatic/TabulateDensityExRTF [ Subroutines ]
[ Top ] [ densityStatic ] [ Subroutines ]
NAME
subroutine TabulateDensityExRTF(nuc)
PURPOSE
Tabulate the density distribution according to Relativistic Thomas-Fermi model code from Horst Lenske.
Tabulates the static density to make it available faster for later use
INPUTS
- type(tNucleus) :: nuc
OUTPUT
- type(tNucleus) :: nuc
NOTES
everything in fm
densityStatic/SearchMaxVals [ Subroutines ]
[ Top ] [ densityStatic ] [ Subroutines ]
NAME
subroutine SearchMaxVals(nuc)
PURPOSE
go through the tabulated distributions to search for the extrema
densityStatic/staticDensity [ Functions ]
[ Top ] [ densityStatic ] [ Functions ]
NAME
type(dichte) function staticDensity(r,nucl)
PURPOSE
gives density in the restframe of the nucleus "nucl" at position "r"
INPUTS
- real, dimension(1:3),intent(in) :: r -- position where density should be calculated
- type(tNucleus),pointer :: nucl -- nucleus which is regarded
USAGE
(dichte)=staticDensity(...)
densityStatic/densityLuis [ Subroutines ]
[ Top ] [ densityStatic ] [ Subroutines ]
NAME
subroutine densityLuis(r,z,a,rhop_mat,rhon_mat,center_in)
PURPOSE
- This routine calculates the density of matter (p and n) and the density of centers ( p and n number densities )
- following J. Nieves et al., Pionic atoms... NPA554
- everything in fm
- returns per default the density of matter, use "center_in" to switch to density of centers
INPUTS
- real, intent(in)::r -- radius
- integer, intent(in)::z -- charge of the nucleus
- integer, intent(in)::a -- atomic number
- logical, optional, intent(in) :: center_in -- if true then density of centers is given in the output
RESULT
- real, intent(out):: rhop,rhon -- Proton and neutron densities at r
- real, intent(out):: rp,ap,rho0p,rn,an,rho0n -- parameters of the density distributions
densityStatic/TabulateDensityBirger [ Subroutines ]
[ Top ] [ densityStatic ] [ Subroutines ]
NAME
subroutine TabulateDensityBirger(nuc)
PURPOSE
Tabulate the density distribution based on a local density approximation first described by Brueckner et al.
INPUTS
- type(tNucleus) :: nuc
OUTPUT
- type(tNucleus) :: nuc
densityStatic/TabulateDensityBirgerWelke [ Subroutines ]
[ Top ] [ densityStatic ] [ Subroutines ]
NAME
subroutine TabulateDensityBirgerWelke(nuc)
PURPOSE
Tabulate the density distribution based on a local density approximation first described by Brueckner et al. and a momentum-dependent potential
INPUTS
- type(tNucleus) :: nuc
OUTPUT
- type(tNucleus) :: nuc
densityStatic/ReAdjust [ Subroutines ]
[ Top ] [ densityStatic ] [ Subroutines ]
NAME
subroutine ReAdjust(nuc, potP, potN)
PURPOSE
This routine recalculates the density distributions for protons and neutrons by considering the given potentials as static and fulfill the condition
sqrt(p_F^2+m_N^2) + U - m_N == E_sep ~ -8MeV
With the Local-Thomas-Fermi, we connect the resulting fermi momentum to a density,
rho = p_F^3/(3pi^2)
Since the potentials are given as function of r, we calculate rho(r).
INPUTS
- type(tNucleus),pointer :: nuc -- the nucleus to consider
- real, dimension(0:) :: potP, potN -- the proton,neutron potentials with p=pF. The dimension has to be identical to nuc%densTab(0: ,1:2).
- real, dimension(0:) :: potC -- The Coulomb potential (>0, in GeV)
OUTPUT
- nuc%densTab(0: ,1:2) is changed