MPI initialized with 16 MPI processes MPI initialized with thread support level 0 AMReX (22.12-2-g1a5b308c8af7) initialized Calling Setup() Calling ReadParameters() reading extern runtime parameters ... Calling VariableSetup() Calling BCSetup() Calling BaseStateGeometry::Init() Calling Init() Calling InitData() initdata model_File = kepler_new_6.25e8.hybrid.hse.320 model file = kepler_new_6.25e8.hybrid.hse.320 reading initial model 320 points found in the initial model 6 variables found in the initial model WARNING: variable not found: ash WARNING: magnesium-24 not provided in inputs file model file mapping (spherical base state) dr of MAESTRO base state = 781250.000000 dr of input file data = 1562500.000000 maximum radius (cell-centered) of input model = 499218750.000000 setting r_cutoff to 239 radius at r_cutoff 187109375 Maximum HSE Error = 0.002311 (after putting initial model into base state arrays, and for density < base_cutoff_density) Writing plotfile wdconvect-amr_pltInitData after InitData inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Time to write plotfile: 0.918793642 inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Doing initial projection Calling nodal solver MLMG: Initial rhs = 214086148.6 MLMG: Initial residual (resid0) = 214086148.6 MLMG: Final Iter. 5 resid, resid/bnorm = 0.01199288666, 5.601897526e-11 MLMG: Timers: Solve = 0.176969371 Iter = 0.170565742 Bottom = 0.000156281 Done calling nodal solver Writing plotfile wdconvect-amr_pltafter_InitProj after InitProj inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Time to write plotfile: 0.915761087 Call to firstdt for level 0 gives dt_lev = 0.005789736987 Multiplying dt_lev by init_shrink; dt_lev = 0.0005789736987 Call to firstdt for level 1 gives dt_lev = 0.002885750885 Multiplying dt_lev by init_shrink; dt_lev = 0.0002885750885 Minimum firstdt over all levels = 0.0002885750885 Doing initial divu iteration #1 Calling nodal solver MLMG: Initial rhs = 98262295.14 MLMG: Initial residual (resid0) = 98262295.14 MLMG: Final Iter. 2 resid, resid/bnorm = 0.3496076327, 3.557902166e-09 MLMG: Timers: Solve = 0.074631179 Iter = 0.068007643 Bottom = 6.3109e-05 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 21.60061365 Call to estdt for level 1 gives dt_lev = 10.10291081 Minimum estdt over all levels = 10.10291081 Call to estdt at end of istep_divu_iter = 1 gives dt = 10.10291081 Multiplying dt by init_shrink; dt = 1.010291081 Ignoring this new dt since it's larger than the previous dt = 0.0002885750885 Doing initial divu iteration #2 Calling nodal solver MLMG: Initial rhs = 57185020 MLMG: Initial residual (resid0) = 57185020 MLMG: Final Iter. 2 resid, resid/bnorm = 0.01799653837, 3.147072147e-10 MLMG: Timers: Solve = 0.07550078 Iter = 0.06913304 Bottom = 6.1327e-05 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 21.60061365 Call to estdt for level 1 gives dt_lev = 10.47579575 Minimum estdt over all levels = 10.47579575 Call to estdt at end of istep_divu_iter = 2 gives dt = 10.47579575 Multiplying dt by init_shrink; dt = 1.047579575 Ignoring this new dt since it's larger than the previous dt = 0.0002885750885 Doing initial divu iteration #3 Calling nodal solver MLMG: Initial rhs = 36948393.74 MLMG: Initial residual (resid0) = 36948393.74 MLMG: Final Iter. 3 resid, resid/bnorm = 0.002165395877, 5.860595435e-11 MLMG: Timers: Solve = 0.11188035 Iter = 0.105709112 Bottom = 9.1785e-05 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 21.60061365 Call to estdt for level 1 gives dt_lev = 10.62248274 Minimum estdt over all levels = 10.62248274 Call to estdt at end of istep_divu_iter = 3 gives dt = 10.62248274 Multiplying dt by init_shrink; dt = 1.062248274 Ignoring this new dt since it's larger than the previous dt = 0.0002885750885 Writing plotfile wdconvect-amr_pltafter_DivuIter after final DivuIter inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Time to write plotfile: 0.888767309 Doing initial pressure iteration #1 Timestep 0 starts with TIME = 0 DT = 0.0002885750885 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 300731026.5 MLMG: Initial residual (resid0) = 300731026.5 MLMG: Final Iter. 22 resid, resid/bnorm = 0.2870682329, 9.545680613e-10 MLMG: Timers: Solve = 0.251433063 Iter = 0.242584499 Bottom = 0.000998055 <<< STEP 4 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 300731026.9 MLMG: Initial residual (resid0) = 12346.73858 MLMG: Final Iter. 2 resid, resid/bnorm = 0.2159021638, 7.179244723e-10 MLMG: Timers: Solve = 0.029831446 Iter = 0.023620483 Bottom = 9.1314e-05 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 16091441.2 MLMG: Initial residual (resid0) = 16091441.2 MLMG: Final Iter. 14 resid, resid/bnorm = 0.0008983202359, 5.582596516e-11 MLMG: Timers: Solve = 0.511449518 Iter = 0.505089636 Bottom = 0.000429239 Done calling nodal solver Timestep 0 ends with TIME = 0.0002885750885 DT = 0.0002885750885 Timing summary: Advection :1.240811585 seconds MAC Proj :0.316100646 seconds Nodal Proj :0.569130006 seconds Reactions :1.177245391 seconds Misc :0.257811112 seconds Base State :0.007294936 seconds Time to advance time step: 3.591864442 Writing plotfile 0 after all initialization inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Time to write plotfile: 0.888922774 Calling Evolve() Timestep 1 starts with TIME = 0 DT = 0.0002885750885 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 300732046.1 MLMG: Initial residual (resid0) = 300732046.1 MLMG: Final Iter. 22 resid, resid/bnorm = 0.2871757597, 9.549223749e-10 MLMG: Timers: Solve = 0.257554227 Iter = 0.251152493 Bottom = 0.001011009 <<< STEP 4 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 300732045.9 MLMG: Initial residual (resid0) = 20.36291351 MLMG: Final Iter. 2 resid, resid/bnorm = 0.1935165525, 6.434849731e-10 MLMG: Timers: Solve = 0.030018147 Iter = 0.02373301 Bottom = 9.4757e-05 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 27933872.46 MLMG: Initial residual (resid0) = 27933872.46 MLMG: Final Iter. 5 resid, resid/bnorm = 0.001084109014, 3.880983617e-11 MLMG: Timers: Solve = 0.17896642 Iter = 0.172557244 Bottom = 0.000150721 Done calling nodal solver Timestep 1 ends with TIME = 0.0002885750885 DT = 0.0002885750885 Timing summary: Advection :1.236705088 seconds MAC Proj :0.322570802 seconds Nodal Proj :0.237485416 seconds Reactions :1.178293027 seconds Misc :0.257767205 seconds Base State :0.006319470001 seconds Time to advance time step: 3.26383691 Writing plotfile 1 inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Time to write plotfile: 0.90556468 Call to estdt for level 0 gives dt_lev = 2.502707215 Call to estdt for level 1 gives dt_lev = 0.5259215932 Minimum estdt over all levels = 0.5259215932 Call to estdt at beginning of step 2 gives dt =0.5259215932 dt_growth factor limits the new dt = 0.0003174325974 Timestep 2 starts with TIME = 0.0002885750885 DT = 0.0003174325974 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 294777754.2 MLMG: Initial residual (resid0) = 294777754.2 MLMG: Final Iter. 5 resid, resid/bnorm = 0.2874117196, 9.750115654e-10 MLMG: Timers: Solve = 0.064418323 Iter = 0.057949737 Bottom = 0.00023579 <<< STEP 4 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 294777754.1 MLMG: Initial residual (resid0) = 721.0340273 MLMG: Final Iter. 2 resid, resid/bnorm = 0.008546076859, 2.899159364e-11 MLMG: Timers: Solve = 0.029879974 Iter = 0.023653793 Bottom = 9.0243e-05 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 21482439.32 MLMG: Initial residual (resid0) = 21482439.32 MLMG: Final Iter. 13 resid, resid/bnorm = 0.002090653579, 9.731918932e-11 MLMG: Timers: Solve = 0.436652526 Iter = 0.430290303 Bottom = 0.000403489 Done calling nodal solver Timestep 2 ends with TIME = 0.0006060076859 DT = 0.0003174325974 Timing summary: Advection :1.271188233 seconds MAC Proj :0.128439825 seconds Nodal Proj :0.494683255 seconds Reactions :1.18172723 seconds Misc :0.260729573 seconds Base State :0.006143661999 seconds Time to advance time step: 3.337281624 Writing plotfile 2 inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Time to write plotfile: 0.951610287 Time to regrid: 0.142459269 Call to estdt for level 0 gives dt_lev = 3.210043458 Call to estdt for level 1 gives dt_lev = 0.6510215137 Minimum estdt over all levels = 0.6510215137 Call to estdt at beginning of step 3 gives dt =0.6510215137 dt_growth factor limits the new dt = 0.0003491758571 Timestep 3 starts with TIME = 0.0006060076859 DT = 0.0003491758571 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 292704981.6 MLMG: Initial residual (resid0) = 292704981.6 MLMG: Final Iter. 5 resid, resid/bnorm = 0.2704412341, 9.239379277e-10 MLMG: Timers: Solve = 0.068558002 Iter = 0.062051043 Bottom = 0.000238952 <<< STEP 4 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 4a: thermal conduct >>> <<< STEP 5 : react state >>> <<< STEP 6 : make new S and new w0 >>> <<< STEP 7 : create MAC velocities >>> MLMG: Initial rhs = 292704981.3 MLMG: Initial residual (resid0) = 0.3555647364 MLMG: Final Iter. 1 resid, resid/bnorm = 0.0220395716, 7.529619583e-11 MLMG: Timers: Solve = 0.018831465 Iter = 0.012521503 Bottom = 4.5557e-05 <<< STEP 8 : advect base >>> : density_advance >>> : tracer_advance >>> : enthalpy_advance >>> <<< STEP 8a: thermal conduct >>> <<< STEP 9 : react state >>> <<< STEP 10: make new S >>> <<< STEP 11: update and project new velocity >>> Calling nodal solver MLMG: Initial rhs = 16800418.13 MLMG: Initial residual (resid0) = 16800418.13 MLMG: Final Iter. 5 resid, resid/bnorm = 0.001573332078, 9.364838815e-11 MLMG: Timers: Solve = 0.173977014 Iter = 0.167583846 Bottom = 0.000162273 Done calling nodal solver Timestep 3 ends with TIME = 0.0009551835431 DT = 0.0003491758571 Timing summary: Advection :1.270755275 seconds MAC Proj :0.121972576 seconds Nodal Proj :0.233046152 seconds Reactions :1.184833293 seconds Misc :0.260421642 seconds Base State :0.006561077999 seconds Time to advance time step: 3.071585932 Writing plotfile 3 inner sponge: r_sp , r_tp : 151953125, 183203125 outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125 Time to write plotfile: 0.916939274 Total Time: 23.80731084 Unused ParmParse Variables: [TOP]::maestro.v(nvals = 1) :: [1] [TOP]::amr.check_file(nvals = 1) :: [wdconvect-amr_chk] [TOP]::amr.checkpoint_files_output(nvals = 1) :: [0] AMReX (22.12-2-g1a5b308c8af7) finalized