[mesa-users] MESA: doing something more than reporting...
Evan O'Connor
evanoc at caltech.edu
Tue Jan 18 20:34:05 EST 2011
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
More information about the Mesa-users
mailing list