[Commits] [svn:einsteintoolkit] AHFinderDirect/trunk/ (Rev. 1550)

[knarf at cct.lsu.edu]

User: knarf
Date: 2010/08/10 09:10 AM

Modified:
 /trunk/
  param.ccl, schedule.ccl

Log:
 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).

File Changes:

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

File [modified]: param.ccl
Delta lines: +3 -0
===================================================================
--- trunk/param.ccl	2010-07-12 14:49:57 UTC (rev 1549)
+++ trunk/param.ccl	2010-08-10 14:10:23 UTC (rev 1550)
@@ -114,6 +114,9 @@
 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?"
 {

File [modified]: schedule.ccl
Delta lines: +50 -0
===================================================================
--- trunk/schedule.ccl	2010-07-12 14:49:57 UTC (rev 1549)
+++ trunk/schedule.ccl	2010-08-10 14:10:23 UTC (rev 1550)
@@ -50,6 +50,14 @@
 		 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 @@
 		 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 \