Fixed a bug that made the re-initialization only work in full mode.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/EHFinder/trunk@186 2a26948c-0e4f-0410-aee8-f1d3e353619c

1 files changed, 33 insertions, 0 deletions

diff --git a/src/EHFinder_SetSym.F90 b/src/EHFinder_SetSym.F90
index 45fe9f5..79dda89 100644
--- a/src/EHFinder_SetSym.F90
+++ b/src/EHFinder_SetSym.F90
@@ -26,6 +26,12 @@ subroutine EHFinder_SetSym(CCTK_ARGUMENTS)
     call CCTK_WARN(1,'Failed to register symmetry for the level set function')
   end if
 
+! Set up symmetries for the temporary source function.
+  call SetCartSymGN ( ierr, cctkGH, sym, 'ehfinder::sftmp' )
+  if ( ierr .gt. 0 ) then
+    call CCTK_WARN(1,'Failed to register symmetry for the temp source function')
+  end if
+
 ! Set up symmetries for the mask function.
   call SetCartSymGN( ierr, cctkGH, sym, 'ehfinder::eh_mask' )
   if ( ierr .gt. 0 ) then
@@ -93,6 +99,33 @@ subroutine EHFinder_ApplySymF(CCTK_ARGUMENTS)
 end subroutine EHFinder_ApplySymF
 
 
+subroutine EHFinder_ApplySymFSFTMP(CCTK_ARGUMENTS)
+
+  use EHFinder_mod
+
+  implicit none
+
+  DECLARE_CCTK_ARGUMENTS
+  DECLARE_CCTK_PARAMETERS
+  DECLARE_CCTK_FUNCTIONS
+
+  character(len=80) :: warn_message
+
+  ierr = Boundary_SelectGroupForBC ( cctkGH, CCTK_ALL_FACES, 1, -1, &
+                                                 'ehfinder::f', 'None' )
+  if ( ierr /= 0 ) then
+    call CCTK_WARN ( 0, 'Could not select f for boundary condition' )
+  end if
+
+  ierr = Boundary_SelectGroupForBC ( cctkGH, CCTK_ALL_FACES, 1, -1, &
+                                                 'ehfinder::sftmp', 'None' )
+  if ( ierr /= 0 ) then
+    call CCTK_WARN ( 0, 'Could not select sftmp for boundary condition' )
+  end if
+
+  return
+end subroutine EHFinder_ApplySymFSFTMP
+
 subroutine EHFinder_ApplySymMask(CCTK_ARGUMENTS)
 
   use EHFinder_mod