MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (23.01-1-g30f2290cf7c3) 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.019069858 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.000259732 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.015993175 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.513921765e-30, 0.0003067501755 MLMG: Timers: Solve = 0.002543753 Iter = 0.002335137 Bottom = 0.000110277 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.268862229e-27 MLMG: Initial residual (resid0) = 2.268862229e-27 MLMG: Final Iter. 3 resid, resid/bnorm = 1.339713804e-34, 5.904782522e-08 MLMG: Timers: Solve = 0.002452305 Iter = 0.002258137 Bottom = 9.6824e-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.37819286e-28 MLMG: Initial residual (resid0) = 7.37819286e-28 MLMG: Final Iter. 4 resid, resid/bnorm = 1.231219843e-37, 1.668728192e-10 MLMG: Timers: Solve = 0.003277642 Iter = 0.002968906 Bottom = 0.000127832 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.014085126 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.757086086e-30, 1.176464708e-10