subroutine initializeFlow(
     *     u,  nX, nY, 
     *     nMax, x, 
     *     y)
      implicit none
      real*8 u, x, y

      real*8 iniUVal, iniVVal

      integer nX, nY, nMax
      integer i,j

      dimension u(2,0:nMax+1,0:nMax+1)
      dimension x(0:nMax+1)
      dimension y(0:nMax+1)

      do j=0,nMax+1
         do i=0,nMax+1
            u(1,i,j) = iniUVal(
     *           X(i), Y(j))
            u(2,i,j) = iniVVal(
     *           X(i), Y(j))
         enddo
      enddo

      return 
      end