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

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

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


Revision 1.4 - (hide annotations) (download)
Thu Apr 14 23:35:42 2011 UTC (14 years, 3 months ago) by roquet
Branch: MAIN
CVS Tags: checkpoint65r, checkpoint65p, checkpoint65q
Changes since 1.3: +1 -1 lines
bug fix for old netcdf toolbox + new function to get a variable.

1 roquet 1.1 function varid = ncdefVar(ncid,varname,xtype,dimlist);
2     % add a variable in a netcdf file.
3    
4     global useNativeMatlabNetcdf;
5     if isempty(useNativeMatlabNetcdf); useNativeMatlabNetcdf = ~isempty(which('netcdf.open')); end;
6    
7     if useNativeMatlabNetcdf;
8 roquet 1.4 if isempty(dimlist), error('ncdefVar 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.defVar(ncid,varname,xtype,iDim);
14     else;%try to use old mex stuff
15     % inverse the order of list dimensions
16     dimlist=fliplr(dimlist);
17 roquet 1.2 switch length(dimlist)
18     case 1,
19 roquet 1.3 eval(sprintf('ncid{''%s''}=nc%s(''%s'');',varname,xtype,dimlist{1}));
20 roquet 1.2 case 2,
21 roquet 1.3 eval(sprintf('ncid{''%s''}=nc%s(''%s'',''%s'');',varname,xtype,dimlist{1},dimlist{2}));
22 roquet 1.2 case 3,
23 roquet 1.3 eval(sprintf('ncid{''%s''}=nc%s(''%s'',''%s'',''%s'');',varname,xtype,dimlist{1},dimlist{2},dimlist{3}));
24 roquet 1.2 otherwise
25     error('ncdefVar: number of dimension > 3');
26     end
27 roquet 1.1 end;
28    
29    

  ViewVC Help
Powered by ViewVC 1.1.22