!==  IN: ptab is an array  ==!
#      define ARRAY_TYPE(i,j,k)    REAL(wp)          , INTENT(in   ) ::   ARRAY_IN(i,j,k)
#      define MASK_TYPE(i,j,k)     REAL(wp)          , INTENT(in   ) ::   MASK_IN(i,j,k)
#   if defined DIM_2d
#      define ARRAY_IN(i,j,k)   ptab(i,j)
#      define MASK_IN(i,j,k)    pmask(i,j)
#      define INDEX_TYPE(k)        INTEGER           , INTENT(  out) ::   kindex(2)
#      define K_SIZE(ptab)      1
#   endif
#   if defined DIM_3d
#      define ARRAY_IN(i,j,k)   ptab(i,j,k)
#      define MASK_IN(i,j,k)    pmask(i,j,k)
#      define INDEX_TYPE(k)        INTEGER           , INTENT(  out) ::   kindex(3)
#      define K_SIZE(ptab)      SIZE(ptab,3)
#   endif
#   if defined OPERATION_MAXLOC
#      define MPI_OPERATION mpi_maxloc
#   endif
#   if defined OPERATION_MINLOC
#      define MPI_OPERATION mpi_minloc
#   endif

   SUBROUTINE ROUTINE_LOC( cdname, ptab, pmask, pmin, kindex )
      !!----------------------------------------------------------------------
      CHARACTER(len=*), INTENT(in   ) ::   cdname  ! name of the calling subroutine
      ARRAY_TYPE(:,:,:)                                ! array on which loctrans operation is applied
      MASK_TYPE(:,:,:)                                 ! local mask
      REAL(wp)        , INTENT(  out) ::   pmin    ! Global minimum of ptab
      INDEX_TYPE(:)                                ! index of minimum in global frame
# if defined key_mpp_mpi
      !
      INTEGER  ::   ierror, ii, idim
      REAL(wp) ::   zmin     ! local minimum
      INTEGER , DIMENSION(:), ALLOCATABLE  ::   ilocs
      REAL(wp), DIMENSION(2,1) ::   zain, zaout
      !!-----------------------------------------------------------------------
      !
      IF( narea == 1 .AND. numcom == -1 ) CALL mpp_report( cdname, ld_glb = .TRUE. )
      !
      idim = SIZE(kindex)
      ALLOCATE ( ilocs(idim) )
      !
      zmin  = MINVAL( ARRAY_IN(:,:,:) , mask= MASK_IN(:,:,:) == 1._wp )
      ilocs = MINLOC( ARRAY_IN(:,:,:) , mask= MASK_IN(:,:,:) == 1._wp )
      !
      DO ii = 1, idim
        IF (ii == 1) kindex(1) = ilocs(1) + nimpp - 1
        IF (ii == 2) kindex(2) = ilocs(2) + njmpp - 1
        IF (ii == 3) kindex(3) = ilocs(3)
      ENDDO
      ! 
      DEALLOCATE (ilocs)
      !
      zain(1,:) = zmin
      zain(2,:) = kindex(1) + 10000.* kindex(2)
      IF ( idim == 3) zain(2,:) = zain(2,:) + 100000000.*kindex(3)
      !
      IF( ln_timing ) CALL tic_tac(.TRUE., ld_global = .TRUE.)
      CALL MPI_ALLREDUCE( zain, zaout, 1, MPI_2DOUBLE_PRECISION, MPI_OPERATION ,MPI_COMM_OCE, ierror)
      IF( ln_timing ) CALL tic_tac(.FALSE., ld_global = .TRUE.)
      !
      pmin = zaout(1,1)
      IF ( idim == 3) THEN
        kindex(3)   = INT( zaout(2,1) / 100000000. )
        kindex(2)   = INT( zaout(2,1) - kindex(3) * 100000000. ) / 10000
        kindex(1)   = INT( zaout(2,1) - kindex(3) * 100000000. -kindex(2) * 10000. )
      ELSE
        kindex(2)   = INT( zaout(2,1) / 10000 )
        kindex(1)   = INT( zaout(2,1) -kindex(2) * 10000. )
      ENDIF
#else
      kindex = 0 ; pmin = 0.
      WRITE(*,*) 'ROUTINE_LOC: You should not have seen this print! error?'
#endif

   END SUBROUTINE ROUTINE_LOC

#undef ARRAY_TYPE
#undef MAX_TYPE
#undef ARRAY_IN
#undef MASK_IN
#undef K_SIZE
#undef MPI_OPERATION
#undef INDEX_TYPE