/[MITgcm]/MITgcm/utils/matlab/rdmnc.m
ViewVC logotype

Diff of /MITgcm/utils/matlab/rdmnc.m

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

revision 1.1 by adcroft, Mon Mar 8 16:45:56 2004 UTC revision 1.2 by mlosch, Wed Feb 2 10:31:22 2005 UTC
# Line 78  end Line 78  end
78    
79  % -----------------------------------------------------------------------------  % -----------------------------------------------------------------------------
80  function [result] = ismncfile(nc);  function [result] = ismncfile(nc);
81  result=~isempty(nc.('MITgcm_mnc_ver'));  result=~isempty(nc.MITgcm_mnc_ver);
82    %MLresult=~isempty(nc.('MITgcm_mnc_ver'));
83  % -----------------------------------------------------------------------------  % -----------------------------------------------------------------------------
84  function [A] = read_att(nc);  function [A] = read_att(nc);
85  global verbose  global verbose
86  allatt=ncnames(att(nc)); if verbose; allatt, end  allatt=ncnames(att(nc)); if verbose; allatt, end
87  A='none';  A='none';
88  for attr=allatt;  for attr=allatt;
89   A.(char(attr))=nc.(char(attr))(:);   tmp = char(attr);
90     eval(['A.' tmp '= nc.' tmp '(:);'])
91    %ML A.(char(attr))=nc.(char(attr))(:);
92  end  end
93  % -----------------------------------------------------------------------------  % -----------------------------------------------------------------------------
94  function [i0,j0,fn] = findTileOffset(S);  function [i0,j0,fn] = findTileOffset(S);
# Line 126  for ivar=1:size(varlist,2); Line 129  for ivar=1:size(varlist,2);
129    disp(['No such variable ''' cvar ''' in netcdf file' name(nc)])    disp(['No such variable ''' cvar ''' in netcdf file' name(nc)])
130   else   else
131    tmpdata=squeeze(permute(tmpdata,[9:-1:1]));    tmpdata=squeeze(permute(tmpdata,[9:-1:1]));
132    S.(cvar)=tmpdata;    eval(['S.' cvar '=tmpdata;'])
133    S.attributes.(cvar)=read_att(nc{cvar});    eval(['S.attributes.' cvar '=read_att(nc{' cvar '});'])
134    %ML  S.(cvar)=tmpdata;
135    %ML  S.attributes.(cvar)=read_att(nc{cvar});
136   end   end
137  end  end
138  % -----------------------------------------------------------------------------  % -----------------------------------------------------------------------------
# Line 188  for ivar=1:size(varlist,2); Line 193  for ivar=1:size(varlist,2);
193    else    else
194     j0=0;     j0=0;
195    end    end
196    S.(cvar)(i0+(1:ni),j0+(1:nj),(1:nk),(1:nm),(1:nn),(1:no),(1:np))=tmpdata;    eval(['S.' cvar ...
197    S.attributes.(cvar)=read_att(nc{cvar});          '(i0+(1:ni),j0+(1:nj),(1:nk),(1:nm),(1:nn),(1:no),(1:np))=tmpdata;'])
198      eval(['S.attributes.' cvar ' =read_att(nc{''' cvar '''});'])
199    %ML  S.(cvar)(i0+(1:ni),j0+(1:nj),(1:nk),(1:nm),(1:nn),(1:no),(1:np))=tmpdata;
200    %ML  S.attributes.(cvar)=read_att(nc{cvar});
201   end   end
202  end  end
203    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22