--- MITgcm_contrib/enderton/Diagnostics/DiagPlotMakePlot.m 2005/09/05 18:50:15 1.5 +++ MITgcm_contrib/enderton/Diagnostics/DiagPlotMakePlot.m 2005/09/12 20:13:21 1.6 @@ -34,13 +34,28 @@ % Make contour (non-filled) plot. Let MATLAB determine the contour % intervals on a differencing plot. elseif isequal(pst,'Con') - try [cs,h] = contour(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint); - catch [cs,h] = contour(xax{inrow}{incol},yax{inrow}{incol},plotdata); end - set(h,'linewidth',linewidth); - if UseConLabel - clabel(cs,h,'fontsize',fs_clabel,'rotation',0); + if isequal(cmap,'solid-dashed') + for ii = 1:3 + cint = contint; lw = linewidth; + if ii == 1, cint = cint(cint>0); lst = '-'; end + if ii == 2, cint = cint(cint<0); lst = '--'; end + if ii == 3, cint = [0,0]; lst = '-'; lw = 2.*lw; end + [cs,h] = contour(xax{inrow}{incol},yax{inrow}{incol},... + plotdata,cint,['k',lst]); + set(h,'linewidth',lw); + if UseConLabel + clabel(cs,h,'fontsize',fs_clabel,'rotation',0); + end + end + else + [cs,h] = contour(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint); + set(h,'linewidth',linewidth); + if UseConLabel + clabel(cs,h,'fontsize',fs_clabel,'rotation',0); + end end + % Make contour (filled) plot. Let MATLAB determine the contour % intervals on a differencing plot. elseif isequal(pst,'Cnf')