Added routines to find out how many surfaces are present. At present

this information is not used. Added beginning support for the spacemask.
Only routine to read in the emask. Nothing is done with it.


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

1 files changed, 28 insertions, 0 deletions

diff --git a/src/EHFinder_ReadData.F90 b/src/EHFinder_ReadData.F90
index 68516b6..75d633f 100644
--- a/src/EHFinder_ReadData.F90
+++ b/src/EHFinder_ReadData.F90
@@ -157,6 +157,34 @@ subroutine EHFinder_Read_Conformal(CCTK_ARGUMENTS)
 end subroutine EHFinder_Read_Conformal
 
 
+! This routine reads in the excision mask (from SpaceMask).
+subroutine EHFinder_Read_Mask(CCTK_ARGUMENTS)
+
+  use EHFinder_mod
+
+  implicit none
+
+  DECLARE_CCTK_PARAMETERS
+  DECLARE_CCTK_ARGUMENTS
+  DECLARE_CCTK_FUNCTIONS
+
+  character(len=128) :: in_files, in_vars
+  character(len=10) :: iteration_string
+  CCTK_INT :: nc, res
+
+  write(iteration_string,'(i10)') last_iteration_number - &
+                                  saved_iteration_every *cctk_iteration
+  iteration_string = adjustl(iteration_string)
+  nc = len_trim(iteration_string)
+
+  in_files = 'emask'
+  in_vars = 'spacemask::emask[cctk_iteration='//iteration_string(1:nc)//']'
+
+  call IOUtil_RecoverVarsFromDatafiles ( res, cctkGH, in_files, in_vars )
+
+end subroutine EHFinder_Read_Mask
+
+
 ! This routine reads in everything.
 subroutine EHFinder_ReadData(CCTK_ARGUMENTS)