diff options
Diffstat (limited to 'src/EHFinder_Sources.F90')
| -rw-r--r-- | src/EHFinder_Sources.F90 | 10 |
1 files changed, 4 insertions, 6 deletions
diff --git a/src/EHFinder_Sources.F90 b/src/EHFinder_Sources.F90 index 86a77ed..c59bf39 100644 --- a/src/EHFinder_Sources.F90 +++ b/src/EHFinder_Sources.F90 @@ -17,7 +17,7 @@ subroutine EHFinder_Sources(CCTK_ARGUMENTS) DECLARE_CCTK_FUNCTIONS CCTK_INT :: i, j, k, l - CCTK_REAL :: idx, idy, idz, mdelta + CCTK_REAL :: idx, idy, idz CCTK_REAL :: a, b, c CCTK_REAL :: gxxc, gxyc, gxzc, gyyc, gyzc, gzzc, psito4 CCTK_REAL :: idetg, alp2, tmp1, tmp2, tmp3 @@ -30,11 +30,9 @@ subroutine EHFinder_Sources(CCTK_ARGUMENTS) #include "include/physical_part.h" ! calculate 1/(2*delta) in each direction - idx = half / cctk_delta_space(1) - idy = half / cctk_delta_space(2) - idz = half / cctk_delta_space(3) - - mdelta = maxval ( cctk_delta_space ) + idx = half / CCTK_DELTA_SPACE(1) + idy = half / CCTK_DELTA_SPACE(2) + idz = half / CCTK_DELTA_SPACE(3) ! Set the sign depending on the surface direction. if ( CCTK_EQUALS ( surface_direction, 'outward' ) ) ssign = one |