diff options
| author | schnetter <schnetter@850bcc8b-0e4f-0410-8c26-8d28fbf1eda9> | 2005-06-14 20:26:02 +0000 |
|---|---|---|
| committer | schnetter <schnetter@850bcc8b-0e4f-0410-8c26-8d28fbf1eda9> | 2005-06-14 20:26:02 +0000 |
| commit | 006083638b7fcad16e5e2c6eed71f72f2c1a2c9c (patch) | |
| tree | 5b65477b75adf2b24b5d21473ef1a730c06f288a /src/apply_dissipation.F77 | |
| parent | efcd06a3138fbffba457960d8992111983695e82 (diff) | |
Add Frank Herrmann's extra dissipation switches.
Two new parameters, Dissipation::extra_dissipation_in_horizons and
Dissipation::extra_dissipation_at_outerbound, increase the disspation
coefficient inside horizons and near the outer boundary.
git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/Dissipation/trunk@10 850bcc8b-0e4f-0410-8c26-8d28fbf1eda9
Diffstat (limited to 'src/apply_dissipation.F77')
| -rw-r--r-- | src/apply_dissipation.F77 | 7 |
1 files changed, 3 insertions, 4 deletions
diff --git a/src/apply_dissipation.F77 b/src/apply_dissipation.F77 index 74cf168..a246774 100644 --- a/src/apply_dissipation.F77 +++ b/src/apply_dissipation.F77 @@ -9,7 +9,7 @@ c $Header$ CCTK_REAL var(ni,nj,nk), rhs(ni,nj,nk) CCTK_REAL dx(3), dt CCTK_INT order - CCTK_REAL epsdis + CCTK_REAL epsdis(ni,nj,nk) integer i, j, k @@ -18,8 +18,7 @@ c $Header$ do k = 2, nk-1 do j = 2, nj-1 do i = 2, ni-1 - - rhs(i,j,k) = rhs(i,j,k) + epsdis + rhs(i,j,k) = rhs(i,j,k) + epsdis(i,j,k) $ * (+ (var(i-1,j,k) - 2*var(i,j,k) + var(i+1,j,k)) / dx(1) $ + (var(i,j-1,k) - 2*var(i,j,k) + var(i,j+1,k)) / dx(2) $ + (var(i,j,k-1) - 2*var(i,j,k) + var(i,j,k+1)) / dx(3)) @@ -34,7 +33,7 @@ c $Header$ do j = 3, nj-2 do i = 3, ni-2 - rhs(i,j,k) = rhs(i,j,k) - epsdis / 16 + rhs(i,j,k) = rhs(i,j,k) - epsdis(i,j,k) / 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)) / dx(1) $ + (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)) / dx(2) $ + (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)) / dx(3)) |