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