!!!     Initialize MPI routines
      comm=MPI_COMM_WORLD
      call MPI_INIT(ierr)
      call MPI_COMM_RANK(comm,myid,ierr)
      call MPI_COMM_SIZE(comm,nprocs,ierr)
      
      if( myid .eq. 0 ) then
       write(0,*) ' MPI initialized with ',nprocs,' live processors'
       if( nprocs .ne. numprocs ) then
        write(0,*) ' Number of live processors: ',nprocs,' does not match'
        write(0,*) ' the number specified in variable_declarations.h: ',numprocs
        write(0,*) ' Check the run.sh script.'
        goto 999     ! shutdown MPI and exit gracefully
       endif
      endif