Initializing AMReX (26.10-76-g347a3a2d2675)... MPI initialized with 4 MPI processes MPI initialized with thread support level 0 AMReX (26.10-76-g347a3a2d2675) 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.1280 model file = kepler_new_6.25e8.hybrid.hse.1280 reading initial model 1280 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 = 390625.000000 maximum radius (cell-centered) of input model = 499804687.500000 setting r_cutoff to 239 radius at r_cutoff 187109375 Maximum HSE Error = 0.003696 (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: 1.346149288 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 = 204587874.9 MLMG: Initial residual (resid0) = 204587874.9 MLMG: Final Iter. 9 resid, resid/bnorm = 0.01924523525, 9.40683081e-11 MLMG: Timers: Solve = 0.34047203 Iter = 0.329610888 Bottom = 0.000168309 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: 1.35025577 Call to firstdt for level 0 gives dt_lev = 0.01172063987 Multiplying dt_lev by init_shrink; dt_lev = 0.001172063987 Call to firstdt for level 1 gives dt_lev = 0.00578879243 Multiplying dt_lev by init_shrink; dt_lev = 0.000578879243 Call to firstdt for level 2 gives dt_lev = 0.002885281332 Multiplying dt_lev by init_shrink; dt_lev = 0.0002885281332 Minimum firstdt over all levels = 0.0002885281332 Doing initial divu iteration #1 Calling nodal solver MLMG: Initial rhs = 93695655.05 MLMG: Initial residual (resid0) = 93695655.05 MLMG: Final Iter. 4 resid, resid/bnorm = 0.005666818703, 6.048112583e-11 MLMG: Timers: Solve = 0.16078974 Iter = 0.149136691 Bottom = 7.475e-05 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 66.92844876 Call to estdt for level 1 gives dt_lev = 30.20205688 Call to estdt for level 2 gives dt_lev = 9.589324705 Minimum estdt over all levels = 9.589324705 Call to estdt at end of istep_divu_iter = 1 gives dt = 9.589324705 Multiplying dt by init_shrink; dt = 0.9589324705 Ignoring this new dt since it's larger than the previous dt = 0.0002885281332 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: 1.351419777 Doing initial pressure iteration #1  Timestep 0 starts with TIME = 0 DT = 0.0002885281332 Cell Count: Level 0, 32768 cells Level 1, 110592 cells Level 2, 262144 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 = 280128357.4 MLMG: Initial residual (resid0) = 280128357.4 MLMG: Final Iter. 5 resid, resid/bnorm = 0.7851953404, 2.802984131e-09 MLMG: Timers: Solve = 0.093864464 Iter = 0.085861373 Bottom = 0.000152652 <<< 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 = 280128358.1 MLMG: Initial residual (resid0) = 4804.307766 MLMG: Final Iter. 2 resid, resid/bnorm = 0.1073008775, 3.830418249e-10 MLMG: Timers: Solve = 0.044196332 Iter = 0.036738267 Bottom = 6.0367e-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 = 13780896.98 MLMG: Initial residual (resid0) = 13780896.98 MLMG: Final Iter. 30 resid, resid/bnorm = 0.001205376582, 8.746720796e-11 MLMG: Timers: Solve = 1.095734014 Iter = 1.084879494 Bottom = 0.000602333 Done calling nodal solver Timestep 0 ends with TIME = 0.0002885281332 DT = 0.0002885281332 Timing summary: Advection :1.928848203 seconds MAC Proj :0.190817946 seconds Nodal Proj :1.202597122 seconds Reactions :1.907875684 seconds Misc :0.378112294 seconds Base State :0.010193743 seconds Time to advance time step: 5.61257484 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: 1.34708221 Beginning main evolution  Timestep 1 starts with TIME = 0 DT = 0.0002885281332 Cell Count: Level 0, 32768 cells Level 1, 110592 cells Level 2, 262144 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 = 280128393.8 MLMG: Initial residual (resid0) = 280128393.8 MLMG: Final Iter. 5 resid, resid/bnorm = 0.7852297546, 2.803106619e-09 MLMG: Timers: Solve = 0.094877515 Iter = 0.08676489 Bottom = 0.000160799 <<< 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 = 280128393.4 MLMG: Initial residual (resid0) = 0.8028606176 MLMG: No iterations needed MLMG: Timers: Solve = 0.007416094 Iter = 0 Bottom = 0 <<< 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 = 50383856.33 MLMG: Initial residual (resid0) = 50383856.33 MLMG: Final Iter. 7 resid, resid/bnorm = 0.004018843421, 7.976450621e-11 MLMG: Timers: Solve = 0.268112502 Iter = 0.257485709 Bottom = 0.000171475 Done calling nodal solver Timestep 1 ends with TIME = 0.0002885281332 DT = 0.0002885281332 Timing summary: Advection :1.936533752 seconds MAC Proj :0.156886256 seconds Nodal Proj :0.374442122 seconds Reactions :1.936758491 seconds Misc :0.370868741 seconds Base State :0.009809689 seconds Time to advance time step: 4.780469431 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: 1.361178937 Call to estdt for level 0 gives dt_lev = 5.556499851 Call to estdt for level 1 gives dt_lev = 1.360999872 Call to estdt for level 2 gives dt_lev = 0.3374737097 Minimum estdt over all levels = 0.3374737097 Call to estdt at beginning of step 2 gives dt =0.3374737097 dt_growth factor limits the new dt = 0.0003173809465  Timestep 2 starts with TIME = 0.0002885281332 DT = 0.0003173809465 Cell Count: Level 0, 32768 cells Level 1, 110592 cells Level 2, 262144 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 = 255394708.2 MLMG: Initial residual (resid0) = 255394708.2 MLMG: Final Iter. 5 resid, resid/bnorm = 0.6347444654, 2.485346974e-09 MLMG: Timers: Solve = 0.093983095 Iter = 0.086437191 Bottom = 0.000151809 <<< 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 = 255394707.9 MLMG: Initial residual (resid0) = 0.6680596396 MLMG: No iterations needed MLMG: Timers: Solve = 0.007459608 Iter = 0 Bottom = 0 <<< 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 = 34141443.43 MLMG: Initial residual (resid0) = 34141443.43 MLMG: Final Iter. 4 resid, resid/bnorm = 0.002735068622, 8.010992939e-11 MLMG: Timers: Solve = 0.160420068 Iter = 0.149855644 Bottom = 8.1238e-05 Done calling nodal solver Timestep 2 ends with TIME = 0.0006059090797 DT = 0.0003173809465 Timing summary: Advection :1.966342004 seconds MAC Proj :0.155494437 seconds Nodal Proj :0.267097886 seconds Reactions :1.933779258 seconds Misc :0.370283033 seconds Base State :0.010242013 seconds Time to advance time step: 4.694295917 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: 1.366354533 Time to regrid: 0.181973595 Call to estdt for level 0 gives dt_lev = 8.534960391 Call to estdt for level 1 gives dt_lev = 2.077552712 Call to estdt for level 2 gives dt_lev = 0.4347551008 Minimum estdt over all levels = 0.4347551008 Call to estdt at beginning of step 3 gives dt =0.4347551008 dt_growth factor limits the new dt = 0.0003491190411  Timestep 3 starts with TIME = 0.0006059090797 DT = 0.0003491190411 Cell Count: Level 0, 32768 cells Level 1, 110592 cells Level 2, 262144 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 = 244508152.9 MLMG: Initial residual (resid0) = 244508152.9 MLMG: Final Iter. 5 resid, resid/bnorm = 0.6593789328, 2.696756427e-09 MLMG: Timers: Solve = 0.093303255 Iter = 0.085615662 Bottom = 0.000156683 <<< 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 = 244508152.5 MLMG: Initial residual (resid0) = 0.667964071 MLMG: No iterations needed MLMG: Timers: Solve = 0.007412592 Iter = 0 Bottom = 0 <<< 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 = 24597178.87 MLMG: Initial residual (resid0) = 24597178.87 MLMG: Final Iter. 5 resid, resid/bnorm = 0.001891462023, 7.689751874e-11 MLMG: Timers: Solve = 0.201717626 Iter = 0.19116102 Bottom = 0.000102044 Done calling nodal solver Timestep 3 ends with TIME = 0.0009550281208 DT = 0.0003491190411 Timing summary: Advection :1.977135478 seconds MAC Proj :0.155373883 seconds Nodal Proj :0.307920918 seconds Reactions :1.936654291 seconds Misc :0.370441868 seconds Base State :0.010640564 seconds Time to advance time step: 4.747840909 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: 1.374341879 Total Time: 32.84057873 Unused ParmParse Variables: [TOP]::amr.check_file(nvals = 1) :: [wdconvect_amr_chk] [TOP]::amr.checkpoint_files_output(nvals = 1) :: [0] [TOP]::amr.plot_file(nvals = 1) :: [wdconvect_amr_plt] AMReX (26.10-76-g347a3a2d2675) finalized