[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