!-list of variables that indicate indices in each spatial direction
!-(and more). Modified by PD to accomodate MPI on 1/26/05
!-nz = total number of points in z-direction (global)
!-1/6/03: Given that the penalty method is discontinuous
!-nz = nzloc*nsubd ...
!-nzloc = number of points in a spectral subdomain.
!-nzm = nzloc-1 (Used in Setup of Legendre interpolant functions etc.)
!-This is essentially the order of polynomial approximation in 
!-each subdomain.
!-nsubd = number of subdomains in vertical direction = M
!
!-nproch = number of processors in horizontal (y-dir. in physical space)
!-         space
!-nprocv = number of processors in vertical. Should be  a multiple
!-         of nsubd
!-nsubdproc = number of subdomains per processor !
!-
!-nzpl = number of grid points in each vertical processor subdomain
!-nypl = number of grid points in each y-dir. horizontal processor subdomain
!-nxpl = number of grid points in each x-dir. horizontal processor subdomain
!-(y-partition in physical space. x-partition in Fourier space)
!-CAREFUL: Processor subdomain is not necessarily the same with a
!-         physical subdomain !!!
!-------------------------------------------------------------------
!-Global grid parameters
      integer, parameter ::   nx =   16, ny =   16
!-Vertical subdomain grid parameters (global vertical resolution is determined by these)
      integer, parameter ::  nzloc = 5, nsubd = 13, nz = nzloc*nsubd
      integer, parameter :: nsdtype = 4
!-Auxilliary parameters
      integer, parameter ::  nxpp = nx + 2,     nxh=nx/2, nxhp = nxh+1,
     &              nxf = 3*nx/2 + 1,  nyf = 2*ny,  nyh=ny/2,
     &              nxppy = nxpp*ny,   ntot = nxppy*nz,  
     &              nzm = nzloc - 1
      integer, parameter :: nxhpy=nxhp*ny    
      integer, parameter :: nzfl=48,nzh = nzfl/2
!-Parameters associated with domain decomposition among processors
      integer, parameter :: nproch=2, nprocv=1
      integer, parameter :: nzpl = nz/nprocv, nypl = ny/nproch,
     &                      nxpl = nx/nproch, nxhpl = nxpl/2, 
     &                      nsubdproc = nsubd/nprocv
!-Number of dimensions (2-D or 3-D code)
      integer, parameter :: ndims=3


!-------------------------------------------------------------------