%  add potential density lines to theta-s diagram
hold on

%xa = get(gca,'XLim');
%ya = get(gca,'YLim');

x = [32:0.05:37];
y = [-2:0.1:28];

[X,Y]=meshgrid(x,y);

hold on

Z=sw_dens(X,Y,0); % NOTE THE REFERENCE PRESSURE
Z=Z-1000;
%cs = contour(x,y,Z,[26:0.02:29],'k:');
cs = contour(x,y,Z,[22:0.5:30],'k-');
clabel(cs,[22:1:29],'fontsize',16);
%clabel(cs,'manual','fontsize',14);
axis([32 37 -2 28])
