!*********************************************************************** 
!  Generates the spatial grid. 
!*********************************************************************** 
subroutine makegrid(params,h,x,y,z) 
implicit none 
real*8                       :: h
real*8,dimension(:)          :: params
real*8,dimension(:)          :: x,y,z
!     local variables 
integer                       :: i


   do i = 1, int(params(1)) 
      x(i) = params(15) + (i-1.0d0)*h 
   end do 
   do i = 1, int(params(2)) 
      y(i) = params(17) + (i-1.0d0)*h 
   end do 
   do i = 1, int(params(3)) 
      z(i) = params(19) + (i-1.0d0)*h 
   end do 
 
return 
end 
