*
* Define the size of the computational domain.  Since the boundary is 
* fixed, Computations are only performed INSIDE the boundary, so
* first and last zone in each direction are not computed on by the 
* CMStab solver.
*

      integer idim,jdim,kdim
      parameter (idim=132+2,jdim=132+2,kdim=132+2)
