From 0c7034c5e863d96b6182b41f77d7eca8b25440e9 Mon Sep 17 00:00:00 2001 From: Grzegorz Kowal Date: Fri, 28 Apr 2017 13:26:10 -0300 Subject: [PATCH] INTERPOLATION: Limit all cells in MLP limiting. Signed-off-by: Grzegorz Kowal --- src/interpolations.F90 | 40 ++++++++++++++++++++-------------------- 1 file changed, 20 insertions(+), 20 deletions(-) diff --git a/src/interpolations.F90 b/src/interpolations.F90 index df099b7..ec2d4fa 100644 --- a/src/interpolations.F90 +++ b/src/interpolations.F90 @@ -1200,32 +1200,32 @@ module interpolations dq(m) = sign(ap(m), dq(m)) end do + end if + ! update the interpolated states ! - dql(1) = qi(i ,j,k,1,1) - q(i,j,k) - dqr(1) = qi(im1,j,k,2,1) - q(i,j,k) - if (max(abs(dql(1)), abs(dqr(1))) > abs(dq(1))) then - qi(i ,j,k,1,1) = q(i,j,k) + dq(1) - qi(im1,j,k,2,1) = q(i,j,k) - dq(1) - end if + dql(1) = qi(i ,j,k,1,1) - q(i,j,k) + dqr(1) = qi(im1,j,k,2,1) - q(i,j,k) + if (max(abs(dql(1)), abs(dqr(1))) > abs(dq(1))) then + qi(i ,j,k,1,1) = q(i,j,k) + dq(1) + qi(im1,j,k,2,1) = q(i,j,k) - dq(1) + end if - dql(2) = qi(i,j ,k,1,2) - q(i,j,k) - dqr(2) = qi(i,jm1,k,2,2) - q(i,j,k) - if (max(abs(dql(2)), abs(dqr(2))) > abs(dq(2))) then - qi(i,j ,k,1,2) = q(i,j,k) + dq(2) - qi(i,jm1,k,2,2) = q(i,j,k) - dq(2) - end if + dql(2) = qi(i,j ,k,1,2) - q(i,j,k) + dqr(2) = qi(i,jm1,k,2,2) - q(i,j,k) + if (max(abs(dql(2)), abs(dqr(2))) > abs(dq(2))) then + qi(i,j ,k,1,2) = q(i,j,k) + dq(2) + qi(i,jm1,k,2,2) = q(i,j,k) - dq(2) + end if #if NDIMS == 3 - dql(3) = qi(i,j,k ,1,3) - q(i,j,k) - dqr(3) = qi(i,j,km1,2,3) - q(i,j,k) - if (max(abs(dql(3)), abs(dqr(3))) > abs(dq(3))) then - qi(i,j,k ,1,3) = q(i,j,k) + dq(3) - qi(i,j,km1,2,3) = q(i,j,k) - dq(3) - end if -#endif /* NDIMS == 3 */ - + dql(3) = qi(i,j,k ,1,3) - q(i,j,k) + dqr(3) = qi(i,j,km1,2,3) - q(i,j,k) + if (max(abs(dql(3)), abs(dqr(3))) > abs(dq(3))) then + qi(i,j,k ,1,3) = q(i,j,k) + dq(3) + qi(i,j,km1,2,3) = q(i,j,k) - dq(3) end if +#endif /* NDIMS == 3 */ end do ! i = ibl, ieu end do ! j = jbl, jeu