/[MITgcm]/MITgcm_contrib/enderton/Diagnostics/DiagPlotMakePlot.m
ViewVC logotype

Diff of /MITgcm_contrib/enderton/Diagnostics/DiagPlotMakePlot.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by enderton, Mon Jan 31 15:43:27 2005 UTC revision 1.6 by enderton, Mon Sep 12 20:13:21 2005 UTC
# Line 16  end Line 16  end
16  % (and in the horizontal for that matter).  % (and in the horizontal for that matter).
17  if isequal(pst,'Grd')  if isequal(pst,'Grd')
18      if isCS      if isCS
19          merccube(xax{inrow}{incol},yax{inrow}{incol},plotdata);          merccube_mod(xax{inrow}{incol},yax{inrow}{incol},plotdata);
20      else      else
21          imagesc(xax{inrow}{incol},yax{inrow}{incol},plotdata);          imagesc(xax{inrow}{incol},yax{inrow}{incol},plotdata);
22      end      end
# Line 24  if isequal(pst,'Grd') Line 24  if isequal(pst,'Grd')
24  % Make interpolated plot.  % Make interpolated plot.
25  elseif isequal(pst,'Int')  elseif isequal(pst,'Int')
26      if isCS      if isCS
27          merccube(xax{inrow}{incol},yax{inrow}{incol},plotdata);          merccube_mod(xax{inrow}{incol},yax{inrow}{incol},plotdata);
28          shading interp;          shading interp;
29      else      else
30          pcolor(xax{inrow}{incol},yax{inrow}{incol},plotdata);          pcolor(xax{inrow}{incol},yax{inrow}{incol},plotdata);
# Line 34  elseif isequal(pst,'Int') Line 34  elseif isequal(pst,'Int')
34  % Make contour (non-filled) plot.  Let MATLAB determine the contour  % Make contour (non-filled) plot.  Let MATLAB determine the contour
35  % intervals on a differencing plot.  % intervals on a differencing plot.
36  elseif isequal(pst,'Con')  elseif isequal(pst,'Con')
37      if isequal(cmp,'Dif')      if isequal(cmap,'solid-dashed')
38          [cs,h] = contour(xax{inrow}{incol},yax{inrow}{incol},plotdata);          for ii = 1:3
39                cint = contint; lw = linewidth;
40                if ii == 1, cint = cint(cint>0); lst = '-';  end
41                if ii == 2, cint = cint(cint<0); lst = '--'; end
42                if ii == 3, cint = [0,0]; lst = '-'; lw = 2.*lw; end
43                [cs,h] = contour(xax{inrow}{incol},yax{inrow}{incol},...
44                                 plotdata,cint,['k',lst]);
45                set(h,'linewidth',lw);
46                if UseConLabel
47                    clabel(cs,h,'fontsize',fs_clabel,'rotation',0);
48                end
49            end
50      else      else
51          [cs,h]=contour(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint);          [cs,h] = contour(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint);
52          set(h,'linewidth',linewidth)          set(h,'linewidth',linewidth);
53      end          if UseConLabel
54      if UseConLabel              clabel(cs,h,'fontsize',fs_clabel,'rotation',0);
55          clabel(cs,h,'fontsize',fs_clabel,'rotation',0);          end
56      end      end
57            
58        
59  % Make contour (filled) plot.  Let MATLAB determine the contour  % Make contour (filled) plot.  Let MATLAB determine the contour
60  % intervals on a differencing plot.  % intervals on a differencing plot.
61  elseif isequal(pst,'Cnf')  elseif isequal(pst,'Cnf')
62      if isequal(cmp,'Dif')      try   [cs,h] = contourf(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint);
63          [cs,h] = contourf(xax{inrow}{incol},yax{inrow}{incol},plotdata);      catch [cs,h] = contourf(xax{inrow}{incol},yax{inrow}{incol},plotdata); end
     else  
         [cs,h] = contourf(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint);  
     end  
64      if UseCnfLabel      if UseCnfLabel
65          clabel(cs,h,'fontsize',fs_clabel,'rotation',0);          clabel(cs,h,'fontsize',fs_clabel,'rotation',0);
66      end      end
# Line 60  elseif isequal(pst,'Cnf') Line 69  elseif isequal(pst,'Cnf')
69  % appropriate data of which to overlay.  % appropriate data of which to overlay.
70  elseif isequal(pst,'Lin')  elseif isequal(pst,'Lin')
71      if ismember(cmp,{'OvC','OvE','OvF'})      if ismember(cmp,{'OvC','OvE','OvF'})
72          for inexp = 1:nexp          for intrl = 1:ntrl
73              plot(xax{inrow}{inexp},data{inrow}{inexp},linecolors{inexp},'linewidth',linewidth);              plot(xax{inrow}{intrl},data{inrow}{intrl},...
74                     linecolors{intrl},'linewidth',linewidth);
75          end          end
76      else      else
77          plot(xax{inrow}{incol},plotdata,linecolors{1},'linewidth',linewidth);          plot(xax{inrow}{incol},plotdata,linecolors{1},'linewidth',linewidth);

Legend:
Removed from v.1.1  
changed lines
  Added in v.1.6

  ViewVC Help
Powered by ViewVC 1.1.22