modd_surfex_mpi.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.
!     ######################
      MODULE modd_surfex_mpi
!     ######################
!
!!****  *MODD_SURFEX_MPI
!!
!!    PURPOSE
!!    -------
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!      None 
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!      S. Faroux   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original       04/06/12
!!      B. Decharme    04/2013  Delete NWG_LAYER_TOT, NWG_SIZE
!!      J.Escobar      10/06/2013: replace DOUBLE PRECISION by REAL to handle problem for promotion of real on IBM SP
!
!-------------------------------------------------------------------------------
!
!*       0.   DECLARATIONS
!             ------------
!
USE modd_surfex_omp, ONLY : nblock
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!-------------------------------------------------------------------------------
!
INTEGER, DIMENSION(:), ALLOCATABLE :: NINDEX
INTEGER, DIMENSION(:), ALLOCATABLE :: NSIZE_TASK
!
INTEGER :: NPROC = 1
INTEGER :: NRANK = 0
INTEGER :: NCOMM = -1
INTEGER :: NPIO  = 0
INTEGER :: NSIZE = 0
INTEGER :: MYPROC = 0
!
INTEGER :: IDX_R = 0
INTEGER :: IDX_W = 0
!
REAL, DIMENSION(5) :: XTIME_WRITE
!$OMP THREADPRIVATE(XTIME_WRITE)
REAL, DIMENSION(6) :: XTIME_CALC
!$OMP THREADPRIVATE(XTIME_CALC)
REAL :: XTIME_NPIO_READ=0
!$OMP THREADPRIVATE(XTIME_NPIO_READ)
REAL :: XTIME_NPIO_WRITE=0
!$OMP THREADPRIVATE(XTIME_NPIO_WRITE)
REAL :: XTIME_COMM_READ=0
!$OMP THREADPRIVATE(XTIME_COMM_READ)
REAL :: XTIME_CALC_READ=0
!$OMP THREADPRIVATE(XTIME_CALC_READ)
REAL :: XTIME_COMM_WRITE=0
!$OMP THREADPRIVATE(XTIME_COMM_WRITE)
REAL :: XTIME_CALC_WRITE=0
!$OMP THREADPRIVATE(XTIME_CALC_WRITE)
REAL :: XTIME_OMP_BARR=0
!$OMP THREADPRIVATE(XTIME_OMP_BARR)
!
REAL :: XTIME_INIT_SEA=0
!$OMP THREADPRIVATE(XTIME_INIT_SEA)
REAL :: XTIME_INIT_WATER=0
!$OMP THREADPRIVATE(XTIME_INIT_WATER)
REAL :: XTIME_INIT_NATURE=0
!$OMP THREADPRIVATE(XTIME_INIT_NATURE)
REAL :: XTIME_INIT_TOWN=0
!$OMP THREADPRIVATE(XTIME_INIT_TOWN)
!
REAL :: XTIME_SEA=0
!$OMP THREADPRIVATE(XTIME_SEA)
REAL :: XTIME_WATER=0
!$OMP THREADPRIVATE(XTIME_WATER)
REAL :: XTIME_NATURE=0
!$OMP THREADPRIVATE(XTIME_NATURE)
REAL :: XTIME_TOWN=0
!$OMP THREADPRIVATE(XTIME_TOWN)
!
 CONTAINS
!
SUBROUTINE wlog_mpi(HLOG,PLOG,KLOG,KLOG2,OLOG)
!
IMPLICIT NONE
!
 CHARACTER(LEN=*), INTENT(IN) :: hlog
REAL, INTENT(IN),OPTIONAL :: plog
INTEGER, INTENT(IN), OPTIONAL :: klog
INTEGER, INTENT(IN), OPTIONAL :: klog2
LOGICAL, INTENT(IN), OPTIONAL :: olog
!
INTEGER :: iunit
!
REAL(KIND=JPRB) :: zhook_handle
!
IF (lhook) CALL dr_hook('MODD_SURFEX_MPI:WLOG_MPI',0,zhook_handle)
!
iunit=100000+nrank
!
!$OMP CRITICAL
!
IF (present(olog)) THEN
  IF (present(plog)) THEN
    IF (present(klog)) THEN
      IF (present(klog2)) THEN
        WRITE(iunit,*) nblock, hlog, klog, klog2, plog, olog
      ELSE
        WRITE(iunit,*) nblock, hlog, klog, plog, olog
      ENDIF
    ELSE
      WRITE(iunit,*) nblock, hlog, plog, olog
    ENDIF
  ELSEIF (present(klog)) THEN
    IF (present(klog2)) THEN
      WRITE(iunit,*) nblock, hlog, klog, klog2, olog
    ELSE
      WRITE(iunit,*) nblock, hlog, klog, olog
    ENDIF
  ELSE
    WRITE(iunit,*) nblock, hlog, olog
  ENDIF
ELSEIF (present(plog)) THEN
  IF (present(klog)) THEN
    IF (present(klog2)) THEN
      WRITE(iunit,*) nblock, hlog, klog, klog2, plog
    ELSE
      WRITE(iunit,*) nblock, hlog, klog, plog
    ENDIF
  ELSE
    WRITE(iunit,*) nblock, hlog, plog
  ENDIF
ELSEIF (present(klog)) THEN
  IF (present(klog2)) THEN
    WRITE(iunit,*) nblock, hlog, klog, klog2
  ELSE
    WRITE(iunit,*) nblock, hlog, klog
  ENDIF
ELSE
  WRITE(iunit,*) nblock, hlog
ENDIF
!
!$OMP END CRITICAL
!
IF (lhook) CALL dr_hook('MODD_SURFEX_MPI:WLOG_MPI',1,zhook_handle)
!
END SUBROUTINE wlog_mpi
!
SUBROUTINE prep_log_mpi
!
 CHARACTER(LEN=100) :: yname
 CHARACTER(LEN=10)  :: yrank
INTEGER :: iunit
!
REAL(KIND=JPRB) :: zhook_handle
!
IF (lhook) CALL dr_hook('MODD_SURFEX_MPI:PREP_LOG_MPI',0,zhook_handle)
!
WRITE(yrank,fmt='(I10)') nrank
yname='log'//adjustl(yrank)
!
iunit=100000+nrank
!
!$OMP SINGLE
OPEN(unit=iunit,file=trim(yname),form='FORMATTED')
WRITE(iunit,*) "Log offline pour proc ",nrank
!$OMP END SINGLE
!
IF (lhook) CALL dr_hook('MODD_SURFEX_MPI:PREP_LOG_MPI',1,zhook_handle)
!
END SUBROUTINE prep_log_mpi
!
SUBROUTINE end_log_mpi
!
INTEGER :: iunit
!
REAL(KIND=JPRB) :: zhook_handle
!
IF (lhook) CALL dr_hook('MODD_SURFEX_MPI:END_LOG_MPI',0,zhook_handle)
!
iunit=100000+nrank
!
!$OMP SINGLE
WRITE(iunit,*) "End log offline pour proc ",nrank
 CLOSE(unit=iunit)
!$OMP END SINGLE
!
IF (lhook) CALL dr_hook('MODD_SURFEX_MPI:END_LOG_MPI',1,zhook_handle)
!
END SUBROUTINE end_log_mpi
!-------------------------------------------------------------------------------
!
END MODULE modd_surfex_mpi