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

  INTEGER,SAVE :: COMM_LMDZ
  INTEGER,SAVE :: MPI_REAL_LMDZ
 

CONTAINS 

  SUBROUTINE init_const_mpi
#ifdef CPP_IOIPSL
    USE ioipsl
#else
! if not using IOIPSL, we still need to use (a local version of) getin
    USE ioipsl_getincom
#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
       CALL prism_init_comp_proto(comp_id, 'lmdz.x', ierr)
       CALL prism_get_localcomm_proto(comm_lmdz,ierr)
!$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
  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
!$OMP END MASTER
#endif

   END SUBROUTINE init_mpi
    
END MODULE mod_const_mpi