MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.11-12-g81e0635ce832) initialized Calling Setup() Calling ReadParameters() reading extern runtime parameters ... Calling VariableSetup() Calling BCSetup() Calling BaseStateGeometry::Init() Calling Init() Calling InitData() initdata model_File = 5peaks_const_grav.dat model file = 5peaks_const_grav.dat reading initial model 1000 points found in the initial model 10 variables found in the initial model WARNING: variable not found: gravity WARNING: variable not found: abar WARNING: variable not found: zbar WARNING: variable not found: gpot WARNING: variable not found: cno model file mapping, level: 0 dr of MAESTRO base state = 93750000.000000 dr of input file data = 10010010.010010 maximum radius (cell-centered) of input model = 10000000000.000000 WARNING: resolution of base state array is not an integer multiple of the initial model's resolution. make sure this is a desired property as this could lead to aliasing when performing the interpolation. modulus = 3659909.909910 setting r_cutoff to 105 radius at r_cutoff 9890625000 Maximum HSE Error = 0.054171 (after putting initial model into base state arrays, and for density < base_cutoff_density) Writing plotfile test_stability_pltInitData after InitData inner sponge: r_sp , r_tp : 9609375000, 9609375000 Time to write plotfile: 0.023641822 inner sponge: r_sp , r_tp : 9609375000, 9609375000 Doing initial projection Calling nodal solver MLMG: Initial rhs = 0 MLMG: Initial residual (resid0) = 0 MLMG: No iterations needed MLMG: Timers: Solve = 0.000306892 Iter = 0 Bottom = 0 Done calling nodal solver Writing plotfile test_stability_pltafter_InitProj after InitProj inner sponge: r_sp , r_tp : 9609375000, 9609375000 Time to write plotfile: 0.021663754 Call to firstdt for level 0 gives dt_lev = 0.01907406555 Multiplying dt_lev by init_shrink; dt_lev = 0.001907406555 Minimum firstdt over all levels = 0.001907406555 Doing initial divu iteration #1 Calling nodal solver MLMG: Initial rhs = 4.942022417e-27 MLMG: Initial residual (resid0) = 4.942022417e-27 MLMG: Final Iter. 3 resid, resid/bnorm = 1.977825976e-30, 0.000400205788 MLMG: Timers: Solve = 0.003146508 Iter = 0.002881421 Bottom = 0.000106616 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 1e+50 Minimum estdt over all levels = 1e+20 Call to estdt at end of istep_divu_iter = 1 gives dt = 1e+20 Multiplying dt by init_shrink; dt = 1e+19 Ignoring this new dt since it's larger than the previous dt = 0.001907406555 Doing initial divu iteration #2 Calling nodal solver MLMG: Initial rhs = 2.288227942e-27 MLMG: Initial residual (resid0) = 2.288227942e-27 MLMG: Final Iter. 3 resid, resid/bnorm = 1.28764231e-34, 5.627246685e-08 MLMG: Timers: Solve = 0.003060335 Iter = 0.002789781 Bottom = 9.488e-05 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 1e+50 Minimum estdt over all levels = 1e+20 Call to estdt at end of istep_divu_iter = 2 gives dt = 1e+20 Multiplying dt by init_shrink; dt = 1e+19 Ignoring this new dt since it's larger than the previous dt = 0.001907406555 Doing initial divu iteration #3 Calling nodal solver MLMG: Initial rhs = 7.441434633e-28 MLMG: Initial residual (resid0) = 7.441434633e-28 MLMG: Final Iter. 4 resid, resid/bnorm = 1.269827073e-37, 1.706427773e-10 MLMG: Timers: Solve = 0.002864227 Iter = 0.002659121 Bottom = 0.000111314 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 1e+50 Minimum estdt over all levels = 1e+20 Call to estdt at end of istep_divu_iter = 3 gives dt = 1e+20 Multiplying dt by init_shrink; dt = 1e+19 Ignoring this new dt since it's larger than the previous dt = 0.001907406555 Writing plotfile test_stability_pltafter_DivuIter after final DivuIter inner sponge: r_sp , r_tp : 9609375000, 9609375000 Time to write plotfile: 0.013908536 Doing initial pressure iteration #1 Timestep 0 starts with TIME = 0 DT = 0.001907406555 Cell Count: Level 0, 16384 cells inner sponge: r_sp , r_tp : 9609375000, 9609375000 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 4.893547633e-20 MLMG: Initial residual (resid0) = 4.893547633e-20 MLMG: Failing to converge after 142 iterations. resid, resid/bnorm = 5.097578092, 1.041693772e+20