re-add update of Tmunu after AtmosphereReset, fix typo in routine name

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@264 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45

1 files changed, 9 insertions, 3 deletions

diff --git a/schedule.ccl b/schedule.ccl
index 152de5f..5d38b42 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -967,8 +967,10 @@ schedule GRHydro_Boundaries IN HydroBase_Select_Boundaries AS GRHydro_Bound
 ### If a point has been marked as requiring resetting to the ###
 ### atmosphere, this is where we do it.                      ###
 ###                                                          ###
-### You didn't see this abuse of MoL. Nothing to see here,   ###
-### move along.                                              ###
+### This is logically still part of the MoL evolution step   ###
+### since it is part of evolving the primitives. The final   ###
+### SetTmunu should be before PseudoEvolution so that eg.    ###
+### constraints use the correct stress energy tensor.        ###
 ################################################################
 
 if (CCTK_Equals(Bvec_evolution_method,"GRHydro"))
@@ -984,10 +986,14 @@ if (CCTK_Equals(Bvec_evolution_method,"GRHydro"))
   } "Reset the atmosphere"
 }
 
-schedule group HydroBase_Boundaries AT EVOL AFTER MoL_Evolution AFTER (GRHydro_AtmosphereRese GRHydro_AtmosphereResetM)
+schedule group HydroBase_Boundaries AT EVOL AFTER MoL_Evolution AFTER (GRHydro_AtmosphereReset GRHydro_AtmosphereResetM)
 {
 } "HydroBase Boundary conditions group"
 
+schedule GROUP SetTmunu AT EVOL AFTER HydroBase_Boundaries BEFORE MoL_PseudoEvolution
+{
+} "Calculate the stress-energy tensor"
+
 } # end of  if (CCTK_EQUALS(hydrobase::evolution_method, "GRHydro"))
   
 if (set_trivial_rp_grid_function)