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