--- MITgcm_contrib/enderton/Diagnostics/DiagPlotMakePlot.m 2005/01/31 15:43:27 1.1 +++ MITgcm_contrib/enderton/Diagnostics/DiagPlotMakePlot.m 2005/08/11 14:40:22 1.4 @@ -16,7 +16,7 @@ % (and in the horizontal for that matter). if isequal(pst,'Grd') if isCS - merccube(xax{inrow}{incol},yax{inrow}{incol},plotdata); + merccube_mod(xax{inrow}{incol},yax{inrow}{incol},plotdata); else imagesc(xax{inrow}{incol},yax{inrow}{incol},plotdata); end @@ -24,7 +24,7 @@ % Make interpolated plot. elseif isequal(pst,'Int') if isCS - merccube(xax{inrow}{incol},yax{inrow}{incol},plotdata); + merccube_mod(xax{inrow}{incol},yax{inrow}{incol},plotdata); shading interp; else pcolor(xax{inrow}{incol},yax{inrow}{incol},plotdata); @@ -34,12 +34,9 @@ % Make contour (non-filled) plot. Let MATLAB determine the contour % intervals on a differencing plot. elseif isequal(pst,'Con') - if isequal(cmp,'Dif') - [cs,h] = contour(xax{inrow}{incol},yax{inrow}{incol},plotdata); - else - [cs,h]=contour(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint); - set(h,'linewidth',linewidth) - end + 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); end @@ -47,11 +44,8 @@ % Make contour (filled) plot. Let MATLAB determine the contour % intervals on a differencing plot. elseif isequal(pst,'Cnf') - if isequal(cmp,'Dif') - [cs,h] = contourf(xax{inrow}{incol},yax{inrow}{incol},plotdata); - else - [cs,h] = contourf(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint); - end + try [cs,h] = contourf(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint); + catch [cs,h] = contourf(xax{inrow}{incol},yax{inrow}{incol},plotdata); end if UseCnfLabel clabel(cs,h,'fontsize',fs_clabel,'rotation',0); end @@ -60,8 +54,8 @@ % appropriate data of which to overlay. elseif isequal(pst,'Lin') if ismember(cmp,{'OvC','OvE','OvF'}) - for inexp = 1:nexp - plot(xax{inrow}{inexp},data{inrow}{inexp},linecolors{inexp},'linewidth',linewidth); + for intrl = 1:ntrl + plot(xax{inrow}{intrl},data{inrow}{intrl},linecolors{intrl},'linewidth',linewidth); end else plot(xax{inrow}{incol},plotdata,linecolors{1},'linewidth',linewidth);