add new option run_at_CCTK_POSTINITIAL

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

2 files changed, 45 insertions, 1 deletions

diff --git a/param.ccl b/param.ccl
index eb712a4..4d31178 100644
--- a/param.ccl
+++ b/param.ccl
@@ -55,7 +55,7 @@ int find_every "how often should we try to find apparent horizons?"
 # These parameters control which schedule bin (or even bins, if you
 # want) AHFinderDirect runs in.  Historically, AHFinderDirect has run
 # at CCTK_ANALYSIS.  However, for some purposes it's more convenient
-# to run at CCTK_POSTSTEP.
+# to run at CCTK_POSTSTEP and/or CCTK_POSTINITIAL
 #
 Boolean run_at_CCTK_ANALYSIS "should we run in CCTK_ANALYSIS schedule bin?"
 {
@@ -63,6 +63,10 @@ Boolean run_at_CCTK_ANALYSIS "should we run in CCTK_ANALYSIS schedule bin?"
 Boolean run_at_CCTK_POSTSTEP "should we run in CCTK_POSTSTEP schedule bin?"
 {
 } false
+Boolean run_at_CCTK_POSTINITIAL \
+  "should we run in CCTK_POSTINITIAL schedule bin?"
+{
+} false
 
 #
 # unless you're trying to debug AHFinderDirect, you should leave
diff --git a/schedule.ccl b/schedule.ccl
index 3cd6cb0..92313c6 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -42,6 +42,20 @@ if (run_at_CCTK_POSTSTEP != 0)
 		 options: global
 		 } "maybe find apparent horizon(s) after this time step"
 	}
+if (run_at_CCTK_POSTINITIAL != 0)
+	{
+	schedule AHFinderDirect_update at CCTK_POSTINITIAL		\
+	before   AHFinderDirect_find_horizons
+		 {
+		 lang: C
+		 options: global
+		 } "update internal data structures from steerable parameters"
+	schedule AHFinderDirect_find_horizons at CCTK_POSTINITIAL
+		 {
+		 lang: C
+		 options: global
+		 } "maybe find apparent horizon(s) after this time step"
+	}
 
 #
 # announce horizon positions to other thorns
@@ -64,6 +78,15 @@ if (run_at_CCTK_POSTSTEP != 0)
 		 lang: C
 		 } "maybe announce horizon position(s) to other thorns"
 	}
+if (run_at_CCTK_POSTINITIAL != 0)
+	{
+	schedule AHFinderDirect_announce at CCTK_POSTINITIAL		\
+	before   DriftCorrect						\
+	after    AHFinderDirect_find_horizons
+		 {
+		 lang: C
+		 } "maybe announce horizon position(s) to other thorns"
+	}
 
 #
 # store horizon info in SphericalSurface variables
@@ -86,6 +109,15 @@ if (run_at_CCTK_POSTSTEP != 0)
 		 options: global
 		 } "maybe store horizon info in SphericalSurface variables"
 	}
+if (run_at_CCTK_POSTINITIAL != 0)
+	{
+	schedule AHFinderDirect_store_SS_info at CCTK_POSTINITIAL	\
+	after    AHFinderDirect_find_horizons
+		 {
+		 lang: C
+		 options: global
+		 } "maybe store horizon info in SphericalSurface variables"
+	}
 
 #
 # set excision mask
@@ -121,6 +153,14 @@ if (run_at_CCTK_POSTSTEP != 0)
 		 lang: C
 		 } "maybe set mask(s) based on apparent horizon position(s)"
 	}
+if (run_at_CCTK_POSTINITIAL != 0)
+	{
+	schedule AHFinderDirect_maybe_do_masks at CCTK_POSTINITIAL	\
+	after    AHFinderDirect_find_horizons
+		 {
+		 lang: C
+		 } "maybe set mask(s) based on apparent horizon position(s)"
+	}
 # if using mesh refinement, reset the mask after regridding
 # if not using mesh refinement, POSTREGRID is ignored so this is a no-op
 schedule AHFinderDirect_maybe_do_masks at POSTREGRID			\