[Mesa-users] Neutron star wind model: Tau Factor
Pero, Jason
jsp5218 at mavs.uta.edu
Tue Oct 25 00:59:09 UTC 2022
Hello MESA users,
MESA version r22.05.1
I'm having trouble simulating the hydrodynamic phase of a neutron star x-ray burst. During the hydrostatic rise the tau_factor is set to a value similar to what is in the ns_he test suite:
set_initial_tau_factor = .true.
set_to_this_tau_factor = 100
During the wind phase I switch from the hydrostatic to the hydrodynamic solver at which time a change the tau_factor to capture regions closer to the photosphere:
relax_initial_tau_factor = .true.
relax_tau_factor = .true.
relax_to_this_tau_factor = 1
Running the hydrodynamic phase then produces the following error:
stopping because of problems dt < min_timestep_limit
failed in relax tau factor
star_relax_tau_factor ierr -1
I've tried different tau factors ranging from 1-99 and still no luck. Based on the debug information it looks like there is an issue with the density becoming unrealistic. I then tried adjusting the control options "min_timestep_limit" and "hydro_mtx_max_allowed_logRho", but still can't get the model to relax to the new tau factor.
Any help/advice would be greatly appreciated.
Here is the full debug output from the terminal:
jpero at Jasons-MacBook-Pro burst22 % ./rn
DATE: 2022-10-24
TIME: 19:11:33
reading user weak rate file /Users/jpero/Apps/mesa-r22.05.1/data/rates_data/rate_tables/S13_r_be7_wk_li7.h5
version_number r22.05.1
read inlist_hydrodynamic
load saved model ns_he_hstatic.mod
WARNING: changing to saved tau_factor = 1.0000000000000000D+02
set_initial_dt (seconds) 1.0000000000000000D-08
set_initial_model_number 1810
set_initial_number_retries 0
relax_tau_factor 1.0000000000000000D+00
current tau_factor 1.0000000000000000D+02
relax to new tau_factor 1.0000000000000000D+00
dlogtau_factor 1.0000000000000001D-01
solver_call_number 1
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
1811 1 coeff 1.0000 avg resid 0.836E-11 max resid equL 2684 0.15028E-06 mix type 00000 avg corr 0.529E-07 max corr lnd 211 -0.16499E-04 mix type 00000 max resid
1811 2 coeff 1.0000 avg resid 0.272E-13 max resid equL 2518 0.41121E-10 mix type 00000 avg corr 0.175E-10 max corr lnT 2686 0.28520E-08 mix type 000xx okay!
1811 9.399050 2.084E+07 0.772050 6.046914 1.400000 0.000000 0.000000 0.000000 1.000000 28.700557 2686 0
-1.5499E+01 8.290527 -4.606047 5.452963 1.355829 -99.000000 1.400000 0.000000 0.000000 0.000000 44.028798 2
1.0000E-08 6.990363 5.018212 5.919235 -99.000000 -3.723369 0.000000 0.000000 0.000000 1.000000 0.000E+00 initial dt
tau_factor target current 1.0000000000000000D+00 1.0000000000000000D+02
solver_call_number 2
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
1812 1 coeff 1.0000 avg resid 0.400E-02 max resid equL 38 0.10603E+02 mix type 00000 avg corr 0.530E-03 max corr lnd 25 0.95381E-01 mix type 00000 avg+max resid
1812 2 coeff 1.0000 avg resid 0.319E+04 max resid equL 28 0.16914E+09 mix type 00000 avg corr 0.491E-03 max corr lnd 27 0.20412E+01 mix type 00000 max residual jumped -- give up
solver rejected trial model
s% model_number 1812
s% solver_call_number 2
nz 2686
s% num_retries 0
dt 1.0000000000000000D-08
log dt/secyer -1.5499103967085226D+01
retry: max residual jumped -- give up in solver 1812
retry log10(dt/yr), log10(dt), retry_factor 1 1812 -16.101164 -8.602060 0.500000
solver_call_number 3
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
1812 1 coeff 0.0100 avg resid 0.181E-04 max resid dv_dt 73 0.22215E+00 mix type 00000 avg corr 0.530E-03 max corr lnd 25 0.95381E-01 mix type 00000 avg+max corr+resid
1812 2 coeff 0.0100 avg resid 0.177E-04 max resid dv_dt 73 0.21993E+00 mix type 00000 avg corr 0.505E-03 max corr lnd 25 -0.15709E+00 mix type 00000 avg+max corr+resid
tol2 residual tolerances: norm, max 9.9999999999999997D+89 9.9999999999999997D+89
1812 3 coeff 0.0400 avg resid 0.180E-04 max resid dv_dt 73 0.21114E+00 mix type 00000 avg corr 0.441E-03 max corr lnd 68 -0.41882E-01 mix type 00000 okay!
1812 9.399050 2.080E+07 0.772050 6.046912 1.400000 0.000000 0.000000 0.000000 1.000000 28.700383 2686 1
-1.5800E+01 8.290523 -4.604521 5.452961 1.355829 -99.000000 1.400000 0.000000 0.000000 0.000000 44.028768 3
1.5000E-08 6.990350 5.017960 5.919233 -99.000000 -3.730244 0.000000 0.000000 0.000000 1.000000 0.000E+00 retry
tau_factor target current 1.0000000000000000D+00 7.9432823472428112D+01
solver_call_number 4
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.7764988665118665D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6647185452688291D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.7764988665118665D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6647185452688291D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D+00 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.280E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.459E-02 max corr lnd 32 -0.20621E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 4
nz 2686
s% num_retries 1
dt 5.0000000000000001D-09
log dt/secyer -1.5800133962749207D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 1 1813 -16.402194 -8.903090 0.500000
solver_call_number 5
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.7762166309271677D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6577176618983742D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.7762166309271677D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6577176618983742D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D+00 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.266E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.458E-02 max corr lnd 32 -0.20620E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 5
nz 2686
s% num_retries 2
dt 2.5000000000000001D-09
log dt/secyer -1.6101163958413188D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 2 1813 -16.703224 -9.204120 0.500000
solver_call_number 6
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.7760663515997678D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6542193595916380D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.7760663515997678D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6542193595916380D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.259E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.457E-02 max corr lnd 32 -0.20619E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 6
nz 2686
s% num_retries 3
dt 1.2500000000000000D-09
log dt/secyer -1.6402193954077170D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 3 1813 -17.004254 -9.505150 0.500000
solver_call_number 7
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.7759749037840308D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6524625347039866D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.7759749037840308D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6524625347039866D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.256E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.457E-02 max corr lnd 32 -0.20619E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 7
nz 2686
s% num_retries 4
dt 6.2500000000000001D-10
log dt/secyer -1.6703223949741151D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 4 1813 -17.305284 -9.806180 0.500000
solver_call_number 8
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.7758996356894414D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6515672981571981D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.7758996356894414D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6515672981571981D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.254E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.457E-02 max corr lnd 32 -0.20619E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 8
nz 2686
s% num_retries 5
dt 3.1250000000000001D-10
log dt/secyer -1.7004253945405132D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 5 1813 -17.606314 -10.107210 0.500000
solver_call_number 9
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.7758061254386817D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6510870936385245D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.7758061254386817D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6510870936385245D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.254E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.457E-02 max corr lnd 32 -0.20619E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 9
nz 2686
s% num_retries 6
dt 1.5625000000000000D-10
log dt/secyer -1.7305283941069113D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 6 1813 -17.907344 -10.408240 0.500000
solver_call_number 10
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.7756504331167193D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6507832619298995D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.7756504331167193D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6507832619298995D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.457E-02 max corr lnd 32 -0.20618E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 10
nz 2686
s% num_retries 7
dt 7.8125000000000002D-11
log dt/secyer -1.7606313936733095D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 7 1813 -18.208374 -10.709270 0.500000
solver_call_number 11
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.7753567016692713D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6505050000936969D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.7753567016692713D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6505050000936969D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.457E-02 max corr lnd 32 -0.20618E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 11
nz 2686
s% num_retries 8
dt 3.9062500000000001D-11
log dt/secyer -1.7907343932397076D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 8 1813 -18.509404 -11.010300 0.500000
solver_call_number 12
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.7747792983475534D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6501138859877848D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.7747792983475534D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6501138859877848D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.457E-02 max corr lnd 32 -0.20616E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 12
nz 2686
s% num_retries 9
dt 1.9531250000000000D-11
log dt/secyer -1.8208373928061057D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 9 1813 -18.810434 -11.311330 0.500000
solver_call_number 13
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.7736308687931725D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6494151478818168D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.7736308687931725D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6494151478818168D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.457E-02 max corr lnd 32 -0.20614E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 13
nz 2686
s% num_retries 10
dt 9.7656250000000002D-12
log dt/secyer -1.8509403923725039D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 10 1813 -19.111464 -11.612360 0.500000
solver_call_number 14
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.7713411341485354D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6480617195286861D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.7713411341485354D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6480617195286861D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.457E-02 max corr lnd 32 -0.20608E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 14
nz 2686
s% num_retries 11
dt 4.8828125000000001D-12
log dt/secyer -1.8810433919389020D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 11 1813 -19.412494 -11.913390 0.500000
solver_call_number 15
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.7667804232278002D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6453857805434993D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.7667804232278002D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6453857805434993D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.456E-02 max corr lnd 32 -0.20598E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 15
nz 2686
s% num_retries 12
dt 2.4414062500000001D-12
log dt/secyer -1.9111463915053001D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 12 1813 -19.713524 -12.214420 0.500000
solver_call_number 16
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.7577284558265234D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6400845222268970D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.7577284558265234D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6400845222268970D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.455E-02 max corr lnd 32 -0.20577E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 16
nz 2686
s% num_retries 13
dt 1.2207031250000000D-12
log dt/secyer -1.9412493910716982D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 13 1813 -20.014554 -12.515450 0.500000
solver_call_number 17
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.7398950304141820D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6296453139892570D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.7398950304141820D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6296453139892570D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.454E-02 max corr lnd 32 -0.20537E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 17
nz 2686
s% num_retries 14
dt 6.1035156250000001D-13
log dt/secyer -1.9713523906380964D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 14 1813 -20.315584 -12.816480 0.500000
solver_call_number 18
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.7052720557328954D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6093803496422780D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.7052720557328954D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.6093803496422780D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.450E-02 max corr lnd 32 -0.20458E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 18
nz 2686
s% num_retries 15
dt 3.0517578125000001D-13
log dt/secyer -2.0014553902044945D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 15 1813 -20.616614 -13.117510 0.500000
solver_call_number 19
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.6399266270713213D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.5711345637503982D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.6399266270713213D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.5711345637503982D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.444E-02 max corr lnd 32 -0.20309E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 19
nz 2686
s% num_retries 16
dt 1.5258789062500000D-13
log dt/secyer -2.0315583897708926D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 16 1813 -20.917644 -13.418540 0.500000
solver_call_number 20
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.5229302378731362D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.5026587129372885D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.5229302378731362D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.5026587129372885D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.432E-02 max corr lnd 32 -0.20042E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 20
nz 2686
s% num_retries 17
dt 7.6293945312500002D-14
log dt/secyer -2.0616613893372907D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 17 1813 -21.218674 -13.719570 0.500000
solver_call_number 21
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.3319639333893221D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.3908898561977708D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.3319639333893221D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.3908898561977708D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.414E-02 max corr lnd 32 -0.19607E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 21
nz 2686
s% num_retries 18
dt 3.8146972656250001D-14
log dt/secyer -2.0917643889036889D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 18 1813 -21.519704 -14.020600 0.500000
solver_call_number 22
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 8.0628342467869913D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.2333735894676629D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 8.0628342467869913D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.2333735894676629D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.387E-02 max corr lnd 32 -0.18994E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 22
nz 2686
s% num_retries 19
dt 1.9073486328125000D-14
log dt/secyer -2.1218673884700870D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 19 1813 -21.820734 -14.321630 0.500000
solver_call_number 23
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 7.7526908398509178D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.0518528397071215D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 7.7526908398509178D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 4.0518528397071215D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.357E-02 max corr lnd 32 -0.18287E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 23
nz 2686
s% num_retries 20
dt 9.5367431640625002D-15
log dt/secyer -2.1519703880364851D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 20 1813 -22.121764 -14.622660 0.500000
solver_call_number 24
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 7.4683370142408990D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 3.8854262283032590D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 7.4683370142408990D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 3.8854262283032590D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.329E-02 max corr lnd 32 -0.17639E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 24
nz 2686
s% num_retries 21
dt 4.7683715820312501D-15
log dt/secyer -2.1820733876028832D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 21 1813 -22.422794 -14.923690 0.500000
solver_call_number 25
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 7.2573303813524603D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 3.7619282661198596D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 7.2573303813524603D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 3.7619282661198596D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.309E-02 max corr lnd 32 -0.17159E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 25
nz 2686
s% num_retries 22
dt 2.3841857910156250D-15
log dt/secyer -2.2121763871692814D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 22 1813 -22.723824 -15.224720 0.500000
solver_call_number 26
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 7.1240479154627408D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 3.6839207052837843D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 7.1240479154627408D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 3.6839207052837843D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.296E-02 max corr lnd 32 -0.16855E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 26
nz 2686
s% num_retries 23
dt 1.1920928955078125D-15
log dt/secyer -2.2422793867356795D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 23 1813 -23.024854 -15.525750 0.500000
solver_call_number 27
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 7.0482306473971676D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 3.6395463753087846D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 7.0482306473971676D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 3.6395463753087846D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.289E-02 max corr lnd 32 -0.16682E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 27
nz 2686
s% num_retries 24
dt 5.9604644775390626D-16
log dt/secyer -2.2723823863020776D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 24 1813 -23.325884 -15.826780 0.500000
solver_call_number 28
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 7.0076508351482289D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 3.6157958232160027D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 7.0076508351482289D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 3.6157958232160027D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.285E-02 max corr lnd 32 -0.16590E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 28
nz 2686
s% num_retries 25
dt 2.9802322387695313D-16
log dt/secyer -2.3024853858684757D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 25 1813 -23.626914 -16.127810 0.500000
solver_call_number 29
gold tolerances level 0
correction tolerances: norm, max 3.0000000000000001D-03 2.9999999999999999D-01
tol1 residual tolerances: norm, max 1.0000000000000000D-10 1.0000000000000001D-09
hydro_mtx: logRho too large 1813 31 1 6.9866375842896275D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 3.6034971879855739D+01 -4.5203457618001934D+00
hydro_mtx: logRho too large 1813 31 1 6.9866375842896275D+01 -4.5216933510852124D+00
hydro_mtx: logRho too large 1813 41 1 3.6034971879855739D+01 -4.5203457618001934D+00
set_vars_for_solver failed: model, nz 1813 2686
eval_equations: set_solver_vars returned ierr -1
adjust_correction: eval_equations returned ierr -1 1.0000000000000000D-02 1.0000000000000000D+00
1813 1 coeff 1.0000 avg resid 0.253E-04 max resid dv_dt 111 0.26024E+00 mix type 00000 avg corr 0.285E-02 max corr lnd 32 -0.16542E+02 mix type 00000 adjust_correction failed in eval_equations -- give up
solver rejected trial model
s% model_number 1813
s% solver_call_number 29
nz 2686
s% num_retries 26
dt 1.4901161193847657D-16
log dt/secyer -2.3325883854348739D+01
retry: logRho > hydro_mtx_max_allowed_logRho 41 1813
retry log10(dt/yr), log10(dt), retry_factor 26 1813 -23.927944 -16.428840 0.500000
dt 7.4505805969238283D-17
min_timestep_limit 9.9999999999999998D-17
stopping because of problems dt < min_timestep_limit
failed in relax tau factor
star_relax_tau_factor ierr -1
do_star_job_controls_after ierr -1
do_before_evolve_loop ierr -1
do_before_evolve_loop ierr -1
DATE: 2022-10-24
TIME: 19:11:48
jpero at Jasons-MacBook-Pro burst22 %
Jason Pero
Graduate Research Assistant
University of Texas at Arlington?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20221025/f71904d9/attachment.htm>
More information about the Mesa-users
mailing list