diff options
| -rw-r--r-- | src/EHFinder_Check.F90 | 4 | ||||
| -rw-r--r-- | src/EHFinder_Init.F90 | 8 | ||||
| -rw-r--r-- | src/EHFinder_Sources.F90 | 10 |
3 files changed, 11 insertions, 11 deletions
diff --git a/src/EHFinder_Check.F90 b/src/EHFinder_Check.F90 index 81d3413..0d391fe 100644 --- a/src/EHFinder_Check.F90 +++ b/src/EHFinder_Check.F90 @@ -83,7 +83,9 @@ subroutine EHFinder_ReInitialize_Check(CCTK_ARGUMENTS) ! If the local extremum indicate and imminent pinch off undo the ! re-initialization if this is requested by the user. if ( re_init_undo .gt. 0 ) then - if ( fmin_extremum .gt. -two * minval(cctk_delta_space) ) then + if ( fmin_extremum .gt. -two * min ( CCTK_DELTA_SPACE(1), & + CCTK_DELTA_SPACE(2), & + CCTK_DELTA_SPACE(3) ) ) then write(info_message,'(a45,i3,a33)') & 'Re-initialization of level set ', l, & ' undone due to imminent pinch-off' diff --git a/src/EHFinder_Init.F90 b/src/EHFinder_Init.F90 index 478839f..10888d8 100644 --- a/src/EHFinder_Init.F90 +++ b/src/EHFinder_Init.F90 @@ -278,11 +278,11 @@ subroutine EHFinder_Init(CCTK_ARGUMENTS) DECLARE_CCTK_ARGUMENTS DECLARE_CCTK_FUNCTIONS -! Set up the value used in interiour inactive cells. - ex_value = - ( one + shell_width ) * maxval(cctk_delta_space) - ! Find the maximal grid spacing. - delta = maxval ( cctk_delta_space ) + delta = max ( CCTK_DELTA_SPACE(1), CCTK_DELTA_SPACE(2), CCTK_DELTA_SPACE(3) ) + +! Set up the value used in interiour inactive cells. + ex_value = - ( one + shell_width ) * delta ! Get handles for various reduction operations. call CCTK_ReductionArrayHandle ( max_handle, 'maximum' ) 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 |