[Mesa-users] ?==?utf-8?q? Timestep issue at the He flash stage
Dreau Guillaume
Guillaume.Dreau at obspm.fr
Fri Feb 5 17:57:05 UTC 2021
Hello Bill,
Ah yes, of course! Attached is the model saved right after the termination of inlist_to_start_he_core_flash that I sent earlier. It is approximately 100 steps before the timestep stopping condition. You can rerun it with inlist_to_end_core_he_burn, I send it to you in this email. I checked that I could reproduce the same problem from this starting model.
Have a good day!
Sincerely,
Guillaume
*******************************************************************************************************************
Guillaume Dréau
PhD student 2nd year – Instructor
Speciality Stellar Physics
LESIA - Observatoire de Paris
Building 14, office 107
5 place Jules Janssen
92195 Meudon Cedex – France
http://www.lesia.obspm.fr/
*******************************************************************************************************************
Le Vendredi, Février 05, 2021 17:26 CET, Bill Paxton <paxton at kitp.ucsb.edu> a écrit:
i should have mentioned that you need to double check that you can start a run from the new saved model from ust before the problem and recreate the same (or similar) problem from that starting model. you probably already thought of that, but i wanted to remind everyone of that important step! b On Feb 5, 2021, at 8:20 AM, Bill Paxton <paxton at kitp.ucsb.edu> wrote: Hi Guillaume, I’ll take a look, but it would be a big help if you’d send a saved model from 100 steps or so before the trouble starts so I don’t have to redo the run up to that point. Cheers,Bill On Feb 5, 2021, at 5:44 AM, Dreau Guillaume <Guillaume.Dreau at obspm.fr> wrote: Dear MESA users,
I am glad to know how flexible MESA is to model stellar evolution up to the asymptotic giant branch! With the latest version of MESA (mesa-r15140), I aim at studying the structure differences between red giants and asymptotic giants. To this end, I started using the test suite case 1M_pre_ms_to_wd and I adjusted the inlists to take into account the micropĥysics and macrophysics adapted to my study. Attached are my inlists. I am facing an issue at the tip of the red giant branch. In my inlists, when I set atm_T_tau_opacity = 'varying' (default is ‘fixed’), the model is no longer able to pass the He flash stage. It seems that the timestep becomes too small, reaching min_timestep_limit that is the smallest timestep allowed before the run stops. I tried to decrease that timestep limit as much as possible, but it does not fix the issue. I also tried to play with the parameters atm_T_tau_errtol and atm_T_tau_max_steps, but without success. Below are the messages I get that seem to be linked to the timestep stopping condition. Other surprising thing is that when I choose a lower mass (0.9 M_sun with the same other inputs), I end up having the same issue but also with default option atm_T_tau_opacity = 'fixed'. Then, it does not seem to be exclusively due to the atmosphere inputs. It seems that the residuals are too high so that the timestep does not stop decreasing. I tried setting max_resid_jump_limit = 1d10 (default is 1d06), but it does not help to fix the issue.
Do you have an idea about how to fix this and how to pass the He flash?
I hope I didn’t forget anything to share with you.
Many thanks in advance!
Sincerely,
Guillaume
*******************************************************************************************************************
Guillaume Dréau
PhD student 2nd year – Instructor
Speciality Stellar Physics
LESIA - Observatoire de Paris
Building 14, office 107
5 place Jules Janssen
92195 Meudon Cedex – France
http://www.lesia.obspm.fr/
*******************************************************************************************************************
------------------------------- Timestep stopping conditions -------------------------------
run inlist_to_end_core_he_burn_header
DATE: 2021-02-04
TIME: 15:58:35
read inlist_to_end_core_he_burn
read inlist_astero_search_controls
nn_redist 7000
call srdist to read redistrb.c.pruned.in
--------------------------------------------------
Entering srdist
nn,icnmsh?
601 0
icase,icvzbn,nsmth,xdisc,dlxdsc,dlgrmx,cacvzb?
0 0 0 0.0000000000000000 1.0000000474974513E-003 10.000000000000000 0.0000000000000000
cg,cx,ca,cdgr,cddgr,cdg,alphsf,adda4,accrmn?
1.0000000000000000 40.000000000000000 9.9999997764825821E-003 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000
nout,cn,irsu,unew,iresa4 ?
100 2.0000000000000000 0 3.0000000000000000 0
nmodel,kmodel,itsaml,ioldex,idsin,idsout?
0 0 0 0 2 3
Using new discontinuity criterion, ndisc = 60
End parameter input
End parameter input
Exiting srdist
--------------------------------------------------
back from 1st call on srdist
call adipls to read adipls.c.pruned.in
*************************************************************
In this version, the maximum number of mesh points is 20000
Maximum number of modes for internal storage is 20000
*************************************************************
--------------------------------------------------
Entering adipls with parameter passing
i_paramset = 1
Files needed:
Model input on unit imlds (default 2)
Trial frequency input on unit itrds (default 10). Optional.
Grand summary on unit idsgsm (default 11).
Short summary on unit idsssm (default 15).
Eigenfunction on unit idsefn (default 4) (optional).
Rotational kernel on unit idsrkr (default 12) (optional).
Gamma1 kernel on unit idsgkr (default 13) (optional).
Printed output on unit istdpr, if not default
Input format: <unit number> <file name>
input ends with EOF or a line containing -1.
if file name is given as 0, /dev/null is used for the file.
files set in s/r ofiles:
2 'amdl.mesa ' @
9 'ttt.adipls.prt ' @
11 'agsm.mesa ' @
15 'ttt.adipls.ssm ' @
16 'ttt.adipls.fsm ' @
-1 '' @
----------------------------------------------
cntrd?
mod.osc.int.out.dgn
ifind,xmod,imlds,in,nprmod?
-1 0.0000000000000000 2 1 0
xtrnct,ntrnsf,imdmod?
0.0000000000000000 0 0
el,nsel,els1,dels,dfsig1,dfsig2,nsig1,nsig2?
1.0000000000000000 0 10.000000000000000 1.0000000000000000 0.0000000000000000 0.0000000000000000 1 1
itrsig,sig1,istsig,inomd1,itrds?
0 10.000000000000000 1 1 10
dfsig,nsig,iscan,sig2?
0.0000000000000000 1 1 0.0000000000000000
eltrw1, eltrw2, sgtrw1, sgtrw2?
0.0000000000000000 -1.0000000000000000 0.0000000000000000 -1.0000000000000000
cgrav?
6.6732000000000005E-008
iplneq,iturpr,icow,alb?
0 0 0 1.0000000000000000
istsbc,fctsbc,ibotbc,fcttbc?
1 0.0000000000000000 0 0.0000000000000000
mdintg,iriche,xfit,fcnorm,eps,epssol,itmax,dsigre?
1 0 0.50000000000000000 0.50000000000000000 2.0000000000000001E-013 1.0000000000000001E-005 8 0.0000000000000000
fsig,dsigmx,irsevn,xmnevn,nftmax,itsord?
1.0000000000000000E-003 0.10000000000000001 2 0.0000000000000000 3 0
istdpr,nout,nprcen,irsord,iekinr?
6 50 0 0 0
iper,ivarf,kvarf,npvarf,nfmode?
0 1 2 0 0
irotkr,nprtkr,igm1kr,npgmkr,ispcpr?
0 0 0 50 0
icaswn, sigwn1, sigwn2, frqwn1, frqwn2, iorwn1, iorwn2, frlwn1, frlwn2?
-1 0.0000000000000000 -1.0000000000000000 0.0000000000000000 -1.0000000000000000 0 -1 0.0000000000000000 -1.0000000000000000
End parameter input
back from 1st call on adipls
search_type == use_first_values
version_number 15140
read inlist_to_end_core_he_burn
enter astero_extras_controls
load saved model start_he_core_flash.mod
read /home/gdreau/Documents/mesa_versions/mesa-r15140/star/test_suite/pre_ms_to_wd_test_reference/history_columns.list
read /home/gdreau/Documents/mesa_versions/mesa-r15140/star/test_suite/pre_ms_to_wd_test_reference/profile_columns.list
net name o18_and_ne22.net
atm_option: T_tau
set_cumulative_energy_error 0.0000000000000000D+00
kap_option a09
kap_CO_option a09_co
kap_lowT_option AESOPUS
OMP_NUM_THREADS 6
__________________________________________________________________________________________________________________________________________________
step lg_Tmax Teff lg_LH lg_Lnuc Mass H_rich H_cntr N_cntr Y_surf eta_cntr zones retry
lg_dt_yr lg_Tcntr lg_R lg_L3a lg_Lneu lg_Mdot He_core He_cntr O_cntr Z_surf gam_cntr iters
age_yr lg_Dcntr lg_L lg_LZ lg_Lphoto lg_Dsurf C_core C_cntr Ne_cntr Z_cntr v_div_cs dt_limit
__________________________________________________________________________________________________________________________________________________
5390 7.963517 2839.875 3.386714 3.389298 0.771023 0.301100 0.000000 0.006654 0.274608 21.289283 1924 0
2.896495 7.872970 2.311806 1.161107 2.221894 -6.980492 0.469923 0.986760 0.001985 0.013359 0.772990 8
1.2692E+10 5.967407 3.391502 -1.340560 -99.000000 -8.810068 0.000000 0.000047 0.001149 0.013240 0.000E+00 max increase
rel_E_err 4.8062360872066043D-10
log_rel_run_E_err -8.3486370097780451
5399 7.973993 2837.996 3.388011 3.392959 0.770304 0.300147 0.000000 0.006654 0.274608 21.293579 1913 0
2.976965 7.872687 2.312959 1.446040 2.223226 -6.978474 0.470156 0.986760 0.001985 0.013359 0.773299 9
1.2692E+10 5.967078 3.392658 -1.159515 -99.000000 -8.813380 0.000000 0.000048 0.001149 0.013240 0.000E+00 dH/H
rel_E_err 5.7977895671571087D-10
log_rel_run_E_err -8.0617472249783830
time_in_oscillation_code and total 3.3681744000000000D-02 3.3681744000000000D-02
save LOGS/profile1.data for model 5399
__________________________________________________________________________________________________________________________________________________
step lg_Tmax Teff lg_LH lg_Lnuc Mass H_rich H_cntr N_cntr Y_surf eta_cntr zones retry
lg_dt_yr lg_Tcntr lg_R lg_L3a lg_Lneu lg_Mdot He_core He_cntr O_cntr Z_surf gam_cntr iters
age_yr lg_Dcntr lg_L lg_LZ lg_Lphoto lg_Dsurf C_core C_cntr Ne_cntr Z_cntr v_div_cs dt_limit
__________________________________________________________________________________________________________________________________________________
5400 7.976457 2838.427 3.387563 3.393294 0.770221 0.300038 0.000000 0.006654 0.274608 21.294125 1911 0
2.895370 7.872618 2.312750 1.509835 2.222781 -6.976925 0.470183 0.986760 0.001985 0.013359 0.773365 8
1.2692E+10 5.966980 3.392504 -1.118943 -99.000000 -8.812733 0.000000 0.000048 0.001149 0.013240 0.000E+00 dH/H
rel_E_err 4.8237413168037073D-10
log_rel_run_E_err -8.0382377949747301
time_in_oscillation_code and total 3.3431141999999997D-02 6.7112885999999997D-02
save LOGS/profile2.data for model 5400
save photos/x400 for model 5400
retry: some xa < min_xa_hard_limit 5404
correction_norm > tol_correction_norm*coeff 5404 4.3623760698828614D-05 3.0000000000000004D-07 1.0000000000000000D-02
max_abs_correction > tol_max_correction*coeff 5404 9.9162462853773594D-04 3.0000000000000001D-05 1.0000000000000000D-02
residual_norm > tol_residual_norm*coeff 5404 3.8411971009334471D-07 1.0000000000000000D-10 1.0000000000000000D-02
max_residual > tol_max_residual*coeff 5404 1.9721214286905079D-03 1.0000000000000001D-07 1.0000000000000000D-02
retry: avg+max corr+resid -- give up in solver 5404
retry: some xa < min_xa_hard_limit 5407
5407 8.008975 2841.907 3.383759 3.418310 0.769937 0.299651 0.000000 0.006654 0.274608 21.296632 1955 3
2.072849 7.871923 2.310219 2.301236 2.219044 -6.981868 0.470286 0.986760 0.001985 0.013359 0.773961 5
1.2692E+10 5.965899 3.389571 -0.562002 -99.000000 -8.806661 0.000000 0.000048 0.001149 0.013240 0.000E+00 retry
rel_E_err 7.2220803081596820D-11
log_rel_run_E_err -7.9659510473464996
time_in_oscillation_code and total 3.4539677999999997D-02 1.0165256400000000D-01
save LOGS/profile3.data for model 5407
time_in_oscillation_code and total 3.5653549999999999D-02 1.3730611400000001D-01
save LOGS/profile4.data for model 5408
retry: some xa < min_xa_hard_limit 5410
5410 8.020523 2843.110 3.382791 3.460149 0.769905 0.299606 0.000000 0.006654 0.274608 21.297346 1973 4
1.851000 7.871524 2.309370 2.672229 2.218163 -6.983908 0.470300 0.986760 0.001985 0.013359 0.774289 5
1.2692E+10 5.965253 3.388607 -0.236227 -99.000000 -8.804573 0.000000 0.000048 0.001149 0.013240 0.000E+00 retry
rel_E_err 4.1792267377913336D-11
log_rel_run_E_err -7.9583783133098596
5420 8.051350 2845.074 3.378820 3.904946 0.769874 0.299561 0.000000 0.006654 0.274608 21.299914 2004 4
1.060546 7.869683 2.307970 3.750811 2.216149 -6.987133 0.470313 0.986760 0.001985 0.013359 0.775788 5
1.2692E+10 5.962252 3.387007 0.903648 -99.000000 -8.801152 0.000000 0.000048 0.001149 0.013240 0.000E+00 Lnuc_He
rel_E_err 1.0003972737816324D-11
log_rel_run_E_err -7.9497963069322486
5430 8.086933 2845.556 3.359695 4.940461 0.769869 0.299553 0.000000 0.006654 0.274608 21.303775 2029 4
0.185043 7.866138 2.307629 4.927887 2.246755 -6.988036 0.470316 0.986760 0.001985 0.013359 0.778664 5
1.2692E+10 5.956438 3.386621 2.298651 -99.000000 -8.800317 0.000000 0.000048 0.001149 0.013240 0.000E+00 Lnuc_He
rel_E_err -2.5731165006873385D-12
log_rel_run_E_err -7.9465224372793299
retry: some xa < min_xa_hard_limit 5438
5440 8.127470 2845.570 3.304937 6.140727 0.769869 0.299553 0.000000 0.006654 0.274608 21.308249 2044 5
-0.962956 7.859996 2.307621 6.138152 2.878432 -6.988062 0.470316 0.986760 0.001985 0.013359 0.783633 5
1.2692E+10 5.946300 3.386614 3.789329 -99.000000 -8.800293 0.000000 0.000048 0.001149 0.013240 0.000E+00 Lnuc_He
rel_E_err 7.0168568104319037D-12
log_rel_run_E_err -7.4989928392373404
5450 8.162695 2845.571 3.233160 7.078313 0.769869 0.299553 0.000000 0.006654 0.274608 21.312123 2042 5
-1.612974 7.852190 2.307621 7.074861 3.976379 -6.988061 0.470316 0.986760 0.001985 0.013359 0.789968 5
1.2692E+10 5.933372 3.386614 4.968948 -99.000000 -8.800292 0.000000 0.000048 0.001149 0.013240 0.000E+00 Lnuc_He
rel_E_err 6.9405974328972019D-09
log_rel_run_E_err -7.2247202882694799
time_in_oscillation_code and total 3.6028908999999998D-02 1.7333502300000000D-01
save LOGS/profile5.data for model 5450
save photos/x450 for model 5450
5460 8.208244 2845.574 3.090441 8.133551 0.769869 0.299553 0.000000 0.006654 0.274608 21.305627 2077 5
-2.388280 7.837498 2.307623 8.126636 5.333810 -6.988055 0.470316 0.986760 0.001985 0.013359 0.801780 6
1.2692E+10 5.908631 3.386620 6.331869 -99.000000 -8.800290 0.000000 0.000048 0.001149 0.013240 0.000E+00 Lnuc_He
rel_E_err -8.3481536991033631D-13
log_rel_run_E_err -7.2100615859016024
5470 8.257764 2845.582 2.803576 9.059055 0.769869 0.299553 0.000000 0.006654 0.274608 21.251561 2123 5
-3.076266 7.811680 2.307629 9.043819 6.596779 -6.988036 0.470316 0.986760 0.001985 0.013359 0.822171 6
1.2692E+10 5.863899 3.386635 7.596557 -99.000000 -8.800286 0.000000 0.000048 0.001149 0.013240 0.000E+00 Lnuc_He
rel_E_err -9.7760767079314158D-12
log_rel_run_E_err -7.1878621529795916
5480 8.298986 3134.810 2.248746 9.563214 0.769869 0.299553 0.000000 0.006654 0.274608 21.066588 2169 5
-3.403045 7.771156 2.294283 9.531103 7.411917 -7.008922 0.470316 0.986760 0.001985 0.013359 0.853323 8
1.2692E+10 5.790780 3.528104 8.416126 -99.000000 -8.557797 0.000000 0.000048 0.001149 0.013240 0.000E+00 Lnuc_He
rel_E_err -2.2182503987318975D-11
log_rel_run_E_err -7.1500966346930088
retry: adjust_correction failed in eval_equations -- give up in solver 5482
correction_norm > tol_correction_norm*coeff 5482 1.3362116697402277D-04 3.0000000000000001D-05 1.0000000000000000D+00
max_abs_correction > tol_max_correction*coeff 5482 9.3310264248890103D-03 3.0000000000000001D-03 1.0000000000000000D+00
residual_norm > tol_residual_norm*coeff 5482 9.9479773990558933D-07 1.0000000000000000D-08 1.0000000000000000D+00
max_residual > tol_max_residual*coeff 5482 9.9714114085482919D-03 1.0000000000000001D-05 1.0000000000000000D+00
retry: avg+max corr+resid -- give up in solver 5482
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
AESOPUS_interp: logR, logT outside of table bounds
retry: get_T_tau -- L <= 0 5483
correction_norm > tol_correction_norm*coeff 5483 1.3831579513077668D-03 3.0000000000000001D-05 1.0000000000000000D+00
max_abs_correction > tol_max_correction*coeff 5483 9.7665522700013427D-02 3.0000000000000001D-03 1.0000000000000000D+00
residual_norm > tol_residual_norm*coeff 5483 3.3947649638033567D-04 1.0000000000000000D-08 1.0000000000000000D+00
max_residual > tol_max_residual*coeff 5483 1.2060561367438729D+00 1.0000000000000001D-05 1.0000000000000000D+00
retry: avg+max corr+resid -- give up in solver 5483
max_residual > tol_max_residual*coeff 5483 2.9124211182551818D-04 1.0000000000000001D-05 1.0000000000000000D+00
retry: max resid -- give up in solver 5483
residual_norm > tol_residual_norm*coeff 5484 2.2695977638159261D-08 1.0000000000000000D-08 1.0000000000000000D+00
max_residual > tol_max_residual*coeff 5484 6.3155097696565576D-04 1.0000000000000001D-05 1.0000000000000000D+00
retry: avg+max resid -- give up in solver 5484
correction_norm > tol_correction_norm*coeff 5484 1.1783796264113738D-06 3.0000000000000004D-07 1.0000000000000000D-02
max_abs_correction > tol_max_correction*coeff 5484 8.5560922867624064D-05 3.0000000000000001D-05 1.0000000000000000D-02
residual_norm > tol_residual_norm*coeff 5484 1.2489872051461098D-09 1.0000000000000000D-10 1.0000000000000000D-02
max_residual > tol_max_residual*coeff 5484 5.9947131981985535D-05 1.0000000000000001D-07 1.0000000000000000D-02
retry: avg+max corr+resid -- give up in solver 5484
residual_norm > tol_residual_norm*coeff 5486 2.1606643034105666D-09 1.0000000000000001D-09 1.0000000000000001D-01
max_residual > tol_max_residual*coeff 5486 9.0748038077848037D-05 1.0000000000000002D-06 1.0000000000000001D-01
retry: avg+max resid -- give up in solver 5486
correction_norm > tol_correction_norm*coeff 5486 6.7063161107737327D-07 3.0000000000000004D-07 1.0000000000000000D-02
max_abs_correction > tol_max_correction*coeff 5486 4.8549014837285764D-05 3.0000000000000001D-05 1.0000000000000000D-02
residual_norm > tol_residual_norm*coeff 5486 8.0046690576254016D-10 1.0000000000000000D-10 1.0000000000000000D-02
max_residual > tol_max_residual*coeff 5486 3.6086651047817497D-05 1.0000000000000001D-07 1.0000000000000000D-02
retry: avg+max corr+resid -- give up in solver 5486
correction_norm > tol_correction_norm*coeff 5488 7.6328901268138305D-05 3.0000000000000001D-06 1.0000000000000001D-01
max_abs_correction > tol_max_correction*coeff 5488 5.4416209314463206D-03 3.0000000000000003D-04 1.0000000000000001D-01
residual_norm > tol_residual_norm*coeff 5488 1.9721955915124521D-08 1.0000000000000001D-09 1.0000000000000001D-01
max_residual > tol_max_residual*coeff 5488 7.6968831160062479D-04 1.0000000000000002D-06 1.0000000000000001D-01
retry: avg+max corr+resid -- give up in solver 5488
residual_norm > tol_residual_norm*coeff 5489 1.1919826478280686D-09 1.0000000000000001D-09 1.0000000000000001D-01
max_residual > tol_max_residual*coeff 5489 5.4975456330375039D-05 1.0000000000000002D-06 1.0000000000000001D-01
retry: avg+max resid -- give up in solver 5489
correction_norm > tol_correction_norm*coeff 5489 3.0997915303175289D-06 3.0000000000000001D-06 1.0000000000000001D-01
residual_norm > tol_residual_norm*coeff 5489 1.9804181565549182D-09 1.0000000000000001D-09 1.0000000000000001D-01
max_residual > tol_max_residual*coeff 5489 9.5162034392970440D-05 1.0000000000000002D-06 1.0000000000000001D-01
retry: avg corr, avg+max resid -- give up in solver 5489
max_residual > tol_max_residual*coeff 5489 2.7350659960845292D-05 1.0000000000000002D-06 1.0000000000000001D-01
retry: max resid -- give up in solver 5489
max_residual > tol_max_residual*coeff 5490 3.0376929410991238D-05 4.3183984527231594D-06 4.3183984527231589D-01
retry: max resid -- give up in solver 5490
5490 8.302142 2709.632 2.171271 9.581968 0.769869 0.299553 0.000000 0.006654 0.274608 21.048060 2160 19
-7.652923 7.766051 2.315676 9.547612 7.458216 -7.049488 0.470316 0.986760 0.001985 0.013359 0.857461 8
1.2692E+10 5.781769 3.317688 8.463110 -99.000000 -8.972460 0.000000 0.000048 0.001149 0.013240 0.000E+00 retry
rel_E_err 3.7694923984974081D-13
log_rel_run_E_err -7.1452607561514121
max_residual > tol_max_residual*coeff 5491 4.2417566794128447D-05 1.0000000000000002D-06 1.0000000000000001D-01
retry: max resid -- give up in solver 5491
max_residual > tol_max_residual*coeff 5491 3.5355960334548797D-05 1.0000000000000002D-06 1.0000000000000001D-01
retry: max resid -- give up in solver 5491
max_residual > tol_max_residual*coeff 5491 3.3282861256367369D-05 1.0000000000000002D-06 1.0000000000000001D-01
retry: max resid -- give up in solver 5491
max_residual > tol_max_residual*coeff 5492 4.3734177451959511D-05 1.0000000000000002D-06 1.0000000000000001D-01
retry: max resid -- give up in solver 5492
max_residual > tol_max_residual*coeff 5492 4.3353651551392147D-05 1.0000000000000002D-06 1.0000000000000001D-01
retry: max resid -- give up in solver 5492
residual_norm > tol_residual_norm*coeff 5492 1.2002374871314211D-09 1.0000000000000001D-09 1.0000000000000001D-01
max_residual > tol_max_residual*coeff 5492 5.6598734313517696D-05 1.0000000000000002D-06 1.0000000000000001D-01
retry: avg+max resid -- give up in solver 5492
correction_norm > tol_correction_norm*coeff 5493 3.1204510306963463D-06 3.0000000000000001D-06 1.0000000000000001D-01
residual_norm > tol_residual_norm*coeff 5493 2.3776171785653902D-09 1.0000000000000001D-09 1.0000000000000001D-01
max_residual > tol_max_residual*coeff 5493 1.1377641317800169D-04 1.0000000000000002D-06 1.0000000000000001D-01
retry: avg corr, avg+max resid -- give up in solver 5493
max_residual > tol_max_residual*coeff 5493 1.9116969975918607D-05 4.2539882362861366D-06 4.2539882362861364D-01
retry: max resid -- give up in solver 5493
max_residual > tol_max_residual*coeff 5493 1.5554131995720866D-05 9.6543155778222150D-07 9.6543155778222145D-02
retry: max resid -- give up in solver 5493
residual_norm > tol_residual_norm*coeff 5493 4.9770169585344934D-10 1.0000000000000000D-10 1.0000000000000000D-02
max_residual > tol_max_residual*coeff 5493 2.2707980907199499D-05 1.0000000000000001D-07 1.0000000000000000D-02
retry: avg+max resid -- give up in solver 5493
correction_norm > tol_correction_norm*coeff 5493 5.6248869191990835D-07 3.0000000000000004D-07 1.0000000000000000D-02
max_abs_correction > tol_max_correction*coeff 5493 3.9627961415967690D-05 3.0000000000000001D-05 1.0000000000000000D-02
residual_norm > tol_residual_norm*coeff 5493 5.9924447051292586D-10 1.0000000000000000D-10 1.0000000000000000D-02
max_residual > tol_max_residual*coeff 5493 2.7668016532089723D-05 1.0000000000000001D-07 1.0000000000000000D-02
retry: avg+max corr+resid -- give up in solver 5493
correction_norm > tol_correction_norm*coeff 5493 4.6697775388045699D-07 3.0000000000000004D-07 1.0000000000000000D-02
max_abs_correction > tol_max_correction*coeff 5493 3.2903649950295061D-05 3.0000000000000001D-05 1.0000000000000000D-02
residual_norm > tol_residual_norm*coeff 5493 5.8950245299098804D-10 1.0000000000000000D-10 1.0000000000000000D-02
max_residual > tol_max_residual*coeff 5493 2.7158341643129802D-05 1.0000000000000001D-07 1.0000000000000000D-02
retry: avg+max corr+resid -- give up in solver 5493
max_residual > tol_max_residual*coeff 5493 2.7018754647833987D-05 8.1327761656072186D-07 8.1327761656072181D-02
retry: max resid -- give up in solver 5493
correction_norm > tol_correction_norm*coeff 5493 6.1912787632503697D-07 3.0000000000000004D-07 1.0000000000000000D-02
max_abs_correction > tol_max_correction*coeff 5493 4.3612094440387417D-05 3.0000000000000001D-05 1.0000000000000000D-02
residual_norm > tol_residual_norm*coeff 5493 6.4070563121962976D-10 1.0000000000000000D-10 1.0000000000000000D-02
max_residual > tol_max_residual*coeff 5493 2.9893030726313441D-05 1.0000000000000001D-07 1.0000000000000000D-02
retry: avg+max corr+resid -- give up in solver 5493
correction_norm > tol_correction_norm*coeff 5493 9.4035871322936061D-07 5.4456226354391813D-07 1.8152075451463936D-02
max_abs_correction > tol_max_correction*coeff 5493 6.6144006245853941D-05 5.4456226354391809D-05 1.8152075451463936D-02
residual_norm > tol_residual_norm*coeff 5493 9.1135064050353081D-10 1.8152075451463935D-10 1.8152075451463936D-02
max_residual > tol_max_residual*coeff 5493 4.2501058226562805D-05 1.8152075451463936D-07 1.8152075451463936D-02
retry: avg+max corr+resid -- give up in solver 5493
max_residual > tol_max_residual*coeff 5493 2.5096750931652423D-05 3.9401542403792785D-06 3.9401542403792783D-01
retry: max resid -- give up in solver 5493
max_residual > tol_max_residual*coeff 5493 3.5387751162296865D-05 2.0000000000000003D-06 2.0000000000000001D-01
retry: max resid -- give up in solver 5493
correction_norm > tol_correction_norm*coeff 5493 5.8823648535986421D-07 5.2683976751124031D-07 1.7561325583708012D-02
residual_norm > tol_residual_norm*coeff 5493 5.8796184510871723D-10 1.7561325583708013D-10 1.7561325583708012D-02
max_residual > tol_max_residual*coeff 5493 2.7135615064288743D-05 1.7561325583708012D-07 1.7561325583708012D-02
retry: avg corr, avg+max resid -- give up in solver 5493
correction_norm > tol_correction_norm*coeff 5494 2.4613432650099367D-05 1.3113380081944059D-05 4.3711266939813526D-01
max_abs_correction > tol_max_correction*coeff 5494 5.4776452848070544D-03 1.3113380081944059D-03 4.3711266939813526D-01
residual_norm > tol_residual_norm*coeff 5494 1.0776637333564953D-06 4.3711266939813526D-09 4.3711266939813526D-01
max_residual > tol_max_residual*coeff 5494 4.2073063953340650D-03 4.3711266939813526D-06 4.3711266939813526D-01
retry: avg+max corr+resid -- give up in solver 5494
correction_norm > tol_correction_norm*coeff 5494 8.6647233797873017D-07 3.0000000000000004D-07 1.0000000000000000D-02
max_abs_correction > tol_max_correction*coeff 5494 6.2254230763755561D-05 3.0000000000000001D-05 1.0000000000000000D-02
residual_norm > tol_residual_norm*coeff 5494 9.0076806917285472D-10 1.0000000000000000D-10 1.0000000000000000D-02
max_residual > tol_max_residual*coeff 5494 4.1808920271613204D-05 1.0000000000000001D-07 1.0000000000000000D-02
retry: avg+max corr+resid -- give up in solver 5494
dt 6.6923891897861795D-07
min_timestep_limit 9.9999999999999995D-07
stopping because of problems dt < min_timestep_limit
time_in_oscillation_code and total 3.8033053999999997D-02 2.1136807699999999D-01
terminated evolution: cannot find acceptable model
termination code: min_timestep_limit
DATE: 2021-02-04
TIME: 16:03:04
<inlist_to_end_core_h_burn><rn><inlist_to_wd><inlist_start><inlist_to_end_core_he_burn_header><inlist_to_start_he_core_flash><inlist_to_wd_header><inlist_to_end_core_h_burn_header><inlist_to_end_agb><inlist_to_end_agb_header><inlist_start_header><inlist_to_end_core_he_burn><inlist_to_start_he_core_flash_header>_______________________________________________
mesa-users at lists.mesastar.org
https://lists.mesastar.org/mailman/listinfo/mesa-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210205/5be0c2bb/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: start_he_core_flash.mod
Type: audio/x-mod
Size: 840687 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210205/5be0c2bb/attachment-0001.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_to_end_core_he_burn_header
Type: application/octet-stream
Size: 737 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210205/5be0c2bb/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_to_end_core_he_burn
Type: application/octet-stream
Size: 11716 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210205/5be0c2bb/attachment-0003.obj>
More information about the Mesa-users
mailing list