diff options
Diffstat (limited to 'src/Dissipation_8_4.F90')
-rw-r--r-- | src/Dissipation_8_4.F90 | 25 |
1 files changed, 20 insertions, 5 deletions
diff --git a/src/Dissipation_8_4.F90 b/src/Dissipation_8_4.F90 index bc225af..ec98407 100644 --- a/src/Dissipation_8_4.F90 +++ b/src/Dissipation_8_4.F90 @@ -1,18 +1,21 @@ ! $Header$ #include "cctk.h" +#include "cctk_Parameters.h" -subroutine dissipation_8_4 (var, ni, nj, nk, bb, gsize, delta, epsdis, rhs) +subroutine dissipation_8_4 (var, ni, nj, nk, bb, gsize, delta, epsilon, rhs) implicit none + DECLARE_CCTK_PARAMETERS + integer :: ni, nj, nk CCTK_REAL, dimension(ni,nj,nk), intent(in) :: var CCTK_REAL, dimension(ni,nj,nk), intent(inout) :: rhs CCTK_INT, dimension(6), intent(in) :: bb CCTK_INT, dimension(3), intent(in) :: gsize CCTK_REAL, dimension(3), intent(in) :: delta - CCTK_REAL, intent(in) :: epsdis + CCTK_REAL, intent(in) :: epsilon CCTK_REAL :: zero = 0.0 integer, parameter :: wp = kind(zero) @@ -23,7 +26,11 @@ subroutine dissipation_8_4 (var, ni, nj, nk, bb, gsize, delta, epsdis, rhs) call set_coeff ( a ) - idel = epsdis / 256 + if ( scale_with_h > 0 ) then + idel = epsilon / ( 256 * delta(1) ) + else + idel = epsilon / 256 + end if if ( bb(1) == 0 ) then il = 1 + gsize(1) @@ -134,7 +141,11 @@ subroutine dissipation_8_4 (var, ni, nj, nk, bb, gsize, delta, epsdis, rhs) call set_coeff ( a ) - idel = epsdis / 256 + if ( scale_with_h > 0 ) then + idel = epsilon / ( 256 * delta(2) ) + else + idel = epsilon / 256 + end if if ( bb(3) == 0 ) then jl = 1 + gsize(2) @@ -245,7 +256,11 @@ subroutine dissipation_8_4 (var, ni, nj, nk, bb, gsize, delta, epsdis, rhs) call set_coeff ( a ) - idel = epsdis / 256 + if ( scale_with_h > 0 ) then + idel = epsilon / ( 256 * delta(3) ) + else + idel = epsilon / 256 + end if if ( bb(5) == 0 ) then kl = 1 + gsize(3) |