make_lcover.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.

!     #########
      SUBROUTINE make_lcover(OCOVER)
!     ##############################################################
!
!!**** *PGD_COVER* monitor for averaging and interpolations of cover fractions
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    10/12/97
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE modd_surfex_mpi, ONLY : nrank, npio, nproc, ncomm
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
#ifdef SFX_MNH
USE mode_fd_ll, ONLY : getfd,fd_ll
USE modd_io_ll, ONLY : isp, isnproc
USE modd_var_ll, ONLY : nmnh_comm_world
#endif
!
IMPLICIT NONE
!
#if defined(SFX_MPI) || defined(SFX_MNH)
include "mpif.h"
#endif
!
!*    0.1    Declaration of arguments
!            ------------------------
!
LOGICAL, DIMENSION(:), INTENT(INOUT) :: OCOVER
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
#ifdef SFX_MNH
TYPE(fd_ll), POINTER                     :: TZFD
#endif
!
INTEGER :: INFOMPI, JPROC, JCOVER
!
INTEGER :: IRANK_SAVE, IPROC_SAVE, IPIO_SAVE, ICOMM_SAVE
!
LOGICAL, DIMENSION(:,:), ALLOCATABLE :: GCOVER_ALL
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!---------------------------------------------------------------
!
!*    1.      Initializations
!             ---------------
!
IF (lhook) CALL dr_hook('MAKE_LCOVER',0,zhook_handle)
!
#ifdef SFX_MNH
tzfd=>getfd(nmnh_comm_world)
!
irank_save = nrank
iproc_save = nproc
ipio_save = npio
icomm_save = ncomm
!
! on met les infos de mésonh
nrank = isp-1
nproc = isnproc
npio = tzfd%OWNER-1
ncomm = tzfd%COMM
#endif
!
ALLOCATE(gcover_all(SIZE(ocover),0:nproc-1))
!
!
IF (nproc>1) THEN
#if defined(SFX_MPI) || defined(SFX_MNH)
  CALL mpi_allgather(ocover,SIZE(ocover),mpi_logical,gcover_all,SIZE(ocover),&
                  mpi_logical,ncomm,infompi)
#endif
ELSE
  gcover_all(:,0) = ocover(:)
ENDIF
!
!
ocover(:) = .false.
DO jproc = 0,nproc-1
  DO jcover=1,SIZE(ocover)
    IF (gcover_all(jcover,jproc)) ocover(jcover) = .true.
  ENDDO
ENDDO
!
DEALLOCATE(gcover_all)
!
!
IF (nproc>1) THEN
#if defined(SFX_MPI) || defined(SFX_MNH)
  CALL mpi_bcast(ocover,SIZE(ocover),mpi_logical,npio,ncomm,infompi)
#endif
ENDIF
!
#ifdef SFX_MNH
nrank = irank_save
nproc = iproc_save
npio = ipio_save
ncomm = icomm_save
#endif
!
IF (lhook) CALL dr_hook('MAKE_LCOVER',1,zhook_handle)
!
END SUBROUTINE make_lcover