diff --git a/HySoP/src/scalesInterface/particles/advec_remesh_Mprime.f90 b/HySoP/src/scalesInterface/particles/advec_remesh_Mprime.f90 index acf7c62b8ccb722da21ed2032414ac1e853c8a1d..60ee52d14a38539a786f462342d4dbb24f3c5505 100644 --- a/HySoP/src/scalesInterface/particles/advec_remesh_Mprime.f90 +++ b/HySoP/src/scalesInterface/particles/advec_remesh_Mprime.f90 @@ -359,7 +359,7 @@ subroutine AC_remesh_Mprime8_array(dir, pos_adim, sca, buffer) bP3=y0*(y0*(y0*(y0*(y0*(y0*(-y0/48. + 3./32.) - 1./12.) & & - 1./32.) - 1./48.) + 1./96.) + 1./60.) + 17./3360. ! bM2=(2-y0)**7/2688.-(y0+2)*(2-y0)**6/640.+(y0+2)**2*(2-y0)**5/960 - ! -xx2**7/336+(y0+2)*xx2**6/80.-(y0+2)**2*xx2**5/120. + ! -(1-y0)**7/336+(y0+2)*(1-y0)**6/80.-(y0+2)**2*(1-y0)**5/120. bM2=y0*(y0*(y0*(y0*(y0*(y0*(y0/48. - 5./96.) - 1./24.) & & + 11./48.) - 1./6.) - 5./48.) + 3./20.) - 17./560. ! bP2=(y0+2)**7/2688.-(2-y0)*(y0+2)**6/640.+(2-y0)**2*(y0+2)**5/960 @@ -446,7 +446,7 @@ subroutine AC_remesh_Mprime8_pter(pos_adim, sca, buffer) bP3=y0*(y0*(y0*(y0*(y0*(y0*(-y0/48. + 3./32.) - 1./12.) & & - 1./32.) - 1./48.) + 1./96.) + 1./60.) + 17./3360. ! bM2=(2-y0)**7/2688.-(y0+2)*(2-y0)**6/640.+(y0+2)**2*(2-y0)**5/960 - ! -xx2**7/336+(y0+2)*xx2**6/80.-(y0+2)**2*xx2**5/120. + ! -(1-y0)**7/336+(y0+2)*(1-y0)**6/80.-(y0+2)**2*(1-y0)**5/120. bM2=y0*(y0*(y0*(y0*(y0*(y0*(y0/48. - 5./96.) - 1./24.) & & + 11./48.) - 1./6.) - 5./48.) + 3./20.) - 17./560. ! bP2=(y0+2)**7/2688.-(2-y0)*(y0+2)**6/640.+(2-y0)**2*(y0+2)**5/960