!&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
subroutine initMPI
!***********************************************************************
!                                                                      !
!     Initialize mpi                                                   !
!                                                                      !
!***********************************************************************
USE mppstaff_comm
!-----------------------------------------------------------------------
implicit none
include "mpif.h"
integer::ierr
!-----------------------------------------------------------------------
!***
!***  Initial MPI calls
!***
      
      call MPI_INIT(ierr)
      call MPI_COMM_RANK(MPI_COMM_WORLD,mype,ierr)
      call MPI_COMM_SIZE(MPI_COMM_WORLD,npes,ierr)

      itype = MPI_REAL
      dgnsnum=npes-1
!      print*,'npes=',npes
      dgnstag=101

!       print *,mype,dgnsnum,dgnstag
!***
!***  Numbers
!***  
      listInVrbls = 10
      listInConst = 11
 !!!!!!!!     listOut     = 21!!!!!!!!!!!!
      listOut     = 98
!***
!***  Chars
!***  
      write(c_mype,1000) mype
!!!!!!! 1000 format(i3.3)
    1000 format(i4.4)    
!-----------------------------------------------------------------------
!
      call MPI_BARRIER(MPI_COMM_WORLD,ierr)
!
!-----------------------------------------------------------------------
end subroutine initMPI