XRB generator for helium burst setup

Historic generator file for APPROX-only XRB simulation of -accretor.

Generator

c This generator only runs APPROX
c
net 1 h1 he3 he4 n14 c12 o16 ne20 mg24
net 1 si28 s32 ar36 ca40 ti44 cr48 fe52
net 1 ni56 fe54 pn1 nt1
m nstar 1.00 fe54
c He star abundances
m acret 1.0 he4 0.0 c12 0.0 o16
c m acret 0.3 he4 0.2 c12 0.5 o16
c m acret 0. h1 0.980 he4 0. c12 0.02 n14 0. o16
c THIS GRID FOR He ACCRETION
g 0   2.0000e25  1 nstar  4.0e+8  1.0e+9
g 1   1.9000e25  1 nstar  4.0e+8  1.0e+9
g 40  1.0000e22  1 nstar  4.0e+8  1.0e+8
g 50  1.0000e21  1 nstar  4.0e+8  1.0e+8
g 51  8.0000e20  1 acret  2.0e+8  1.0e+8
g 54  2.0000e20  1 acret  1.0e+8  1.0e+6
g 55  0.         1 acret  5.0e+7  1.0e+4
dstat
c genburn  rphebg
p dtnew 1.e-4
p maxit 40
p dtcr .05
p dtct .05
p dtcd .10
p dtcq .1
p dtcdt .99
p nstop 1000000
p nedit 100000
p ndump 10
p ipup 2
p dyemult 50.
p dyqmult 2.
p dtcp .15
p yfloorx 3.e-3
p cenu 1.
p t7peek 1.e+50
p dtcut .1
p tfcrbu 2.
p tfctbu 10.
p fracneut .05
p radius0 1.0e+06
p summ0 2.8e+33
p xlum0 1.6e+34
p tnucmin 1.0e+99
p fclmax 1.e+99
p fclbu 1.e+99
p dtsmult 1.e+99
p dnratmax .25
p rnmin 1.e+6
p tnmin 1.e+4
p dnmin 2.e-5
p izonef 0
p idzonef 0
p jshell0 51
p rnmax 1.e+14
p tqsemin 3.e+9
p irtype 6
p fracrz1 .33
p fracrz2 .5
p abarratm 1.3
p frcsound .0
p convlim .0
p woversht 0.
p fmax0 .01
p fmax1 .015
p fmax2 .03
p nsdump 10
p iautoout 5
p iflgabar 0
p siqselim .02
p yfloorbx .003
p iwinsize 14001000
p wilsonmt -1.
p iold 1
p kaptab 0
p btempmin 1.1e8
p ncycqq 100000
p mazful -1
c p alpth 0.
p npixedit 1

p h1hdep -1.d0
p he4hedep -1.d0
p si28dep -1.d0
p npixedit 1
p lburn 1
p iplotb 2
p bmasslow 2.80000019999895D33
p convlim 1.
p frcsound 1.
p btempmin 1.d7
p tnucmin 1.d7
p accrate 1.75d-9
p minzone 51
p jp0 40
p irtype 4
p zonemmin 1.5d19
p zonemmax 1.d20
p lcout 10
p ncnvout 1
p fcrmax 1.d-8
p fctmax 1.d-8
p abunlim 1.d-4

p nadapb 1
c=======================================================================
c Now follows the command file
c=======================================================================
//*
c ------- 4U 0614+091 - --------
c base heating (MeV/nuclen)
o xqheat 0.45 def
c accretion rate relative to Eddington
o xledd 0.1 def
c hydrogen mass fraction
o xhfrac 0. def
c .... accretion rate 1.75D-8 (L/Ledd) * 1.7/(X + 1)
p accrate 1.75D-8
c .... substrate luminosity - accrate * 6.0816737e+43 * (Q/MeV)
c .... 1.0642929e+36 (L/Ledd) * (Q/MeV)
p xlum0 1.0642929e+36
c -------------------------
c ..... SCALE to He/C/O L_Edd accretion
o xeddf {1.7/(1.+xhfrac)} def
p accrate {xeddf} *
p xlum0 {xeddf} *
c .......
c substrate L, Q/MeV
p xlum0 {xqheat} *
c .......
c set fraction of Eddington accretion rate
p accrate {xledd} *
p xlum0 {xledd} *

c -------------------------
c get model in equilibrium
p ncnvout 0
p nstop 1000000000
p tnucmin 1.d10
p tnumin 1.d7
p accmass 1.d13
p optconv 0.67
p iaccadv 0
p iacceadv 1
c for APPROX ONLY
p jp0 0
p irtype 4
plot
c .........................
c MATCH TO accmass
@xm(jm)<1.01d18
p maxbak 20
p accdepth 1.d99
p iterbarm 999999
c .........................
@time>1.d17
zerotime
p toffset 0.
setcycle 0
cutbin
resetacc
p lburn 0
p dtnew 1.
p iaccadv 1
c use accdepth 5.d20 for He
c use accdepth 1.d20 for H
p accdepth 0.5d21
p tnucmin 1.d7
p izonef 1
p idzonef 1
p iazonef 0
p zonermax 10.
p zonemmax 1.d99
p ddmin 1.d4
c decretion
p decrate -1.D0
p idecmode 1
p jshell0 0
c use y coordinate for zoning
c p fmax0 1.
c p fmax1 1.
c p fmax2 1.
c p zonemmax 1.d99
c p zonermax 1.d99
c p zoneymax 0.02
c
c some other stuff
p pbound 5.d18
c
c rotation
c p centmult .667
c 600 Hz rotation?
c p angjacc 2.51327d15
c p nangdis 1
c p magnet 8
c p ipixtype 31700
p irtype 11
c p xmagfbr 1.d-10
c p xmagfbt 1.d-5