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)