KIS Solar Evolution Code
The solar evolution code at the KIS is an evolved version of the Hofmeister
et.al. stellar evolution code (see Hofmeister, E., Kippenhahn, R.,
Weigert, A.: Z.f.Aph.59, p.215-241,1964; Kippenhahn, R., Weigert, A.,
Hofmeister, E., Methods in Comp. Phys.7,p.129-191, 1967). It consists
of approximately fifty
Fortran 90
or Fortran 77 subroutines. The code is currently able to follow low mass evolution
(< 3 solar masses, < 14 solar radii) up to (but not covering) the helium flash ,
but has several restrictions.
The following physical ingredients are implemented:
- Opacity Tables:
- Equation of State:
- partial Ionization, Degeneracy, Electrostatic Correction
(see
Stix & Skaley)
- MHD tables
- OPAL 95 tables
- Convection:
- local mixing length / artificial overshoot
- local Canuto-Mazzitelli formulation (1991/92 and 1996)
- nonlocal integration with overshoot
(see
Skaley & Stix,
Shaviv & Salpeter)
- Energy and Element Production
- gravitational energy release
- hydrogen burning in equilibrium
- equilibrium / nonequilibrium CNO cycle
- lithium / beryllium reaction chains
- neutrino fluxes (see also John Bahcall's excellent
neutrino viewgraphs)
- Diffusion
Contact me to get the source code
and documentation.
A map of the program.
Home
last update 19. August 2004, U.Grabowski