Initializing CUDA... CUDA initialized with 1 GPU AMReX (19.11-284-g9290bd95f8e9) initialized Calling Setup() Calling ReadParameters() WARNING: burning_cutoff_density not supplied in the inputs file WARNING: setting burning_cutoff_density = base_cutoff_density reading extern runtime parameters ... ERROR: problem in the namelist Calling VariableSetup() Calling set_method_params() Initializing Helmholtz EOS and using Coulomb corrections. Calling BCSetup() Calling init_base_state_geometry() Calling Init() Calling InitData() ------------------------------------------------------------------------------ reading initial model 320 points found in the initial model file 6 variables found in the initial model file ------------------------------------------------------------------------------ model file mapping (spherical base state) dr of MAESTRO base state = 781250.0000 dr of input file data = 1562500.000 maximum radius (cell-centered) of input model = 499218750.0 setting r_cutoff to 239 radius at r_cutoff 187109375.0000000 Maximum HSE Error = 0.2310976572E-02 (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 : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Time to write plotfile: 0.804862867 inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Doing initial projection Calling nodal solver MLMG: Initial rhs = 188240648.4 MLMG: Initial residual (resid0) = 188240648.4 MLMG: Final Iter. 12 resid, resid/bnorm = 0.01633389294, 8.677133807e-11 MLMG: Timers: Solve = 2.514499838 Iter = 2.477853229 Bottom = 0.000835186 Done calling nodal solver Writing plotfile wdconvect-amr_pltafter_InitProj after InitProj inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Time to write plotfile: 0.701517395 Call to firstdt for level 0 gives dt_lev = 0.005789678157 Multiplying dt_lev by init_shrink; dt_lev = 0.0005789678157 Call to firstdt for level 1 gives dt_lev = 0.0005789678157 Multiplying dt_lev by init_shrink; dt_lev = 5.789678157e-05 Minimum firstdt over all levels = 5.789678157e-05 Doing initial divu iteration #1 Calling nodal solver MLMG: Initial rhs = 80008358.59 MLMG: Initial residual (resid0) = 80008358.59 MLMG: Final Iter. 2 resid, resid/bnorm = 0.482713986, 6.033294452e-09 MLMG: Timers: Solve = 0.460680178 Iter = 0.425987588 Bottom = 0.000145757 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 5.916319527 Call to estdt for level 1 gives dt_lev = 2.341804478 Minimum estdt over all levels = 2.341804478 Call to estdt at end of istep_divu_iter = 1 gives dt = 2.341804478 Multiplying dt by init_shrink; dt = 0.2341804478 Ignoring this new dt since it's larger than the previous dt = 5.789678157e-05 Doing initial divu iteration #2 Calling nodal solver MLMG: Initial rhs = 44170163.99 MLMG: Initial residual (resid0) = 44170163.99 MLMG: Final Iter. 2 resid, resid/bnorm = 0.02610203723, 5.909427285e-10 MLMG: Timers: Solve = 0.464621426 Iter = 0.425819118 Bottom = 0.000147652 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 5.916319527 Call to estdt for level 1 gives dt_lev = 2.341804478 Minimum estdt over all levels = 2.341804478 Call to estdt at end of istep_divu_iter = 2 gives dt = 2.341804478 Multiplying dt by init_shrink; dt = 0.2341804478 Ignoring this new dt since it's larger than the previous dt = 5.789678157e-05 Doing initial divu iteration #3 Calling nodal solver MLMG: Initial rhs = 26504865.08 MLMG: Initial residual (resid0) = 26504865.08 MLMG: Final Iter. 3 resid, resid/bnorm = 0.001814276911, 6.845071294e-11 MLMG: Timers: Solve = 0.669511181 Iter = 0.63043422 Bottom = 0.000209082 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 5.916319527 Call to estdt for level 1 gives dt_lev = 2.341804478 Minimum estdt over all levels = 2.341804478 Call to estdt at end of istep_divu_iter = 3 gives dt = 2.341804478 Multiplying dt by init_shrink; dt = 0.2341804478 Ignoring this new dt since it's larger than the previous dt = 5.789678157e-05 Writing plotfile wdconvect-amr_pltafter_DivuIter after final DivuIter inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Time to write plotfile: 0.726997612 Doing initial pressure iteration #1 Timestep 0 starts with TIME = 0 DT = 5.789678157e-05 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 173628411.4 MLMG: Initial residual (resid0) = 173628411.4 MLMG: Final Iter. 6 resid, resid/bnorm = 0.02237719041, 1.288797739e-10 MLMG: Timers: Solve = 0.124011877 Iter = 0.10264127 Bottom = 0.006688599 <<< 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 = 173628411.8 MLMG: Initial residual (resid0) = 2158.123377 MLMG: Final Iter. 2 resid, resid/bnorm = 0.04143690597, 2.386527962e-10 MLMG: Timers: Solve = 0.052923318 Iter = 0.036006226 Bottom = 0.002233826 <<< 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 = 35718516.1 MLMG: Initial residual (resid0) = 35718516.1 MLMG: Final Iter. 65 resid, resid/bnorm = 0.003011575318, 8.431412182e-11 MLMG: Timers: Solve = 0.793111123 Iter = 0.78382051 Bottom = 0.022872982 Done calling nodal solver Timestep 0 ends with TIME = 5.789678157e-05 DT = 5.789678157e-05 Timing summary: Advection :4.899143491 seconds MAC Proj :0.399626792 seconds Nodal Proj :0.911960788 seconds Reactions :0.218516179 seconds Misc :0.110656184 seconds Base State :0.04322494 seconds Time to advance time step: 6.540330275 Writing plotfile 0 after all initialization inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Time to write plotfile: 0.549561428 Calling Evolve() Timestep 1 starts with TIME = 0 DT = 5.789678157e-05 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 173628606.2 MLMG: Initial residual (resid0) = 173628606.2 MLMG: Final Iter. 6 resid, resid/bnorm = 0.02241812111, 1.291153664e-10 MLMG: Timers: Solve = 0.112153719 Iter = 0.103639358 Bottom = 0.006656832 <<< 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 = 173628606.2 MLMG: Initial residual (resid0) = 0.3535211849 MLMG: Final Iter. 1 resid, resid/bnorm = 0.02251003103, 1.296447142e-10 MLMG: Timers: Solve = 0.028141772 Iter = 0.019618124 Bottom = 0.001089041 <<< 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 = 16938159.33 MLMG: Initial residual (resid0) = 16938159.33 MLMG: Final Iter. 18 resid, resid/bnorm = 0.001470665447, 8.682557642e-11 MLMG: Timers: Solve = 0.225011839 Iter = 0.220175262 Bottom = 0.006345171 Done calling nodal solver Timestep 1 ends with TIME = 5.789678157e-05 DT = 5.789678157e-05 Timing summary: Advection :4.247419661 seconds MAC Proj :0.304615728 seconds Nodal Proj :0.340996536 seconds Reactions :0.19685321 seconds Misc :0.11761082 seconds Base State :0.030810744 seconds Time to advance time step: 5.20795374 Writing plotfile 1 inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Time to write plotfile: 0.550977562 Call to estdt for level 0 gives dt_lev = 1.112476812 Call to estdt for level 1 gives dt_lev = 0.2935011254 Minimum estdt over all levels = 0.2935011254 Timestep 2 starts with TIME = 5.789678157e-05 DT = 6.368645973e-05 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 166900886 MLMG: Initial residual (resid0) = 166900886 MLMG: Final Iter. 6 resid, resid/bnorm = 0.02281807736, 1.367163344e-10 MLMG: Timers: Solve = 0.11223242 Iter = 0.103707937 Bottom = 0.006720338 <<< 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 = 166900886 MLMG: Initial residual (resid0) = 142.3037968 MLMG: Final Iter. 1 resid, resid/bnorm = 0.09430622882, 5.650433084e-10 MLMG: Timers: Solve = 0.027973689 Iter = 0.01953689 Bottom = 0.001113755 <<< 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 = 11608423.25 MLMG: Initial residual (resid0) = 11608423.25 MLMG: Final Iter. 17 resid, resid/bnorm = 0.001011844724, 8.716469951e-11 MLMG: Timers: Solve = 0.211828237 Iter = 0.206819839 Bottom = 0.00604216 Done calling nodal solver Timestep 2 ends with TIME = 0.0001215832413 DT = 6.368645973e-05 Timing summary: Advection :4.21425959 seconds MAC Proj :0.300524913 seconds Nodal Proj :0.319152195 seconds Reactions :0.195139172 seconds Misc :0.111347162 seconds Base State :0.032035426 seconds Time to advance time step: 5.140817503 Writing plotfile 2 inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Time to write plotfile: 0.54306818 Time to regrid: 0.106728597 Call to estdt for level 0 gives dt_lev = 1.539392536 Call to estdt for level 1 gives dt_lev = 0.3670776924 Minimum estdt over all levels = 0.3670776924 Timestep 3 starts with TIME = 0.0001215832413 DT = 7.00551057e-05 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 166879353.2 MLMG: Initial residual (resid0) = 166879353.2 MLMG: Final Iter. 6 resid, resid/bnorm = 0.02258890867, 1.353607156e-10 MLMG: Timers: Solve = 0.11254106 Iter = 0.103983986 Bottom = 0.006685518 <<< 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 = 166879353.1 MLMG: Initial residual (resid0) = 0.3059132241 MLMG: Final Iter. 1 resid, resid/bnorm = 0.01941337054, 1.163317701e-10 MLMG: Timers: Solve = 0.028200605 Iter = 0.01956644 Bottom = 0.001110238 <<< 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 = 9504180.582 MLMG: Initial residual (resid0) = 9504180.582 MLMG: Final Iter. 18 resid, resid/bnorm = 0.0008338112384, 8.773099703e-11 MLMG: Timers: Solve = 0.224035851 Iter = 0.219109449 Bottom = 0.006331856 Done calling nodal solver Timestep 3 ends with TIME = 0.000191638347 DT = 7.00551057e-05 Timing summary: Advection :4.211049063 seconds MAC Proj :0.308062551 seconds Nodal Proj :0.330341535 seconds Reactions :0.19277138 seconds Misc :0.11205568 seconds Base State :0.032182478 seconds Time to advance time step: 5.154681206 Writing plotfile 3 inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Time to write plotfile: 0.568035865 Total Time: 58.71954537 Unused ParmParse Variables: [TOP]::maestro.v(nvals = 1) :: [1] [TOP]::system.regtest_reduction(nvals = 1) :: [1] [TOP]::amr.check_file(nvals = 1) :: [wdconvect-amr_chk] [TOP]::amr.checkpoint_files_output(nvals = 1) :: [0] Total GPU global memory (MB): 12066 Free GPU global memory (MB): 9480 [The Arena] space (MB): 9049 [The Device Arena] space (MB): 8 [The Managed Arena] space (MB): 8 [The Pinned Arena] space (MB): 8 AMReX (19.11-284-g9290bd95f8e9) finalized