mod_const_mpi.F90

Go to the documentation of this file.

! 
! $Id: mod_const_para.F90 1279 2009-12-10 09:02:56Z fairhead $
!
MODULE mod_const_mpi
  IMPLICIT NONE
  INTEGER,SAVE :: comm_lmdz
  INTEGER,SAVE :: mpi_real_lmdz
 

CONTAINS 

  SUBROUTINE init_const_mpi
#ifdef CPP_IOIPSL
    USE ioipsl, ONLY: getin
#else
! if not using IOIPSL, we still need to use (a local version of) getin
    USE ioipsl_getincom, only: getin
#endif
! Use of Oasis-MCT coupler 
#ifdef CPP_OMCT
    USE mod_prism
#endif
#ifdef CPP_XIOS
    USE wxios, only: wxios_init
#endif
    IMPLICIT NONE
#ifdef CPP_MPI
    include 'mpif.h'
#endif

    INTEGER             :: ierr
    INTEGER             :: comp_id
    INTEGER             :: thread_required
    INTEGER             :: thread_provided
    CHARACTER(len = 6)  :: type_ocean

!$OMP MASTER
    type_ocean = 'force '
    CALL getin('type_ocean', type_ocean)
!$OMP END MASTER
!$OMP BARRIER

    IF (type_ocean=='couple') THEN
#ifdef CPP_COUPLE
!$OMP MASTER
#ifdef CPP_XIOS
      CALL wxios_init("LMDZ", outcom=comm_lmdz, type_ocean=type_ocean)
#else
       CALL prism_init_comp_proto (comp_id, 'LMDZ', ierr)
       CALL prism_get_localcomm_proto(comm_lmdz,ierr)
#endif
!$OMP END MASTER
#endif
#ifdef CPP_MPI
      mpi_real_lmdz=mpi_real8
#endif
    ELSE
      CALL init_mpi
    ENDIF

  END SUBROUTINE init_const_mpi
  
  SUBROUTINE init_mpi
#ifdef CPP_XIOS
    USE wxios, only: wxios_init
#endif
  IMPLICIT NONE
#ifdef CPP_MPI
     include 'mpif.h'
#endif
    INTEGER             :: ierr
    INTEGER             :: thread_required
    INTEGER             :: thread_provided

#ifdef CPP_MPI
!$OMP MASTER
      thread_required=mpi_thread_serialized

      CALL mpi_init_thread(thread_required,thread_provided,ierr)
      IF (thread_provided < thread_required) THEN
        print *,'Warning : The multithreaded level of MPI librairy do not provide the requiered level',  &
                ' in mod_const_mpi::Init_const_mpi'
      ENDIF
      comm_lmdz=mpi_comm_world
      mpi_real_lmdz=mpi_real8
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Initialisation de XIOS
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
#ifdef CPP_XIOS
      CALL wxios_init("LMDZ", outcom=comm_lmdz)
#endif
!$OMP END MASTER
#else
#ifdef CPP_XIOS
!$OMP MASTER
      CALL wxios_init("LMDZ")
!$OMP END MASTER
#endif
#endif

   END SUBROUTINE init_mpi
    
END MODULE mod_const_mpi