Changing the way the ahf_persits option works.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@75 89daf98e-ef62-4674-b946-b8ff9de2216c

1 files changed, 15 insertions, 5 deletions

diff --git a/schedule.ccl b/schedule.ccl
index 0d7018c..6521882 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -11,24 +11,34 @@ schedule AHFinder_SetSym at CCTK_BASEGRID
 } "Set symmetries for AHFinder grid functions"
 
 
-if (ahf_active && CCTK_Equals(ahf_mask,"off"))
+if (ahf_active && ahf_persists)
+{
+   STORAGE: ahfindergrid,ahfinderexp,ahfgradient,ahfinder_gauss,find3grid,ahfmask
+
+   schedule ahfinder at CCTK_POSTSTEP
+   {
+    LANG: Fortran
+   } "Call apparent horizon finder"
+}
+
+else if (ahf_active && CCTK_Equals(ahf_mask,"off"))
 {
    schedule ahfinder at CCTK_POSTSTEP
    {
     LANG: Fortran
-    STORAGE:ahfgradient,ahfinder_gauss,find3grid,ahfindergrid,ahfmask
+    STORAGE: ahfindergrid,ahfinderexp,ahfgradient,ahfinder_gauss,find3grid
    } "Call apparent horizon finder"
 }
 
 else if (ahf_active)
 {
-   STORAGE: ahfindergrid,ahfmask
+   STORAGE: ahfmask
 
    schedule ahfinder at CCTK_POSTSTEP
    {
     LANG: Fortran
-    STORAGE: ahfgradient,ahfinder_gauss,find3grid
-   } "Call apparent horizon finder with persisting grid"
+    STORAGE: ahfindergrid,ahfinderexp,ahfgradient,ahfinder_gauss,find3grid
+   } "Call apparent horizon finder with persisting mask"
 }