parallel_mpi_mod.f

Go to the documentation of this file.

        module parallel_mpi

!       A module to carry out init, finalize and check for any parallel errors

        use geometry
        use compar
        implicit none

        contains

        subroutine parallel_init()
          implicit none
          integer :: ierr

          numpes = 1
          mype = 0

#ifdef MPI
          call mpi_init(ierr)
          call mpi_check( 'parallel_init:MPI_Init ', ierr)

          call mpi_comm_size( mpi_comm_world, numpes, ierr )
          call mpi_check( 'parallel_init:MPI_Comm_size ', ierr )

          call mpi_comm_rank( mpi_comm_world, mype, ierr )
          call mpi_check( 'parallel_init:MPI_Comm_size ', ierr )
#endif
          return
        end subroutine parallel_init

        subroutine parallel_fin()
          implicit none

#ifdef MPI
          integer :: ierr

          call mpi_finalize(ierr)
          call mpi_check( 'parallel_init:MPI_Finalize ', ierr)
#endif

          return
        end subroutine parallel_fin

        subroutine mpi_check( msg, ierr )
#ifdef MPI
          use mpi, only: mpi_success ! ignore-depcomp
#endif
          implicit none
          character(len=*),intent(in) :: msg
          integer, intent(in) :: ierr

#ifdef MPI
          character(len=512) :: errmsg
          integer :: resultlen, ierror

          if (ierr .ne. mpi_success ) then
             call mpi_error_string( ierr, errmsg, resultlen, ierror )
             print*, 'Error: ', msg
             print*, errmsg(1:resultlen)
             stop '** ERROR ** '
          endif
#endif

          return
        end subroutine mpi_check


        end module parallel_mpi