module Communication for Ensemble system More...

Functions/Subroutines
subroutine, public	comm_ensemble_setup
	Setup. More...

subroutine, public	comm_ensemble_finalize
	finalize More...

Variables
integer, public	comm_ensemble_world

integer, public	comm_ensemble_myrank

integer, public	comm_ensemble_nprocs

Detailed Description

module Communication for Ensemble system

Description: MPI Communication module for Ensemble system

Author: Team SCALE

Function/Subroutine Documentation

◆ comm_ensemble_setup()

subroutine, public scale_comm_ensemble::comm_ensemble_setup

Setup.

Definition at line 52 of file scale_comm_ensemble.F90.

     use scale_prc, only: &
        prc_universal_comm_world, &
        prc_global_myrank,        &
        prc_global_nprocs,        &
        prc_local_comm_world,     &
        prc_nprocs,               &
        prc_myrank,               &
        prc_abort
     implicit none
  
     integer, allocatable :: ensemble2color(:)
     integer, allocatable :: ensemble2key  (:)
  
     integer :: n
     integer :: myrank, nprocs
     integer :: ierr
     !---------------------------------------------------------------------------
  
     log_newline
     log_info("COMM_ENSEMBLE_setup",*) 'Setup'
  
     !!--- read namelist
     !rewind(IO_FID_CONF)
     !read(IO_FID_CONF,nml=PARAM_COMM_ENSEMBLE,iostat=ierr)
     !if( ierr < 0 ) then !--- missing
     !   LOG_INFO("COMM_ENSEMBLE_setup",*) 'Not found namelist. Default used.'
     !elseif( ierr > 0 ) then !--- fatal error
     !   LOG_ERROR("COMM_ENSEMBLE_setup",*) 'Not appropriate names in namelist PARAM_COMM_ENSEMBLE. Check!'
     !   call PRC_abort
     !endif
     !LOG_NML(PARAM_COMM_ENSEMBLE)
  
     ! [global] -> [ensemble]
     if( prc_global_nprocs < 2 ) then
        log_error("COMM_ENSEMBLE_setup",*) 'PRC_GLOBAL_nprocs must be greater than or equal to 2:', prc_global_nprocs
        call prc_abort
     end if
  
     allocate( ensemble2color(0:prc_global_nprocs-1) )
     allocate( ensemble2key(0:prc_global_nprocs-1) )
  
     do n = 0, prc_global_nprocs-1
        ensemble2color(n) = n
        ensemble2key(n) = n
     end do
  
     call mpi_comm_split( prc_universal_comm_world,          & ! [IN]
                          ensemble2color(prc_global_myrank), & ! [IN]
                          ensemble2key(prc_global_myrank), & ! [IN]
                          comm_ensemble_world,               & ! [OUT]
                          ierr                               ) ! [OUT]
  
     call mpi_comm_rank( comm_ensemble_world, comm_ensemble_myrank, ierr )
     call mpi_comm_size( comm_ensemble_world, comm_ensemble_nprocs, ierr )
  
     deallocate( ensemble2color )
     deallocate( ensemble2key   )
  
     return

References comm_ensemble_myrank, comm_ensemble_nprocs, comm_ensemble_world, scale_prc::prc_abort(), scale_prc::prc_global_myrank, scale_prc::prc_global_nprocs, scale_prc::prc_local_comm_world, scale_prc::prc_myrank, scale_prc::prc_nprocs, and scale_prc::prc_universal_comm_world.

Referenced by mod_da_driver::da_driver_setup().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ comm_ensemble_finalize()

subroutine, public scale_comm_ensemble::comm_ensemble_finalize

finalize

Definition at line 116 of file scale_comm_ensemble.F90.

     implicit none
  
     integer :: ierr
     !---------------------------------------------------------------------------
  
     call mpi_comm_free( comm_ensemble_world, ierr )
  
     return