/[MITgcm]/MITgcm_contrib/gael/profilesMatlabProcessing/profiles_misc/ncaddVar.m
ViewVC logotype

Diff of /MITgcm_contrib/gael/profilesMatlabProcessing/profiles_misc/ncaddVar.m

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

revision 1.1 by roquet, Wed Apr 13 20:49:14 2011 UTC revision 1.4 by roquet, Fri Apr 22 21:31:02 2011 UTC
# Line 1  Line 1 
1  function ncaddVar(ncid,varname,xtype,dimlist);  function ncaddVar(ncid,varname,xtype,dimlist);
2  % add a variable in an existing netcdf file.  % function ncaddVar(ncid,varname,xtype,dimlist);
3    %   add a variable in an existing netcdf file.
4    %   The netcdf file must be open in 'write' mode.
5    
6  global useNativeMatlabNetcdf;  global useNativeMatlabNetcdf;
7  if isempty(useNativeMatlabNetcdf); useNativeMatlabNetcdf = ~isempty(which('netcdf.open')); end;  if isempty(useNativeMatlabNetcdf); useNativeMatlabNetcdf = ~isempty(which('netcdf.open')); end;
8    
9  if useNativeMatlabNetcdf;  if useNativeMatlabNetcdf;
10      if isempty(dimlist), error('nddefVar error: no dimensions allocated'); end      if isempty(dimlist), error('ncaddVar error: no dimension allocated'); end
11      iDim=[];      iDim=[];
12      for ii=1:length(dimlist),      for ii=1:length(dimlist),
13          iDim(ii)=netcdf.inqDimID(ncid,dimlist{ii});          iDim(ii)=netcdf.inqDimID(ncid,dimlist{ii});
# Line 18  else;%try to use old mex stuff Line 20  else;%try to use old mex stuff
20      dimlist=fliplr(dimlist);      dimlist=fliplr(dimlist);
21      switch length(dimlist)      switch length(dimlist)
22          case 1,          case 1,
23              sprintf('ncid{''%s''}=nc%s(''%s'');',varname,xtype,dimlist{1});              eval(sprintf('ncid{''%s''}=nc%s(''%s'');',varname,xtype,dimlist{1}));
24          case 2,          case 2,
25              sprintf('ncid{''%s''}=nc%s(''%s'',''%s'');',varname,xtype,dimlist{1},dimlist{2});              eval(sprintf('ncid{''%s''}=nc%s(''%s'',''%s'');',varname,xtype,dimlist{1},dimlist{2}));
26          case 3,          case 3,
27              sprintf('ncid{''%s''}=nc%s(''%s'',''%s'',''%s'');',varname,xtype,dimlist{1},dimlist{2},dimlist{3});              eval(sprintf('ncid{''%s''}=nc%s(''%s'',''%s'',''%s'');',varname,xtype,dimlist{1},dimlist{2},dimlist{3}));
28          otherwise          otherwise
29              error('ncdefVar: number of dimension > 3');              error('ncaddVar: number of dimension > 3');
30      end      end
31  end;  end;
32    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22