diff options
Diffstat (limited to 'src/apply_dissipation.F77')
-rw-r--r-- | src/apply_dissipation.F77 | 27 |
1 files changed, 27 insertions, 0 deletions
diff --git a/src/apply_dissipation.F77 b/src/apply_dissipation.F77 new file mode 100644 index 0000000..9363824 --- /dev/null +++ b/src/apply_dissipation.F77 @@ -0,0 +1,27 @@ +c $Header$ + +#include "cctk.h" + + subroutine apply_dissipation (var, rhs, ni, nj, nk, epsdis) + implicit none + + integer ni, nj, nk + CCTK_REAL var(ni,nj,nk), rhs(ni,nj,nk) + CCTK_REAL epsdis + + integer i, j, k + + do k = 3, nk-2 + do j = 3, nj-2 + do i = 3, ni-2 + + rhs(i,j,k) = rhs(i,j,k) - epsdis / 16 + $ * ( var(i-2,j,k) - 4*var(i-1,j,k) + 6*var(i,j,k) - 4*var(i+1,j,k) + var(i+2,j,k) + $ + var(i,j-2,k) - 4*var(i,j-1,k) + 6*var(i,j,k) - 4*var(i,j+1,k) + var(i,j+2,k) + $ + var(i,j,k-2) - 4*var(i,j,k-1) + 6*var(i,j,k) - 4*var(i,j,k+1) + var(i,j,k+2)) + + end do + end do + end do + + end |