MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (23.01-29-g225c605680e8) 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.017126934 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.00024018 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.014545277 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.935357453e-27 MLMG: Initial residual (resid0) = 4.935357453e-27 MLMG: Final Iter. 3 resid, resid/bnorm = 1.943958012e-30, 0.0003938839346 MLMG: Timers: Solve = 0.002507659 Iter = 0.002326658 Bottom = 0.000116579 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.26897467e-27 MLMG: Initial residual (resid0) = 2.26897467e-27 MLMG: Final Iter. 3 resid, resid/bnorm = 1.768457277e-34, 7.794081178e-08 MLMG: Timers: Solve = 0.002435835 Iter = 0.002217712 Bottom = 0.000103176 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.378118428e-28 MLMG: Initial residual (resid0) = 7.378118428e-28 MLMG: Final Iter. 4 resid, resid/bnorm = 7.786789467e-38, 1.055389601e-10 MLMG: Timers: Solve = 0.003114912 Iter = 0.002912159 Bottom = 0.000126817 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.013463476 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.893547634e-20 MLMG: Initial residual (resid0) = 4.893547634e-20 MLMG: Failed to converge after 200 iterations. resid, resid/bnorm = 5.757092105e-30, 1.176465938e-10