[mesa-users] MESA: doing something more than reporting...

[Evan O'Connor]

Hello all,

I just ran into something weird.  For the situation I am looking at, when I switch on 

report_all_dt_limits = .true.

the code crashes, if it is false then it is fine.  some of the error message is below.  

>>>>>>>>>>>>>>>
       1510   8.180391  1.244E+04   4.636411   4.755512  14.674456  11.547502   0.000000   0.012976   0.280000   0.477836  -3.242907    705      1
   2.248222   3.082165   1.683802   3.145216   3.469588  -6.929051   3.126954   0.979459   0.000502   0.020000   0.502238   0.041347      2      1
 1.1563E+07  19.104688   4.700862   4.088597   2.141492 -10.288952   0.000000   0.000728   0.004313   0.000396  1.993E-02  0.260E-04    varcontrol
 
                        dX n14  623    0.2006471728E-02    0.2000000000E-02    0.8955172923E+00    0.1445299224E-02
                             s% D_mix(k)         623    0.0000000000000000E+00
                                   X_old         623    1.4982242077306349E-02
                              s% xa(j,k)         623    1.2975770349680954E-02
                                      dX         623   -2.0064717276253950E-03
                                   D_mix         618    0.0000000000000000E+00    4.7017707893630550E-03
                                   D_mix         619    0.0000000000000000E+00    3.9011863067364982E-03
                                   D_mix         620    0.0000000000000000E+00    3.1006018241099415E-03
                                   D_mix         621    0.0000000000000000E+00    1.4994328588572581E-03
                                   D_mix         622    0.0000000000000000E+00    7.4971642942862904E-04
                                   D_mix         623    0.0000000000000000E+00    0.0000000000000000E+00
                                   D_mix         624    1.2954886313713439E+13   -1.4994328588571471E-03
                                   D_mix         625    1.6061308510682049E+13   -2.6642224508062681E-03
                                   D_mix         626    1.8407963009267238E+13   -3.8290120427557223E-03
                                   D_mix         627    2.0443522267985812E+13   -4.9938016347049544E-03
                                   D_mix         628    2.2259202775794098E+13   -6.1585912266542975E-03
 
check_dX
>>>>>>>>>>>>>>>>>

For me this is occurring right after helium burning in a 15Msun star (when the Ne20 and N14 levels are changing).  For reference, the inlist_15_z2m2 is identical to mesa/start/test_suite/15_si_burn/inlist_15_z2m2 except this flag is changed and I do not load the test_suite model but rather start from pre-main sequence).  I've also seen this occur in the same place in stars up to 19Msun.

It seems like the extra reporting causes mesa to fail.  I'm not sure why, thoughts?

Evan