diff options
Diffstat (limited to 'src/Dissipation_6_5_min_err_coeff.F90')
| -rw-r--r-- | src/Dissipation_6_5_min_err_coeff.F90 | 22 |
1 files changed, 17 insertions, 5 deletions
diff --git a/src/Dissipation_6_5_min_err_coeff.F90 b/src/Dissipation_6_5_min_err_coeff.F90 index a88d3ad..312ecef 100644 --- a/src/Dissipation_6_5_min_err_coeff.F90 +++ b/src/Dissipation_6_5_min_err_coeff.F90 @@ -6,7 +6,7 @@ subroutine dissipation_6_5_opt (var, lsh, gsh, lbnd, bb, gsize, & - delta, epsdisl, dfl, rhs) + delta, epsilon, dfl, rhs) implicit none @@ -20,7 +20,7 @@ subroutine dissipation_6_5_opt (var, lsh, gsh, lbnd, bb, gsize, & 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) :: epsdisl, dfl + CCTK_REAL, intent(in) :: epsilon, dfl CCTK_REAL :: zero = 0.0 integer, parameter :: wp = kind(zero) @@ -79,7 +79,11 @@ subroutine dissipation_6_5_opt (var, lsh, gsh, lbnd, bb, gsize, & end if - idel = epsdisl / 64.0_wp + if ( scale_with_h > 0 ) then + idel = epsilon / ( 64 * delta(1) ) + else + idel = epsilon / 64 + end if if ( bb(1) == 0 ) then @@ -275,7 +279,11 @@ subroutine dissipation_6_5_opt (var, lsh, gsh, lbnd, bb, gsize, & end if - idel = epsdisl / 64.0_wp + if ( scale_with_h > 0 ) then + idel = epsilon / ( 64 * delta(2) ) + else + idel = epsilon / 64 + end if if ( bb(3) == 0 ) then @@ -471,7 +479,11 @@ subroutine dissipation_6_5_opt (var, lsh, gsh, lbnd, bb, gsize, & end if - idel = epsdisl / 64.0_wp + if ( scale_with_h > 0 ) then + idel = epsilon / ( 64 * delta(3) ) + else + idel = epsilon / 64 + end if if ( bb(5) == 0 ) then |