C 
C  This file contains a routine for producing a decomposition of a 1-d array 
C  when given a number of processors.  It may be used in "direct" product 
C  decomposition.  The values returned assume a "global" domain in [1:n] 
C 
      subroutine MPE_DECOMP1D( n, numprocs, myid, s, e ) 
      integer n, numprocs, myid, s, e 
      integer nlocal 
      integer deficit 
C 
      nlocal  = n / numprocs 
      s	      = myid * nlocal + 1 
      deficit = mod(n,numprocs) 
      s	      = s + min(myid,deficit) 
      if (myid .lt. deficit) then 
          nlocal = nlocal + 1 
      endif 
      e = s + nlocal - 1 
      if (e .gt. n .or. myid .eq. numprocs-1) e = n 
      return 
      end