[Mesa-users] crashes at Psurf=1 Ba
Sho Shibata
s.shibata423 at gmail.com
Sat Nov 22 02:19:08 UTC 2025
Hi,
The attached file shows the first 1000 lines of the output.
In extras_start_step, I calculate the next step's Psurf and output them like
---------------------
Set boundary conditions:
Times: 1000025.0142534898 1000028.3153213252
3.3010678353955778
Surface temperature: 93.672200228723625
Surface pressure: 0.69632028405282687
---------------------
In this case, I try to change Psurf linearly from 0.3 Ba to 10 Ba from t=0
to t=1000 yr. Once a retry happens and the timestep gets shorter, I
recalculate the next Psurf. The retries happen at model=17, 32, 39...
Best,
Sho
2025年11月21日(金) 19:57 Farag, Ebraheem <ebraheem.farag at yale.edu>:
> Can you add "report_ierr = .true.", and check if something gets spit out
> at the moment the model begins to timestep dive, rather than the last
> timestep? The last timestep is rarely informative.
>
> -EbF
> ------------------------------
> *From:* Mesa-users <mesa-users-bounces at lists.mesastar.org> on behalf of
> Sho Shibata via Mesa-users <mesa-users at lists.mesastar.org>
> *Sent:* Friday, November 21, 2025 8:50 PM
> *To:* Philip Mocz <pmocz at flatironinstitute.org>
> *Cc:* Jeremy Sakstein via Mesa-users <mesa-users at lists.mesastar.org>
> *Subject:* Re: [Mesa-users] crashes at Psurf=1 Ba
>
> Hi Philip,
>
> I wrote Psurf = 1 bar in the email, but that was a typo. The crashes occur
> at Psurf = 1 Ba. The figure is correct. Sorry for the confusion.
>
> Yes, the timestep is becoming really small and in most cases MESA stops
> because of the too small timestep.
> The bellow is one of the error message.
>
> ----------------
> solver rejected trial model
> s% model_number 135
> s% solver_call_number 169
> nz 354
> s% num_retries 34
> dt
> 1.4937230988914816D+07
> log dt/secyer
> -3.2483387012067205D-01
>
> retry: max residual jumped -- give up in solver 135
> retry log10(dt/yr), log10(dt), retry_factor 1 135
> -0.926894 6.572210 0.500000
> residual_norm > tol_residual_norm 135
> 9.4041696866144896D-03 1.0000000000000000D-08
> max_residual > tol_max_residual 135
> 4.8329438443848197D+01 1.0000000000000001D-05
> correction_norm > tol_correction_norm*coeff 135
> 2.1671809254193686D-04 3.0000000000000004D-07 1.0000000000000000D-02
> max_abs_correction > tol_max_correction*coeff 135
> 2.9765389377772104D-03 3.0000000000000001D-05 1.0000000000000000D-02
>
>
> ..... repeat this messages .....
>
>
> stopping because of problems dt < min_timestep_limit
> terminated evolution: cannot find acceptable model
> termination code: min_timestep_limit
> ----------------
>
> Best,
> Sho
>
> 2025年11月21日(金) 18:59 Philip Mocz <pmocz at flatironinstitute.org>:
>
> Hi Sho,
>
> It does look like something odd is going on right at 1 Ba
>
> Can you share what error messages you get when you crash? Is timestep
> becoming really small?
>
> Best,
> Philip
>
> On Fri, Nov 21, 2025, 7:54 PM Sho Shibata via Mesa-users <
> mesa-users at lists.mesastar.org> wrote:
>
> Hi all,
>
> I'm running test simulations where I smoothly change the surface pressure
> while keeping the surface temperature fixed, using
> atm_option='fixed_Tsurf_and_Psurf'. I'm using MESA Version: mesa-24.08.1,
> and I use run_star_extras.f90 to implement the smooth change of Psurf, but
> no other modification.
>
> I am encountering unexpected crashes: MESA can evolve with a varying
> surface pressure, but the code consistently crashes exactly when the model
> attempts to cross P_surf = 1 bar. The attached figures show what I’m doing.
> Changing the fixed temperature does not shift the crash location—the
> failure always occurs at P_surf = 1 bar.
>
> I see a similar issue when varying T_surf with fixed P_surf, specifically
> when crossing T_surf = 1000 K.
>
> Initially, I suspected the EOS tables, so I reran the experiment using the
> ideal EOS (by turning off other EOS tables), but the crash still occurs.
>
> Do you have any idea what might be happening internally in MESA in this
> situation? And is there a way to resolve or work around this issue?
>
> For context, these tests involve a planetary-sized object—an atmosphere on
> a Moon-sized protoplanet—so the envelope is extremely thin. The imposed
> change is rapid, but the timescale itself should not be the root of the
> problem.
>
> Any advice or suggestions would be greatly appreciated. Thanks in advance
> for your help.
>
> Best regards,
> Sho Shibata
> _______________________________________________
> mesa-users at lists.mesastar.org
>
> https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!DSb-azq1wVFtOg!Wig7tLQuDHlHxgddN8l0R6wshsctpY6atBNv6NpTKNt3dZ3b_Xnu2flgmR3Imum2iG_VS5Kvz8eGGEX8_xyg0alFQr6LHJU-$
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20251121/fec46a6d/attachment.htm>
-------------- next part --------------
reading user weak rate file /home/shibatasho/codes/MESA/mesa-24.08.1/data/rates_data/rate_tables/S13_r_be7_wk_li7.h5
version_number r24.08.1
load saved model model_P03_Tfree.mod
WARNING: changing to saved tau_factor = 1.5000000000000000D+00
WARNING -- inlist initial_z ignored 2.0000000000000000D-02
using saved initial_z instead 2.9999999999999999D-02
set_initial_age 1.0000000000000000D+06
set_age 1.0000000000000000D+06
set_initial_dt (years) 1.0000000000000000D+00
set_initial_model_number 0
set_initial_number_retries 0
net name basic.net
tau_factor 1.5000000000000000D+00
atm_option: fixed_Psurf_and_Tsurf
xmstar/mstar 4.9858253073330688D-06
xmstar (g) 2.9754858387044125D+21
M_center (g) 5.9678605246668854D+26
xmstar/Msun 1.4964147395123519D-12
M_center/Msun 3.0013231238747628D-07
R_center (cm) 3.1621095078850460D+08
R_center/Rsun 4.5452199337143104D-03
core density 4.5060845594279284D+00
L_center/Lsun 5.9550779532496703D-09
set_cumulative_energy_error 0.0000000000000000D+00
kap_option gs98
kap_CO_option gs98_co
kap_lowT_option lowT_Freedman11
OMP_NUM_THREADS 2
Set boundary conditions:
Times: 1000000.0000000000 1032622.7001355864 32622.700135586398
Surface temperature: 93.672200228723625
Surface pressure: 10.000000000000000
__________________________________________________________________________________________________________________________________________________
step lg_Tmax Teff lg_LH lg_Lnuc Mass H_rich H_cntr N_cntr Y_surf eta_cntr zones retry
lg_dt_yrs lg_Tcntr lg_R lg_L3a lg_Lneu lg_Mdot He_core He_cntr O_cntr Z_surf gam_cntr iters
age_yrs lg_Dcntr lg_L lg_LZ lg_Lphoto lg_Dsurf CO_core C_cntr Ne_cntr Z_cntr v_div_cs dt_limit
__________________________________________________________________________________________________________________________________________________
save LOGS/profile1.data for model 1
save photos/x001 for model 1
Set boundary conditions:
Times: 1000001.0000000000 1000002.0000000000 1.0000000000000000
Surface temperature: 93.672200228723625
Surface pressure: 0.41545114958073726
save LOGS/profile2.data for model 2
save photos/x002 for model 2
Set boundary conditions:
Times: 1000002.2000000000 1000003.3999999999 1.2000000000000000
Surface temperature: 93.672200228723625
Surface pressure: 0.43039365335545254
save LOGS/profile3.data for model 3
save photos/x003 for model 3
Set boundary conditions:
Times: 1000003.6400000000 1000005.0800000000 1.4399999999999999
Surface temperature: 93.672200228723625
Surface pressure: 0.44832465788685050
save LOGS/profile4.data for model 4
save photos/x004 for model 4
Set boundary conditions:
Times: 1000005.3680000000 1000007.0960000000 1.7280000000000000
Surface temperature: 93.672200228723625
Surface pressure: 0.46984186332452793
save LOGS/profile5.data for model 5
save photos/x005 for model 5
Set boundary conditions:
Times: 1000007.3902036850 1000009.4124073702 2.0222036850823848
Surface temperature: 93.672200228723625
Surface pressure: 0.49456538180632847
save LOGS/profile6.data for model 6
save photos/x006 for model 6
Set boundary conditions:
Times: 1000009.6533576389 1000011.9165115926 2.2631539537339203
Surface temperature: 93.672200228723625
Surface pressure: 0.52129222951921350
save LOGS/profile7.data for model 7
save photos/x007 for model 7
Set boundary conditions:
Times: 1000012.1816895372 1000014.7100214357 2.5283318984542404
Surface temperature: 93.672200228723625
Surface pressure: 0.55110796621861058
save LOGS/profile8.data for model 8
save photos/x008 for model 8
Set boundary conditions:
Times: 1000015.0539444768 1000017.9261994163 2.8722549394790398
Surface temperature: 93.672200228723625
Surface pressure: 0.58543493165964744
save LOGS/profile9.data for model 9
save photos/x009 for model 9
Set boundary conditions:
Times: 1000018.3160588202 1000021.5781731637 3.2621143434493596
Surface temperature: 93.672200228723625
Surface pressure: 0.62441323988075248
save LOGS/profile10.data for model 10
save photos/x010 for model 10
Set boundary conditions:
Times: 1000021.7131856544 1000025.1103124885 3.3971268340943803
Surface temperature: 93.672200228723625
Surface pressure: 0.66211252930746212
save LOGS/profile11.data for model 11
save photos/x011 for model 11
Set boundary conditions:
Times: 1000025.0142534898 1000028.3153213252 3.3010678353955778
Surface temperature: 93.672200228723625
Surface pressure: 0.69632028405282687
save LOGS/profile12.data for model 12
save photos/x012 for model 12
Set boundary conditions:
Times: 1000028.8871178986 1000032.7599823073 3.8728644088252562
Surface temperature: 93.672200228723625
Surface pressure: 0.74375911512927828
save LOGS/profile13.data for model 13
save photos/x013 for model 13
Set boundary conditions:
Times: 1000033.2380615333 1000037.5890051681 4.3509436347927410
Surface temperature: 93.672200228723625
Surface pressure: 0.79530032393609473
save LOGS/profile14.data for model 14
save photos/x014 for model 14
Set boundary conditions:
Times: 1000037.4298800703 1000041.6216986072 4.1918185369198513
Surface temperature: 93.672200228723625
Surface pressure: 0.83834213603679142
save LOGS/profile15.data for model 15
save photos/x015 for model 15
Set boundary conditions:
Times: 1000042.4600623145 1000047.4902445588 5.0301822443038215
Surface temperature: 93.672200228723625
Surface pressure: 0.90097840035038246
save LOGS/profile16.data for model 16
save photos/x016 for model 16
Set boundary conditions:
Times: 1000048.3026921805 1000054.1453220465 5.8426298659203963
Surface temperature: 93.672200228723625
Surface pressure: 0.97200948641290752
solver rejected trial model
s% model_number 17
s% solver_call_number 17
nz 460
s% num_retries 0
dt 1.7468532295716450D+08
log dt/secyer 7.4315245005291386D-01
retry: max residual jumped -- give up in solver 17
retry log10(dt/yr), log10(dt), retry_factor 1 17 0.141092 7.640196 0.500000
get_eos ierr -1
k 95
xa(j) equ_h1 1 7.1938774594784971D-01
xa(j) equ_he3 2 1.7415464504652846D-34
xa(j) equ_he4 3 2.8061225405215035D-01
xa(j) equ_c12 4 1.0000000000000002D-99
xa(j) equ_n14 5 1.0000000000000002D-99
xa(j) equ_o16 6 1.0000000000000002D-99
xa(j) equ_ne20 7 1.0000000000000002D-99
xa(j) equ_mg24 8 1.0000000000000002D-99
log10Rho 9.7696467619556160D-01
log10T 1.9716543405538405D+00
do_eos_for_cell: get_eos ierr -1
get_eos ierr -1
k 94
xa(j) equ_h1 1 7.1938774594784971D-01
xa(j) equ_he3 2 1.7414634979004675D-34
xa(j) equ_he4 3 2.8061225405215035D-01
xa(j) equ_c12 4 1.0000000000000002D-99
xa(j) equ_n14 5 1.0000000000000002D-99
xa(j) equ_o16 6 1.0000000000000002D-99
xa(j) equ_ne20 7 1.0000000000000002D-99
xa(j) equ_mg24 8 1.0000000000000002D-99
log10Rho 8.5670835038755055D-01
log10T 1.9716537601056767D+00
do_eos_for_cell: get_eos ierr -1
do_eos returned ierr -1
set_hydro_vars failed in call to set_micro_vars
set_hydro_vars returned ierr -1
eval_equations: set_solver_vars returned ierr -1
solver failed for hydro
hydro_solver_step returned ierr -1
s% model_number 17
nz 460
s% num_retries 1
retry: get_eos failed 17
retry log10(dt/yr), log10(dt), retry_factor 2 17 -0.159938 7.339166 0.500000
residual_norm > tol_residual_norm 17 6.8862853837799534D-04 1.0000000000000000D-08
max_residual > tol_max_residual 17 3.5349532394868932D-01 1.0000000000000001D-05
correction_norm > tol_correction_norm*coeff 17 5.6344059483153170D-04 3.0000000000000004D-07 1.0000000000000000D-02
max_abs_correction > tol_max_correction*coeff 17 9.4794713595840879D-03 3.0000000000000001D-05 1.0000000000000000D-02
solver rejected trial model
s% model_number 17
s% solver_call_number 19
nz 460
s% num_retries 2
dt 4.3671330739291124D+07
log dt/secyer 1.4109245872495155D-01
retry: avg+max corr+resid -- give up in solver 17
retry log10(dt/yr), log10(dt), retry_factor 3 17 -0.460968 7.038136 0.500000
residual_norm > tol_residual_norm 17 3.5178378355546236D-04 1.0000000000000000D-08
max_residual > tol_max_residual 17 1.5742498167789143D-01 1.0000000000000001D-05
correction_norm > tol_correction_norm*coeff 17 2.2595807397593118D-04 3.0000000000000004D-07 1.0000000000000000D-02
max_abs_correction > tol_max_correction*coeff 17 3.7216506880185592D-03 3.0000000000000001D-05 1.0000000000000000D-02
solver rejected trial model
s% model_number 17
s% solver_call_number 20
nz 460
s% num_retries 3
dt 2.1835665369645562D+07
log dt/secyer -1.5993753693902962D-01
retry: avg+max corr+resid -- give up in solver 17
retry log10(dt/yr), log10(dt), retry_factor 4 17 -0.761998 6.737106 0.500000
residual_norm > tol_residual_norm 17 1.8610651800032513D-04 1.0000000000000000D-08
max_residual > tol_max_residual 17 7.4333772915789353D-02 1.0000000000000001D-05
correction_norm > tol_correction_norm*coeff 17 1.0624434053518751D-04 3.0000000000000004D-07 1.0000000000000000D-02
max_abs_correction > tol_max_correction*coeff 17 2.3240675204075972D-03 3.0000000000000001D-05 1.0000000000000000D-02
solver rejected trial model
s% model_number 17
s% solver_call_number 21
nz 460
s% num_retries 4
dt 1.0917832684822781D+07
log dt/secyer -4.6096753260301082D-01
retry: avg+max corr+resid -- give up in solver 17
retry log10(dt/yr), log10(dt), retry_factor 5 17 -1.063028 6.436076 0.500000
residual_norm > tol_residual_norm 17 8.2657417075605652D-05 1.0000000000000000D-08
max_residual > tol_max_residual 17 4.9701360467261091D-02 1.0000000000000001D-05
correction_norm > tol_correction_norm*coeff 17 3.7494084569368033D-05 3.0000000000000004D-07 1.0000000000000000D-02
max_abs_correction > tol_max_correction*coeff 17 1.4548972412914936D-03 3.0000000000000001D-05 1.0000000000000000D-02
solver rejected trial model
s% model_number 17
s% solver_call_number 22
nz 460
s% num_retries 5
dt 5.4589163424113905D+06
log dt/secyer -7.6199752826699196D-01
retry: avg+max corr+resid -- give up in solver 17
retry log10(dt/yr), log10(dt), retry_factor 6 17 -1.364058 6.135046 0.500000
save LOGS/profile17.data for model 17
save photos/x017 for model 17
Set boundary conditions:
Times: 1000048.3891834907 1000048.4756748009 8.6491310213884939E-002
Surface temperature: 93.672200228723625
Surface pressure: 0.91149611114605567
save LOGS/profile18.data for model 18
save photos/x018 for model 18
Set boundary conditions:
Times: 1000048.4756748009 1000048.5621661112 8.6491310213884939E-002
Surface temperature: 93.672200228723625
Surface pressure: 0.91241925166768634
save LOGS/profile19.data for model 19
save photos/x019 for model 19
Set boundary conditions:
Times: 1000048.5794643732 1000048.6832539454 0.10378957225666192
Surface temperature: 93.672200228723625
Surface pressure: 0.91371164839647834
20 2.702683 187.889 -36.511016 -36.511016 0.000000 0.000000 0.719388 0.000000 0.280612*********** 455 6
-9.0467E-01 2.702683 -1.137682 -99.000000 -37.908884 -99.000000 0.000000 0.280612 0.000000 0.000000 9.916461 10
1.0000E+06 -5.138183 -8.225061 -99.000000 -99.000000 -9.827905 0.000000 0.000000 0.000000 0.000000 0.000E+00 max increase
rel_E_err 1.6614947623035167D-13
log_rel_run_E_err -10.9479127006996020
save LOGS/profile20.data for model 20
save photos/x020 for model 20
Set boundary conditions:
Times: 1000048.7040118598 1000048.8285593465 0.12454748670799430
Surface temperature: 93.672200228723625
Surface pressure: 0.91526252447078027
save LOGS/profile21.data for model 21
save photos/x021 for model 21
Set boundary conditions:
Times: 1000048.8534688439 1000049.0029258279 0.14945698404959315
Surface temperature: 93.672200228723625
Surface pressure: 0.91712357576093662
save LOGS/profile22.data for model 22
save photos/x022 for model 22
Set boundary conditions:
Times: 1000049.0328172247 1000049.2121656056 0.17934838085951177
Surface temperature: 93.672200228723625
Surface pressure: 0.91935683731036666
save LOGS/profile23.data for model 23
save photos/x023 for model 23
Set boundary conditions:
Times: 1000049.2480352818 1000049.4632533388 0.21521805703141411
Surface temperature: 93.672200228723625
Surface pressure: 0.92203675116844031
save LOGS/profile24.data for model 24
save photos/x024 for model 24
Set boundary conditions:
Times: 1000049.5062969503 1000049.7645586188 0.25826166843769693
Surface temperature: 93.672200228723625
Surface pressure: 0.92525264779937111
save LOGS/profile25.data for model 25
save photos/x025 for model 25
Set boundary conditions:
Times: 1000049.8162109524 1000050.1261249544 0.30991400212523629
Surface temperature: 93.672200228723625
Surface pressure: 0.92911172375350604
save LOGS/profile26.data for model 26
save photos/x026 for model 26
Set boundary conditions:
Times: 1000050.1881077549 1000050.5600045575 0.37189680255028351
Surface temperature: 93.672200228723625
Surface pressure: 0.93374261490120147
save LOGS/profile27.data for model 27
save photos/x027 for model 27
Set boundary conditions:
Times: 1000050.6343839180 1000051.0806600811 0.44627616306034018
Surface temperature: 93.672200228723625
Surface pressure: 0.93929968427818755
save LOGS/profile28.data for model 28
save photos/x028 for model 28
Set boundary conditions:
Times: 1000051.1699153136 1000051.7054467093 0.53553139567240826
Surface temperature: 93.672200228723625
Surface pressure: 0.94596816752858282
save LOGS/profile29.data for model 29
save photos/x029 for model 29
Set boundary conditions:
Times: 1000051.8125529884 1000052.4551906632 0.64263767480688982
Surface temperature: 93.672200228723625
Surface pressure: 0.95397034743029963
save LOGS/profile30.data for model 30
save photos/x030 for model 30
Set boundary conditions:
Times: 1000052.5837181981 1000053.3548834079 0.77116520976826775
Surface temperature: 93.672200228723625
Surface pressure: 0.96357296331235975
save LOGS/profile31.data for model 31
save photos/x031 for model 31
Set boundary conditions:
Times: 1000053.5091164499 1000054.4345147016 0.92539825172192125
Surface temperature: 93.672200228723625
Surface pressure: 0.97509610237083200
save LOGS/profile32.data for model 32
save photos/x032 for model 32
Set boundary conditions:
Times: 1000054.6195943519 1000055.7300722540 1.1104779020663056
Surface temperature: 93.672200228723625
Surface pressure: 0.98892386924075015
hydro_mtx: logRho too large 33 51 4 1.2006700827182824D+01 -9.7829015857949866D+00
hydro_mtx: logRho too large 33 52 4 1.2741390534446255D+01 -9.7757670512116608D+00
hydro_mtx: logRho too large 33 53 4 1.3984299543328341D+01 -9.7635339815523299D+00
hydro_mtx: logRho too large 33 54 4 1.5917153428200809D+01 -9.7449456268113366D+00
hydro_mtx: logRho too large 33 57 4 2.2187327510861376D+01 -9.6815831789055355D+00
hydro_mtx: logRho too large 33 58 4 2.4320067725921110D+01 -9.6570048969240236D+00
hydro_mtx: logRho too large 33 59 4 2.5895111063923849D+01 -9.6368264552427849D+00
hydro_mtx: logRho too large 33 60 4 2.6860048606136068D+01 -9.6229054435824946D+00
hydro_mtx: logRho too large 33 61 4 2.8165192641543882D+01 -9.6012962069240757D+00
hydro_mtx: logRho too large 33 62 4 2.9309582180815958D+01 -9.5773294446903598D+00
hydro_mtx: logRho too large 33 63 4 2.9928314493168546D+01 -9.5607963666923261D+00
hydro_mtx: logRho too large 33 64 4 3.0420862410611612D+01 -9.5435249634841881D+00
hydro_mtx: logRho too large 33 65 4 3.0774190369626432D+01 -9.5254791655472477D+00
hydro_mtx: logRho too large 33 66 4 3.1003535646108261D+01 -9.5040894962063867D+00
hydro_mtx: logRho too large 33 67 4 3.1009687938495667D+01 -9.4786507083089759D+00
hydro_mtx: logRho too large 33 68 4 3.0760236220170988D+01 -9.4517476056569336D+00
hydro_mtx: logRho too large 33 69 4 3.0376385016063868D+01 -9.4302436103447924D+00
hydro_mtx: logRho too large 33 70 4 3.0074280593200296D+01 -9.4156289877314929D+00
hydro_mtx: logRho too large 33 71 4 2.9561820626725687D+01 -9.3932884673814065D+00
hydro_mtx: logRho too large 33 72 4 2.8665032144608865D+01 -9.3615650678169793D+00
hydro_mtx: logRho too large 33 73 4 2.7854926966128943D+01 -9.3362377993126451D+00
hydro_mtx: logRho too large 33 74 4 2.7310776100351852D+01 -9.3190781414318291D+00
hydro_mtx: logRho too large 33 75 4 2.6745317315718619D+01 -9.3014739919476099D+00
hydro_mtx: logRho too large 33 76 4 2.6163992262311670D+01 -9.2834218241675508D+00
hydro_mtx: logRho too large 33 77 4 2.5571982272779962D+01 -9.2649192014307324D+00
hydro_mtx: logRho too large 33 78 4 2.4974117590286482D+01 -9.2459648762918238D+00
hydro_mtx: logRho too large 33 79 4 2.4374807076499668D+01 -9.2265588865171022D+00
hydro_mtx: logRho too large 33 80 4 2.3777989248384671D+01 -9.2067026459881554D+00
hydro_mtx: logRho too large 33 81 4 2.3187104015312357D+01 -9.1863990284861572D+00
hydro_mtx: logRho too large 33 82 4 2.2605083533648241D+01 -9.1656524422475005D+00
hydro_mtx: logRho too large 33 83 4 2.2034359878176847D+01 -9.1444688931770415D+00
hydro_mtx: logRho too large 33 84 4 2.1476886551027057D+01 -9.1228560346396073D+00
hydro_mtx: logRho too large 33 85 4 2.0934170613093130D+01 -9.1008232018756310D+00
hydro_mtx: logRho too large 33 86 4 2.0407312339861438D+01 -9.0783814292871643D+00
hydro_mtx: logRho too large 33 87 4 1.9897049594437540D+01 -9.0555434491106386D+00
hydro_mtx: logRho too large 33 88 4 1.9403804288901405D+01 -9.0323236703399861D+00
hydro_mtx: logRho too large 33 89 4 1.9040202926234890D+01 -9.0145155285822014D+00
hydro_mtx: logRho too large 33 90 4 1.8806681217057381D+01 -9.0026368254366140D+00
hydro_mtx: logRho too large 33 91 4 1.8577415395013009D+01 -8.9906723088090246D+00
hydro_mtx: logRho too large 33 92 4 1.8352373238042798D+01 -8.9786243667916938D+00
hydro_mtx: logRho too large 33 93 4 1.8026789907553304D+01 -8.9605557350893097D+00
hydro_mtx: logRho too large 33 94 4 1.7601101818244260D+01 -8.9359846747939322D+00
hydro_mtx: logRho too large 33 95 4 1.7191342312313061D+01 -8.9111268296044823D+00
hydro_mtx: logRho too large 33 96 4 1.6890742344363861D+01 -8.8921506089814990D+00
hydro_mtx: logRho too large 33 97 4 1.6697393075137441D+01 -8.8795309053321620D+00
hydro_mtx: logRho too large 33 98 4 1.6417202190983687D+01 -8.8606498553815705D+00
hydro_mtx: logRho too large 33 99 4 1.6051468977384729D+01 -8.8350728491737414D+00
hydro_mtx: logRho too large 33 100 4 1.5699035244649412D+01 -8.8093058239317177D+00
hydro_mtx: logRho too large 33 101 4 1.5440525344868673D+01 -8.7897136776418439D+00
hydro_mtx: logRho too large 33 102 4 1.5273722383476251D+01 -8.7767148162437856D+00
hydro_mtx: logRho too large 33 103 4 1.5109846379674183D+01 -8.7636830986563012D+00
hydro_mtx: logRho too large 33 104 4 1.4948812455546561D+01 -8.7506217956867225D+00
hydro_mtx: logRho too large 33 55 4 1.7950631369920789D+01 -9.7251975656601815D+00
hydro_mtx: logRho too large 33 105 4 1.4790537093954367D+01 -8.7375341762923462D+00
hydro_mtx: logRho too large 33 106 4 1.4634938308324770D+01 -8.7244235021699161D+00
hydro_mtx: logRho too large 33 107 4 1.4481935785380735D+01 -8.7112930224583618D+00
hydro_mtx: logRho too large 33 108 4 1.4331451003319035D+01 -8.6981459685753837D+00
hydro_mtx: logRho too large 33 109 4 1.4112100530824153D+01 -8.6785519307650762D+00
hydro_mtx: logRho too large 33 110 4 1.3825212783365249D+01 -8.6521904902068343D+00
hydro_mtx: logRho too large 33 111 4 1.3547091089526736D+01 -8.6258194076346246D+00
hydro_mtx: logRho too large 33 112 4 1.3342432122905173D+01 -8.6059036005169887D+00
hydro_mtx: logRho too large 33 113 4 1.3209496450263408D+01 -8.5927364662237036D+00
hydro_mtx: logRho too large 33 114 4 1.3078462574730455D+01 -8.5795802827211940D+00
hydro_mtx: logRho too large 33 115 4 1.2949272485962378D+01 -8.5664378608804164D+00
hydro_mtx: logRho too large 33 116 4 1.2760119770130947D+01 -8.5468948263857847D+00
hydro_mtx: logRho too large 33 56 4 2.0042296275517792D+01 -9.7042504122220183D+00
hydro_mtx: logRho too large 33 119 4 1.2209904365818119D+01 -8.4880141740689794D+00
hydro_mtx: logRho too large 33 120 4 1.2092160970622118D+01 -8.4750412601117446D+00
hydro_mtx: logRho too large 33 117 4 1.2512080010027885D+01 -8.5207190073288075D+00
hydro_mtx: logRho too large 33 118 4 1.2329140442979110D+01 -8.5010186797995768D+00
hydro_mtx: logRho too large 33 51 4 1.2006700827182824D+01 -9.7829015857949866D+00
hydro_mtx: logRho too large 33 52 4 1.2741390534446255D+01 -9.7757670512116608D+00
hydro_mtx: logRho too large 33 53 4 1.3984299543328341D+01 -9.7635339815523299D+00
hydro_mtx: logRho too large 33 54 4 1.5917153428200809D+01 -9.7449456268113366D+00
hydro_mtx: logRho too large 33 55 4 1.7950631369920789D+01 -9.7251975656601815D+00
hydro_mtx: logRho too large 33 56 4 2.0042296275517792D+01 -9.7042504122220183D+00
hydro_mtx: logRho too large 33 57 4 2.2187327510861376D+01 -9.6815831789055355D+00
hydro_mtx: logRho too large 33 58 4 2.4320067725921110D+01 -9.6570048969240236D+00
hydro_mtx: logRho too large 33 59 4 2.5895111063923849D+01 -9.6368264552427849D+00
hydro_mtx: logRho too large 33 60 4 2.6860048606136068D+01 -9.6229054435824946D+00
hydro_mtx: logRho too large 33 61 4 2.8165192641543882D+01 -9.6012962069240757D+00
hydro_mtx: logRho too large 33 62 4 2.9309582180815958D+01 -9.5773294446903598D+00
hydro_mtx: logRho too large 33 63 4 2.9928314493168546D+01 -9.5607963666923261D+00
hydro_mtx: logRho too large 33 64 4 3.0420862410611612D+01 -9.5435249634841881D+00
hydro_mtx: logRho too large 33 65 4 3.0774190369626432D+01 -9.5254791655472477D+00
hydro_mtx: logRho too large 33 66 4 3.1003535646108261D+01 -9.5040894962063867D+00
hydro_mtx: logRho too large 33 67 4 3.1009687938495667D+01 -9.4786507083089759D+00
hydro_mtx: logRho too large 33 68 4 3.0760236220170988D+01 -9.4517476056569336D+00
hydro_mtx: logRho too large 33 69 4 3.0376385016063868D+01 -9.4302436103447924D+00
hydro_mtx: logRho too large 33 70 4 3.0074280593200296D+01 -9.4156289877314929D+00
hydro_mtx: logRho too large 33 71 4 2.9561820626725687D+01 -9.3932884673814065D+00
hydro_mtx: logRho too large 33 72 4 2.8665032144608865D+01 -9.3615650678169793D+00
hydro_mtx: logRho too large 33 73 4 2.7854926966128943D+01 -9.3362377993126451D+00
hydro_mtx: logRho too large 33 74 4 2.7310776100351852D+01 -9.3190781414318291D+00
hydro_mtx: logRho too large 33 75 4 2.6745317315718619D+01 -9.3014739919476099D+00
hydro_mtx: logRho too large 33 76 4 2.6163992262311670D+01 -9.2834218241675508D+00
hydro_mtx: logRho too large 33 77 4 2.5571982272779962D+01 -9.2649192014307324D+00
hydro_mtx: logRho too large 33 78 4 2.4974117590286482D+01 -9.2459648762918238D+00
hydro_mtx: logRho too large 33 79 4 2.4374807076499668D+01 -9.2265588865171022D+00
hydro_mtx: logRho too large 33 80 4 2.3777989248384671D+01 -9.2067026459881554D+00
hydro_mtx: logRho too large 33 81 4 2.3187104015312357D+01 -9.1863990284861572D+00
hydro_mtx: logRho too large 33 82 4 2.2605083533648241D+01 -9.1656524422475005D+00
hydro_mtx: logRho too large 33 83 4 2.2034359878176847D+01 -9.1444688931770415D+00
hydro_mtx: logRho too large 33 84 4 2.1476886551027057D+01 -9.1228560346396073D+00
hydro_mtx: logRho too large 33 85 4 2.0934170613093130D+01 -9.1008232018756310D+00
hydro_mtx: logRho too large 33 86 4 2.0407312339861438D+01 -9.0783814292871643D+00
hydro_mtx: logRho too large 33 87 4 1.9897049594437540D+01 -9.0555434491106386D+00
hydro_mtx: logRho too large 33 88 4 1.9403804288901405D+01 -9.0323236703399861D+00
hydro_mtx: logRho too large 33 89 4 1.9040202926234890D+01 -9.0145155285822014D+00
hydro_mtx: logRho too large 33 90 4 1.8806681217057381D+01 -9.0026368254366140D+00
hydro_mtx: logRho too large 33 91 4 1.8577415395013009D+01 -8.9906723088090246D+00
hydro_mtx: logRho too large 33 92 4 1.8352373238042798D+01 -8.9786243667916938D+00
hydro_mtx: logRho too large 33 93 4 1.8026789907553304D+01 -8.9605557350893097D+00
hydro_mtx: logRho too large 33 94 4 1.7601101818244260D+01 -8.9359846747939322D+00
hydro_mtx: logRho too large 33 95 4 1.7191342312313061D+01 -8.9111268296044823D+00
hydro_mtx: logRho too large 33 96 4 1.6890742344363861D+01 -8.8921506089814990D+00
hydro_mtx: logRho too large 33 97 4 1.6697393075137441D+01 -8.8795309053321620D+00
hydro_mtx: logRho too large 33 98 4 1.6417202190983687D+01 -8.8606498553815705D+00
hydro_mtx: logRho too large 33 99 4 1.6051468977384729D+01 -8.8350728491737414D+00
hydro_mtx: logRho too large 33 100 4 1.5699035244649412D+01 -8.8093058239317177D+00
hydro_mtx: logRho too large 33 101 4 1.5440525344868673D+01 -8.7897136776418439D+00
hydro_mtx: logRho too large 33 102 4 1.5273722383476251D+01 -8.7767148162437856D+00
hydro_mtx: logRho too large 33 103 4 1.5109846379674183D+01 -8.7636830986563012D+00
hydro_mtx: logRho too large 33 104 4 1.4948812455546561D+01 -8.7506217956867225D+00
hydro_mtx: logRho too large 33 105 4 1.4790537093954367D+01 -8.7375341762923462D+00
hydro_mtx: logRho too large 33 106 4 1.4634938308324770D+01 -8.7244235021699161D+00
hydro_mtx: logRho too large 33 107 4 1.4481935785380735D+01 -8.7112930224583618D+00
hydro_mtx: logRho too large 33 108 4 1.4331451003319035D+01 -8.6981459685753837D+00
hydro_mtx: logRho too large 33 109 4 1.4112100530824153D+01 -8.6785519307650762D+00
hydro_mtx: logRho too large 33 110 4 1.3825212783365249D+01 -8.6521904902068343D+00
hydro_mtx: logRho too large 33 111 4 1.3547091089526736D+01 -8.6258194076346246D+00
hydro_mtx: logRho too large 33 112 4 1.3342432122905173D+01 -8.6059036005169887D+00
hydro_mtx: logRho too large 33 113 4 1.3209496450263408D+01 -8.5927364662237036D+00
hydro_mtx: logRho too large 33 114 4 1.3078462574730455D+01 -8.5795802827211940D+00
hydro_mtx: logRho too large 33 115 4 1.2949272485962378D+01 -8.5664378608804164D+00
hydro_mtx: logRho too large 33 116 4 1.2760119770130947D+01 -8.5468948263857847D+00
hydro_mtx: logRho too large 33 117 4 1.2512080010027885D+01 -8.5207190073288075D+00
hydro_mtx: logRho too large 33 118 4 1.2329140442979110D+01 -8.5010186797995768D+00
hydro_mtx: logRho too large 33 119 4 1.2209904365818119D+01 -8.4880141740689794D+00
hydro_mtx: logRho too large 33 120 4 1.2092160970622118D+01 -8.4750412601117446D+00
set_vars_for_solver failed: model, nz 33 449
eval_equations: set_solver_vars returned ierr -1
solver failed for hydro
hydro_solver_step returned ierr -1
s% model_number 33
nz 449
s% num_retries 6
retry: logRho > hydro_mtx_max_allowed_logRho 120 33
retry log10(dt/yr), log10(dt), retry_factor 1 33 -0.477369 7.021735 0.500000
residual_norm > tol_residual_norm 33 3.5160186415156459D-04 1.0000000000000000D-08
max_residual > tol_max_residual 33 6.6153665354484037D-01 1.0000000000000001D-05
correction_norm > tol_correction_norm*coeff 33 1.1508576751279624D-04 3.0000000000000004D-07 1.0000000000000000D-02
max_abs_correction > tol_max_correction*coeff 33 1.9642743794676286D-03 3.0000000000000001D-05 1.0000000000000000D-02
solver rejected trial model
s% model_number 33
s% solver_call_number 40
nz 449
s% num_retries 7
dt 2.1026410465348586D+07
log dt/secyer -1.7633882888058225D-01
retry: avg+max corr+resid -- give up in solver 33
retry log10(dt/yr), log10(dt), retry_factor 2 33 -0.778399 6.720705 0.500000
residual_norm > tol_residual_norm 33 3.9053758383352808D-06 1.0000000000000000D-08
max_residual > tol_max_residual 33 2.6514498513849748D-02 1.0000000000000001D-05
solver rejected trial model
s% model_number 33
s% solver_call_number 41
nz 449
s% num_retries 8
dt 1.0513205232674293D+07
log dt/secyer -4.7736882454456342D-01
retry: avg+max resid -- give up in solver 33
retry log10(dt/yr), log10(dt), retry_factor 3 33 -1.079429 6.419675 0.500000
save LOGS/profile33.data for model 33
save photos/x033 for model 33
Set boundary conditions:
Times: 1000054.7861660372 1000054.9527377225 0.16657168530994584
Surface temperature: 93.672200228723625
Surface pressure: 0.98062720911879930
save LOGS/profile34.data for model 34
save photos/x034 for model 34
Set boundary conditions:
Times: 1000054.9527377225 1000055.1193094079 0.16657168530994584
Surface temperature: 93.672200228723625
Surface pressure: 0.98240506485966106
save LOGS/profile35.data for model 35
save photos/x035 for model 35
Set boundary conditions:
Times: 1000055.1526237449 1000055.3525097673 0.19988602237193501
Surface temperature: 93.672200228723625
Surface pressure: 0.98489406289587356
save LOGS/profile36.data for model 36
save photos/x036 for model 36
Set boundary conditions:
Times: 1000055.3924869717 1000055.6323501986 0.23986322684632197
Surface temperature: 93.672200228723625
Surface pressure: 0.98788086053957713
save LOGS/profile37.data for model 37
save photos/x037 for model 37
Set boundary conditions:
Times: 1000055.6803228440 1000055.9681587162 0.28783587221558637
Surface temperature: 93.672200228723625
Surface pressure: 0.99146501771301532
save LOGS/profile38.data for model 38
save photos/x038 for model 38
Set boundary conditions:
Times: 1000056.0257258905 1000056.3711289372 0.34540304665870364
Surface temperature: 93.672200228723625
Surface pressure: 0.99576600631940171
hydro_mtx: logRho too large 39 71 5 1.2922955493478241D+01 -9.3749606302539821D+00
hydro_mtx: logRho too large 39 72 5 1.3564341104064011D+01 -9.3423472778145538D+00
hydro_mtx: logRho too large 39 73 5 1.3998128808215320D+01 -9.3163254446659458D+00
hydro_mtx: logRho too large 39 74 5 1.4257644626936228D+01 -9.2987087895673817D+00
hydro_mtx: logRho too large 39 75 5 1.4494301605947124D+01 -9.2806455977997935D+00
hydro_mtx: logRho too large 39 76 5 1.4707026989179440D+01 -9.2621335669911478D+00
hydro_mtx: logRho too large 39 77 5 1.4894939293409280D+01 -9.2431715845855908D+00
hydro_mtx: logRho too large 39 78 5 1.5057363624375032D+01 -9.2237598224537862D+00
hydro_mtx: logRho too large 39 79 5 1.5193842526347005D+01 -9.2038998262713321D+00
hydro_mtx: logRho too large 39 80 5 1.5304141934062729D+01 -9.1835945976546220D+00
hydro_mtx: logRho too large 39 81 5 1.5388252000283222D+01 -9.1628486669775899D+00
hydro_mtx: logRho too large 39 82 5 1.5446382794983267D+01 -9.1416681547763758D+00
hydro_mtx: logRho too large 39 83 5 1.5478955095727992D+01 -9.1200608197081205D+00
hydro_mtx: logRho too large 39 84 5 1.5486586698182958D+01 -9.0980360911456408D+00
hydro_mtx: logRho too large 39 85 5 1.5470074852931614D+01 -9.0756050847061758D+00
hydro_mtx: logRho too large 39 86 5 1.5430375583261744D+01 -9.0527805992709407D+00
hydro_mtx: logRho too large 39 87 5 1.5368580731513115D+01 -9.0295770944178884D+00
hydro_mtx: logRho too large 39 88 5 1.5285893611151971D+01 -9.0060106475922321D+00
hydro_mtx: logRho too large 39 89 5 1.5207782245962441D+01 -8.9879545392030451D+00
hydro_mtx: logRho too large 39 90 5 1.5152041072063840D+01 -8.9759182076267710D+00
hydro_mtx: logRho too large 39 91 5 1.5091827670525946D+01 -8.9638016328189458D+00
hydro_mtx: logRho too large 39 92 5 1.5027315772156133D+01 -8.9516073732336565D+00
hydro_mtx: logRho too large 39 93 5 1.4926133713469577D+01 -8.9333309960463900D+00
hydro_mtx: logRho too large 39 94 5 1.4773288209569328D+01 -8.9085034569487540D+00
hydro_mtx: logRho too large 39 95 5 1.4606392426764762D+01 -8.8834147174814326D+00
hydro_mtx: logRho too large 39 96 5 1.4469464303486635D+01 -8.8642829931334646D+00
hydro_mtx: logRho too large 39 97 5 1.4376846147966283D+01 -8.8515679036658383D+00
hydro_mtx: logRho too large 39 98 5 1.4236128354317659D+01 -8.8325559157621392D+00
hydro_mtx: logRho too large 39 99 5 1.4035225668907735D+01 -8.8068294134509699D+00
hydro_mtx: logRho too large 39 100 5 1.3825541733894442D+01 -8.7809409472852611D+00
hydro_mtx: logRho too large 39 101 5 1.3659992057550845D+01 -8.7612778359728036D+00
hydro_mtx: logRho too large 39 102 5 1.3549737795984045D+01 -8.7482394222399318D+00
hydro_mtx: logRho too large 39 103 5 1.3437926287030967D+01 -8.7351751736073968D+00
hydro_mtx: logRho too large 39 104 5 1.3324686463119995D+01 -8.7220883310758701D+00
hydro_mtx: logRho too large 39 105 5 1.3210142618457075D+01 -8.7089821225981048D+00
hydro_mtx: logRho too large 39 106 5 1.3094414391687353D+01 -8.6958597579838859D+00
hydro_mtx: logRho too large 39 107 5 1.2977616771347630D+01 -8.6827244239553423D+00
hydro_mtx: logRho too large 39 108 5 1.2859860120771865D+01 -8.6695792793707636D+00
hydro_mtx: logRho too large 39 109 5 1.2684560802453909D+01 -8.6499979382059582D+00
hydro_mtx: logRho too large 39 110 5 1.2444370004415873D+01 -8.6236798706288216D+00
hydro_mtx: logRho too large 39 111 5 1.2202105956380821D+01 -8.5973774720987208D+00
hydro_mtx: logRho too large 39 112 5 1.2016763147230202D+01 -8.5775324594731899D+00
hydro_mtx: logRho too large 39 69 5 1.2040863715206330D+01 -9.4129774027563329D+00
hydro_mtx: logRho too large 39 70 5 1.2406992910323755D+01 -9.3979380906960532D+00
hydro_mtx: logRho too large 39 69 5 1.2040863715206330D+01 -9.4129774027563329D+00
hydro_mtx: logRho too large 39 70 5 1.2406992910323755D+01 -9.3979380906960532D+00
hydro_mtx: logRho too large 39 71 5 1.2922955493478241D+01 -9.3749606302539821D+00
hydro_mtx: logRho too large 39 72 5 1.3564341104064011D+01 -9.3423472778145538D+00
hydro_mtx: logRho too large 39 73 5 1.3998128808215320D+01 -9.3163254446659458D+00
hydro_mtx: logRho too large 39 74 5 1.4257644626936228D+01 -9.2987087895673817D+00
hydro_mtx: logRho too large 39 75 5 1.4494301605947124D+01 -9.2806455977997935D+00
hydro_mtx: logRho too large 39 76 5 1.4707026989179440D+01 -9.2621335669911478D+00
hydro_mtx: logRho too large 39 77 5 1.4894939293409280D+01 -9.2431715845855908D+00
hydro_mtx: logRho too large 39 78 5 1.5057363624375032D+01 -9.2237598224537862D+00
hydro_mtx: logRho too large 39 79 5 1.5193842526347005D+01 -9.2038998262713321D+00
hydro_mtx: logRho too large 39 80 5 1.5304141934062729D+01 -9.1835945976546220D+00
hydro_mtx: logRho too large 39 81 5 1.5388252000283222D+01 -9.1628486669775899D+00
hydro_mtx: logRho too large 39 82 5 1.5446382794983267D+01 -9.1416681547763758D+00
hydro_mtx: logRho too large 39 83 5 1.5478955095727992D+01 -9.1200608197081205D+00
hydro_mtx: logRho too large 39 84 5 1.5486586698182958D+01 -9.0980360911456408D+00
hydro_mtx: logRho too large 39 85 5 1.5470074852931614D+01 -9.0756050847061758D+00
hydro_mtx: logRho too large 39 86 5 1.5430375583261744D+01 -9.0527805992709407D+00
hydro_mtx: logRho too large 39 87 5 1.5368580731513115D+01 -9.0295770944178884D+00
hydro_mtx: logRho too large 39 88 5 1.5285893611151971D+01 -9.0060106475922321D+00
hydro_mtx: logRho too large 39 89 5 1.5207782245962441D+01 -8.9879545392030451D+00
hydro_mtx: logRho too large 39 90 5 1.5152041072063840D+01 -8.9759182076267710D+00
hydro_mtx: logRho too large 39 91 5 1.5091827670525946D+01 -8.9638016328189458D+00
hydro_mtx: logRho too large 39 92 5 1.5027315772156133D+01 -8.9516073732336565D+00
hydro_mtx: logRho too large 39 93 5 1.4926133713469577D+01 -8.9333309960463900D+00
hydro_mtx: logRho too large 39 94 5 1.4773288209569328D+01 -8.9085034569487540D+00
hydro_mtx: logRho too large 39 95 5 1.4606392426764762D+01 -8.8834147174814326D+00
hydro_mtx: logRho too large 39 96 5 1.4469464303486635D+01 -8.8642829931334646D+00
hydro_mtx: logRho too large 39 97 5 1.4376846147966283D+01 -8.8515679036658383D+00
hydro_mtx: logRho too large 39 98 5 1.4236128354317659D+01 -8.8325559157621392D+00
hydro_mtx: logRho too large 39 99 5 1.4035225668907735D+01 -8.8068294134509699D+00
hydro_mtx: logRho too large 39 100 5 1.3825541733894442D+01 -8.7809409472852611D+00
hydro_mtx: logRho too large 39 101 5 1.3659992057550845D+01 -8.7612778359728036D+00
hydro_mtx: logRho too large 39 102 5 1.3549737795984045D+01 -8.7482394222399318D+00
hydro_mtx: logRho too large 39 103 5 1.3437926287030967D+01 -8.7351751736073968D+00
hydro_mtx: logRho too large 39 104 5 1.3324686463119995D+01 -8.7220883310758701D+00
hydro_mtx: logRho too large 39 105 5 1.3210142618457075D+01 -8.7089821225981048D+00
hydro_mtx: logRho too large 39 106 5 1.3094414391687353D+01 -8.6958597579838859D+00
hydro_mtx: logRho too large 39 107 5 1.2977616771347630D+01 -8.6827244239553423D+00
hydro_mtx: logRho too large 39 108 5 1.2859860120771865D+01 -8.6695792793707636D+00
hydro_mtx: logRho too large 39 109 5 1.2684560802453909D+01 -8.6499979382059582D+00
hydro_mtx: logRho too large 39 110 5 1.2444370004415873D+01 -8.6236798706288216D+00
hydro_mtx: logRho too large 39 111 5 1.2202105956380821D+01 -8.5973774720987208D+00
hydro_mtx: logRho too large 39 112 5 1.2016763147230202D+01 -8.5775324594731899D+00
set_vars_for_solver failed: model, nz 39 447
eval_equations: set_solver_vars returned ierr -1
solver failed for hydro
hydro_solver_step returned ierr -1
s% model_number 39
nz 447
s% num_retries 9
retry: logRho > hydro_mtx_max_allowed_logRho 112 39
retry log10(dt/yr), log10(dt), retry_factor 1 39 -0.984553 6.514551 0.500000
save LOGS/profile39.data for model 39
save photos/x039 for model 39
Set boundary conditions:
Times: 1000056.2329677185 1000056.4402095465 0.20724182799522217
Surface temperature: 93.672200228723625
Surface pressure: 0.99650331865117958
40 2.703070 193.310 -36.509258 -36.509258 0.000000 0.000000 0.719388 0.000000 0.280612*********** 443 10
-6.8352E-01 2.703070 -1.162431 -99.000000 -37.907126 -99.000000 0.000000 0.280612 0.000000 0.000000 9.910936 10
1.0001E+06 -5.137746 -8.225153 -99.000000 -99.000000 -9.790341 0.000000 0.000000 0.000000 0.000000 0.000E+00 hold
rel_E_err 1.7857052753862711D-13
log_rel_run_E_err -10.8305508690837957
save LOGS/profile40.data for model 40
save photos/x040 for model 40
Set boundary conditions:
Times: 1000056.4402095465 1000056.6474513744 0.20724182799522217
Surface temperature: 93.672200228723625
Surface pressure: 0.99871525564899821
hydro_mtx: logRho too large 41 73 3 1.3092062935218790D+01 -9.3104899419068854D+00
hydro_mtx: logRho too large 41 74 3 1.3505126895467571D+01 -9.2927410764679248D+00
hydro_mtx: logRho too large 41 75 3 1.3907321967815998D+01 -9.2745453284060915D+00
hydro_mtx: logRho too large 41 76 3 1.4296948450982418D+01 -9.2559007677133351D+00
hydro_mtx: logRho too large 41 77 3 1.4672320911718066D+01 -9.2368066820087975D+00
hydro_mtx: logRho too large 41 78 3 1.5031789488943152D+01 -9.2172636694274832D+00
hydro_mtx: logRho too large 41 72 3 1.2443100376023448D+01 -9.3367131333533919D+00
hydro_mtx: logRho too large 41 81 3 1.5999266601422233D+01 -9.1559684803767016D+00
hydro_mtx: logRho too large 41 82 3 1.6280095950028478D+01 -9.1346648188560025D+00
hydro_mtx: logRho too large 41 83 3 1.6538061503102661D+01 -9.1129376065780896D+00
hydro_mtx: logRho too large 41 84 3 1.6772167660636295D+01 -9.0907967851820253D+00
hydro_mtx: logRho too large 41 79 3 1.5373761796450360D+01 -9.1972737256697545D+00
hydro_mtx: logRho too large 41 80 3 1.5696724874817082D+01 -9.1768403231802136D+00
hydro_mtx: logRho too large 41 85 3 1.6981588057995175D+01 -9.0682539798032327D+00
hydro_mtx: logRho too large 41 86 3 1.7165675799450394D+01 -9.0453224900733034D+00
hydro_mtx: logRho too large 41 89 3 1.7536802029069705D+01 -8.9802304635954204D+00
hydro_mtx: logRho too large 41 87 3 1.7323970229268269D+01 -9.0220172613353480D+00
hydro_mtx: logRho too large 41 88 3 1.7456200089802230D+01 -8.9983548355591498D+00
hydro_mtx: logRho too large 41 91 3 1.7623293732716448D+01 -8.9559927253237159D+00
hydro_mtx: logRho too large 41 92 3 1.7656833658702816D+01 -8.9437586003838074D+00
hydro_mtx: logRho too large 41 90 3 1.7583294910582520D+01 -8.9681508810781221D+00
hydro_mtx: logRho too large 41 95 3 1.7730278069814108D+01 -8.8753791427786428D+00
hydro_mtx: logRho too large 41 96 3 1.7715112238280497D+01 -8.8562065851941600D+00
hydro_mtx: logRho too large 41 97 3 1.7699434886553494D+01 -8.8434666509388489D+00
hydro_mtx: logRho too large 41 98 3 1.7667972397669956D+01 -8.8244208735967060D+00
hydro_mtx: logRho too large 41 99 3 1.7603014913334615D+01 -8.7986566867487923D+00
hydro_mtx: logRho too large 41 100 3 1.7516986711428629D+01 -8.7727386448359397D+00
hydro_mtx: logRho too large 41 101 3 1.7436173403635305D+01 -8.7530591976118615D+00
hydro_mtx: logRho too large 41 102 3 1.7378474197121811D+01 -8.7400120773152565D+00
hydro_mtx: logRho too large 41 103 3 1.7316166700342347D+01 -8.7269411307781581D+00
hydro_mtx: logRho too large 41 104 3 1.7249403890204196D+01 -8.7138495860994407D+00
hydro_mtx: logRho too large 41 105 3 1.7178340336935609D+01 -8.7007406552399491D+00
hydro_mtx: logRho too large 41 106 3 1.7103131739680723D+01 -8.6876175290305007D+00
hydro_mtx: logRho too large 41 107 3 1.7023934483014774D+01 -8.6744833723441648D+00
hydro_mtx: logRho too large 41 108 3 1.6940905215478924D+01 -8.6613413194432649D+00
hydro_mtx: logRho too large 41 109 3 1.6812866151751706D+01 -8.6417673368059464D+00
hydro_mtx: logRho too large 41 110 3 1.6625551640380102D+01 -8.6154666055016680D+00
hydro_mtx: logRho too large 41 111 3 1.6425691424102961D+01 -8.5891886414642311D+00
hydro_mtx: logRho too large 41 112 3 1.6264889517986315D+01 -8.5693673981605674D+00
hydro_mtx: logRho too large 41 113 3 1.6156749779501418D+01 -8.5562698469937857D+00
hydro_mtx: logRho too large 41 114 3 1.6046141847798353D+01 -8.5431905706337865D+00
hydro_mtx: logRho too large 41 115 3 1.5933206955831613D+01 -8.5301321886116437D+00
hydro_mtx: logRho too large 41 116 3 1.5763041220352452D+01 -8.5107236228427787D+00
hydro_mtx: logRho too large 41 117 3 1.5524691669954432D+01 -8.4847564240013345D+00
hydro_mtx: logRho too large 41 118 3 1.5338085445883371D+01 -8.4652337889263336D+00
hydro_mtx: logRho too large 41 119 3 1.5213983983238636D+01 -8.4523518390264716D+00
hydro_mtx: logRho too large 41 120 3 1.5088457961133962D+01 -8.4395071546693128D+00
hydro_mtx: logRho too large 41 121 3 1.4961622828587222D+01 -8.4267017791246062D+00
hydro_mtx: logRho too large 41 122 3 1.4772240190425887D+01 -8.4076936096193382D+00
hydro_mtx: logRho too large 41 123 3 1.4511877699451604D+01 -8.3823362541490472D+00
hydro_mtx: logRho too large 41 124 3 1.4311393333582089D+01 -8.3633255560293112D+00
hydro_mtx: logRho too large 41 125 3 1.4178927689549223D+01 -8.3507946503922845D+00
hydro_mtx: logRho too large 41 126 3 1.4045862162825491D+01 -8.3383152444763606D+00
hydro_mtx: logRho too large 41 127 3 1.3848152747929928D+01 -8.3198063945546590D+00
hydro_mtx: logRho too large 41 128 3 1.3643840065590476D+01 -8.3011980434405039D+00
hydro_mtx: logRho too large 41 129 3 1.3509213452656043D+01 -8.2889383336404148D+00
hydro_mtx: logRho too large 41 130 3 1.3374390088971607D+01 -8.2767366321841163D+00
hydro_mtx: logRho too large 41 131 3 1.3176585290226699D+01 -8.2588341633798059D+00
hydro_mtx: logRho too large 41 132 3 1.2910716102753032D+01 -8.2351626280522030D+00
hydro_mtx: logRho too large 41 133 3 1.2645217410945184D+01 -8.2117328732672430D+00
hydro_mtx: logRho too large 41 134 3 1.2444255822121345D+01 -8.1942313948724710D+00
hydro_mtx: logRho too large 41 135 3 1.2312239465003120D+01 -8.1827053567021242D+00
hydro_mtx: logRho too large 41 136 3 1.2117059388902737D+01 -8.1656230213940706D+00
hydro_mtx: logRho too large 41 93 3 1.7697066917662042D+01 -8.9254258493431049D+00
hydro_mtx: logRho too large 41 94 3 1.7726029546050356D+01 -8.9005293232163307D+00
hydro_mtx: logRho too large 41 72 3 1.2443100376023448D+01 -9.3367131333533919D+00
hydro_mtx: logRho too large 41 73 3 1.3092062935218790D+01 -9.3104899419068854D+00
hydro_mtx: logRho too large 41 74 3 1.3505126895467571D+01 -9.2927410764679248D+00
hydro_mtx: logRho too large 41 75 3 1.3907321967815998D+01 -9.2745453284060915D+00
hydro_mtx: logRho too large 41 76 3 1.4296948450982418D+01 -9.2559007677133351D+00
hydro_mtx: logRho too large 41 77 3 1.4672320911718066D+01 -9.2368066820087975D+00
hydro_mtx: logRho too large 41 78 3 1.5031789488943152D+01 -9.2172636694274832D+00
hydro_mtx: logRho too large 41 79 3 1.5373761796450360D+01 -9.1972737256697545D+00
hydro_mtx: logRho too large 41 80 3 1.5696724874817082D+01 -9.1768403231802136D+00
hydro_mtx: logRho too large 41 81 3 1.5999266601422233D+01 -9.1559684803767016D+00
hydro_mtx: logRho too large 41 82 3 1.6280095950028478D+01 -9.1346648188560025D+00
hydro_mtx: logRho too large 41 83 3 1.6538061503102661D+01 -9.1129376065780896D+00
hydro_mtx: logRho too large 41 84 3 1.6772167660636295D+01 -9.0907967851820253D+00
hydro_mtx: logRho too large 41 85 3 1.6981588057995175D+01 -9.0682539798032327D+00
hydro_mtx: logRho too large 41 86 3 1.7165675799450394D+01 -9.0453224900733034D+00
hydro_mtx: logRho too large 41 87 3 1.7323970229268269D+01 -9.0220172613353480D+00
hydro_mtx: logRho too large 41 88 3 1.7456200089802230D+01 -8.9983548355591498D+00
hydro_mtx: logRho too large 41 89 3 1.7536802029069705D+01 -8.9802304635954204D+00
hydro_mtx: logRho too large 41 90 3 1.7583294910582520D+01 -8.9681508810781221D+00
hydro_mtx: logRho too large 41 91 3 1.7623293732716448D+01 -8.9559927253237159D+00
hydro_mtx: logRho too large 41 92 3 1.7656833658702816D+01 -8.9437586003838074D+00
hydro_mtx: logRho too large 41 93 3 1.7697066917662042D+01 -8.9254258493431049D+00
hydro_mtx: logRho too large 41 94 3 1.7726029546050356D+01 -8.9005293232163307D+00
hydro_mtx: logRho too large 41 95 3 1.7730278069814108D+01 -8.8753791427786428D+00
hydro_mtx: logRho too large 41 96 3 1.7715112238280497D+01 -8.8562065851941600D+00
hydro_mtx: logRho too large 41 97 3 1.7699434886553494D+01 -8.8434666509388489D+00
hydro_mtx: logRho too large 41 98 3 1.7667972397669956D+01 -8.8244208735967060D+00
hydro_mtx: logRho too large 41 99 3 1.7603014913334615D+01 -8.7986566867487923D+00
hydro_mtx: logRho too large 41 100 3 1.7516986711428629D+01 -8.7727386448359397D+00
hydro_mtx: logRho too large 41 101 3 1.7436173403635305D+01 -8.7530591976118615D+00
hydro_mtx: logRho too large 41 102 3 1.7378474197121811D+01 -8.7400120773152565D+00
hydro_mtx: logRho too large 41 103 3 1.7316166700342347D+01 -8.7269411307781581D+00
hydro_mtx: logRho too large 41 104 3 1.7249403890204196D+01 -8.7138495860994407D+00
hydro_mtx: logRho too large 41 105 3 1.7178340336935609D+01 -8.7007406552399491D+00
hydro_mtx: logRho too large 41 106 3 1.7103131739680723D+01 -8.6876175290305007D+00
hydro_mtx: logRho too large 41 107 3 1.7023934483014774D+01 -8.6744833723441648D+00
hydro_mtx: logRho too large 41 108 3 1.6940905215478924D+01 -8.6613413194432649D+00
hydro_mtx: logRho too large 41 109 3 1.6812866151751706D+01 -8.6417673368059464D+00
hydro_mtx: logRho too large 41 110 3 1.6625551640380102D+01 -8.6154666055016680D+00
hydro_mtx: logRho too large 41 111 3 1.6425691424102961D+01 -8.5891886414642311D+00
hydro_mtx: logRho too large 41 112 3 1.6264889517986315D+01 -8.5693673981605674D+00
hydro_mtx: logRho too large 41 113 3 1.6156749779501418D+01 -8.5562698469937857D+00
hydro_mtx: logRho too large 41 114 3 1.6046141847798353D+01 -8.5431905706337865D+00
hydro_mtx: logRho too large 41 115 3 1.5933206955831613D+01 -8.5301321886116437D+00
hydro_mtx: logRho too large 41 116 3 1.5763041220352452D+01 -8.5107236228427787D+00
hydro_mtx: logRho too large 41 117 3 1.5524691669954432D+01 -8.4847564240013345D+00
hydro_mtx: logRho too large 41 118 3 1.5338085445883371D+01 -8.4652337889263336D+00
hydro_mtx: logRho too large 41 119 3 1.5213983983238636D+01 -8.4523518390264716D+00
hydro_mtx: logRho too large 41 120 3 1.5088457961133962D+01 -8.4395071546693128D+00
hydro_mtx: logRho too large 41 121 3 1.4961622828587222D+01 -8.4267017791246062D+00
hydro_mtx: logRho too large 41 122 3 1.4772240190425887D+01 -8.4076936096193382D+00
hydro_mtx: logRho too large 41 123 3 1.4511877699451604D+01 -8.3823362541490472D+00
hydro_mtx: logRho too large 41 124 3 1.4311393333582089D+01 -8.3633255560293112D+00
hydro_mtx: logRho too large 41 125 3 1.4178927689549223D+01 -8.3507946503922845D+00
hydro_mtx: logRho too large 41 126 3 1.4045862162825491D+01 -8.3383152444763606D+00
hydro_mtx: logRho too large 41 127 3 1.3848152747929928D+01 -8.3198063945546590D+00
hydro_mtx: logRho too large 41 128 3 1.3643840065590476D+01 -8.3011980434405039D+00
hydro_mtx: logRho too large 41 129 3 1.3509213452656043D+01 -8.2889383336404148D+00
hydro_mtx: logRho too large 41 130 3 1.3374390088971607D+01 -8.2767366321841163D+00
hydro_mtx: logRho too large 41 131 3 1.3176585290226699D+01 -8.2588341633798059D+00
hydro_mtx: logRho too large 41 132 3 1.2910716102753032D+01 -8.2351626280522030D+00
hydro_mtx: logRho too large 41 133 3 1.2645217410945184D+01 -8.2117328732672430D+00
hydro_mtx: logRho too large 41 134 3 1.2444255822121345D+01 -8.1942313948724710D+00
hydro_mtx: logRho too large 41 135 3 1.2312239465003120D+01 -8.1827053567021242D+00
hydro_mtx: logRho too large 41 136 3 1.2117059388902737D+01 -8.1656230213940706D+00
set_vars_for_solver failed: model, nz 41 447
eval_equations: set_solver_vars returned ierr -1
solver failed for hydro
hydro_solver_step returned ierr -1
s% model_number 41
nz 447
s% num_retries 10
retry: logRho > hydro_mtx_max_allowed_logRho 136 41
retry log10(dt/yr), log10(dt), retry_factor 1 41 -1.206401 6.292703 0.500000
solver rejected trial model
s% model_number 41
s% solver_call_number 52
nz 447
s% num_retries 11
dt 3.9240328266852135D+06
log dt/secyer -9.0537133525075819D-01
retry: avg corr too large -- give up in solver 41
retry log10(dt/yr), log10(dt), retry_factor 2 41 -1.507431 5.991673 0.500000
save LOGS/profile41.data for model 41
save photos/x041 for model 41
Set boundary conditions:
Times: 1000056.5023820950 1000056.5645546434 6.2172548398566646E-002
Surface temperature: 93.672200228723625
Surface pressure: 0.99783048085086479
save LOGS/profile42.data for model 42
save photos/x042 for model 42
Set boundary conditions:
Times: 1000056.5645546434 1000056.6267271917 6.2172548398566646E-002
Surface temperature: 93.672200228723625
Surface pressure: 0.99849406195008605
get_eos ierr -1
k 60
xa(j) equ_h1 1 7.1938774594784971D-01
xa(j) equ_he3 2 1.5309352696426230D-34
xa(j) equ_he4 3 2.8061225405215035D-01
xa(j) equ_c12 4 1.0000000000000002D-99
xa(j) equ_n14 5 1.0000000000000002D-99
xa(j) equ_o16 6 1.0000000000000002D-99
xa(j) equ_ne20 7 1.0000000000000002D-99
xa(j) equ_mg24 8 1.0000000000000002D-99
log10Rho 9.3036235439489523D-01
log10T 1.9716119621324397D+00
do_eos_for_cell: get_eos ierr -1
get_eos ierr -1
k 61
xa(j) equ_h1 1 7.1938774594784971D-01
xa(j) equ_he3 2 1.5415796923546096D-34
xa(j) equ_he4 3 2.8061225405215035D-01
xa(j) equ_c12 4 1.0000000000000002D-99
xa(j) equ_n14 5 1.0000000000000002D-99
xa(j) equ_o16 6 1.0000000000000002D-99
xa(j) equ_ne20 7 1.0000000000000002D-99
xa(j) equ_mg24 8 1.0000000000000002D-99
log10Rho 1.5334100993451505D+00
log10T 1.9716123204854852D+00
do_eos_for_cell: get_eos ierr -1
do_eos returned ierr -1
set_hydro_vars failed in call to set_micro_vars
set_hydro_vars returned ierr -1
eval_equations: set_solver_vars returned ierr -1
solver failed for hydro
hydro_solver_step returned ierr -1
s% model_number 43
nz 447
s% num_retries 12
retry: get_eos failed 43
retry log10(dt/yr), log10(dt), retry_factor 1 43 -1.729280 5.769824 0.500000
save LOGS/profile43.data for model 43
save photos/x043 for model 43
Set boundary conditions:
Times: 1000056.6018581723 1000056.6391617013 3.7303529039139985E-002
Surface temperature: 93.672200228723625
Surface pressure: 0.99862677816893630
save LOGS/profile44.data for model 44
save photos/x044 for model 44
Set boundary conditions:
Times: 1000056.6391617014 1000056.6764652304 3.7303529039139985E-002
Surface temperature: 93.672200228723625
Surface pressure: 0.99902492682921462
save LOGS/profile45.data for model 45
save photos/x045 for model 45
Set boundary conditions:
Times: 1000056.6839259362 1000056.7286901710 4.4764234846967975E-002
Surface temperature: 93.672200228723625
Surface pressure: 0.99958233495236182
solver rejected trial model
s% model_number 46
s% solver_call_number 59
nz 443
s% num_retries 13
dt 1.6951821811280118D+06
log dt/secyer -1.2698875884358460D+00
retry: avg corr too large -- give up in solver 46
retry log10(dt/yr), log10(dt), retry_factor 1 46 -1.871948 5.627156 0.500000
save LOGS/profile46.data for model 46
save photos/x046 for model 46
Set boundary conditions:
Times: 1000056.7107844771 1000056.7376430180 2.6858540908180784E-002
Surface temperature: 93.672200228723625
Surface pressure: 0.99967789063082857
solver rejected trial model
s% model_number 47
s% solver_call_number 61
nz 447
s% num_retries 14
dt 8.4759109056400589D+05
log dt/secyer -1.5709175840998273D+00
retry: max residual jumped -- give up in solver 47
retry log10(dt/yr), log10(dt), retry_factor 1 47 -2.172978 5.326126 0.500000
save LOGS/profile47.data for model 47
save photos/x047 for model 47
Set boundary conditions:
Times: 1000056.7242137477 1000056.7376430181 1.3429270454090392E-002
Surface temperature: 93.672200228723625
Surface pressure: 0.99967789063207113
save LOGS/profile48.data for model 48
save photos/x048 for model 48
Set boundary conditions:
Times: 1000056.7376430181 1000056.7510722886 1.3429270454090392E-002
Surface temperature: 93.672200228723625
Surface pressure: 0.99982122414977126
save LOGS/profile49.data for model 49
save photos/x049 for model 49
Set boundary conditions:
Times: 1000056.7537581427 1000056.7698732673 1.6115124544908467E-002
Surface temperature: 93.672200228723625
Surface pressure: 1.0000218910745515
solver rejected trial model
s% model_number 50
s% solver_call_number 65
nz 443
s% num_retries 15
dt 6.1026558520608419D+05
log dt/secyer -1.7135850876685588D+00
retry: avg corr too large -- give up in solver 50
retry log10(dt/yr), log10(dt), retry_factor 1 50 -2.315645 5.183459 0.500000
do_kap_for_cell: get_kap ierr -1
eos info
More information about the Mesa-users
mailing list