function makeGraph2(dOlambda)
dx    = 0.01;
zOlam = (-dOlambda:dx:dx);
xOlam = (0:dx:1);
xi    = 2*pi*xOlam;
zeta  = 2*pi*(zOlam+dOlambda)';
phi   = cosh(zeta)*cos(xi);

figure(1)  % the contour graph
plot(xOlam,-dx*sin(xi)) % plot surface wave
hold on
contour(xOlam,zOlam,phi,20) % plot contour
texText('$z/\lambda$','y')  % add lables
texText('$x/\lambda$','x')
axis equal
hold off
name = sprintf('contour%03d.pdf',100*dOlambda);
graphToPDF(name,1+4,1+4*dOlambda)

figure(2) % the contour graph with velocity field.
dx    = 0.05;
plot(xOlam,-0.01*sin(xi),'k') % plot surface wave
hold on
contour(xOlam,zOlam,phi,20) % plot contour

zOlam = (-dOlambda:dx:0); % compute velocity field
xOlam = (0:dx:1);
xi    = 2*pi*xOlam;
zeta  = 2*pi*(zOlam+dOlambda)';
ux = -cosh(zeta)*sin(xi);
uz =  sinh(zeta)*cos(xi);
quiver(xOlam,zOlam,ux,uz)  % plot the velocity field
ylim([-dOlambda,0.01])
axis equal
hold off
name = sprintf('contourQuiver%03d.pdf',100*dOlambda);
graphToPDF(name,1+4,1+4*dOlambda)
end