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

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

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

revision 1.5 by molod, Fri Apr 29 17:34:39 2005 UTC revision 1.6 by enderton, Wed Aug 10 19:58:03 2005 UTC
# Line 42  for inrow = 1:nrow Line 42  for inrow = 1:nrow
42            
43      ntrl = length(page{inrow}); if ntrl ~= 1, ntrl = ntrl - 1; end      ntrl = length(page{inrow}); if ntrl ~= 1, ntrl = ntrl - 1; end
44      if ntrl == 1, cmp = 'Sep'; else, cmp = page{inrow}{end}; end      if ntrl == 1, cmp = 'Sep'; else, cmp = page{inrow}{end}; end
45      if ntrl == 1, ncol = 1; elseif cmp == 'Sbs', ncol = ntrl; else ncol = 1; end      if ntrl == 1, ncol = 1; elseif isequal(cmp,'Sbs'), ncol = ntrl; else ncol = 1; end
46            
47      dx = (1-dxl-dxr-(ncol-1)*dxm)/ncol;      dx = (1-dxl-dxr-(ncol-1)*dxm)/ncol;
48      dy = (1-dyb-dyt-(nrow-1)*dym)/nrow;      dy = (1-dyb-dyt-(nrow-1)*dym)/nrow;
49            
50      for incol = 1:ncol      for incol = 1:ncol
51                    
52          if size(data{inrow}{incol}) == [6*hres,hres], isCS = 1; else, isCS = 0; end          % Very crude test to see if cube sphere, must be fixed!
53            test = size(data{inrow}{incol});
54            if test(1)./test(2) == 6, isCS = 1; else, isCS = 0; end
55          fln = page{inrow}{incol}{ifln};          fln = page{inrow}{incol}{ifln};
56          pst = page{inrow}{incol}{ipst};          pst = page{inrow}{incol}{ipst};
57          flu = page{inrow}{incol}{iflu};          flu = page{inrow}{incol}{iflu};
# Line 75  for inrow = 1:nrow Line 77  for inrow = 1:nrow
77                '; Cmp: ',cmp]);                '; Cmp: ',cmp]);
78                    
79          % Load contour intervals, units (loaded in 'DiagFieldParam[A,O]').          % Load contour intervals, units (loaded in 'DiagFieldParam[A,O]').
80          try, contint; catch          if ~isequal(cmp,'Dif')
81              try, eval(['contint = ',fln,'contour',flu,';']);              try, contint; catch
82              catch, disp(['***Warning***  No contour information for ',fln]);                  try, eval(['contint = ',fln,'contour',flu,';']);
83                     disp(['               Using 10 generic contour levels.']);                  catch, disp(['***Warning***  No contour information for ',fln]);
84                     contint = 10; end, end                         disp(['               Using 10 generic contour levels.']);
85                           contint = 10; end, end
86            end
87          try, units; catch          try, units; catch
88              try, eval(['units = ',fln,'units',flu,';']);              try, eval(['units = ',fln,'units',flu,';']);
89              catch, disp(['***Warning***  No unit information found for ',fln]);              catch, disp(['***Warning***  No unit information found for ',fln]);
90                     disp(['               Using question mark.']);                     disp(['               Using question mark.']);
91                     units = '?'; end, end                     units = '?'; end, end
92                  
93          xi = dxl + (incol-1)*(dx+dxm);          xi = dxl + (incol-1)*(dx+dxm);
94          yi = 1-dyt-inrow*dy-(inrow-1)*dym;          yi = 1-dyt-inrow*dy-(inrow-1)*dym;
95                    
# Line 109  for inrow = 1:nrow Line 113  for inrow = 1:nrow
113          DiagPlotResetAxes;          DiagPlotResetAxes;
114          DiagPlotMisc;          DiagPlotMisc;
115          DiagPlotTitles;          DiagPlotTitles;
116          clear contint units          clear contint units crange
117      end      end
118  end  end
119    
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22