      dimension z(100,100)
      ix=100
      iy=100
      x0=14.0
      y0=20.0
      dx=x0/(ix-1)
      dy=y0/(iy-1)
      
      do 10 i=1,ix
      do 10 j=1,iy
      x=(i-1)*dx-4.0
      y=(j-1)*dy-10.0
      z(i,j)=exp(0.2*x)*sin(y)
   10 continue
      
      open(99,file="2dlat_sample2.dat")
      write(99,*) ix,iy
      do 20 j=1,iy
      write(99,*) (z(i,j),i=1,ix)
   20 continue
      
      stop
      end
      