--- MITgcm/verification/front_relax/input/gendata.m	2000/04/10 14:38:01	1.1
+++ MITgcm/verification/front_relax/input/gendata.m	2000/04/10 16:14:07	1.2
@@ -20,9 +20,17 @@
 Lx=dx*nx;
 Ly=dy*(ny-1);  % Solid wall in North
 
+% Variable resolution
+y=(1:ny)/(ny-0)-0.5;
+dy=1-0.3*exp( -(5*y).^2 );
+%dy=ones(1,ny);              % Constant resolution
+dy=dy/sum(dy(1:ny-1))*Ly;
+fid=fopen('dy.bin','w',ieee); fwrite(fid,dy,prec); fclose(fid);
+
 % Coordinates
 xc=((1:nx)-0.5)*dx;
-yc=((1:ny)-0.5)*dy-Ly/2;
+yf=-Ly/2+[0 cumsum(dy)];
+yc=(yf(1:end-1)+yf(2:end))/2;
 zf=[0 -cumsum(dz)];
 zc=(zf(1:end-1)+zf(2:end))/2;
 [X,Y]=ndgrid(xc,yc);