/[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.2 - (hide annotations) (download)
Thu Apr 14 23:35:42 2011 UTC (14 years, 3 months ago) by roquet
Branch: MAIN
Changes since 1.1: +2 -2 lines
bug fix for old netcdf toolbox + new function to get a variable.

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

  ViewVC Help
Powered by ViewVC 1.1.22