Changeset 10292 for NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/OCE/LBC/lib_mpp.F90

Timestamp:

2018-11-09T16:35:08+01:00 (5 years ago)

Author:

smasson

Message:

dev_r10164_HPC09_ESIWACE_PREP_MERGE: action 4b: reduce communications in si3, see #2133

File:

• NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/OCE/LBC/lib_mpp.F90 (modified) (3 diffs)

NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/OCE/LBC/lib_mpp.F90

@@ -83 +83 @@
    PUBLIC   mpp_lbc_north_icb
    PUBLIC   mpp_min, mpp_max, mpp_sum, mpp_minloc, mpp_maxloc
+   PUBLIC   mpp_ilor
    PUBLIC   mpp_max_multiple
    PUBLIC   mppscatter, mppgather
@@ -633 +634 @@
       ptab = zwork
    END SUBROUTINE mppmax_real
+   !!
+   SUBROUTINE mpp_ilor( ld_switch, ldlast, kcom )
+      ! WARNING: must be used only once (by ice_dyn_adv_umx) because ll_switch and ireq are SAVE
+      !!----------------------------------------------------------------------
+      LOGICAL, INTENT(inout), DIMENSION(2) ::   ld_switch
+      LOGICAL, INTENT(in   ), OPTIONAL     ::   ldlast
+      INTEGER, INTENT(in   ), OPTIONAL     ::   kcom 
+      INTEGER  ::   ierror, ilocalcomm
+      LOGICAL, SAVE ::   ll_switch 
+      INTEGER, SAVE ::   ireq = -1
+      !!----------------------------------------------------------------------
+      ilocalcomm = mpi_comm_oce
+      IF( PRESENT(kcom) )   ilocalcomm = kcom
+      
+      IF ( ireq /= -1 ) THEN   ! get ld_switch(2) from ll_switch (from previous call)
+         CALL mpi_wait(ireq, MPI_STATUS_IGNORE, ierror )
+         ld_switch(2) = ll_switch
+      ENDIF
+      IF( .NOT. ldlast ) &     ! send ll_switch to be received on next call
+         CALL mpi_iallreduce( ld_switch(1), ll_switch, 1, MPI_LOGICAL, mpi_lor, ilocalcomm, ireq, ierror )
+
+   END SUBROUTINE mpp_ilor
 
 
@@ -1720 +1743 @@
    END SUBROUTINE mpp_maxloc3d
 
+   SUBROUTINE mpp_ilor( ld_switch, ldlast, kcom )
+      LOGICAL, INTENT(in   ), DIMENSION(2) ::   ld_switch
+      LOGICAL, INTENT(in   ), OPTIONAL     ::   ldlast
+      INTEGER, INTENT(in   ), OPTIONAL     ::   kcom    ! ???
+      WRITE(*,*) 'mpp_ilor: You should not have seen this print! error?', ld_switch
+   END SUBROUTINE mpp_ilor
+
    SUBROUTINE mppstop
       STOP      ! non MPP case, just stop the run