nx=56;
ny=110;
dx=15;
dy=15;
for j=1:ny+1;
    y(1:nx+1,j)=(j-1)*dy;
end
for i=1:nx+1
    x(i,1:ny+1)=(i-1)*dx;
end
zb=-13.6+13.6/(765-90)*(x-90);
for j=56:60;
   zb(:,j)=-13.6+12.3/(765-90)*(x(:,j)-90);
   zb(:,j)=min(zb(:,j),-1.3);
end
zb=max(-13.6,zb);
zb(:,2:ny)=.25*zb(:,1:ny-1)+.5*zb(:,2:ny)+.25*zb(:,3:ny+1);
figure(1);
pcolor(x,y,zb);shading interp;colorbar;axis equal
hold on
v=[-15:2:1];
contour(x,y,zb,v,'k')
fid=fopen('outflow.dep','wt')
for j=1:ny+1
    fprintf(fid,'%7.3f ',zb(:,j));fprintf(fid,'\n')
end
fclose(fid)