diff options
| -rw-r--r-- | schedule.ccl | 5 | ||||
| -rw-r--r-- | src/CheckGridSizes.F90 | 54 | ||||
| -rw-r--r-- | src/make.code.defn | 2 |
3 files changed, 60 insertions, 1 deletions
diff --git a/schedule.ccl b/schedule.ccl index e1bc52a..7575d69 100644 --- a/schedule.ccl +++ b/schedule.ccl @@ -8,3 +8,8 @@ SCHEDULE SBP_SetNormMask AT initial { LANG: Fortran } "Setup the mask for the calulcation of the norm" + +SCHEDULE SBP_CheckGridSizes AT postinitial +{ + LANG: Fortran +} "Check grid sizes and ghost zones" diff --git a/src/CheckGridSizes.F90 b/src/CheckGridSizes.F90 new file mode 100644 index 0000000..3240f6a --- /dev/null +++ b/src/CheckGridSizes.F90 @@ -0,0 +1,54 @@ +! Check if there are enough grid cells and ghost zones to compute the finite +! differences. +! $Header$ + +#include "cctk.h" +#include "cctk_Parameters.h" +#include "cctk_Arguments.h" +#include "cctk_Functions.h" + +subroutine SBP_CheckGridSizes(CCTK_ARGUMENTS) + + implicit none + + DECLARE_CCTK_PARAMETERS + DECLARE_CCTK_ARGUMENTS + DECLARE_CCTK_FUNCTIONS + + integer :: i + integer, dimension(3) :: bc + integer, dimension(3,3) :: np + character(len=128) :: info_message + character(len=1), dimension(3) :: dir = (/ 'x', 'y', 'z' /) + + bc(1) = cctk_bbox(1)+cctk_bbox(2)+1 ! 0 processor boundaries -> bc = 3 + bc(2) = cctk_bbox(3)+cctk_bbox(4)+1 ! 1 processor boundary -> bc = 2 + bc(3) = cctk_bbox(5)+cctk_bbox(6)+1 ! 2 processor boundaries -> bc = 1 + + if ( any (cctk_nghostzones < order / 2) .and. CCTK_nProcs(cctkGH) > 1 ) then + write(info_message,'(a22,i1,a13,i1,a25)') 'You need ghostsize >= ', & + order/2, ' to run with ', order, & + ' order finite differences' + call CCTK_WARN(0,info_message) + end if + + if ( CCTK_EQUALS(norm_type,"Diagonal") ) then + np(1,:) = 2*cctk_nghostzones + np(2,:) = 3*order/2 + np(3,:) = 2*order + else + np(1,:) = 2*cctk_nghostzones + np(2,:) = 3*order/2 + 1 + np(3,:) = 2*(order+1) + end if + + do i = 1, 3 + if ( cctk_lsh(i) < np(bc(i),i) ) then + write(info_message,'(a27,a1,a35,i1,a25)') & + 'Your local gridsize in the ', dir(i), & + '-direction is to small to run with ', order, & + ' order finite differences' + call CCTK_WARN(0,info_message) + end if + end do +end subroutine SBP_CheckGridSizes diff --git a/src/make.code.defn b/src/make.code.defn index 893ee3f..5810dc0 100644 --- a/src/make.code.defn +++ b/src/make.code.defn @@ -2,7 +2,7 @@ # $Header$ # Source files in this directory -SRCS = Derivatives_2_1.F90 Derivatives_4_2.F90 Derivatives_6_3.F90 Derivatives_8_4.F90 call_derivs.c Derivatives_4_3.F90 call_derivs_name.c Get_Coeff.F90 call_norm_mask.c set_norm_mask.F90 +SRCS = Derivatives_2_1.F90 Derivatives_4_2.F90 Derivatives_6_3.F90 Derivatives_8_4.F90 call_derivs.c Derivatives_4_3.F90 call_derivs_name.c Get_Coeff.F90 call_norm_mask.c set_norm_mask.F90 CheckGridSizes.F90 # Subdirectories containing source files SUBDIRS = |