diff options
Diffstat (limited to 'src/Dissipation_2_1.F90')
| -rw-r--r-- | src/Dissipation_2_1.F90 | 7 |
1 files changed, 6 insertions, 1 deletions
diff --git a/src/Dissipation_2_1.F90 b/src/Dissipation_2_1.F90 index 0b625f4..d6a3d5e 100644 --- a/src/Dissipation_2_1.F90 +++ b/src/Dissipation_2_1.F90 @@ -19,7 +19,6 @@ subroutine dissipation_2_1 (var, ni, nj, nk, bb, gsize, delta, epsilon, rhs) CCTK_REAL :: zero = 0.0 integer, parameter :: wp = kind(zero) - integer :: i, j, k CCTK_REAL, dimension(3,2) :: a CCTK_REAL :: idel @@ -55,11 +54,13 @@ subroutine dissipation_2_1 (var, ni, nj, nk, bb, gsize, delta, epsilon, rhs) ir = ni - 2 end if +!$omp parallel workshare rhs(il:ir,:,:) = rhs(il:ir,:,:) + & ( -2.0_wp * var(il:ir,:,:) + & ( var(il-1:ir-1,:,:) + & var(il+1:ir+1,:,:) ) ) * idel +!$omp end parallel workshare if ( zero_derivs_y == 0 ) then call set_coeff ( a ) @@ -91,10 +92,12 @@ subroutine dissipation_2_1 (var, ni, nj, nk, bb, gsize, delta, epsilon, rhs) jr = nj - 2 end if +!$omp parallel workshare rhs(:,jl:jr,:) = rhs(:,jl:jr,:) + & ( -2.0_wp * var(:,jl:jr,:) + & ( var(:,jl-1:jr-1,:) + & var(:,jl+1:jr+1,:) ) ) * idel +!$omp end parallel workshare end if if ( zero_derivs_z == 0 ) then @@ -128,10 +131,12 @@ subroutine dissipation_2_1 (var, ni, nj, nk, bb, gsize, delta, epsilon, rhs) kr = nk - 2 end if +!$omp parallel workshare rhs(:,:,kl:kr) = rhs(:,:,kl:kr) + & ( -2.0_wp * var(:,:,kl:kr) + & ( var(:,:,kl-1:kr-1) + & var(:,:,kl+1:kr+1) ) ) * idel +!$omp end parallel workshare end if contains |