gibuu is hosted by Hepforge, IPPP Durham
GiBUU

TABLE OF CONTENTS


/thermoDyn [ Modules ]

[ Top ] [ 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.