[Mesa-users] How to Solve "hdyro_newton_step returned ierr > Retry cnt, step, log 10(dt/yr), retry_factor x xxxx" MESA r7624 ?

Rob Farmer robert.j.farmer37 at gmail.com
Mon Sep 26 12:06:54 UTC 2022


The error is telling you that the MESA solver is running into difficulties
and some part of the code is returning an error.

You may get more information by setting report_ierr = .true.

However, retries are a normal part of a MESA run, so as long as your model
is still evolving on a reasonable time step for your problem then they
don't necessarily need to be fixed. Only if they stop your model from
progressing forward do they need attention.

Rob

On Wed, 21 Sept 2022 at 16:39, Agung Prawira Negara via Mesa-users <
mesa-users at lists.mesastar.org> wrote:

> Dear MESA Users,
>
> Good Evening from Indonesia,
> Right now, I am trying to reproduce some results from paper "A Systematic
> Survey of the Effects of Wind Mass Loss Algorithms on the Evolution of
> Single Massive Stars" (https://arxiv.org/abs/1703.09705) and fortunately
> the author of the paper has shared the *mesa inlist* and customized
> *run_start_extras.f* that he used in that paper that we can access here (
> https://stellarcollapse.org/renzo2017.html). Anyway we use the *r7624*
> (running via NuDocker).
>
> Long story short, this time what I have done and what I am trying to do is
> I want to do some modifications from the paper above that is to enhance the
> stellar winds mass loss coefficient (\eta ) let say \eta = 5 and \eta=10.
> Actually this idea to enhance is from this paper (
> https://www.aanda.org/articles/aa/full_html/2015/03/aa24671-14/aa24671-14.html
> and several other papers).
>
> *This simulation works well *in the 12Msun with \eta=5 and \eta=10, but
> the *problem *happens when I do the same thing with the higher masses,
> i.e. 20Msun with \eta=10 and 25Msun with \eta=10. Here I attach the error
> message from the terminal that I got when running the 25Msun \eta=10 (the
> error code happens the same with the 20Msun \eta=10) . it says
> "hdyro_newton_step returned ierr
>        Retry cnt, step, log 10(dt/yr), retry_factor x xxxx -x.xxxxxx
> x.xxxxxx".
> (idk if this is a error code or not to be honest, but since this line code
> happen some properties does not change along with the increasing time step,
> namely central temperature and mass of the star)
>
> I don't know if this additional information will be useful, but I do
> little modification on the code, just to reduce the computational
> resolution, viz.
> *inlist_massive_defaults*
> mesh_delta_coeff = 2.5 (previously 1.5)
> mesh_delta_coeff_for_highT = 3.5 (previously 2.5)
>
> *inlist project*
> mesh_delta_coeff_for_highT = 2.5 (previously 1.0)
> mesh_delta_coeff = 2.5 (previously 1.0)
> max_allowed_nz=100000000 (previously 300000000)
>
> the model evolve until Tc = 10^9
>
> So here my question is, *what is the meaning of this "error code" *? How
> could I solve this problem? Thank you in advance, any assistance will be
> very helpful for me.
>
> here is the fully error code while mesa was running
>
> __________________________________________________________________________________________________________________________________________________
>
>       14930   8.647701  1.357E+05   1.156180   4.034533   1.090658
> 0.000000   0.000000   0.000000   0.341239   0.000000   5.651026  32347
>  0
>    0.678192   6.341081  -0.658587   3.607572   4.051800  -5.510072
> 1.090658   0.000000   0.516614   0.658761   0.022589   1.341948      3
>  0
>  1.0188E+07  23.008543   4.166206   3.830088   5.931632  -8.122070
> 0.993534   0.458404   0.001656   0.000519  9.774E-01  0.248E-05
> dX_nuc_drop
>
> save LOGS/profile36.data for model 14930
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14940
>  0.074935        0.500000
>  using: N&L, eta=   10.000000000000000
>
>
> __________________________________________________________________________________________________________________________________________________
>
>        step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass
> H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones retry
>    lg_dt_yr    lg_Dcntr    lg_R       lg_L3a    lg_Lneu     lg_Mdot
>  He_core    He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
>      age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf
> C_core     C_cntr     Ne_cntr    Si_cntr    Z_avg     v_div_cs     dt_limit
>
> __________________________________________________________________________________________________________________________________________________
>
>       14940   8.647589  1.358E+05   1.161274   4.034659   1.090518
> 0.000000   0.000000   0.000000   0.341186   0.000000   5.664049  32347
>  1
>    0.375965   6.342524  -0.659394   3.607692   4.053447  -5.509744
> 1.090518   0.000000   0.516614   0.658814   0.022521   1.343781      3
>  0
>  1.0188E+07  23.010534   4.166561   3.830208   5.933603  -8.122451
> 0.993652   0.458404   0.001657   0.000519  9.775E-01  0.447E-05
> dX_nuc_drop
>
> __________________________________________________________________________________________________________________________________________________
>
>       14930   8.647701  1.357E+05   1.156180   4.034533   1.090658
> 0.000000   0.000000   0.000000   0.341239   0.000000   5.651026  32347
>  0
>    0.678192   6.341081  -0.658587   3.607572   4.051800  -5.510072
> 1.090658   0.000000   0.516614   0.658761   0.022589   1.341948      3
>  0
>  1.0188E+07  23.008543   4.166206   3.830088   5.931632  -8.122070
> 0.993534   0.458404   0.001656   0.000519  9.774E-01  0.248E-05
> dX_nuc_drop
>
> save LOGS/profile36.data for model 14930
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14940
>  0.074935        0.500000
>  using: N&L, eta=   10.000000000000000
>
>
> __________________________________________________________________________________________________________________________________________________
>
>        step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass
> H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones retry
>    lg_dt_yr    lg_Dcntr    lg_R       lg_L3a    lg_Lneu     lg_Mdot
>  He_core    He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
>      age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf
> C_core     C_cntr     Ne_cntr    Si_cntr    Z_avg     v_div_cs     dt_limit
>
> __________________________________________________________________________________________________________________________________________________
>
>       14940   8.647589  1.358E+05   1.161274   4.034659   1.090518
> 0.000000   0.000000   0.000000   0.341186   0.000000   5.664049  32347
>  1
>    0.375965   6.342524  -0.659394   3.607692   4.053447  -5.509744
> 1.090518   0.000000   0.516614   0.658814   0.022521   1.343781      3
>  0
>  1.0188E+07  23.010534   4.166561   3.830208   5.933603  -8.122451
> 0.993652   0.458404   0.001657   0.000519  9.775E-01  0.447E-05
> dX_nuc_drop
> save LOGS/profile40.data for model 14940
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14941
> -0.146914        0.500000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     2 14941
> -0.447944        0.500000
>  using: N&L, eta=   10.000000000000000
>  have reached retry limit so now backup
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14940
> -0.748974        0.500000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     2 14940
> -1.050004        0.500000
>  using: N&L, eta=   10.000000000000000
>  have reached retry limit so now backup
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14940
> -1.953094        0.500000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     2 14940
> -2.254124        0.500000
>  using: N&L, eta=   10.000000000000000
>  have reached retry limit so now backup
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14940
> -3.759274        0.500000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     2 14940
> -4.060304        0.500000
>  using: N&L, eta=   10.000000000000000
>  have reached retry limit so now backup
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14940
> -6.167514        0.500000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     2 14940
> -6.468544        0.500000
>  using: N&L, eta=   10.000000000000000
>  have reached retry limit so now backup
>  using: N&L, eta=   10.000000000000000
>
>       14940   8.647595  1.354E+05   1.161006   4.034652   1.090525
> 0.000000   0.000000   0.000000   0.341189   0.000000   5.663363  32347
> 11
>   -8.575754   6.342448  -0.659266   3.607685   4.053361  -5.509744
> 1.090525   0.000000   0.516614   0.658811   0.022525   1.343684      3
>  5
>  1.0188E+07  23.010430   4.160783   3.830201   5.928113  -8.103213
> 0.993652   0.458404   0.001657   0.000519  9.775E-01  0.000E+00  max
> increase
>
> save LOGS/profile40.data for model 14940
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>
>
> __________________________________________________________________________________________________________________________________________________
>
>        step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass
> H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones retry
>    lg_dt_yr    lg_Dcntr    lg_R       lg_L3a    lg_Lneu     lg_Mdot
>  He_core    He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
>      age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf
> C_core     C_cntr     Ne_cntr    Si_cntr    Z_avg     v_div_cs     dt_limit
>
> __________________________________________________________________________________________________________________________________________________
>
>       14950   8.647595  1.355E+05   1.161006   4.034652   1.090525
> 0.000000   0.000000   0.000000   0.341189   0.000000   5.663363  32347
> 11
>   -8.021485   6.342448  -0.659394   3.607685   4.053361  -5.517288
> 1.090525   0.000000   0.516614   0.658811   0.022525   1.343684      3
>  5
>  1.0188E+07  23.010430   4.162168   3.830201   5.929559  -8.108337
> 0.993652   0.458404   0.001657   0.000519  9.775E-01  0.000E+00  max
> increase
>
> save LOGS/profile42.data for model 14950
> save photos/x950 for model 14950
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14956
> -8.148458        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>
>
> __________________________________________________________________________________________________________________________________________________
>
>        step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass
> H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones retry
>    lg_dt_yr    lg_Dcntr    lg_R       lg_L3a    lg_Lneu     lg_Mdot
>  He_core    He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
>      age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf
> C_core     C_cntr     Ne_cntr    Si_cntr    Z_avg     v_div_cs     dt_limit
>
> __________________________________________________________________________________________________________________________________________________
>     14960   8.647595  1.358E+05   1.161006   4.034652   1.090525
> 0.000000   0.000000   0.000000   0.341189   0.000000   5.663363  32347
> 12
>   -7.530703   6.342448  -0.659173   3.607685   4.053361  -5.511110
> 1.090525   0.000000   0.516614   0.658811   0.022525   1.343684      3
>  5
>  1.0188E+07  23.010430   4.166278   3.830201   5.932877  -8.121889
> 0.993652   0.458404   0.001657   0.000519  9.775E-01  0.000E+00  max
> increase
>
> save LOGS/profile43.data for model 14960
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14963
> -7.895219        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14964
> -8.117068        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14970
> -7.943010        0.500000
>  using: N&L, eta=   10.000000000000000
>
>
> __________________________________________________________________________________________________________________________________________________
>
>        step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass
> H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones retry
>    lg_dt_yr    lg_Dcntr    lg_R       lg_L3a    lg_Lneu     lg_Mdot
>  He_core    He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
>      age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf
> C_core     C_cntr     Ne_cntr    Si_cntr    Z_avg     v_div_cs     dt_limit
>
> __________________________________________________________________________________________________________________________________________________
>
>       14970   8.647595  1.357E+05   1.161006   4.034652   1.090525
> 0.000000   0.000000   0.000000   0.341189   0.000000   5.663363  32347
> 15
>   -7.641980   6.342448  -0.659153   3.607685   4.053361  -5.510231
> 1.090525   0.000000   0.516614   0.658811   0.022525   1.343684      3
>  5
>  1.0188E+07  23.010430   4.166108   3.830201   5.932667  -8.121468
> 0.993652   0.458404   0.001657   0.000519  9.775E-01  0.000E+00  max
> increase
>
> save LOGS/profile46.data for model 14970
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14973
> -8.006496        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14977
> -7.990801        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>
>
> __________________________________________________________________________________________________________________________________________________
>        step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass
> H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones retry
>    lg_dt_yr    lg_Dcntr    lg_R       lg_L3a    lg_Lneu     lg_Mdot
>  He_core    He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
>      age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf
> C_core     C_cntr     Ne_cntr    Si_cntr    Z_avg     v_div_cs     dt_limit
>
> __________________________________________________________________________________________________________________________________________________
>
>       14980   8.647595  1.358E+05   1.161006   4.034653   1.090525
> 0.000000   0.000000   0.000000   0.341189   0.000000   5.663363  32347
> 17
>   -7.452228   6.342448  -0.659147   3.607685   4.053361  -5.510620
> 1.090525   0.000000   0.516614   0.658811   0.022525   1.343684      3
>  5
>  1.0188E+07  23.010430   4.166492   3.830201   5.933038  -8.122370
> 0.993652   0.458404   0.001657   0.000519  9.775E-01  0.000E+00  max
> increase
>
> save LOGS/profile48.data for model 14980
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14983
> -7.816744        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14984
> -8.038593        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14989
> -7.943716        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>
>
> __________________________________________________________________________________________________________________________________________________
>
>        step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass
> H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones retry
>    lg_dt_yr    lg_Dcntr    lg_R       lg_L3a    lg_Lneu     lg_Mdot
>  He_core    He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
>      age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf
> C_core     C_cntr     Ne_cntr    Si_cntr    Z_avg     v_div_cs     dt_limit
>
> __________________________________________________________________________________________________________________________________________________
>
>       14990   8.647595  1.358E+05   1.161006   4.034652   1.090525
> 0.000000   0.000000   0.000000   0.341189   0.000000   5.663363  32347
> 20
>   -7.563505   6.342448  -0.659133   3.607685   4.053361  -5.510752
> 1.090525   0.000000   0.516614   0.658811   0.022525   1.343684      3
>  5
>  1.0188E+07  23.010430   4.166395   3.830201   5.932913  -8.122115
> 0.993652   0.458404   0.001657   0.000519  9.775E-01  0.000E+00  max
> increase
>
> save LOGS/profile50.data for model 14990
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14993
> -7.928021        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 14994
> -8.149870        0.500000
> using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 15000
> -7.975813        0.500000
>  using: N&L, eta=   10.000000000000000
>
>
> __________________________________________________________________________________________________________________________________________________
>
>        step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass
> H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones retry
>    lg_dt_yr    lg_Dcntr    lg_R       lg_L3a    lg_Lneu     lg_Mdot
>  He_core    He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
>      age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf
> C_core     C_cntr     Ne_cntr    Si_cntr    Z_avg     v_div_cs     dt_limit
>
> __________________________________________________________________________________________________________________________________________________
>
>       15000   8.647595  1.357E+05   1.161006   4.034652   1.090525
> 0.000000   0.000000   0.000000   0.341189   0.000000   5.663363  32347
> 23
>   -7.674783   6.342448  -0.659127   3.607685   4.053361  -5.510138
> 1.090525   0.000000   0.516614   0.658811   0.022525   1.343684      3
>  5
>  1.0188E+07  23.010430   4.166176   3.830201   5.932683  -8.121570
> 0.993652   0.458404   0.001657   0.000519  9.775E-01  0.000E+00  max
> increase
>
> save LOGS/profile2.data for model 15000
> save photos/15000 for model 15000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 15001
> -8.197661        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 15007
> -8.023604        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>
> __________________________________________________________________________________________________________________________________________________
>
>        step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass
> H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones retry
>    lg_dt_yr    lg_Dcntr    lg_R       lg_L3a    lg_Lneu     lg_Mdot
>  He_core    He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
>      age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf
> C_core     C_cntr     Ne_cntr    Si_cntr    Z_avg     v_div_cs     dt_limit
>
> __________________________________________________________________________________________________________________________________________________
>
>       15010   8.647595  1.357E+05   1.161006   4.034652   1.090525
> 0.000000   0.000000   0.000000   0.341189   0.000000   5.663363  32347
> 25
>   -7.485030   6.342448  -0.659111   3.607685   4.053361  -5.512164
> 1.090525   0.000000   0.516614   0.658811   0.022525   1.343684      3
>  5
>  1.0188E+07  23.010430   4.166093   3.830201   5.932568  -8.121363
> 0.993652   0.458404   0.001657   0.000519  9.775E-01  0.000E+00  max
> increase
>
> save LOGS/profile4.data for model 15010
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 15011
> -8.007909        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 15014
> -8.071395        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 15018
> -8.055700        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> ydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 15020
> -8.198368        0.500000
>  using: N&L, eta=   10.000000000000000
>
>
> __________________________________________________________________________________________________________________________________________________
>
>        step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass
> H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones retry
>    lg_dt_yr    lg_Dcntr    lg_R       lg_L3a    lg_Lneu     lg_Mdot
>  He_core    He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
>      age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf
> C_core     C_cntr     Ne_cntr    Si_cntr    Z_avg     v_div_cs     dt_limit
>
> __________________________________________________________________________________________________________________________________________________
>
>       15020   8.647595  1.356E+05   1.161006   4.034652   1.090525
> 0.000000   0.000000   0.000000   0.341189   0.000000   5.663363  32347
> 29
>   -7.897338   6.342448  -0.659098   3.607685   4.053361  -5.512807
> 1.090525   0.000000   0.516614   0.658811   0.022525   1.343684      3
>  5
>  1.0188E+07  23.010430   4.164390   3.830201   5.930898  -8.116656
> 0.993652   0.458404   0.001657   0.000519  9.775E-01  0.000E+00  max
> increase
>
> save LOGS/profile6.data for model 15020
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 15025
> -8.103491        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>
>
> __________________________________________________________________________________________________________________________________________________
>        step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass
> H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones retry
>    lg_dt_yr    lg_Dcntr    lg_R       lg_L3a    lg_Lneu     lg_Mdot
>  He_core    He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
>      age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf
> C_core     C_cntr     Ne_cntr    Si_cntr    Z_avg     v_div_cs     dt_limit
>
> __________________________________________________________________________________________________________________________________________________
>
>       15030   8.647595  1.357E+05   1.161006   4.034652   1.090525
> 0.000000   0.000000   0.000000   0.341189   0.000000   5.663363  32347
> 30
>   -7.406555   6.342448  -0.659051   3.607685   4.053361  -5.511683
> 1.090525   0.000000   0.516614   0.658811   0.022525   1.343684      3
>  5
>  1.0188E+07  23.010430   4.166313   3.830201   5.932667  -8.121888
> 0.993652   0.458404   0.001657   0.000519  9.775E-01  0.000E+00  max
> increase
>
> save LOGS/profile8.data for model 15030
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 15031
> -7.929434        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 15032
> -8.151283        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 15038
> -7.977225        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 15040
> -8.119893        0.500000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     2 15040
> -8.420923        0.500000
>  using: N&L, eta=   10.000000000000000
>  have reached retry limit so now backup
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>
>
> __________________________________________________________________________________________________________________________________________________
>
>        step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass
> H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones retry
>    lg_dt_yr    lg_Dcntr    lg_R       lg_L3a    lg_Lneu     lg_Mdot
>  He_core    He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
>      age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf
> C_core     C_cntr     Ne_cntr    Si_cntr    Z_avg     v_div_cs     dt_limit
>
> __________________________________________________________________________________________________________________________________________________
>
>       15040   8.647595  1.354E+05   1.161006   4.034652   1.090525
> 0.000000   0.000000   0.000000   0.341189   0.000000   5.663363  32347
> 35
>   -8.119893   6.342448  -0.659237   3.607685   4.053361  -5.510984
> 1.090525   0.000000   0.516614   0.658811   0.022525   1.343684      3
>  6
>  1.0188E+07  23.010430   4.161498   3.830201   5.928680  -8.105585
> 0.993652   0.458404   0.001657   0.000519  9.775E-01  0.000E+00
>  hold
>
> save LOGS/profile12.data for model 15040
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>
>
> __________________________________________________________________________________________________________________________________________________
>
>        step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass
> H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones retry
>    lg_dt_yr    lg_Dcntr    lg_R       lg_L3a    lg_Lneu     lg_Mdot
>  He_core    He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
>      age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf
> C_core     C_cntr     Ne_cntr    Si_cntr    Z_avg     v_div_cs     dt_limit
>
> __________________________________________________________________________________________________________________________________________________
>
>       15050   8.647595  1.358E+05   1.161006   4.034652   1.090525
> 0.000000   0.000000   0.000000   0.341189   0.000000   5.663363  32347
> 35
>   -7.565624   6.342448  -0.659200   3.607685   4.053361  -5.511205
> 1.090525   0.000000   0.516614   0.658811   0.022525   1.343684      3
>  6
>  1.0188E+07  23.010430   4.166352   3.830201   5.933005  -8.121837
> 0.993652   0.458404   0.001657   0.000519  9.775E-01  0.000E+00  max
> increase
>
> save LOGS/profile13.data for model 15050
> save photos/x050 for model 15050
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 15053
> -7.930140        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 15056
> -7.993627        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
> hydro_newton_step returned ierr
>        retry cnt, step, log10(dt/yr), retry_factor     1 15058
> -8.136294        0.500000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>  using: N&L, eta=   10.000000000000000
>
>
> __________________________________________________________________________________________________________________________________________________
>
>
>
>
>
> Regards,
> Agung Prawira
> _______________________________________________
> 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/20220926/e0c58474/attachment.htm>


More information about the Mesa-users mailing list