1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
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
|