diff options
| author | Anton Khirnov <anton@khirnov.net> | 2016-02-26 17:13:59 +0100 |
|---|---|---|
| committer | Anton Khirnov <anton@khirnov.net> | 2016-02-26 17:13:59 +0100 |
| commit | 8c0efc3ebf2e8f5e8f55b52e7a977ffb39cd7520 (patch) | |
| tree | 7b74958d40dcecc50dc551eeba2890c85f57e2b2 | |
| parent | 653b6bf82b056616db102b7eb1486af4cf7afc8b (diff) | |
cartoon: disable the y part of the dissipation operator
| -rw-r--r-- | src/apply_dissipation.F77 | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/src/apply_dissipation.F77 b/src/apply_dissipation.F77 index a54b80c..222d10c 100644 --- a/src/apply_dissipation.F77 +++ b/src/apply_dissipation.F77 @@ -25,8 +25,8 @@ c$omp parallel do private(i,j,k) do i = 2, ni-1 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)) * idx(1) - $ + (var(i,j-1,k) - 2*var(i,j,k) + var(i,j+1,k)) * idx(2) $ + (var(i,j,k-1) - 2*var(i,j,k) + var(i,j,k+1)) * idx(3)) +c $ + (var(i,j-1,k) - 2*var(i,j,k) + var(i,j+1,k)) * idx(2) end do end do @@ -41,8 +41,8 @@ c$omp parallel do private(i,j,k) 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)) * idx(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)) * idx(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)) * idx(3)) +c $ + (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)) * idx(2) end do end do @@ -57,8 +57,8 @@ c$omp parallel do private(i,j,k) rhs(i,j,k) = rhs(i,j,k) + epsdis(i,j,k) / 64 $ * (+ (var(i-3,j,k) - 6*var(i-2,j,k) + 15*var(i-1,j,k) - 20*var(i,j,k) + 15*var(i+1,j,k) - 6*var(i+2,j,k) + var(i+3,j,k)) * idx(1) - $ + (var(i,j-3,k) - 6*var(i,j-2,k) + 15*var(i,j-1,k) - 20*var(i,j,k) + 15*var(i,j+1,k) - 6*var(i,j+2,k) + var(i,j+3,k) ) * idx(2) $ + (var(i,j,k-3) - 6*var(i,j,k-2) + 15*var(i,j,k-1) - 20*var(i,j,k) + 15*var(i,j,k+1) - 6*var(i,j,k+2) + var(i,j,k+3) ) * idx(3)) +c $ + (var(i,j-3,k) - 6*var(i,j-2,k) + 15*var(i,j-1,k) - 20*var(i,j,k) + 15*var(i,j+1,k) - 6*var(i,j+2,k) + var(i,j+3,k) ) * idx(2) end do end do @@ -73,8 +73,8 @@ c$omp parallel do private(i,j,k) rhs(i,j,k) = rhs(i,j,k) - epsdis(i,j,k) / 256 $ * (+ (var(i-4,j,k) - 8*var(i-3,j,k) + 28*var(i-2,j,k) - 56*var(i-1,j,k) + 70*var(i,j,k) - 56*var(i+1,j,k) + 28*var(i+2,j,k) - 8*var(i+3,j,k) + var(i+4,j,k)) * idx(1) - $ + (var(i,j-4,k) - 8*var(i,j-3,k) + 28*var(i,j-2,k) - 56*var(i,j-1,k) + 70*var(i,j,k) - 56*var(i,j+1,k) + 28*var(i,j+2,k) - 8*var(i,j+3,k) + var(i,j+4,k)) * idx(2) $ + (var(i,j,k-4) - 8*var(i,j,k-3) + 28*var(i,j,k-2) - 56*var(i,j,k-1) + 70*var(i,j,k) - 56*var(i,j,k+1) + 28*var(i,j,k+2) - 8*var(i,j,k+3) + var(i,j,k+4)) * idx(3)) +c $ + (var(i,j-4,k) - 8*var(i,j-3,k) + 28*var(i,j-2,k) - 56*var(i,j-1,k) + 70*var(i,j,k) - 56*var(i,j+1,k) + 28*var(i,j+2,k) - 8*var(i,j+3,k) + var(i,j+4,k)) * idx(2) end do end do @@ -89,8 +89,8 @@ c$omp parallel do private(i,j,k) rhs(i,j,k) = rhs(i,j,k) + epsdis(i,j,k) / 1024 $ * (+ (var(i-5,j,k) - 10*var(i-4,j,k) + 45*var(i-3,j,k) - 120*var(i-2,j,k) + 210*var(i-1,j,k) - 252*var(i,j,k) + 210*var(i+1,j,k) - 120*var(i+2,j,k) + 45*var(i+3,j,k) - 10*var(i+4,j,k) + var(i+5,j,k)) * idx(1) - $ + (var(i,j-5,k) - 10*var(i,j-4,k) + 45*var(i,j-3,k) - 120*var(i,j-2,k) + 210*var(i,j-1,k) - 252*var(i,j,k) + 210*var(i,j+1,k) - 120*var(i,j+2,k) + 45*var(i,j+3,k) - 10*var(i,j+4,k) + var(i,j+5,k)) * idx(2) $ + (var(i,j,k-5) - 10*var(i,j,k-4) + 45*var(i,j,k-3) - 120*var(i,j,k-2) + 210*var(i,j,k-1) - 252*var(i,j,k) + 210*var(i,j,k+1) - 120*var(i,j,k+2) + 45*var(i,j,k+3) - 10*var(i,j,k+4) + var(i,j,k+5)) * idx(3)) +c $ + (var(i,j-5,k) - 10*var(i,j-4,k) + 45*var(i,j-3,k) - 120*var(i,j-2,k) + 210*var(i,j-1,k) - 252*var(i,j,k) + 210*var(i,j+1,k) - 120*var(i,j+2,k) + 45*var(i,j+3,k) - 10*var(i,j+4,k) + var(i,j+5,k)) * idx(2) end do end do |