MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (23.01-15-gf972714bd111) 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.018039874 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.000303937 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.015665316 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.938689935e-27 MLMG: Initial residual (resid0) = 4.938689935e-27 MLMG: Final Iter. 3 resid, resid/bnorm = 1.96290527e-30, 0.000397454648 MLMG: Timers: Solve = 0.002487909 Iter = 0.002299654 Bottom = 9.5181e-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 = 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.283312526e-27 MLMG: Initial residual (resid0) = 2.283312526e-27 MLMG: Final Iter. 3 resid, resid/bnorm = 1.801183061e-34, 7.888464852e-08 MLMG: Timers: Solve = 0.00226481 Iter = 0.002097582 Bottom = 9.6632e-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.42747586e-28 MLMG: Initial residual (resid0) = 7.42747586e-28 MLMG: Final Iter. 4 resid, resid/bnorm = 1.247408091e-37, 1.679450886e-10 MLMG: Timers: Solve = 0.003026288 Iter = 0.002857301 Bottom = 0.000116014 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.014007217 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: Failed to converge after 200 iterations. resid, resid/bnorm = 1.19485416e-29, 2.441693122e-10