Index: param.ccl =================================================================== --- param.ccl (revision 1548) +++ param.ccl (working copy) @@ -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?" { Index: schedule.ccl =================================================================== --- schedule.ccl (revision 1548) +++ schedule.ccl (working copy) @@ -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 \