!----------------------------------------------------------------------
!  This routine initializes the following grid functions
!  u ures 
!----------------------------------------------------------------------
      subroutine initializer0(u_n,ures_n,g1_Nx,g1_Ny,x,y,amp,xc,xwid,yc,
     &  ywid)
        implicit none

        include 'globals.inc'

        integer g1_Nx,g1_Ny
        real*8  u_n(1:g1_Nx,1:g1_Ny)
        real*8  ures_n(1:g1_Nx,1:g1_Ny)
        real*8  x(*)
        real*8  y(*)
        real*8  amp
        real*8  xc
        real*8  xwid
        real*8  yc
        real*8  ywid
        integer i,j,k

          do j=1, g1_Ny, 1
          do i=1, g1_Nx, 1
          u_n(i,j)=amp*exp(-((x(i)-xc)**2/xwid**2+(y(j)-yc)**2/ywid**2))
          end do
          end do
          do j=1, g1_Ny, 1
          do i=1, g1_Nx, 1
          ures_n(i,j)=0
          end do
          end do
        return
      end