Changeset 9971 for NEMO/trunk/src/OCE/LBC/mpp_nfd_generic.h90

Timestamp:

2018-07-19T14:39:57+02:00 (6 years ago)

Author:

mathiot

Message:

initialization of ztabl, ztabr and ztab to 0 only in case of land suppression as suggested by Sebastien #2118

File:

• NEMO/trunk/src/OCE/LBC/mpp_nfd_generic.h90 (modified) (2 diffs)

NEMO/trunk/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
+         ! when some processors of the north fold are suppressed, 
+         ! values of ztab* arrays corresponding to these suppressed domain won't be defined 
+         ! and we need a default definition to 0.
+         ! a better test should be: a testing if "suppressed land-processors" belongs to the north-pole folding
+         IF ( jpni*jpnj /= jpnij ) THEN
+            ztabr(:,:,:,:,:) = 0._wp
+            ztabl(:,:,:,:,:) = 0._wp
+         END IF
          !
          DO jf = 1, ipf
@@ -183 +189 @@
          ALLOCATE( znorthgloio(jpimax,4,ipk,ipl,ipf,jpni) )
          !
-         ztab(:,:,:,:,:)=0._wp
+         ! when some processors of the north fold are suppressed,
+         ! values of ztab* arrays corresponding to these suppressed domain won't be defined
+         ! and we need a default definition to 0.
+         ! a better test should be: a testing if "suppressed land-processors" belongs to the north-pole folding
+         IF ( jpni*jpnj /= jpnij ) ztab(:,:,:,:,:)=0._wp
          !
          CALL MPI_ALLGATHER( znorthloc  , itaille, MPI_DOUBLE_PRECISION,                &