Michael Zingale
Michael Zingale
this seems to work, but I need to 1. figure out the T source term 2. figure out the T geometry term
this is part of #781
one thing that would be useful would be to make a plot of the average de/dt over a timestep and another with the instantaneous de/dt -- we can hack the...
tests pass: http://groot.astro.sunysb.edu/Castro/test-suite/gfortran/2022-04-28-005/index.html
I think one first step for this would be to break the main `Castro_F.H` into separate files for the different physics. I've been wanting to do that for a while.
we should make sure that all the *_F.H and .H headers in each Source directory match the contents of those directories.
it surely depends on resolution though, right?
there seem to be a lot of CFL violations and retries
this discussion may help: https://github.com/AMReX-Codes/amrex/discussions/2939#discussioncomment-3579931
I think the way to do this is to add a new runtime parameter, `radial_sponge_dir`. If it is 0, then it is radius, like now, computed from the center. If...