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

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

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

revision 1.2 by enderton, Mon Jan 31 19:58:15 2005 UTC revision 1.3 by molod, Fri Apr 29 17:34:39 2005 UTC
# Line 5  function [data,xax,yax,zax,months,time,d Line 5  function [data,xax,yax,zax,months,time,d
5  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
6  %                          Parse table file                               %  %                          Parse table file                               %
7  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
8    ShiftData = 0;
9  format = 'NONSEQUENTIAL';  format = 'NONSEQUENTIAL';
10    
11  tablfile = [dad,'/',Grads,'.tabl'];  tablfile = [dad,'/',Grads,'.tabl'];
# Line 27  for iline = 1:length(file) Line 27  for iline = 1:length(file)
27          inc = str2num(tokens{5});          inc = str2num(tokens{5});
28          xax = [ini:inc:ini+(num-1)*inc];          xax = [ini:inc:ini+(num-1)*inc];
29          if min(xax) >= 0 && max(xax) > 180          if min(xax) >= 0 && max(xax) > 180
30              xax = xax - 180;              ShiftData = 1;
31          end          end
32          nx = length(xax);          nx = length(xax);
33      end      end
# Line 147  if isequal(format,'SEQUENTIAL') Line 147  if isequal(format,'SEQUENTIAL')
147  end  end
148  data = reshape(data,[nx,ny,nz,nv,nt]);  data = reshape(data,[nx,ny,nz,nv,nt]);
149  data = squeeze(data(:,:,:,ivar,:));  data = squeeze(data(:,:,:,ivar,:));
 data( abs((data-undef)/undef) < tol ) = NaN;  
150    data( abs((data-undef)/undef) < tol ) = NaN;
151    if ShiftData
152        indexWestHemi = xax>=180;
153        indexEastHemi = xax<180;
154        data = cat(1,data(indexWestHemi,:,:),data(indexEastHemi,:,:));
155        xax = cat(2,xax(indexWestHemi)-360,xax(indexEastHemi));
156    end
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22