From 645af56ed21e9b59234569cf56aa9e9f81fe69cc Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Franck=20P=C3=A9rignon?= <franck.perignon@imag.fr> Date: Thu, 28 Mar 2013 10:42:50 +0000 Subject: [PATCH] fix a small bug in fortran scales. Comp/tests ok. --- .../particles/advec_remesh_lambda.f90 | 14 +++++++++----- 1 file changed, 9 insertions(+), 5 deletions(-) diff --git a/HySoP/src/scalesInterface/particles/advec_remesh_lambda.f90 b/HySoP/src/scalesInterface/particles/advec_remesh_lambda.f90 index 066394072..8f68902e2 100644 --- a/HySoP/src/scalesInterface/particles/advec_remesh_lambda.f90 +++ b/HySoP/src/scalesInterface/particles/advec_remesh_lambda.f90 @@ -331,11 +331,13 @@ subroutine AC_remesh_lambda2corrected_array(direction, p_pos_adim, scal1d, bl_ty !end if ! XXX Debug - end if (bl_type(bl_ind)) then - ! tagged, the first particle belong to a centered block and the last to left block. - call AC_remesh_tag_CL(direction, p_pos_adim(p_ind), scal1D(p_ind), p_pos_adim(p_ind+1), scal1D(p_ind+1), remesh_buffer) + ! tagged, the first particle belong to a centered block and the last to left block. + call AC_remesh_tag_CL(direction, p_pos_adim(p_ind), scal1D(p_ind), p_pos_adim(p_ind+1), & + scal1D(p_ind+1), remesh_buffer) else ! tagged, the first particle belong to a left block and the last to centered block. - call AC_remesh_tag_LC(direction, p_pos_adim(p_ind), scal1D(p_ind), p_pos_adim(p_ind+1), scal1D(p_ind+1), remesh_buffer) + call AC_remesh_tag_LC(direction, p_pos_adim(p_ind), scal1D(p_ind), p_pos_adim(p_ind+1), & + scal1D(p_ind+1), remesh_buffer) end if else ! First particle @@ -629,9 +631,11 @@ subroutine AC_remesh_lambda2limited_array(direction, p_pos_adim, scal1d, bl_type end if else ! First particle - call AC_remesh_limitO2(direction, p_pos_adim(p_ind),scal1D(p_ind), bl_type(bl_ind), limit(p_ind:p_ind+1), remesh_buffer) + call AC_remesh_limitO2(direction, p_pos_adim(p_ind),scal1D(p_ind), bl_type(bl_ind), & + limit(p_ind:p_ind+1), remesh_buffer) ! Second particle is remeshed with left formula - call AC_remesh_limitO2(direction, p_pos_adim(p_ind+1),scal1D(p_ind+1), bl_type(bl_ind+1), limit(p_ind+1:p_ind+2), remesh_buffer) + call AC_remesh_limitO2(direction, p_pos_adim(p_ind+1),scal1D(p_ind+1), bl_type(bl_ind+1), & + limit(p_ind+1:p_ind+2), remesh_buffer) end if end do -- GitLab