profilesMatlabProcessing/profiles_misc/ncgetFillVal.m

function [FillVal]=ncgetFillVal(ncid,varname);
% [FillVal]=ncgetFillVal(ncid,varname)
%           return the missing_value or _FillValue of varname
%  return an error if varname does not exist

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

FillVal=[];
if useNativeMatlabNetcdf;
    varid = netcdf.inqVarID(ncid,varname);
    [varname,xtype,dimids,natts] = netcdf.inqVar(ncid,varid);
    [atts]=ncatts(ncid,varid);
    if any(ismember(atts,'missing_value'))
        FillVal = netcdf.getAtt(ncid,varid,'missing_value');
    elseif any(ismember(atts,'_FillValue'))
        FillVal = netcdf.getAtt(ncid,varid,'_FillValue');
    end
    if strcmp(xtype,'single') | strcmp(xtype,'double')
        FillVal=double(FillVal);
    end
else
    eval(['FillVal = ncid{''' varname '''}.missing_value(:);']);
    if isempty(FillVal);
        eval(['FillVal = ncid{''' varname '''}.FillValue_(:);']);
    end
end

1	roquet	1.1	function [FillVal]=ncgetFillVal(ncid,varname);
2			% [FillVal]=ncgetFillVal(ncid,varname)
3			% return the missing_value or _FillValue of varname
4	roquet	1.4	% return an error if varname does not exist
5	roquet	1.1
6			global useNativeMatlabNetcdf;
7			if isempty(useNativeMatlabNetcdf); useNativeMatlabNetcdf = ~isempty(which('netcdf.open')); end;
8
9			FillVal=[];
10			if useNativeMatlabNetcdf;
11			varid = netcdf.inqVarID(ncid,varname);
12	roquet	1.2	[varname,xtype,dimids,natts] = netcdf.inqVar(ncid,varid);
13	roquet	1.1	[atts]=ncatts(ncid,varid);
14			if any(ismember(atts,'missing_value'))
15	roquet	1.2	FillVal = netcdf.getAtt(ncid,varid,'missing_value');
16	roquet	1.1	elseif any(ismember(atts,'_FillValue'))
17	roquet	1.2	FillVal = netcdf.getAtt(ncid,varid,'_FillValue');
18			end
19			if strcmp(xtype,'single') \| strcmp(xtype,'double')
20			FillVal=double(FillVal);
21			end
22	roquet	1.1	else
23	roquet	1.3	eval(['FillVal = ncid{''' varname '''}.missing_value(:);']);
24	roquet	1.1	if isempty(FillVal);
25	roquet	1.3	eval(['FillVal = ncid{''' varname '''}.FillValue_(:);']);
26	roquet	1.1	end
27			end
28