User: rhaas Date: 2011/08/05 01:07 PM

Modified: /trunk/ schedule.ccl

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

File Changes:

Directory: /trunk/ ==================

File [modified]: schedule.ccl Delta lines: +9 -3 =================================================================== --- trunk/schedule.ccl 2011-08-04 17:05:23 UTC (rev 263) +++ trunk/schedule.ccl 2011-08-05 18:07:18 UTC (rev 264) @@ -967,8 +967,10 @@ ### 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 @@ } "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)