[Mesa-users] Timestep issue at the He flash stage
Bill Paxton
paxton at kitp.ucsb.edu
Fri Feb 5 16:20:17 UTC 2021
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/ <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/ecd65e61/attachment.htm>
More information about the Mesa-users
mailing list