diff --git a/src/OCE/DYN/dynvor.F90 b/src/OCE/DYN/dynvor.F90 index 0f43b4bd275ceeeb91202bb4c162a9c104604bed..b3fe561915a6ec73bb66996d59697f0e4e7b2888 100644 --- a/src/OCE/DYN/dynvor.F90 +++ b/src/OCE/DYN/dynvor.F90 @@ -955,10 +955,10 @@ CONTAINS ztsw_jp1 = zwz(ji ,jj ) + zwz(ji-1,jj ) + zwz(ji-1,jj+1) ztsw_ip1 = zwz(ji+1,jj-1) + zwz(ji ,jj-1) + zwz(ji ,jj ) ! - zua = + r1_12 * r1_e1u(ji,jj) * ( ztne * zwy + ztnw_ip1 * zwy_ip1 & - & + ztse * zwy_jm1 + ztsw_ip1 * zwy_ip1_jm1 ) - zva = - r1_12 * r1_e2v(ji,jj) * ( ztsw_jp1 * zwx_im1_jp1 + ztse_jp1 * zwx_jp1 & - & + ztnw * zwx_im1 + ztne * zwx ) + zua = + r1_12 * r1_e1u(ji,jj) * ( ( ztne * zwy + ztnw_ip1 * zwy_ip1 ) & ! add () for + & + ( ztse * zwy_jm1 + ztsw_ip1 * zwy_ip1_jm1 ) ) ! NP repro + zva = - r1_12 * r1_e2v(ji,jj) * ( ( ztsw_jp1 * zwx_im1_jp1 + ztse_jp1 * zwx_jp1 ) & + & + ( ztnw * zwx_im1 + ztne * zwx ) ) pu_rhs(ji,jj,jk) = pu_rhs(ji,jj,jk) + zua pv_rhs(ji,jj,jk) = pv_rhs(ji,jj,jk) + zva END_2D diff --git a/src/OCE/LDF/ldfslp.F90 b/src/OCE/LDF/ldfslp.F90 index f25ec900e5d8a13b8beefeb338ac1526f39fb5a5..e6247b49f6124844a48ede07ed6c755638df4ffe 100644 --- a/src/OCE/LDF/ldfslp.F90 +++ b/src/OCE/LDF/ldfslp.F90 @@ -243,7 +243,6 @@ CONTAINS & + 2.*( ( zww(ji ,jj-1,jk) + zww(ji-1,jj ,jk) ) & & + ( zww(ji+1,jj ,jk) + zww(ji ,jj+1,jk) ) ) & & + 4.* zww(ji,jj ,jk) ) - END_2D ! !* decrease along coastal boundaries uslp(ji,jj,jk) = zuslp * ( umask(ji,jj+1,jk) + umask(ji,jj-1,jk ) ) * 0.5_wp & & * ( umask(ji,jj ,jk) + umask(ji,jj ,jk+1) ) * 0.5_wp