Changeset 15299 for NEMO/trunk/src/OCE/LBC/mppini.F90

Timestamp:

2021-09-28T17:28:00+02:00 (21 months ago)

Author:

smasson

Message:

trunk: important reproducibility bugfix (coastal F points can be used in the code)

File:

• NEMO/trunk/src/OCE/LBC/mppini.F90 (modified) (1 diff)

NEMO/trunk/src/OCE/LBC/mppini.F90

@@ -1170 +1170 @@
          zmsk(jh+1:jh+Ni_0,jh+1:jh+Nj_0) = REAL(COUNT(lloce, dim = 3), wp)   ! define inner domain -> need REAL to use lbclnk
          CALL lbc_lnk('mppini', zmsk, 'T', 1._wp, khls = jh)                 ! fill halos
+         ! Beware, coastal F points can be used in the code -> we may need communications for these points F points even if tmask = 0
+         ! -> the mask we must use here is equal to 1 as soon as one of the 4 neighbours is oce (sum of the mask, not multiplication)
+         zmsk(jh+1:jh+Ni_0,jh+1:jh+Nj_0) = zmsk(jh+1:jh+Ni_0,jh+1  :jh+Nj_0  ) + zmsk(jh+1+1:jh+Ni_0+1,jh+1  :jh+Nj_0  )   &
+            &                            + zmsk(jh+1:jh+Ni_0,jh+1+1:jh+Nj_0+1) + zmsk(jh+1+1:jh+Ni_0+1,jh+1+1:jh+Nj_0+1)
+         CALL lbc_lnk('mppini', zmsk, 'T', 1._wp, khls = jh)                 ! fill halos again!
          !        
          iiwe = jh   ;   iiea = Ni_0   ! bottom-left corner - 1 of the sent data