TABLE OF CONTENTS
- 1. /thermoDyn
- 1.1. thermoDyn/temperatureSwitch
- 1.2. thermoDyn/linearExtrapolation
- 1.3. thermoDyn/initThermoDynamics
- 1.4. thermoDyn/upDateTemperature
- 1.5. thermoDyn/temperatureAt
- 1.6. thermoDyn/tempe
- 1.7. thermoDyn/muAt
- 1.8. thermoDyn/integral
/thermoDyn [ Modules ]
NAME
module thermoDyn
PURPOSE
Includes routines for the determination of temperature T and chemical potential mu.
thermoDyn/temperatureSwitch [ Global module-variables ]
[ Top ] [ thermoDyn ] [ Global module-variables ]
SOURCE
integer,save :: temperatureSwitch = 1
PURPOSE
- 1=groundstate calculations (T=0,mu=E_F)
- 2=the full procedure
thermoDyn/linearExtrapolation [ Global module-variables ]
[ Top ] [ thermoDyn ] [ Global module-variables ]
SOURCE
logical,save :: linearExtrapolation=.true.
PURPOSE
- .true.= Use linear extrapolation for temperature between gridPoints
- .false.= Do not use it
thermoDyn/initThermoDynamics [ Namelists ]
[ Top ] [ thermoDyn ] [ Namelists ]
NAME
NAMELIST initThermoDynamics
PURPOSE
Includes the input switches:
thermoDyn/upDateTemperature [ Subroutines ]
[ Top ] [ thermoDyn ] [ 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.
thermoDyn/temperatureAt [ Functions ]
[ Top ] [ thermoDyn ] [ 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)).
thermoDyn/tempe [ Functions ]
[ Top ] [ thermoDyn ] [ 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".
thermoDyn/muAt [ Functions ]
[ Top ] [ thermoDyn ] [ 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
thermoDyn/integral [ Functions ]
[ Top ] [ thermoDyn ] [ 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. Switch=3 : \partial\rho/\partial\mu = (4/T)*Integral exp((E(p)-mu)/T)/(1+exp((E(p)-mu)/T))**2 (dp)**3)/(2pi)**3 over the full p-space. Switch=4 : rho=4*Integral exp((-E(p)+mu)/T) (dp)**3)/(2pi)**3 over the full p-space. Boltzmann limit. The first Factor of 4 is due to spin&isospin degeneracy.