let GRHydro set the atmosphere of initial data, so that ID thorns don't have to worry about it. This changes the data slitely for EOSG testsuites because this now uses EOSG calls instead of the otherwise fixed calculations of the TOVSolver, thus the change of two testsuites

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

1 files changed, 13 insertions, 5 deletions

diff --git a/schedule.ccl b/schedule.ccl
index 5b565e2..798bb05 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -230,19 +230,27 @@ if (rho_abs_min_after_recovery > 0.0)
 
 if (CCTK_IsThornActive("PUGH" ))
 {
-   schedule GRHydro_Rho_Minima_Setup_Final_PUGH  AT CCTK_PostInitial BEFORE MoL_PostStep
+   schedule GRHydro_Rho_Minima_Setup_Final_PUGH AT CCTK_PostInitial BEFORE MoL_PostStep
    {
      LANG: C
    } "Set the value of the rest-mass density of the atmosphere which will be used during the evolution (PUGH)"
-}  
+   schedule GRHydro_InitialAtmosphereReset AT CCTK_PostInitial BEFORE MoL_PostStep AFTER GRHydro_Rho_Minima_Setup_Final_PUGH
+   {
+     LANG: Fortran
+   } "Use mask to enforce atmosphere at initial time"
+}
 else if (CCTK_IsThornActive("Carpet"))
- {
-#   schedule GRHydro_Rho_Minima_Setup_Final  AT CCTK_PostInitial BEFORE MoL_PostStep
+{
    schedule GRHydro_Rho_Minima_Setup_Final AT CCTK_PostPostInitial BEFORE Con2Prim
    {
      LANG: C
    } "Set the value of the rest-mass density of the atmosphere which will be used during the evolution"
- }    
+   schedule GRHydro_InitialAtmosphereReset AT CCTK_PostPostInitial BEFORE Con2Prim AFTER GRHydro_Rho_Minima_Setup_Final
+   {
+     LANG: Fortran
+   } "Use mask to enforce atmosphere at initial time"
+}
+
 
 # this is a temporary hack, because carpet does allow only one loop over the levels for Basegrid-Iniital-Postinitial
 # so we must compute GRHydro_rho_min at postpostinitial and the constraints after this (and afer C2P)