diff options
| -rw-r--r-- | src/Dissipation_6_5_min_err_coeff.F90 | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/src/Dissipation_6_5_min_err_coeff.F90 b/src/Dissipation_6_5_min_err_coeff.F90 index 312ecef..d4e15a4 100644 --- a/src/Dissipation_6_5_min_err_coeff.F90 +++ b/src/Dissipation_6_5_min_err_coeff.F90 @@ -265,7 +265,7 @@ subroutine dissipation_6_5_opt (var, lsh, gsh, lbnd, bb, gsize, & a(1:jr,1:center) = tmp(1:jr,1:center) - jl = center + 1 - gsize(1) + jl = center + 1 - gsize(2) d = zero call set_dmatrix_half ( d(jl:nj,jl:nj), 2 ) tmp(jl:nj,jl:nj) = -transpose ( & @@ -285,11 +285,11 @@ subroutine dissipation_6_5_opt (var, lsh, gsh, lbnd, bb, gsize, & idel = epsilon / 64 end if - if ( bb(3) == 0 ) then + if ( bb(3) == 0 ) then - jl = 1 + gsize(2) + jl = 1 + gsize(2) - else + else rhs(:,1,:) = rhs(:,1,:) + & ( a(1,1) * var(:,1,:) + a(2,1) * var(:,2,:) + & |