diff options
Diffstat (limited to 'src/apply_dissipation.F77')
| -rw-r--r-- | src/apply_dissipation.F77 | 8 |
1 files changed, 6 insertions, 2 deletions
diff --git a/src/apply_dissipation.F77 b/src/apply_dissipation.F77 index d67a02e..b4e3545 100644 --- a/src/apply_dissipation.F77 +++ b/src/apply_dissipation.F77 @@ -2,12 +2,12 @@ c $Header$ #include "cctk.h" - subroutine apply_dissipation (var, rhs, ni, nj, nk, dx, dt, order, epsdis) + subroutine apply_dissipation (var, rhs, ni, nj, nk, dx, order, epsdis) implicit none integer ni, nj, nk CCTK_REAL var(ni,nj,nk), rhs(ni,nj,nk) - CCTK_REAL dx(3), dt + CCTK_REAL dx(3) CCTK_INT order CCTK_REAL epsdis(ni,nj,nk) @@ -15,6 +15,7 @@ c $Header$ if (order .eq. 1) then +c$omp parallel do do k = 2, nk-1 do j = 2, nj-1 do i = 2, ni-1 @@ -29,6 +30,7 @@ c $Header$ else if (order .eq. 3) then +c$omp parallel do do k = 3, nk-2 do j = 3, nj-2 do i = 3, ni-2 @@ -44,6 +46,7 @@ c $Header$ else if (order .eq. 5) then +c$omp parallel do do k = 4, nk-3 do j = 4, nj-3 do i = 4, ni-3 @@ -59,6 +62,7 @@ c $Header$ else if (order .eq. 7) then +c$omp parallel do do k = 5, nk-4 do j = 5, nj-4 do i = 5, ni-4 |