Changeset 2587 for branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/mpp_tam.F90

Timestamp:

2011-02-15T12:58:59+01:00 (13 years ago)

Author:

vidard

Message:

refer to ticket #798

File:

• branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/mpp_tam.F90 (modified) (3 diffs)

branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/mpp_tam.F90

@@ -47 +47 @@
       & mpp_sum_reals,      &
       & mpp_global_or,      &
-      & mpp_global_max_real
+      & mpp_global_max_real, &
+      & mpp_global_max_real2
 
 CONTAINS
@@ -151 +152 @@
    SUBROUTINE mpp_global_max_real( zin, zout )
       !!----------------------------------------------------------------------
-      !!               ***  ROUTINE mpp_global_or ***
+      !!               ***  ROUTINE mpp_global_max_real ***
       !!          
       !! ** Purpose : Copy a local zin array to a global array and
@@ -213 +214 @@
    END SUBROUTINE mpp_global_max_real
 
+
+   SUBROUTINE mpp_global_max_real2( zin, zout )
+      !!----------------------------------------------------------------------
+      !!               ***  ROUTINE mpp_global_max_real2 ***
+      !!          
+      !! ** Purpose : Copy a local zin array to a global array and
+      !!              apply the "max" operation for all elements in
+      !!              a global (jpiglo,jpjglo) array across processors
+      !!
+      !! ** Method  : MPI allreduce
+      !!
+      !! ** Action  : This does only work for MPI. 
+      !!              It does not work for SHMEM.
+      !!
+      !! References : http://www.mpi-forum.org
+      !!
+      !! History :
+      !!        !  08-01  (K. Mogensen)  Original code
+      !!----------------------------------------------------------------------
+
+      !! * Arguments
+      REAL(wp), DIMENSION(jpi,jpj), INTENT(IN) :: &
+         & zin
+      REAL(wp), DIMENSION(jpiglo,jpjglo), INTENT(OUT) :: &
+         & zout
+      !! * Local declarations
+      INTEGER :: &
+         & ierr
+#if defined key_mpp_mpi
+#include <mpif.h>
+      REAL(wp), DIMENSION(:,:), ALLOCATABLE :: &
+         & zcp
+      INTEGER :: &
+         & ji, &
+         & jj
+      
+      ! Copy data for input to MPI
+
+      ALLOCATE( &
+         & zcp(jpiglo,jpjglo) &
+         & )
+      zcp(:,:) = -1e+38
+      DO jj = nldj, nlej
+         DO ji = nldi, nlei
+            zcp(mig(ji),mjg(jj)) = zin(ji,jj)
+         ENDDO
+      ENDDO
+
+      ! Call the MPI library to find the coast lines globally
+
+      CALL mpi_allreduce( zcp, zout, jpiglo*jpjglo, mpivar, &
+         &                mpi_max, mpi_comm_opa, ierr )
+
+      DEALLOCATE( &
+         & zcp &
+         & )
+
+#elif defined key_mpp_shmem
+#error "Only MPI support for MPP in NEMOVAR"
+#else
+      zout(:,:) = zin(:,:)
+#endif
+      
+   END SUBROUTINE mpp_global_max_real2
+   
 END MODULE mpp_tam