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

Diff of /MITgcm_contrib/enderton/Diagnostics/DiagSlice.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.4 by molod, Tue Jun 28 21:33:52 2005 UTC
# Line 1  Line 1 
1  function [data,xax,yax,pltslc] = ...  function [data,xax,yax,pltslc] = ...
2      DiagSlice(data,fln,exp,dat,dad,grd,itr,tst,...      DiagSlice(data,fln,trl,dat,dad,grd,itr,tst,flu,ddf,gdf,...
3                flu,ddf,gdf,avg,slc,pst,LoadGridData,GridSuffix,ZcordFile);                avg,slc,pst,LoadGridData,GridSuffix,ZcordFile);
4    
5  % Function: DiagSlice  % Function: DiagSlice
6  % Author:   Daniel Enderton  % Author:   Daniel Enderton
# Line 102  elseif isequal(slc,'Zon') Line 102  elseif isequal(slc,'Zon')
102              data=data(:,:,1)'; xax=Ylat;   yax=ZC;   pltslc='lathgt';              data=data(:,:,1)'; xax=Ylat;   yax=ZC;   pltslc='lathgt';
103          end          end
104      elseif isequal(datasize(1:2),[length(XL),length(YL)])      elseif isequal(datasize(1:2),[length(XL),length(YL)])
105          if ismember(fln,{'U','V','uVel','vVel'})          if ~isequal(pst,'Lin')
106              data = squeeze(mean(data,1))'; xax=YL; yax=ZC;  pltslc='lathgt';              data = squeeze(mean(data,1))'; xax=YL; yax=ZC;  pltslc='lathgt';
         elseif ismember(fln,{'TX','TY','USTR','VSTR'})  
             data = squeeze(mean(data,1))'; xax=YL; yax=NaN; pltslc='latfld';  
107          else          else
108              error(['Unknown field for U,V type zonal average data:  ',fln]);              data = squeeze(mean(data,1))'; xax=YL; yax=NaN; pltslc='latfld';
109          end          end
110      else      else
111          error('Incorrect dimensions for slc = ''Zon''');          error('Incorrect dimensions for slc = ''Zon''');
# Line 150  elseif isequal(slc(1:2),'j=') Line 148  elseif isequal(slc(1:2),'j=')
148  elseif isequal(slc(1:2),'k=')  elseif isequal(slc(1:2),'k=')
149      kk = str2num(slc(3:end));      kk = str2num(slc(3:end));
150      data = squeeze(data(:,:,kk));      data = squeeze(data(:,:,kk));
151      if ismember(fln,{'U','V','uVel','vVel'})      if ismember(fln,{'U','V','uVel','vVel','fizhi_U','fizhi_V'})
152          data = data'; xax = XL; yax = YL;          data = data'; xax = XL; yax = YL;
153      elseif ismember(pst,{'Grd','Int'})      elseif ismember(pst,{'Grd','Int'})
154          if isequal(pst,'Grd'),          if isequal(pst,'Grd'),
# Line 187  try Line 185  try
185  catch  catch
186      disp(['***Warning***  No range information found for ',fln]);      disp(['***Warning***  No range information found for ',fln]);
187      disp(['               Data range:  ',mat2str(datarange)]);      disp(['               Data range:  ',mat2str(datarange)]);
 end  
188    end
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22