MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.12-24-gea7e8a507e77) 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.017173047 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.000284866 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.014979414 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 = 6.462616052e-31, 0.0001309452479 MLMG: Timers: Solve = 0.002461671 Iter = 0.002280804 Bottom = 0.000119986 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.25137793e-27 MLMG: Initial residual (resid0) = 2.25137793e-27 MLMG: Final Iter. 3 resid, resid/bnorm = 1.41877687e-34, 6.301815662e-08 MLMG: Timers: Solve = 0.00228725 Iter = 0.002124671 Bottom = 0.000104582 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.32729272e-28 MLMG: Initial residual (resid0) = 7.32729272e-28 MLMG: Final Iter. 4 resid, resid/bnorm = 1.249213268e-37, 1.70487698e-10 MLMG: Timers: Solve = 0.00307682 Iter = 0.00288381 Bottom = 0.000138908 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.01373547 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 176 iterations. resid, resid/bnorm = 5.836468853, 1.192686634e+20