change the code so AHFinderDirect can tell DriftCorrect about a

specified AH centroid -- this uses the new function-aliasing DriftCorrect
interface developed by Peter Diener and myself (with much assistance
from Ian Hawke)


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

1 files changed, 29 insertions, 0 deletions

diff --git a/schedule.ccl b/schedule.ccl
index 2abf1f6..59d6430 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -16,8 +16,37 @@ if (find_every != 0)
 		 triggers: always
 		 } "find apparent horizon(s) after this time step"
 
+	if (which_horizon_to_announce_centroid != 0)
+		{
+		schedule AHFinderDirect_announce at CCTK_ANALYSIS \
+                   before DriftCorrect \
+		   after AHFinderDirect_find_horizons
+			 {
+			 lang: C
+			 triggers: always
+			 } "announce horizon position(s) to other thorns"
+		}
+
 	if (set_mask != 0)
 		{
+		#
+		# *** KLUDGE ***
+		#
+		# We would really like to give this routine the
+		#	options:global
+		# attribute, so it only runs on those time levels where
+		# we've found (or at least tried to fnid) horizons.  But
+		# this doesn't work, because we need a GH -- and the other
+		# thorns' routines we call need valid grid variables in it
+		# -- in the announcing, and  options:global  routines don't
+		# have a "normal" GH, but instead a special one without
+		# usual grid variables defined in the usual manner.
+		#
+		# So, what to do?  The present kludge is to have this
+		# routine called on every grid, then have it explicitly
+		# test if it's at a horizon-finding time and on a base
+		# grid.
+		#
 		schedule AHFinderDirect_do_masks at CCTK_ANALYSIS \
 		   after AHFinderDirect_find_horizons
 			 {