diff options
author | diener <diener@f69c4107-0314-4c4f-9ad4-17e986b73f4a> | 2005-12-15 03:54:19 +0000 |
---|---|---|
committer | diener <diener@f69c4107-0314-4c4f-9ad4-17e986b73f4a> | 2005-12-15 03:54:19 +0000 |
commit | 4c7144d8e695475f86fec845e0730179a8959fc2 (patch) | |
tree | 5bc120016e74900734d0cd7967c1861f94a2fbad /src/CheckGridSizes.F90 | |
parent | 74de5846f7ddeb3bf181b3205f6539b4f0081705 (diff) |
Introduce a new kind of dissipation operators. Kreiss-Oliger type dissipation
operator suitably modified near the boundary. Since the dissipation operator
in the interior uses a stencil that is one point wider than the derivative
operator, an additional check on the number of ghostzones in multi processor
runs have been added. Only implemented for the 8-4 and optimized 6-5
operators.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@53 f69c4107-0314-4c4f-9ad4-17e986b73f4a
Diffstat (limited to 'src/CheckGridSizes.F90')
-rw-r--r-- | src/CheckGridSizes.F90 | 10 |
1 files changed, 10 insertions, 0 deletions
diff --git a/src/CheckGridSizes.F90 b/src/CheckGridSizes.F90 index 3240f6a..b6af08e 100644 --- a/src/CheckGridSizes.F90 +++ b/src/CheckGridSizes.F90 @@ -25,6 +25,16 @@ subroutine SBP_CheckGridSizes(CCTK_ARGUMENTS) 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 ( ( use_dissipation > 0 ) .and. & + ( CCTK_EQUALS(dissipation_type,"Kreiss-Oliger" ) ) .and. & + ( any (cctk_nghostzones < order / 2 + 1) ) .and. & + ( CCTK_nProcs(cctkGH) > 1 ) ) then + write(info_message,'(a22,i1,a13,i2,a32)') 'You need ghostsize >= ', & + order/2+1, ' to run with ', order+2, & + ' order Kreiss-Oliger dissipation' + call CCTK_WARN(0,info_message) + end if + 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, & |