subroutine mpii_end

    ! local variables
    integer :: ierr
    
    call mpi_finalize( ierr )
    
    
  end subroutine mpii_end

  subroutine mpii_init

    ! local variables
    integer :: ierr
    
    call mpi_init( ierr )
    call mpi_comm_size( mpi_comm_world, nprocs, ierr )
    call mpi_comm_rank( mpi_comm_world, myrank, ierr )
    
!    write(*,*) mpi_comm_world, nprocs, myrank, ierr

  end subroutine mpii_init

  subroutine mpii_irecv_dble_2d( idep, im, jm, buf, ireq )

    integer, intent(in)  :: idep
    integer, intent(in)  :: im, jm
    real(8), intent(out) :: buf( im, jm )
    integer, intent(out) :: ireq
    integer :: ierr
    integer :: isize
    
    isize = size( buf )
    call mpi_irecv( buf, isize, mpi_double_precision, idep, 1, mpi_comm_world, ireq, ierr )
    
  end subroutine mpii_irecv_dble_2d

  subroutine mpii_irecv_dble_3d( idep, im, jm, km, buf, ireq )

    integer, intent(in)  :: idep
    integer, intent(in)  :: im, jm, km
    real(8), intent(out) :: buf( im, jm, km )
    integer, intent(out) :: ireq
    integer :: ierr
    integer :: isize
    
    isize = size( buf )
    call mpi_irecv( buf, isize, mpi_double_precision, idep, 1, mpi_comm_world, ireq, ierr )
    
  end subroutine mpii_irecv_dble_3d

  subroutine mpii_irecv_dble_4d( idep, im, jm, km, lm, buf, ireq )

    integer, intent(in)  :: idep
    integer, intent(in)  :: im, jm, km, lm
    real(8), intent(out) :: buf( im, jm, km, lm )
    integer, intent(out) :: ireq
    integer :: ierr
    integer :: isize
    
    isize = size( buf )
    call mpi_irecv( buf, isize, mpi_double_precision, idep, 1, mpi_comm_world, ireq, ierr )
    
  end subroutine mpii_irecv_dble_4d

  subroutine mpii_isend_dble_2d( idest, im, jm, buf, ireq )
    
    integer, intent(in)  :: idest
    integer, intent(in)  :: im, jm
    real(8), intent(in)  :: buf( im, jm )
    integer, intent(out) :: ireq
    integer :: ierr
    integer :: isize

    isize = size( buf )
    call mpi_isend( buf, isize, mpi_double_precision, idest, 1, mpi_comm_world, ireq, ierr )

  end subroutine mpii_isend_dble_2d

  subroutine mpii_isend_dble_3d( idest, im, jm, km, buf, ireq )

    integer, intent(in)  :: idest
    integer, intent(in)  :: im, jm, km
    real(8), intent(in)  :: buf( im, jm, km )
    integer, intent(out) :: ireq
    integer :: ierr
    integer :: isize

    isize = size( buf )
    call mpi_isend( buf, isize, mpi_double_precision, idest, 1, mpi_comm_world, ireq, ierr )
    
  end subroutine mpii_isend_dble_3d

  subroutine mpii_isend_dble_4d( idest, im, jm, km, lm, buf, ireq )

    integer, intent(in)  :: idest
    integer, intent(in)  :: im, jm, km, lm
    real(8), intent(in)  :: buf( im, jm, km, lm )
    integer, intent(out) :: ireq
    integer :: ierr
    integer :: isize

    isize = size( buf )
    call mpi_isend( buf, isize, mpi_double_precision, idest, 1, mpi_comm_world, ireq, ierr )
    
  end subroutine mpii_isend_dble_4d

  subroutine mpii_recv_dble_3d( idep, im, jm, km, buf )

    integer, intent(in)  :: idep
    integer, intent(in)  :: im, jm, km
    real(8), intent(out) :: buf( im, jm, km )
    integer :: ierr
    integer :: isize
    
    isize = size( buf )
    call mpi_irecv( buf, isize, mpi_double_precision, idep, 1, mpi_comm_world, ierr )
    
  end subroutine mpii_recv_dble_3d

  subroutine mpii_send_dble_3d( idest, im, jm, km, buf)

    integer, intent(in)  :: idest
    integer, intent(in)  :: im, jm, km
    real(8), intent(in)  :: buf( im, jm, km )
    integer :: ierr
    integer :: isize

    isize = size( buf )
    call mpi_send( buf, isize, mpi_double_precision, idest, 1, mpi_comm_world, ierr )
    
  end subroutine mpii_send_dble_3d

  subroutine mpii_wait( ireq )
    
    integer, intent(in) :: ireq
    integer :: ierr
    integer :: istatus( MPI_STATUS_SIZE )
    
    call mpi_wait( ireq, istatus, ierr )
    
  end subroutine mpii_wait