GRHD: Add `atm_evolve` parameter
Adds the following new parameter:
-
GRHD.atm_evolve
: Iffalse
then the RHSs of all state variables are set to zero at any grid point where we identified the data to describe atmosphere.
This is parameter is implemented to work with all formulations and schemes, hopefully.