TABLE OF CONTENTS
- 1. /checks
- 1.1. checks/Do_CheckDensity
- 1.2. checks/Do_CheckCoulomb
- 1.3. checks/Do_CheckFermiSurface
- 1.4. checks/Do_CheckRadius
- 1.5. checks/Do_CheckMomentumDensity
- 1.6. checks/Do_CheckEnergyLDA
- 1.7. checks/Do_OccupiedReal
- 1.8. checks/Do_OccupiedPert
- 1.9. checks/Do_CallEnergy
- 1.10. checks/Do_TachyonsReal
- 1.11. checks/Do_TachyonsPert
- 1.12. checks/TachyonIsBlocking
- 1.13. checks/ChecksSetDefaulSwitches
- 1.14. checks/ChecksSwitchRealRun
- 1.15. checks/ChecksSwitchGS
- 1.16. checks/Init
- 1.17. checks/Checks
- 1.18. checks/ChecksCallEnergy
- 1.19. checks/ChecksCallOccupied
- 1.20. checks/ChecksCallAll
- 1.21. checks/evaluateFermiSurface
- 1.22. checks/particleTrack
- 1.23. checks/particleTrack_byNumber
- 1.24. checks/occupied
- 1.25. checks/Occupation.dat
- 1.26. checks/checkPauli
- 1.27. checks/checkMomentumDensity
- 1.28. checks/checkDensity
- 1.29. checks/checkCoulomb
- 1.30. checks/checkPot
- 1.31. checks/checkRadius
- 1.32. checks/checkRadius_getAmplitude
- 1.33. checks/checkRadius_resetMinMax
- 1.34. checks/saveEpot
- 1.35. checks/checkEnergyLDA
- 1.36. checks/checkEnergyLDAWelke
- 1.37. checks/evaluateEnergy
- 1.38. checks/evaluateBindingEnergy_theiss
- 1.39. checks/checkTrajectories
- 1.40. checks/evaluateTimeStep
- 1.41. checks/evaluateTotal4Momentum_RMF
- 1.42. checks/CheckGridSize
/checks [ Modules ]
NAME
module check
PURPOSE
Here some routines are collected, which e.g. calculated the binding energy from the particle vector
checks/Do_CheckDensity [ Global module-variables ]
[ Top ] [ checks ] [ Global module-variables ]
SOURCE
logical, save :: Do_CheckDensity = .false.
PURPOSE
Flag to indicate, whether the density check routine should be called
checks/Do_CheckCoulomb [ Global module-variables ]
[ Top ] [ checks ] [ Global module-variables ]
SOURCE
logical, save :: Do_CheckCoulomb = .false.
PURPOSE
Flag to indicate, whether the Coulomb check routine should be called
checks/Do_CheckFermiSurface [ Global module-variables ]
[ Top ] [ checks ] [ Global module-variables ]
SOURCE
logical, save :: Do_CheckFermiSurface = .false.
PURPOSE
Flag to indicate, whether the fermi surface check routine should be called
checks/Do_CheckRadius [ Global module-variables ]
[ Top ] [ checks ] [ Global module-variables ]
SOURCE
logical, save :: Do_CheckRadius = .false.
PURPOSE
Flag to indicate, whether the radius check routine should be called
checks/Do_CheckMomentumDensity [ Global module-variables ]
[ Top ] [ checks ] [ Global module-variables ]
SOURCE
logical, save :: Do_CheckMomentumDensity = .false.
PURPOSE
Flag to indicate, whether the momentum density check routine should be called
checks/Do_CheckEnergyLDA [ Global module-variables ]
[ Top ] [ checks ] [ Global module-variables ]
SOURCE
logical, save :: Do_CheckEnergyLDA = .false.
PURPOSE
Flag to indicate, whether the local density approximation check routine should be called
checks/Do_OccupiedReal [ Global module-variables ]
[ Top ] [ checks ] [ Global module-variables ]
SOURCE
logical, save :: Do_OccupiedReal = .false.
PURPOSE
Flag to indicate, whether the occupation chek routine for the real particle vector should be called
checks/Do_OccupiedPert [ Global module-variables ]
[ Top ] [ checks ] [ Global module-variables ]
SOURCE
logical, save :: Do_OccupiedPert = .false.
PURPOSE
Flag to indicate, whether the occupation chek routine for the perturbative particle vector should be called
checks/Do_CallEnergy [ Global module-variables ]
[ Top ] [ checks ] [ Global module-variables ]
SOURCE
logical, save :: Do_CallEnergy = .false.
PURPOSE
Flag to indicate, whether the energy check routine should be called
checks/Do_TachyonsReal [ Global module-variables ]
[ Top ] [ checks ] [ Global module-variables ]
SOURCE
logical, save :: Do_TachyonsReal = .false.
PURPOSE
Flag to indicate, whether the tachyon check routine for the real particle vector should be called
checks/Do_TachyonsPert [ Global module-variables ]
[ Top ] [ checks ] [ Global module-variables ]
SOURCE
logical, save :: Do_TachyonsPert = .false.
PURPOSE
Flag to indicate, whether the tachyon check routine for the perturbative particle vector should be called
checks/TachyonIsBlocking [ Global module-variables ]
[ Top ] [ checks ] [ Global module-variables ]
SOURCE
logical, save :: TachyonIsBlocking = .false.
PURPOSE
select, whether the occurence of a 'tachyon' in the check routines will stop the code or not (Error messages will hopefully occur later in the code)
checks/ChecksSetDefaulSwitches [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine ChecksSetDefaulSwitches(eventtype)
PURPOSE
Set the default values of the switches "Do_Check???" according the given eventtype
NOTES
We repeat the specific codes snippet from GiBUU.f90, since not all conditions have been transported. Delete the commented lines after fulfilment.
checks/ChecksSwitchRealRun [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine ChecksSwitchRealRun(isReal)
PURPOSE
Set Do_OccupiedReal,Do_OccupiedPert
checks/ChecksSwitchGS [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine ChecksSwitchGS(isFrozen)
PURPOSE
Set the Flag, whether Ground State checks are only done at timestep 0 or for all timesteps.
checks/Init [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine init()
PURPOSE
Read the namelist "Checks".
checks/Checks [ Namelists ]
[ Top ] [ checks ] [ Namelists ]
NAME
NAMElist /Checks/
PURPOSE
Includes the following switches:
- Do_CheckDensity
- Do_CheckCoulomb
- Do_CheckFermiSurface
- Do_CheckRadius
- Do_CheckMomentumDensity
- Do_CheckEnergyLDA
- Do_OccupiedReal
- Do_OccupiedPert
- Do_CallEnergy
- Do_TachyonsReal
- Do_TachyonsPert
- TachyonIsBlocking
checks/ChecksCallEnergy [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine ChecksCallEnergy(time,realParticles)
PURPOSE
This is one of the upper level routines, called to check the energies.
Major switches: 'Do_CallEnergy',
checks/ChecksCallOccupied [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine ChecksCallOccupied(realParticles,pertParticles,text)
PURPOSE
This is one of the upper level routines, called to check the occupansies of the real and perturbative particle vectors.
Major switches: 'Do_OccupiedReal', 'Do_OccupiedPert'
checks/ChecksCallAll [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine ChecksCallAll(timestep,time,realParticles,pertParticles)
PURPOSE
This is one of the upper level routines, called to perform major checks as indicated by the flags 'Do_Check???'.
NOTES
Some test are done by other routines, as e.g. 'ChecksCallOccupied'
checks/evaluateFermiSurface [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine evaluateFermiSurface(fileName)
PURPOSE
checks/particleTrack [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine particleTrack(p,name)
PURPOSE
Prints the track of a particle
checks/particleTrack_byNumber [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine particleTrack_byNumber(parts,number,time,DoNew)
PURPOSE
Prints the track of a particle with given number
checks/occupied [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine occupied(particleVector,name)
PURPOSE
Prints information, how many entries and of every ensemble is actually occupied by particles and the percentage of maximum occupation
checks/Occupation.dat [ Output files ]
[ Top ] [ checks ] [ Output files ]
NAME
file Occupation.dat
PURPOSE
Show how many entries of every ensemble is actually occupied by particles. If it is close to the maximal possible (lengthPert for perturbative runs), one should increase parameter length_Perturbative in the namelist "input". Columns:
- 1: ensemble number
- 2: occupation
checks/checkPauli [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine checkPauli(realparticles)
PURPOSE
checks/checkMomentumDensity [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine checkMomentumDensity(realParticles,filename)
PURPOSE
checks/checkDensity [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine checkDensity(realParticles,filename)
PURPOSE
checks/checkCoulomb [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine CheckCoulomb(filename)
PURPOSE
checks/checkPot [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine CheckPot
PURPOSE
checks/checkRadius [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine checkRadius(teilchen, time)
PURPOSE
Average radius over all particle vector
checks/checkRadius_getAmplitude [ Functions ]
[ Top ] [ checks ] [ Functions ]
NAME
subroutine checkRadius_getAmplitude()
PURPOSE
return the value (MaximalRadius-MinimalRadius)
checks/checkRadius_resetMinMax [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine checkRadius_resetMinMax()
PURPOSE
reset the values of MaximalRadius and MinimalRadius
checks/saveEpot [ Functions ]
[ Top ] [ checks ] [ Functions ]
NAME
real function saveEpot(E,rho,i)
PURPOSE
Sums the potential Enrgy of all nucleons
INPUTS
- real :: E -- energy
- real :: rho -- density
- logical :: i -- .true.: store; .false.: retrieve and delete
OUTPUT
Energy in units of MeV
checks/checkEnergyLDA [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine checkEnergyLDA(teilchen)
PURPOSE
checks/checkEnergyLDAWelke [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine checkEnergyLDAWelke(teilchen)
PURPOSE
checks/evaluateEnergy [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine evaluateEnergy(teilchen,time)
PURPOSE
print the averaged kinetic energy sqrt(p^2+m^2) and also p0 per time step
checks/evaluateBindingEnergy_theiss [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine evaluateBindingEnergy_theiss(teilchen,time)
PURPOSE
Print binding energy of all nucleons, no Coulomb effects included
checks/checkTrajectories [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
NAME
subroutine checkTrajectories(teilchen)
PURPOSE
checks/evaluateTimeStep [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
PURPOSE
Compute time step using the frequency of two-body collisions
INPUTS
- integer, intent(in) :: iflag=1 ! use real-real counted collisions =2 ! use real-perturbative counted collisions
- real, intent(in) :: coll_num ! allow .leq.coll_num collisions/particle
- real, intent(in) :: delta_T_max ! maximal allowed time step (fm/c)
- real, intent(in) :: time ! current time (fm/c)
OUTPUT
- real, intent(out) :: delta_T_new ! calculated time step
NOTES
If gridSpacing(3) < delta_T_max, the upper limit on the time step is chosen as gridSpacing(3). This is important for high energy HIC due to Lorentz contraction along z-axis.
checks/evaluateTotal4Momentum_RMF [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
PURPOSE
Compute the total energy and momentum of the system within the relativistic mean field model.
INPUTS
- type(particle), dimension(:,:), intent(in) :: teilchen ! particle array
- real, intent(in) :: time ! current time (fm/c)
OUTPUT
Prints out the total 4-momentum vs time.
checks/CheckGridSize [ Subroutines ]
[ Top ] [ checks ] [ Subroutines ]
PURPOSE
Checks if real particles start to escape out of grid.
INPUTS
- type(particle), dimension(:,:) :: teilchen -- particle array
- real, :: time -- current time (fm/c)
- real, :: time_max -- max. time of simulation (fm/c)
- integer :: TheEventType -- Eventclass
- real :: gridSpacing -- size of cells
- real :: gridSize -- Size of the box