Changeset 5857 for branches/2015/dev_r5776_UKMO2_OBS_efficiency_improvs/NEMOGCM/NEMO/OPA_SRC/OBS/obs_mpp.F90

Timestamp:

2015-11-02T18:28:38+01:00 (8 years ago)

Author:

jenniewaters

Message:

Removing global arrays when ln_grid_global=.false.

File:

• branches/2015/dev_r5776_UKMO2_OBS_efficiency_improvs/NEMOGCM/NEMO/OPA_SRC/OBS/obs_mpp.F90 (modified) (3 diffs)

branches/2015/dev_r5776_UKMO2_OBS_efficiency_improvs/NEMOGCM/NEMO/OPA_SRC/OBS/obs_mpp.F90

@@ -13 +13 @@
    !! obs_mpp_max_integer   : Find maximum on all processors of each value in an integer on all processors
    !! obs_mpp_find_obs_proc : Find processors which should hold the observations
+   !! obs_mpp_find_obs_proc_local : Find processors which should hold the observations, avoiding global arrays
    !! obs_mpp_sum_integers  : Sum an integer array from all processors
    !! obs_mpp_sum_integer   : Sum an integer from all processors
@@ -28 +29 @@
       &   obs_mpp_max_integer,   & !: Find maximum across processors in an integer array
       &   obs_mpp_find_obs_proc, & !: Find processors which should hold the observations
+      &   obs_mpp_find_obs_proc_local, & !: Find processors which should hold the observations
       &   obs_mpp_sum_integers,  & !: Sum an integer array from all processors
       &   obs_mpp_sum_integer,   & !: Sum an integer from all processors
@@ -227 +229 @@
    END SUBROUTINE obs_mpp_find_obs_proc
 
+   SUBROUTINE obs_mpp_find_obs_proc_local( kobsp,kno )
+      !!----------------------------------------------------------------------
+      !!               ***  ROUTINE obs_mpp_find_obs_proc_local ***
+      !!         
+      !! ** Purpose : From the array kobsp containing the results of the grid
+      !!              grid search on each processor the processor return a
+      !!              decision of which processors should hold the observation.
+      !!
+      !! ** Method  : Synchronize the processor number for each obs using
+      !!              obs_mpp_max_integer. If an observation exists on two 
+      !!              processors it will be allocated to the lower numbered
+      !!              processor.
+      !!
+      !! ** Action  : This does only work for MPI.
+      !!              It does not work for SHMEM.
+      !!
+      !! References : http://www.mpi-forum.org
+      !!----------------------------------------------------------------------
+      INTEGER                , INTENT(in   ) ::   kno
+      INTEGER, DIMENSION(kno), INTENT(inout) ::   kobsp
+      !
+#if defined key_mpp_mpi
+      !
+      !
+      INTEGER :: ji, isum
+      INTEGER, DIMENSION(kno) ::   iobsp
+INCLUDE 'mpif.h'
+      !!
+
+      iobsp=kobsp
+
+      WHERE( iobsp(:) == -1 )
+         iobsp(:) = 9999999
+      END WHERE
+
+      iobsp=-1*iobsp
+
+      CALL obs_mpp_max_integer( iobsp, kno )
+
+      kobsp=-1*iobsp
+
+      isum=0
+      DO ji = 1, kno
+         IF ( kobsp(ji) == 9999999 ) THEN
+            isum=isum+1
+            kobsp(ji)=-1
+         ENDIF
+      ENDDO
+
+
+      IF ( isum > 0 ) THEN
+         IF (lwp) WRITE(numout,*) isum, ' observations failed the grid search.'
+         IF (lwp) WRITE(numout,*)'If ln_grid_search_lookup=.TRUE., try reducing grid_search_res'
+      ENDIF
+
+#else
+      ! no MPI: empty routine
+#endif     
+      
+   END SUBROUTINE obs_mpp_find_obs_proc_local
+
 
    SUBROUTINE obs_mpp_sum_integers( kvalsin, kvalsout, kno )