Changeset 9850

Timestamp:

2018-06-27T17:12:25+02:00 (7 years ago)

Author:

smasson

Message:

dev_r9759_HPC09_ESIWACE: report bugfix from NEMOtrunk on not-masked NP, see #2101

File:

• NEMO/branches/2018/dev_r9759_HPC09_ESIWACE/src/OCE/LBC/mpp_nfd_generic.h90 (modified) (3 diffs)

NEMO/branches/2018/dev_r9759_HPC09_ESIWACE/src/OCE/LBC/mpp_nfd_generic.h90

@@ -95 +95 @@
          ALLOCATE( ztabl(jpimax   ,4,ipk,ipl,ipf) , ztabr(jpimax*jpmaxngh,4,ipk,ipl,ipf) ) 
          !
-         ztabr(:,:,:,:,:) = 0._wp
-         ztabl(:,:,:,:,:) = 0._wp
-         !
          DO jf = 1, ipf
             DO jl = 1, ipl
@@ -123 +120 @@
             iproc = nfipproc(isendto(jr),jpnj)
             IF(iproc /= -1) THEN
+               iilb = nimppt(iproc+1)
+               ildi = nldit (iproc+1)
                ilei = nleit (iproc+1)
-               ildi = nldit (iproc+1)
+               IF( iilb           ==      1 )   ildi = 1     ! e-w boundary already done -> force to take 1st column
+               IF( iilb + jpi - 1 == jpiglo )   ilei = jpi   ! e-w boundary already done -> force to take last column
                iilb = nfiimpp(isendto(jr),jpnj) - nfiimpp(isendto(1),jpnj)
             ENDIF
@@ -199 +199 @@
          ENDIF
          !
-         ztab(:,:,:,:,:) = 0._wp
          DO jr = 1, ndim_rank_north         ! recover the global north array
             iproc = nrank_north(jr) + 1
+            iilb  = nimppt(iproc)
             ildi  = nldit (iproc)
             ilei  = nleit (iproc)
-            iilb  = nimppt(iproc)
+            IF( iilb           ==      1 )   ildi = 1     ! e-w boundary already done -> force to take 1st column
+            IF( iilb + jpi - 1 == jpiglo )   ilei = jpi   ! e-w boundary already done -> force to take last column
             DO jf = 1, ipf
                DO jl = 1, ipl