matlab_class/gcmfaces_IO/ncgetvar.m

function [VARvalue]=ncgetvar(nc,VARname,varargin)
% function [VARvalue]=ncgetvar(nc,VARname, [index_prof, index_depth] )
%   get data to MITprof netcdf file

global useNativeMatlabNetcdf; 
if isempty(useNativeMatlabNetcdf); useNativeMatlabNetcdf = ~isempty(which('netcdf.open')); end;

if useNativeMatlabNetcdf
    
    %get variable id:
    vv = netcdf.inqVarID(nc,VARname);
    if nargin>2;
        %get and flip position vectors:
        VARpos=fliplr(varargin);
        %convert VARpos to start,count:
        start=[]; count=[];
        for ii=1:length(VARpos);
            start=[start VARpos{ii}(1)-1];
            count=[count VARpos{ii}(end)-VARpos{ii}(1)+1];
        end;
        %write to file:
        VARvalue=netcdf.getVar(nc,vv,start,count);
    else;
        %write to file:
        VARvalue=netcdf.getVar(nc,vv);
    end;
    %flip order of dimensions:
    bb=length(size(VARvalue)); VARvalue=permute(VARvalue,[bb:-1:1]);
    
    
else;%try to use old mex stuff
    
    if nargin==3;
        eval(['VARvalue=nc{''' VARname '''}([' num2str(varargin{1}) ']);']);
    elseif nargin==4,
        eval(['VARvalue=nc{''' VARname '''}([' num2str(varargin{1}) '],[' num2str(varargin{2}) ']);']);
    else;
        eval(['VARvalue=nc{''' VARname '''}(:);']);
    end;
    
end;
1	gforget	1.1	function [VARvalue]=ncgetvar(nc,VARname,varargin)
2			% function [VARvalue]=ncgetvar(nc,VARname, [index_prof, index_depth] )
3			% get data to MITprof netcdf file
4
5			global useNativeMatlabNetcdf;
6			if isempty(useNativeMatlabNetcdf); useNativeMatlabNetcdf = ~isempty(which('netcdf.open')); end;
7
8			if useNativeMatlabNetcdf
9
10			%get variable id:
11			vv = netcdf.inqVarID(nc,VARname);
12			if nargin>2;
13			%get and flip position vectors:
14			VARpos=fliplr(varargin);
15			%convert VARpos to start,count:
16			start=[]; count=[];
17			for ii=1:length(VARpos);
18			start=[start VARpos{ii}(1)-1];
19			count=[count VARpos{ii}(end)-VARpos{ii}(1)+1];
20			end;
21			%write to file:
22			VARvalue=netcdf.getVar(nc,vv,start,count);
23			else;
24			%write to file:
25			VARvalue=netcdf.getVar(nc,vv);
26			end;
27			%flip order of dimensions:
28			bb=length(size(VARvalue)); VARvalue=permute(VARvalue,[bb:-1:1]);
29
30
31			else;%try to use old mex stuff
32
33			if nargin==3;
34			eval(['VARvalue=nc{''' VARname '''}([' num2str(varargin{1}) ']);']);
35			elseif nargin==4,
36			eval(['VARvalue=nc{''' VARname '''}([' num2str(varargin{1}) '],[' num2str(varargin{2}) ']);']);
37			else;
38			eval(['VARvalue=nc{''' VARname '''}(:);']);
39			end;
40
41			end;