Historic KEPLER Generator Files

25 solar mass star of solar composition

c s25sha--25 solar mass star of solar composition
c abundances from anders and grevesse (geochem et cosmochem 1988)
c semi-convection and convective overshoot mixing turned on everywhere
c outer region of the star (abar.ge.abarsemi (p 324) = 4.)
c note drmult (p 24) = 0.l, woversht (p 148) = 0.01
c note drmultlo (p 325) = 0.l, woverslo (p 326) = 0.01
c note c12(ag) rate is 0.5 times the cf88 current value
c
c box and id information:
box v98 weaver
c
c specify burn-genera tor-file name to turn on iso-topic co-processing
genburn sol160bg
c
c input for approx network (network #1)
net 1   h1  he3  he4  n14  c12  o16 ne20 mg24 si28 s32
net 1 ar36 ca40 ti44 cr48 fe52 ni56 fe54  pn1  nt1
c
c include ise and nse networks (networks #2 and #3)
isenet
c
c solar abundances (weight %):
m sol    70.57   h1 2.92ge-3  he3    27.52  he4  0.3069  c12
m sol   0.1109  n14   0.9618  o16   0.1753 ne20 0.06935 mg24
m sol 7.688e-2 si28 4.259e-2  s32 9.61ge-3 ar36
m sol 6.571e-3 ca40 2.947e-4 cr48 1.813e-3 fe52
m sol   0.1285 fe54
c
c initial grid (zone #, exterior mass(g), network #, temp(K), rho(g/cc)):
g   0 5.0e+34    1 sol 6.0e+6 1.e-1
g   1 4.9996e+34 1 sol 5.e+6  8.e-2
g   2 4.9990e+34 1 sol 5.e+6  8.e-2
g   3 4.9980e+34 1 sol 5.e+6  8.e-2
g   4 4.9962e+34 1 sol 5.e+6  8.e-2
g   5 4.9930e+34 1 sol 5.e+6  8.e-2
g   6 4.9900e+34 1 sol 5.e+6  8.e-2
g   7 4.9850e+34 1 sol 5.e+6  8.e-2
g  20 4.ge+34    1 sol 1.e+6  7.e-2
g  30 4.8e+34    1 sol 1.e+6  6.e-2
g  45 4.5e+34    1 sol 1.e+6  5.e-2
g  65 4.0e+34    1 sol 1.e+6  4.e-2
g  85 3.Oe+34    1 sol 3.e+5  3.5e-2
g 115 1.0e+34    1 sol 1.e+5  3.e-2
g 130 1.0e+33    1 sol 1.e+4  2.e-2
g 140 1.0e+32    1 sol 1.e+4  1.e-2
g 150 1.e+31     1 sol 1.e+4  6.e-3
g 160 0.0        1 sol 5.e+3  1.e-3
c
c adjust initial temperature to yield hydrostatic equilibrium
hstat
c
c reset default parameter values:
c
c time-step and back-up controls
p   6 .07
P   7 .05
p   8 .15
p   9 .07
P  46 .15
p  47 1.e-3
c
c problem termination criteria
p 158 999999
p 306 9.e+7
p 307 .05
c
c linear artificial viscosity coefficient (reset to 0.1 at zero-age ms)
p  13 1000.
c
c edit and dump controls
p  16 20
p  18 10
p 156 25
p 268 53
c
c semi convection and overshoot mixing coefficients
p  24 0.1
p 148 0.01
p 324 4.
P 325 0.1
P 326 0.01
c
c graphics parameters
p  42 10240750
p  64 50
p 113 31
c
c minimum nuclear burning temperature (K)
p  65 1.e+7
c
c rezoning criteria
p  83 1.e+4
p  84 1.e-4
p  86 0
P  87 1
P 138 .15
P 139 .5
p 150 .008
p 151 .01
p 152 .02
p 193 .1
p 195 .005
c
c ise control parameters
p 185 1.5e+9
p 185 .04
p 203 1.e+5
c
c c12(a,g) rate multipliers
p 227 0.5
p 228 0.5
c
c post-processor-dump control parameters
p 303 0.5
c
c set the time at which to make zero-age-main-sequence parameter changes
p 308 1.e+12
c
c turn on rezoner at the zero-age main sequence by reseting p 86
c to the value of p 309
p 309 1
c
c turn down the linear artificial viscosity at the zero-age main
c sequence by reseting p 13 to the value of p 310
p 310 .1
c
c set the core temperature at which to make pre-carbon-burning
c parameter changes
p 311 5.e+8
c
c raise floor on abundances considered in calculating the time-step
C just before carbon ignition by reseting p47 to the value of p312
p 312 .003
c
c finely zone the central 2-4 solar masses just before carbon ignition
c by reseting p195 to the value of p313 and p150 to the value of p314
p 313 .0015
p 314 .004
c
c list of dump variables (name, dump ratio, dezone ratio, adzone ratio):
c
c thermodynamic quantities
dump rn      .01   .01   .025
dump tn      .01   .01   .025
dump dn      .03   .03   .1
dump un      .03   .03   .1
dump pn      .03   .03   .1
dump sgain   .05   .05   .15
dump sloss   .05   .05   .15
dump sneut   .05   .05   .15
c
c convection sentinels
dump convect .1    -1.   0.
c
c entropies
dump sig     .01   .01   .025
dump sigi    .01   .01   .025
dump sige    .01   .01   .025
dump sigr    .01   .01   .025
dump sigp    .01   .01   .025
dump sigion  .01   .01   .025
c
c approx and ise ions
dump ionnt1  .02   .02   .06
dump ionh1   .02   .02   .06
dump ionpn1  .02   .02   .06
dump ionhe3  .02   .02   .06
dump ionhe4  .02   .02   .06
dump ionc12  .02   .02   .06
dump ionn14  .02   .02   .06
dump iono16  .02   .02   .06
dump ionne20 .02   .02   .06
dump ionmg24 .02   .02   .06
dump ionsi28 .02   .02   .06
dump ions32  .02   .02   .06
dump ionar36 .02   .02   .06
dump ionca40 .02   .02   .06
dump ionti44 .02   .02   .06
dump ioncr48 .02   .02   .06
dump ionfe52 .02   .02   .06
dump ionfe54 .02   .02   .06
dump ionni56 .02   .02   .06
dump ionfe56 .02   .02   .06
dump ion'fe' .02   .02   .06
dump ionye   .001  .001  .0025
dump ionyq   .003  .003  .01
dump ionebO  .001  .001  .0025
dump ionyf   .003  .003  .01
c
c parameters
dump qparm   1.    1.    1.
dump parm    1.    1.    1.
c
c misc. zonal edit quantities
dump netnum  .1    .1    .3
dump en      .03   .03   .1
dump zn      .03   .03   .1
dump etan    .01   .01   .025
dump abar    .01   .01   .025
dump zbar    .01   .01   .025
dump xkn     .02   .02   .06
dump gamma   .001  .001  .0025
dump xln     .03   .03   .1
dump taumix  .05   .05   .15
dump mixlen  .05   .05   .15
c
c weak rates and neutrino losses
dump wrate   .05   .05 .15
dump rectot  .05   .05   .15
dump rpdtot  .05   .05   .15
dump redtot  .05   .05   .15
dump snuw    .05   .05   .15
dump snubps  .05   .05   .15
c
c burn variables
dump yeburn  .001 .001 .0025
dump sburn   .05 .05 .15
c
c burn isotopes
dump isont1   .02   .02   .06
dump isohl    .02   .02   .06
dump isoh2    .02   .02   .06
dump isohe3   .02   .02   .06
dump isohe4   .02   .02   .06
dump isoli7   .02   .02   .06
dump isobe7   .02   .02   .06
dump isoc12   .02   .02   .06
dump isoc13   .02   .02   .06
dump ison14   .02   .02   .06
dump isoo16   .02   .02   .06
dump isoo18   .02   .02   .06
dump isone20  .02   .02   .06
dump isone22  .02   .02   .06
dump isona22  .02   .02   .06
dump isomg24  .02   .02   .06
dump isoa126  .02   .02   .06
dump isosi28  .02   .02   .06
dump isosi30  .02   .02   .06
dump isos32   .02   .02   .06
dump isos34   .02   .02   .06
dump isoar36  .02   .02   .06
dump isoca40  .02   .02   .06
dump isoti44  .02   .02   .06
dump isofe52  .02   .02   .06
dump isofe54  .02   .02   .06
dump isofe56  .02   .02   .06
dump isofe58  .02   .02   .06
c
c
c
c special ascii edit requests:
c zedit 1 10 sig sige sigi sigr sigion sigp
c
c
c parameter changes before silicon ignition:
c (when o16.lt. 5 weight % and tn(l).gt.tqselim (p184) --use p307)
c
c tighten up time-step controls on radius, temp, and density changes
c p 6 .02
c p 7 .02
c p 8 .02
c
c
c relax convergence criteria for radius and temperature
c p 11 1.e-5
c p 12 1. e- 5
c
c tolerate more convergence error in temperature and density
c before backing up
c p 54 10.
c p 55 10.
c
c make code insensitive to convergence errors in luminosity
c p 70 1.d+99
c p 73 1.d+99.
c :
c tolerate larger abundance changes before backing up
c p 206 3.e-3
c
c turn off stop at oxygen depletion
c p 307 -1.e+99
c
c zero the problem time to achieve more precision
c zerotime

Sample KEPLER Generator Variatons

Listed below are typical variations in the sample generator given above. Only (non-comment) cards, or sets of cards, that differ from the original generator are listed.

LMC Abundances

c lmc abundances (weight %):
m lmc    74.60   h1 2.929e-3  he3    25.00 he4  0.0522 c12
m lmc   0.0122  n14   0.2500  o16   0.0561 ne20 0.0354 mg24
m lmc 4.300e-2 si28 1.830e-2  s32 5.290e-3 ar36
m Imc 3.020e-3 ca40 1.240e-4 cr48 7.610e-4 fe52
m lmc   0.0797 fe54

Restricted Semiconvection

c semi convection and overshoot mixing coefficients
p  24 1.e-4
p 148 0.0
p 324 4.
p 325 1.e-4
p 326 0.0

Nominal deJager Mass Loss Rate

c mass loss parameters
p 220 1.
p 225 .01

8 Solar Mass Star

c initial grid (zone #, exterior mass(g), network #, ternp(K), rho(g/cc)):
g   0 1.6e+34    1 sol 6.0e+6 1.e-1
g   1 1.5996e+34 1 sol 5.e+6  8.e-2
g   2 1.5990e+34 1 sol 5.e+6  8.e-2
g   3 1.5980e+34 1 sol 5.e+6  8.e-2
g   4 1.5962e+34 1 sol 5.e+6  8.e-2
g   5 1.5930e+34 1 sol 5.e+6  8.e-2
g   6 1.5900e+34 1 sol 5.e+6  8.e-2
g   7 1.5850e+34 1 sol 5.e+6  8.e-2
g  20 1.5e+34    1 sol 1.e+6  7.e-2
g  30 1.4e+34    1 sol 1.e+6  6.e-2
g  50 1.2e+34    1 sol 1.e+6  5.e-2
g  70 1.0e+34    1 sol 1.e+6  4.e-2
g 100 5.0e+33    1 sol 1.e+5  3.e-2
g 130 1.0e+33    1 sol 1.e+4  2.e-2
g 140 1.0e+32    1 sol 1.e+4  1.e-2
g 150 1.e+31     1 sol 1.e+4  6.e-3
g 160 0.0        1 sol 5.e+3  1.e-3
c
c rezoning criteria
p  83 1. e+4
p  84 1.e-4
p  86 0
p  87 1
p 138 .4
p 139 .5
p 150 .008
p 151 .01
p 152 .02
p 193 .27
p 195 .005
c
c finely zone the central 2-4 solar masses just before carbon ignition
c by reseting p195 to the value of p313 and p150 to the value of p314
p 313 .003
p 314 .007

15 Solar Mass Star

c initial grid (zone #, exterior mass(g), network #, temp(K), rho(g/cc)):
g   0 3.0e+34    1 sol 6.0e+6 1.e-1
g   1 2.9996e+34 1 sol 5.e+6  8.e-2
g   2 2.9990e+34 1 sol 5.e+6  8.e-2
g   3 2.9980e+34 1 sol 5.e+6  8.e-2
g   4 2.9962e+34 1 sol 5.e+6  8.e-2
g   5 2.9930e+34 1 sol 5.e+6  8.e-2
g   6 2.9900e+34 1 sol 5.e+6  8.e-2
g   7 2.9850e+34 1 sol 5.e+6  8.e-2
g  20 2.ge+34    1 sol 1.e+6  7.e-2
g  30 2.8e+34    1 sol 1.e+6  6.e-2
g  45 2.5e+34    1 sol 1.e+6  5.e-2
g  70 2.0e+34    1 sol 1.e+6  4.e-2
g 115 1.0e+34    1 sol 1.e+5  3.e-2
g 130 1.0e+33    1 sol 1.e+4  2.e-2
g 140 1.0e+32    1 sol 1.e+4  1.e-2
g 150 1.e+31     1 sol 1.e+4  6.e-3
g 160 0.0        1 sol 5.e+3  1.e-3
c
c rezoning criteria
p  83 1.e+4
p  84 1.e-3
p  86 0
p  87 1
p 138 .24
P 139 .5
p 150 .008
p 151 .01
P 152 .02
P 193 .16
p 195 .005
c
c finely zone the central 2-4 solar masses just before carbon ignition
c by reseting p195 to the value of p313 and p150 to the value of p314
p 313 .002
p 314 .005

40 Solar Mass Star

c initial grid (zone #, exterior mass(g), network #, temp(K), rho(g/cc):
g   0 8.0000e+34 1 sol 6.0e+6 l.e-l
g   1 7.9996e+34 1 sol 5.e+6  8.e-2
g   2 7.9990e+34 1 sol 5.e+6  8.e-2
g   3 7.9980e+34 1 sol 5.e+6  8.e-2
g   4 7.9962e+34 1 sol 5.e+6  8.e-2
g   5 7.9930e+34 1 sol 5.e+6  8.e-2
g   6 7.9900e+34 1 sol 5.e+6  8.e-2
g   7 7.9850e+34 1 sol 5.e+6  8.e-2
g  20 7.9e+34    1 sol 1.e+6  7.e-2
g  30 7.8e+34    1 sol 1.e+6  6.e-2
g  45 7.5e+34    1 sol 1.e+6  5.e-2
g  57 7.0e+34    1 sol 1.e+6  4.e-2
g  70 6.0e+34    1 sol 1.e+6  4.e-2
g 115 1.0e+34    1 sol 1.e+5  3.e-2
g 130 1.0e+33    1 sol 1.e+4  2.e-2
g 140 1.0e+32    1 sol 1.e+4  1.e-2
g 150 1.e+31     1 sol 1.e+4  6.e-3
g 160 0.0        1 sol 5.e+3  1.e-3
c
c set very small surface boundary pressure to discourage mass loss
p  69 3.
c
c rezoning criteria
p  83 1.e+4
p  84 1.e-4
p  86 0
p  87 1
p 138 .1
p 139 .5
p 150 .005
p 151 .01
p 152 .02
p 193 .065
p 195 .003
c
c finely zone the central 2-4 solar masses just before carbon ignition
c by reseting p195 to the value of p313 and p150 to the value of p314
p 313 .00085
p 314 .00225