
      open(unit=11,file='mesh.dat',status='unknown')
c maryam

      nim1=169 
      njm1=193 ! for grid 3
C      njm1=265 ! for grid 5
      nkm1=65

c maryam
      do j=1,njm1
      do i=1,nim1
        read (11,*) xc(i,j,1),yc(i,j,1)
      end do
      end do

      close(11)
