Splitted the input up so that the metric, conformal factor, lapse and shift

are read in independently. Also scheduled these routines at CCTK_INITIAL
extending the ADMBase parameter for initial_data, initial_lapse and
initial_shift with "read from file".


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

1 files changed, 128 insertions, 0 deletions

diff --git a/src/EHFinder_ReadData.F90 b/src/EHFinder_ReadData.F90
index 40a9bec..fa31ac4 100644
--- a/src/EHFinder_ReadData.F90
+++ b/src/EHFinder_ReadData.F90
@@ -5,6 +5,134 @@
 #include "cctk_Parameters.h"
 #include "cctk_Arguments.h"
 
+subroutine EHFinder_Read_Metric(CCTK_ARGUMENTS)
+
+  use EHFinder_mod
+
+  implicit none
+
+  DECLARE_CCTK_PARAMETERS
+  DECLARE_CCTK_ARGUMENTS
+  DECLARE_CCTK_FUNCTIONS
+
+  character(len=512) :: in_files, in_vars
+  character(len=64), dimension(6)  :: var_names
+  character(len=10) :: iteration_string
+  CCTK_INT :: i, nc, ntot, res
+
+  write(iteration_string,'(i10)') last_iteration_number - cctk_iteration
+  iteration_string = adjustl(iteration_string)
+  nc = len_trim(iteration_string)
+
+  in_files = 'gxx gxy gxz gyy gyz gzz'
+  var_names(1) = 'admbase::gxx[cctk_iteration='//iteration_string(1:nc)//']'
+  var_names(2) = 'admbase::gxy[cctk_iteration='//iteration_string(1:nc)//']'
+  var_names(3) = 'admbase::gxz[cctk_iteration='//iteration_string(1:nc)//']'
+  var_names(4) = 'admbase::gyy[cctk_iteration='//iteration_string(1:nc)//']'
+  var_names(5) = 'admbase::gyz[cctk_iteration='//iteration_string(1:nc)//']'
+  var_names(6) = 'admbase::gzz[cctk_iteration='//iteration_string(1:nc)//']'
+
+  in_vars = ' '
+  ntot = 0
+  do i = 1, 6
+    nc = len_trim(var_names(i))
+    in_vars(ntot+1:ntot+1+nc+1) = var_names(i)(1:nc+1)
+    ntot = ntot + nc + 1
+  end do
+   
+  call IOUtil_RecoverVarsFromDatafiles ( res, cctkGH, in_files, in_vars )
+
+end subroutine EHFinder_Read_Metric
+
+
+subroutine EHFinder_Read_Lapse(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 - cctk_iteration
+  iteration_string = adjustl(iteration_string)
+  nc = len_trim(iteration_string)
+
+  in_files = 'alp'
+  in_vars = 'admbase::alp[cctk_iteration='//iteration_string(1:nc)//']'
+
+  call IOUtil_RecoverVarsFromDatafiles ( res, cctkGH, in_files, in_vars )
+
+end subroutine EHFinder_Read_Lapse
+
+
+subroutine EHFinder_Read_Shift(CCTK_ARGUMENTS)
+
+  use EHFinder_mod
+
+  implicit none
+
+  DECLARE_CCTK_PARAMETERS
+  DECLARE_CCTK_ARGUMENTS
+  DECLARE_CCTK_FUNCTIONS
+
+  character(len=256) :: in_files, in_vars
+  character(len=64), dimension(3)  :: var_names
+  character(len=10) :: iteration_string
+  CCTK_INT :: i, nc, ntot, res
+
+  write(iteration_string,'(i10)') last_iteration_number - cctk_iteration
+  iteration_string = adjustl(iteration_string)
+  nc = len_trim(iteration_string)
+
+  in_files = 'betax betay betaz'
+  var_names(1) = 'admbase::betax[cctk_iteration='//iteration_string(1:nc)//']'
+  var_names(2) = 'admbase::betay[cctk_iteration='//iteration_string(1:nc)//']'
+  var_names(3) = 'admbase::betaz[cctk_iteration='//iteration_string(1:nc)//']'
+  in_vars = ' '
+  ntot = 0
+  do i = 1, 3
+    nc = len_trim(var_names(i))
+    in_vars(ntot+1:ntot+1+nc+1) = var_names(i)(1:nc+1)
+    ntot = ntot + nc + 1
+  end do
+ 
+  call IOUtil_RecoverVarsFromDatafiles ( res, cctkGH, in_files, in_vars )
+
+end subroutine EHFinder_Read_Shift
+
+
+subroutine EHFinder_Read_Conformal(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 - cctk_iteration
+  iteration_string = adjustl(iteration_string)
+  nc = len_trim(iteration_string)
+
+  in_files = 'psi'
+  in_vars = 'staticconformal::psi[cctk_iteration='//iteration_string(1:nc)//']'
+
+  call IOUtil_RecoverVarsFromDatafiles ( res, cctkGH, in_files, in_vars )
+
+end subroutine EHFinder_Read_Conformal
+
+
 subroutine EHFinder_ReadData(CCTK_ARGUMENTS)
 
   use EHFinder_mod