Patch from Bruno Giacomazzo:

  the patch allows the possibility to schedule
  AHFinderDirect also at POSTPOSTINITIAL by adding the parameter
  run_at_CCTK_POSTPOSTINITIAL (default equals to false). I have used it to
  be sure that an AH was found after initial data computation and before
  the last call to the Con2Prim routine of Whisky before evolution was
  started (in order to have hydro variables properly excised before the
  beginning of the evolution).



git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1550 f88db872-0e4f-0410-b76b-b9085cfa78c5

2 files changed, 53 insertions, 0 deletions

diff --git a/param.ccl b/param.ccl
index 6282504..dbb29b6 100644
--- a/param.ccl
+++ b/param.ccl
@@ -114,6 +114,9 @@ Boolean run_at_CCTK_POSTSTEP "should we run at CCTK_POSTSTEP?"
 Boolean run_at_CCTK_POSTINITIAL "should we run at CCTK_POSTINITIAL?"
 {
 } false
+Boolean run_at_CCTK_POSTPOSTINITIAL "should we run at CCTK_POSTPOSTINITIAL?"
+{
+} false
 Boolean run_at_CCTK_POST_RECOVER_VARIABLES				\
   "should we run at CCTK_POST_RECOVER_VARIABLES?"
 {
diff --git a/schedule.ccl b/schedule.ccl
index ea7bbfa..094444f 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -50,6 +50,14 @@ if (run_at_CCTK_POSTINITIAL != 0)
 		 options: global
 		 } "find apparent horizon(s) after this time step"
 	}
+if (run_at_CCTK_POSTPOSTINITIAL != 0)
+	{
+	schedule AHFinderDirect_find_horizons at CCTK_POSTPOSTINITIAL
+		 {
+		 lang: C
+		 options: global
+		 } "find apparent horizon(s) after this time step"
+	}
 if (run_at_CCTK_POST_RECOVER_VARIABLES != 0)
 	{
 	schedule AHFinderDirect_find_horizons at CCTK_POST_RECOVER_VARIABLES \
@@ -201,6 +209,48 @@ if (run_at_CCTK_POSTINITIAL != 0)
 		 options: global loop-local
 		 } "set mask(s) based on apparent horizon position(s)"
 	}
+if (run_at_CCTK_POSTPOSTINITIAL != 0)
+	{
+	schedule AHFinderDirect_import_mask at CCTK_POSTPOSTINITIAL \
+	   before AHFinderDirect_find_horizons
+		 {
+		 lang: C
+		 options: global loop-local
+		 } "import the excision mask"
+
+	schedule AHFinderDirect_store at CCTK_POSTPOSTINITIAL \
+	   after AHFinderDirect_find_horizons \
+	   before SphericalSurface_HasBeenSet
+		 {
+		 lang: C
+		 options: global
+		 } "store apparent horizon(s) into spherical surface(s)"
+
+	schedule AHFinderDirect_save at CCTK_POSTPOSTINITIAL \
+	   after AHFinderDirect_find_horizons \
+		 {
+		 lang: C
+		 options: global
+		 } "save apparent horizon(s) into Cactus variables"
+
+	if (which_horizon_to_announce_centroid != 0)
+		{
+		schedule AHFinderDirect_announce at CCTK_POSTPOSTINITIAL \
+		   before DriftCorrect \
+		   after AHFinderDirect_find_horizons
+			 {
+			 lang: C
+			 options: global
+			 } "announce horizon position(s) to other thorns"
+		}
+
+	schedule AHFinderDirect_maybe_do_masks at CCTK_POSTPOSTINITIAL \
+	   after AHFinderDirect_find_horizons
+		 {
+		 lang: C
+		 options: global loop-local
+		 } "set mask(s) based on apparent horizon position(s)"
+	}
 if (run_at_CCTK_POST_RECOVER_VARIABLES != 0)
 	{
 	schedule AHFinderDirect_import_mask at CCTK_POST_RECOVER_VARIABLES \