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

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

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

revision 1.10 by enderton, Wed Aug 10 20:00:35 2005 UTC revision 1.13 by enderton, Mon Sep 5 18:48:27 2005 UTC
# Line 53  if ~isequal(Grads,0) Line 53  if ~isequal(Grads,0)
53      if DiagDebug, disp(['  Debug -- Loading GRADS field.']); end      if DiagDebug, disp(['  Debug -- Loading GRADS field.']); end
54      [data,xax,yax,zax,months,time,dim] = ...      [data,xax,yax,zax,months,time,dim] = ...
55          DiagLoadGradsData(Grads,dad,fln,DiagDebug);          DiagLoadGradsData(Grads,dad,fln,DiagDebug);
56      if DiagDebug, disp(['  Debug -- ''data'' size after load:  ',mat2str(size(data))]); end      if DiagDebug, disp(['  Debug -- ''data'' size after load:  ',...
57                            mat2str(size(data))]); end
58      data = DiagAverage(data,fln,avg,months,ddf,dim);      data = DiagAverage(data,fln,avg,months,ddf,dim);
59      if DiagDebug, disp(['  Debug -- ''data'' size after averaging:  ',mat2str(size(data))]); end      if DiagDebug, disp(['  Debug -- ''data'' size after averaging:  ',...
60      if isequal(slc,'Sur')                           mat2str(size(data))]); end
61          data = data';      [data,xax,yax,pltslc] = ...
62          pltslc='lonlat';          DiagSliceGradsData(fln,flu,slc,data,xax,yax,zax);
63      else      if DiagDebug, disp(['  Debug -- ''data'' size after slicing:  ',...
64          error('Loading grads data only set to handle surface plots!');                           mat2str(size(data))]); end
     end  
65    
66  % Load monitor data.  % Load monitor data.
67  elseif isequal(dat(1:2),'Mn')  elseif isequal(dat(1:2),'Mn')
# Line 189  elseif ismember(fln,{'U','V','uVel','vVe Line 189  elseif ismember(fln,{'U','V','uVel','vVe
189      V = DiagAverage(V,fln,avg,months,ddf,Dim);      V = DiagAverage(V,fln,avg,months,ddf,Dim);
190          [XC,XG,YC,YG,Ylat,ZC,ZG,RAC,drC,drF,HFacC,HFacW,HFacS,dxG,dyG,dxC,dyC] = ...          [XC,XG,YC,YG,Ylat,ZC,ZG,RAC,drC,drF,HFacC,HFacW,HFacS,dxG,dyG,dxC,dyC] = ...
191          DiagLoadGridData(LoadGridData,grd,gdf,flu,GridSuffix,ZcordFile);          DiagLoadGridData(LoadGridData,grd,gdf,flu,GridSuffix,ZcordFile);
192          %[U,V]=uvcube2latlon(XC,YC,U,V,XL,YL);          [U,V]=uvcube2latlon(XC,YC,U,V,XL,YL);
193      if ismember(fln,{'U','uVel','fizhi_U'})      if ismember(fln,{'U','uVel','fizhi_U'})
194          data = U;          data = U;
195      elseif ismember(fln,{'V','vVel','fizhi_V'})      elseif ismember(fln,{'V','vVel','fizhi_V'})
# Line 370  else Line 370  else
370       if DiagDebug, disp(['  Debug -- ''data'' size after vector manipulation:  ',mat2str(size(data))]); end       if DiagDebug, disp(['  Debug -- ''data'' size after vector manipulation:  ',mat2str(size(data))]); end
371      else      else
372  % Sequence here for Vector = 0  % Sequence here for Vector = 0
373       if ~isequal(Index,0)       if ~isequal(Index,0) & isequal(size(Index),[1 1]) & isequal(sign(Index),1)
374          if     isequal(Dim,2), data = squeeze(data(:,:,Index,:));          if     isequal(Dim,2), data = squeeze(data(:,:,Index,:));
375          elseif isequal(Dim,3), data = squeeze(data(:,:,:,Index,:)); end          elseif isequal(Dim,3), data = squeeze(data(:,:,:,Index,:)); end
376         elseif isequal(size(Index),[1 1]) & ~isequal(sign(Index),1)
377            if     isequal(Dim,2), data = sign(Index(1)).*squeeze(data(:,:,abs(Index(1)),:));
378            elseif isequal(Dim,3), data = sign(Index(1)).*squeeze(data(:,:,:,abs(Index(1)),:)); end
379         elseif isequal(size(Index),[1 2])
380            if     isequal(Dim,2), data = sign(Index(1)).*squeeze(data(:,:,abs(Index(1)),:)) + sign(Index(2)).*squeeze(data(:,:,abs(Index(2)),:));
381            elseif isequal(Dim,3), data = sign(Index(1)).*squeeze(data(:,:,:,abs(Index(1)),:)) + sign(Index(2)).*squeeze(data(:,:,:,abs(Index(2)),:)); end
382       end       end
383       if DiagDebug, disp(['  Debug -- ''data'' size after indexing:  ',mat2str(size(data))]); end       if DiagDebug, disp(['  Debug -- ''data'' size after indexing:  ',mat2str(size(data))]); end
384       data = DiagAverage(data,fln,avg,months,ddf,Dim);       data = DiagAverage(data,fln,avg,months,ddf,Dim);
# Line 411  function data = DiagLoad_Local(fln,dat,d Line 417  function data = DiagLoad_Local(fln,dat,d
417          data = rdmds([dad,'/',fln,filesuffix],itr);          data = rdmds([dad,'/',fln,filesuffix],itr);
418          elseif isequal(dfm,'MNC')          elseif isequal(dfm,'MNC')
419          data = rdmnc_mod2([dad,mnchandle],[fln,filesuffix],'iter','T',itr);          data = rdmnc_mod2([dad,mnchandle],[fln,filesuffix],'iter','T',itr);
420          if ~isequal(itr,data.iter'), error('Missing iterations in data!'); end          %if ~isequal(itr,data.iter'), error('Missing iterations in data!'); end
421          eval(['data = data.',fln,filesuffix,';']);          eval(['data = data.',fln,filesuffix,';']);
422          else          else
423          error(['Unrecognized data type:  ',dfm]);          error(['Unrecognized data type:  ',dfm]);
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22