/[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.7 by enderton, Mon Sep 12 21:11: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                try [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                catch, end
50            end
51      else      else
52          [cs,h]=contour(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint);          [cs,h] = contour(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint);
53          set(h,'linewidth',linewidth)          set(h,'linewidth',linewidth);
54      end          if UseConLabel
55      if UseConLabel              clabel(cs,h,'fontsize',fs_clabel,'rotation',0);
56          clabel(cs,h,'fontsize',fs_clabel,'rotation',0);          end
57      end      end
58            
59        
60  % Make contour (filled) plot.  Let MATLAB determine the contour  % Make contour (filled) plot.  Let MATLAB determine the contour
61  % intervals on a differencing plot.  % intervals on a differencing plot.
62  elseif isequal(pst,'Cnf')  elseif isequal(pst,'Cnf')
63      if isequal(cmp,'Dif')      try   [cs,h] = contourf(xax{inrow}{incol},yax{inrow}{incol},plotdata,contint);
64          [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  
65      if UseCnfLabel      if UseCnfLabel
66          clabel(cs,h,'fontsize',fs_clabel,'rotation',0);          clabel(cs,h,'fontsize',fs_clabel,'rotation',0);
67      end      end
# Line 60  elseif isequal(pst,'Cnf') Line 70  elseif isequal(pst,'Cnf')
70  % appropriate data of which to overlay.  % appropriate data of which to overlay.
71  elseif isequal(pst,'Lin')  elseif isequal(pst,'Lin')
72      if ismember(cmp,{'OvC','OvE','OvF'})      if ismember(cmp,{'OvC','OvE','OvF'})
73          for inexp = 1:nexp          for intrl = 1:ntrl
74              plot(xax{inrow}{inexp},data{inrow}{inexp},linecolors{inexp},'linewidth',linewidth);              plot(xax{inrow}{intrl},data{inrow}{intrl},...
75                     linecolors{intrl},'linewidth',linewidth);
76          end          end
77      else      else
78          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.7
  ViewVC Help
Powered by ViewVC 1.1.22