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

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

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


Revision 1.4 - (hide annotations) (download)
Fri Apr 22 21:31:02 2011 UTC (14 years, 2 months ago) by roquet
Branch: MAIN
CVS Tags: checkpoint65r, checkpoint65p, checkpoint65q
Changes since 1.3: +4 -3 lines
bug fix: old netcdf procedure

1 roquet 1.1 function ncaddVar(ncid,varname,xtype,dimlist);
2 roquet 1.3 % function ncaddVar(ncid,varname,xtype,dimlist);
3     % add a variable in an existing netcdf file.
4 roquet 1.4 % The netcdf file must be open in 'write' mode.
5 roquet 1.1
6     global useNativeMatlabNetcdf;
7     if isempty(useNativeMatlabNetcdf); useNativeMatlabNetcdf = ~isempty(which('netcdf.open')); end;
8    
9     if useNativeMatlabNetcdf;
10 roquet 1.2 if isempty(dimlist), error('ncaddVar error: no dimension allocated'); end
11 roquet 1.1 iDim=[];
12     for ii=1:length(dimlist),
13     iDim(ii)=netcdf.inqDimID(ncid,dimlist{ii});
14     end
15     netcdf.reDef(ncid);
16     netcdf.defVar(ncid,varname,xtype,iDim);
17     netcdf.endDef(ncid);
18     else;%try to use old mex stuff
19     % inverse the order of list dimensions
20     dimlist=fliplr(dimlist);
21     switch length(dimlist)
22     case 1,
23 roquet 1.4 eval(sprintf('ncid{''%s''}=nc%s(''%s'');',varname,xtype,dimlist{1}));
24 roquet 1.1 case 2,
25 roquet 1.4 eval(sprintf('ncid{''%s''}=nc%s(''%s'',''%s'');',varname,xtype,dimlist{1},dimlist{2}));
26 roquet 1.1 case 3,
27 roquet 1.4 eval(sprintf('ncid{''%s''}=nc%s(''%s'',''%s'',''%s'');',varname,xtype,dimlist{1},dimlist{2},dimlist{3}));
28 roquet 1.1 otherwise
29 roquet 1.2 error('ncaddVar: number of dimension > 3');
30 roquet 1.1 end
31     end;
32    
33    

  ViewVC Help
Powered by ViewVC 1.1.22