From d95549fda707bce381713748e0e5929b60a904c6 Mon Sep 17 00:00:00 2001 From: Grzegorz Kowal Date: Mon, 24 Oct 2022 15:55:46 -0300 Subject: [PATCH] USER_PROBLEMS: Fix user boundary conditions. Signed-off-by: Grzegorz Kowal --- sources/user_problem.F90 | 18 +++++++++++------- 1 file changed, 11 insertions(+), 7 deletions(-) diff --git a/sources/user_problem.F90 b/sources/user_problem.F90 index 939b075..b43dc01 100644 --- a/sources/user_problem.F90 +++ b/sources/user_problem.F90 @@ -413,24 +413,28 @@ module user_problem real(kind=8), dimension(:) , intent(in) :: x, y, z real(kind=8), dimension(:,:,:,:), intent(inout) :: qn - integer :: i, j + integer :: i, j, jr !------------------------------------------------------------------------------- ! if (js == 1) then do j = nbl, 1, -1 + jr = nb + (nbl - j) do i = il, iu - qn( : ,i,j,:) = qn( : ,i,nb,:) - qn(ivy,i,j,:) = 0.0d+00 - qn(iby,i,j,:) = 0.0d+00 + qn( : ,i,j,:) = qn( : ,i,jr,:) + + qn(ivy,i,j,:) = - qn(ivy,i,jr,:) + qn(iby,i,j,:) = - qn(iby,i,jr,:) end do end do else do j = neu, nn + jr = ne - (j - neu) do i = il, iu - qn( : ,i,j,:) = qn( : ,i,ne,:) - qn(ivy,i,j,:) = 0.0d+00 - qn(iby,i,j,:) = 0.0d+00 + qn( : ,i,j,:) = qn( : ,i,jr,:) + + qn(ivy,i,j,:) = - qn(ivy,i,jr,:) + qn(iby,i,j,:) = - qn(iby,i,jr,:) end do end do end if