[mesa-users] Issue with ns_he in test_suite

[Joe Hughto]

Hello,

I am trying to start a project with mesa to look at x-ray bursts on
accreting neutron stars.  I assumed that if I was to run the the system in
test_suite/ns_he then I would get the same output that was in the
instrument paper.  However, after 27 steps, the code quits and gives the
following:

DATE: 2012-03-12
TIME: 13:01:57
 read inlist_ns_he
load saved model ns_1.4M.mod

 change to cno_extras_plus_fe56.net
 number of species =          22
 read log_columns.list
 read profile_columns.list
 net name cno_extras_plus_fe56.net
          which_atm_option: grey_and_kap
                            xmstar/mstar    4.3562806036287768D-12
                              xmstar (g)    1.2131718727433707D+22
                            M_center (g)    2.7848799999878681D+33
                             xmstar/Msun    6.0987928450802873D-12
                           M_center/Msun    1.3999999999939012D+00
                           R_center (cm)    1.0000000000000000D+06
                           R_center/Rsun    1.4368228972096900D-05
                            core density    6.6484112687374000D+14
 chem_isotopes_filename isotopes.data_real
 kappa_file_prefix gs98
 kappa_lowT_prefix lowT_fa05_gs98
   eos_file_prefix mesa
                         OMP_NUM_THREADS           1


__________________________________________________________________________________________________________________________________________________

       step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass
H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr   pts  retry
   lg_dt_yr    lg_Dcntr    lg_R       lg_L3a    lg_Lneu     lg_Mdot
 H_poor     He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
     age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf
He_poor    C_cntr     Ne_cntr    Z_cntr     Z_avg     v_div_cs     dt_limit
__________________________________________________________________________________________________________________________________________________

         10   8.069069  2.049E+06 -99.000000 -68.808059   1.400000
0.000000   0.000000   0.000000   1.000000   0.000000  32.492604    893
 0
  -8.511350  -1.972256  -4.842337 -68.808059  -9.429144  -9.769551
1.400000   0.000000   0.000000   0.000000   0.000000  57.244719      5
 0
 1.5985E-08  23.369543   0.514568 -99.000000  13.720536  -0.378808
1.400000   0.000000   0.000000  1.000E+00  1.000E+00 -0.647E-07
 varcontrol

         20   8.074257  1.890E+06 -99.000000 -57.371656   1.400000
0.000000   0.000000   0.000000   1.000000   0.000000  32.102577    902
 0
  -7.719538  -1.972256  -4.842338 -57.371656  -9.423835  -9.769551
1.400000   0.000000   0.000000   0.000000   0.000000  56.562059      6
 0
 1.1195E-07  23.369537   0.374310 -99.000000  13.638372  -0.424759
1.400000   0.000000   0.000000  1.000E+00  1.000E+00 -0.217E-07
 varcontrol


stop because envelope mass < envelope_mass_limit
   -0.2220446049E-15    0.0000000000E+00    0.1000000000+100
 0.1000000000+100
         27   8.081876  1.760E+06 -99.000000 -50.748326   1.400000
 -0.000000   0.000000   0.000000   1.000000   0.000000  31.538523    908
   0
  -7.165269  -1.671227  -4.842338 -50.748326  -9.408800  -9.769551
1.400000   0.000000   0.000000   0.000000   0.000000  55.574405      6
 0
 4.0759E-07  23.369535   0.250906 -99.000000  13.563874  -0.467850
1.400000   0.000000   0.000000  1.000E+00  1.000E+00  0.175E-08
 varcontrol


save profile LOGS/log1.data     27

save photos/x027

runtime (minutes), retries, ba          0.447383         12.000000
0         0

         27   8.081876  1.760E+06 -99.000000 -50.748326   1.400000
 -0.000000   0.000000   0.000000   1.000000   0.000000  31.538523    908
   0
  -7.165269  -1.671227  -4.842338 -50.748326  -9.408800  -9.769551
1.400000   0.000000   0.000000   0.000000   0.000000  55.574405      6
 0
 4.0759E-07  23.369535   0.250906 -99.000000  13.563874  -0.467850
1.400000   0.000000   0.000000  1.000E+00  1.000E+00  0.175E-08
 varcontrol

DATE: 2012-03-12
TIME: 13:02:38

What am I doing wrong?

-Joe
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120312/34b783af/attachment.html>