beaufort/input/mk_domain.m

cd /home/dimitri/cube/grid/cube66
d=read_cs_face('Depth.data',3);
yc=read_cs_face('YC.data',3);     

clf, mypcolor(d'), caxis([0 1]), grid, hold on
contour(yc,50:10:90,'k')
plot([326 326 365 365 326],[311 350 350 311 311],'k','linewidth',2) 

print -djpeg tmp2
1	cd /home/dimitri/cube/grid/cube66
2	d=read_cs_face('Depth.data',3);
3	yc=read_cs_face('YC.data',3);
4
5	clf, mypcolor(d'), caxis([0 1]), grid, hold on
6	contour(yc,50:10:90,'k')
7	plot([326 326 365 365 326],[311 350 350 311 311],'k','linewidth',2)
8
9	print -djpeg tmp2