diff options
| author | diener <diener@2a26948c-0e4f-0410-aee8-f1d3e353619c> | 2003-06-13 09:06:12 +0000 |
|---|---|---|
| committer | diener <diener@2a26948c-0e4f-0410-aee8-f1d3e353619c> | 2003-06-13 09:06:12 +0000 |
| commit | 270d4606bcdd581c5251a51179f96c8c6b01bd4c (patch) | |
| tree | d4e2dcd41e8dcfc0511a5393f036612cc22f5a4c /src/EHFinder_ReadData.F90 | |
| parent | f7b5353dcdc301c459f61831d87857dfd67bbd00 (diff) | |
Added a parameter read_conformal_factor_once to support reading the
conformal factor only once (since it's not changed during the evolution).
The user now only has to output staticconformal::psi at iteration 0, thereby
saving some diskspace.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/EHFinder/trunk@117 2a26948c-0e4f-0410-aee8-f1d3e353619c
Diffstat (limited to 'src/EHFinder_ReadData.F90')
| -rw-r--r-- | src/EHFinder_ReadData.F90 | 18 |
1 files changed, 11 insertions, 7 deletions
diff --git a/src/EHFinder_ReadData.F90 b/src/EHFinder_ReadData.F90 index 75d633f..08e3ce8 100644 --- a/src/EHFinder_ReadData.F90 +++ b/src/EHFinder_ReadData.F90 @@ -27,7 +27,7 @@ subroutine EHFinder_Read_Metric(CCTK_ARGUMENTS) ! saved) and the current iteration and save it in a string variable. write(iteration_string,'(i10)') last_iteration_number - & - saved_iteration_every *cctk_iteration + saved_iteration_every * cctk_iteration ! Trim the string variable. iteration_string = adjustl(iteration_string) @@ -79,7 +79,7 @@ subroutine EHFinder_Read_Lapse(CCTK_ARGUMENTS) CCTK_INT :: nc, res write(iteration_string,'(i10)') last_iteration_number - & - saved_iteration_every *cctk_iteration + saved_iteration_every * cctk_iteration iteration_string = adjustl(iteration_string) nc = len_trim(iteration_string) @@ -108,7 +108,7 @@ subroutine EHFinder_Read_Shift(CCTK_ARGUMENTS) CCTK_INT :: i, nc, ntot, res write(iteration_string,'(i10)') last_iteration_number - & - saved_iteration_every *cctk_iteration + saved_iteration_every * cctk_iteration iteration_string = adjustl(iteration_string) nc = len_trim(iteration_string) @@ -144,8 +144,12 @@ subroutine EHFinder_Read_Conformal(CCTK_ARGUMENTS) character(len=10) :: iteration_string CCTK_INT :: nc, res - write(iteration_string,'(i10)') last_iteration_number - & - saved_iteration_every *cctk_iteration + if ( read_conformal_factor_once .gt. 0 ) then + write(iteration_string,'(i10)') 0 + else + write(iteration_string,'(i10)') last_iteration_number - & + saved_iteration_every * cctk_iteration + end if iteration_string = adjustl(iteration_string) nc = len_trim(iteration_string) @@ -173,7 +177,7 @@ subroutine EHFinder_Read_Mask(CCTK_ARGUMENTS) CCTK_INT :: nc, res write(iteration_string,'(i10)') last_iteration_number - & - saved_iteration_every *cctk_iteration + saved_iteration_every * cctk_iteration iteration_string = adjustl(iteration_string) nc = len_trim(iteration_string) @@ -202,7 +206,7 @@ subroutine EHFinder_ReadData(CCTK_ARGUMENTS) CCTK_INT :: i, nc, ntot, res write(iteration_string,'(i10)') last_iteration_number - & - saved_iteration_every *cctk_iteration + saved_iteration_every * cctk_iteration iteration_string = adjustl(iteration_string) nc = len_trim(iteration_string) |