% depthmap.m   Makes a nice map of the depth, for overheading
clf

ax = [-2500 20000 5000 25000]/1000;  
ax = [14 26 5 15];
ax = [ 4 70 -26 39]

fieldmap(ax,meshele(:,:),meshnod,bathy(:,2),[0 20])
xlabel('km')
title('Wassaw Sound Depths, NOS data')

% shading interp
colorbar('vert')

kids = get(gcf,'children');
pos = get(kids(1),'position');
pos=pos*.8;
pos(1) = .2;
pos(2) = .2;
set(kids(1),'position',pos);
set(kids(1),'LineWidth',1.0)
set(kids(1),'FontWeight','bold')
set(kids(1),'fontsize',14)
h=get(kids(1),'xlabel');  
set(h,'fontsize',20)
set(h,'FontWeight','bold')
h=get(kids(1),'title');  
set(h,'fontsize',20)
set(h,'FontWeight','bold')


pos = get(kids(2),'position');
pos=pos*.8;
pos(1) = .8;
pos(2) = .2;
set(kids(2),'position',pos);

set(kids(2),'LineWidth',1.0)
set(kids(2),'FontWeight','bold')
set(kids(2),'fontsize',14)

h=get(kids(2),'ylabel');  
set(h,'string','Depth meters')
set(h,'fontsize',20)
set(h,'FontWeight','bold')
