TABLE OF CONTENTS
- 1. /thermoDynamics
/thermoDynamics [ Modules ]
NAME
module thermoDynamics
PURPOSE
Includes routines for the determination of temperature T and chemical potential mu.
thermoDynamics/temperatureSwitch [ Global module-variables ]
[ Top ] [ thermoDynamics ] [ Global module-variables ]
SOURCE
integer,save :: temperatureSwitch = 1
PURPOSE
- 1=groundstate calculations (T=0,mu=E_F)
- 2=the full procedure
thermoDynamics/linearExtrapolation [ Global module-variables ]
[ Top ] [ thermoDynamics ] [ Global module-variables ]
SOURCE
logical,save :: linearExtrapolation=.true.
PURPOSE
- .true.= Use linear extrapolation for temperature between gridPoints
- .false.= Do not use it
thermoDynamics/initThermoDynamics [ Namelists ]
[ Top ] [ thermoDynamics ] [ Namelists ]
NAME
NAMELIST initThermoDynamics
PURPOSE
Includes the input switches:
thermoDynamics/upDateTemperature [ Subroutines ]
[ Top ] [ thermoDynamics ] [ Subroutines ]
PURPOSE
Calculation of local 'temperature'.
upDateTemperature/evaluatePSquared [ Subroutines ]
[ Top ] [ upDateTemperature ] [ Subroutines ]
PURPOSE
Calculation of pSquared(x,y,z). This is stored in the field psquared(-gridPoints(1:3):gridPoints(1:3)). Usage of the smearing weigths to smear particles' momentum over gridpoints in their neighborhood.
upDateTemperature/evaluateTemperature [ Subroutines ]
[ Top ] [ upDateTemperature ] [ Subroutines ]
PURPOSE
Calculation of temperature(x,y,z). This is stored in the field temperature(-gridPoints(1:3):gridPoints(1:3)). Use pSquared(x,y,z) and rho(x,y,z) to determine the temperature.
thermoDynamics/temperatureAt [ Functions ]
[ Top ] [ thermoDynamics ] [ Functions ]
PURPOSE
Evaluates temperature at some space point r. Therefore it uses the values which are stored in the field temperature(-gridPoints(1:3):gridPoints(1:3)).
thermoDynamics/tempe [ Functions ]
[ Top ] [ thermoDynamics ] [ Functions ]
NAME
real function tempe (p2av, rho)
INPUTS
- real, intent(in) :: p2av ! average momentum squared <p**2> in GeV**2
- real, intent(in) :: rho ! baryon density in fm^-3
OUTPUT
- temperature in GeV
PURPOSE
Evaluates temperature as a function of rho(0) and <p**2>. At first call it initializes a field "temSave" which holds this information. Therefore no further calculation is necessary after initializing this field. It is generated in the subroutine "initTempe".
thermoDynamics/muAt [ Functions ]
[ Top ] [ thermoDynamics ] [ Functions ]
PURPOSE
Determine chemical potential as function of temperature and rho. Everything in GeV!!!
INPUTS
real :: rho , density in GEV**3 real :: temp , temperature in GeV
thermoDynamics/integral [ Functions ]
[ Top ] [ thermoDynamics ] [ Functions ]
NAME
real function integral(temperature,mu,Switch)
PURPOSE
Evaluates for a gas of degenerate neutrons and protons of given "temperature" and given chemical potential "mu" the following: Switch=1 : rho=<1>=4*Integral(1/1+exp((E(p)-mu)/T) (dp)**3)/(2pi)**3 over the full p-space. Switch=2 : <p**2>=4*Integral(p**2/1+exp((E(p)-mu)/T) (dp)**3)/(2pi)**3 over the full p-space. The first Factor of 4 is due to spin&isospin degeneracy.