/[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.4 by molod, Tue Jun 28 21:33:52 2005 UTC revision 1.5 by molod, Tue Jul 5 18:57:48 2005 UTC
# Line 1  Line 1 
1  function [data,xax,yax,pltslc] = ...  function [data,xax,yax,pltslc] = ...
2      DiagSlice(data,fln,trl,dat,dad,grd,itr,tst,flu,ddf,gdf,...      DiagSlice(data,fln,trl,dat,dad,grd,itr,tst,flu,ddf,gdf,...
3                avg,slc,pst,LoadGridData,GridSuffix,ZcordFile);                avg,slc,pst,LoadGridData,GridSuffix,ZcordFile,Vector,FieldName);
4    
5  % Function: DiagSlice  % Function: DiagSlice
6  % Author:   Daniel Enderton  % Author:   Daniel Enderton
# Line 53  datasize = size(data); Line 53  datasize = size(data);
53  % since those handle the cubed sphere data.  If the 'pst' is 'Con' or 'Cnf',  % since those handle the cubed sphere data.  If the 'pst' is 'Con' or 'Cnf',
54  % convert data to a lat-lon grid so that it can easily be contours.  % convert data to a lat-lon grid so that it can easily be contours.
55  if isequal(slc,'Sur')  if isequal(slc,'Sur')
56      if ismember(fln,{'TX','TY','USTR','VSTR'})      if ismember(fln,{'TX','TY','USTR','VSTR'}) | ~isequal(Vector,0)
57          data = data'; xax = XL; yax = YL;          data = data'; xax = XL; yax = YL;
58      elseif ~isequal(datasize,[faces*hres,hres])      elseif ~isequal(datasize,[faces*hres,hres])
59          error('Incorrect dimensions for slc:  ',slc);          error('Incorrect dimensions for slc:  ',slc);
# Line 90  if isequal(slc,'Sur') Line 90  if isequal(slc,'Sur')
90  elseif isequal(slc,'Zon')  elseif isequal(slc,'Zon')
91      if isequal(datasize(1:2),[faces*hres,hres])      if isequal(datasize(1:2),[faces*hres,hres])
92          if isequal(flu,'O'), nBas = 0; end          if isequal(flu,'O'), nBas = 0; end
93            if isequal(flu,'A'), nBas = 0; end
94              [data,dump1,dump2] = ...              [data,dump1,dump2] = ...
95                  calc_ZonAv_CS(data,kpr,kwr,nBas,XC,YC,XG,YG,RAC,dad,HFacC);                  calc_ZonAv_CS(data,kpr,kwr,nBas,XC,YC,XG,YG,RAC,dad,HFacC);
96          if isequal(avg,'Tse')          if isequal(avg,'Tse')
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22