eos_idealplusrad

Module

Description

Ideal gas equation of state plus radiation pressure, assumes

inputs are in cgs units

References:

Stellar Structure and Evolution (2nd Edition) (Kippenhahn, Weigert, Weiss)

Owner:

Mike Lau

Runtime parameters:

None

Dependencies:

physcon

Quick access

Variables:

tolerance

Routines:

egas_from_rhot(), erad_from_rhot(), get_idealgasplusrad_tempfrompres(), get_idealplusrad_enfromtemp(), get_idealplusrad_pres(), get_idealplusrad_rhofromprest(), get_idealplusrad_spsoundi(), get_idealplusrad_temp()

Needed modules

• physcon (rg(), radconst()): Physical and mathematical constants (as in sphNG)

Variables

• eos_idealplusrad/tolerance [real,private/parameter/optional/default=1e-15]

Subroutines and functions

subroutine  eos_idealplusrad/get_idealplusrad_temp(rhoi, eni, mu, tempi, ierr)

Parameters:

• rhoi [real,in]

• eni [real,in]

• mu [real,in]

• tempi [real,inout]

• ierr [integer,out]

subroutine  eos_idealplusrad/get_idealplusrad_pres(rhoi, tempi, mu, presi)

Parameters:

• rhoi [real,in]

• tempi [real,in] :: Eq 13.2 (Kippenhahn et al.)

• mu [real,in]

• presi [real,out]

subroutine  eos_idealplusrad/get_idealplusrad_spsoundi(rhoi, presi, eni, spsoundi, gammai)

Parameters:

• rhoi [real,in]

• presi [real,in]

• eni [real,in]

• spsoundi [real,out]

• gammai [real,out]

subroutine  eos_idealplusrad/get_idealgasplusrad_tempfrompres(presi, rhoi, mu, tempi)

Parameters:

• presi [real,in]

• rhoi [real,in]

• mu [real,in]

• tempi [real,inout]

subroutine  eos_idealplusrad/get_idealplusrad_enfromtemp(densi, tempi, mu, eni)

Parameters:

• densi [real,in]

• tempi [real,in]

• mu [real,in]

• eni [real,out]

Call to:

egas_from_rhot(), erad_from_rhot()

function  eos_idealplusrad/egas_from_rhot(tempi, mu)

Parameters:

• tempi [real,in]

• mu [real,in]

Return:

egasi [real]

Called from:

get_idealplusrad_enfromtemp()

function  eos_idealplusrad/erad_from_rhot(densi, tempi, mu)

Parameters:

• densi [real,in]

• tempi [real,in]

• mu [real,in]

Return:

eradi [real]

Called from:

get_idealplusrad_enfromtemp()

subroutine  eos_idealplusrad/get_idealplusrad_rhofromprest(presi, tempi, mu, densi)

Parameters:

• presi [real,in]

• tempi [real,in]

• mu [real,in]

• densi [real,out]