profilesMatlabProcessing/profiles_misc/ncputAtt.m

function [nc] = ncputAtt(ncid,varname,attrname,attrvalue);
% add an attribute to a variable in a netcdf file.

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

if isempty(varname),
    if useNativeMatlabNetcdf;
        netcdf.putAtt(ncid,-1,attrname,attrvalue);
    else;%try to use old mex stuff
        if ischar(attrvalue)
            attrvalue(find(double(attrvalue)==10))=[];
            eval(['ncid.' attrname '=''' attrvalue ''';']);
        else
            eval(['ncid.' attrname '=' num2str(attrvalue) ';']);
        end
    end;
else
    if useNativeMatlabNetcdf;
        varid=netcdf.inqVarID(ncid,varname);
        netcdf.putAtt(ncid,varid,attrname,attrvalue);
    else;%try to use old mex stuff
        if strcmp(attrname,'_FillValue'),
            attrname='FillValue_';
        end
        if ischar(attrvalue)
            attrvalue(find(double(attrvalue)==10))=[];
            eval(['ncid{''' varname '''}.' attrname '=''' attrvalue ''';']);
        else
            eval(['ncid{''' varname '''}.' attrname '=' num2str(attrvalue) ';']);
        end
    end;
end
    


1	function [nc] = ncputAtt(ncid,varname,attrname,attrvalue);
2	% add an attribute to a variable in a netcdf file.
3
4	global useNativeMatlabNetcdf;
5	if isempty(useNativeMatlabNetcdf); useNativeMatlabNetcdf = ~isempty(which('netcdf.open')); end;
6
7	if isempty(varname),
8	if useNativeMatlabNetcdf;
9	netcdf.putAtt(ncid,-1,attrname,attrvalue);
10	else;%try to use old mex stuff
11	if ischar(attrvalue)
12	attrvalue(find(double(attrvalue)==10))=[];
13	eval(['ncid.' attrname '=''' attrvalue ''';']);
14	else
15	eval(['ncid.' attrname '=' num2str(attrvalue) ';']);
16	end
17	end;
18	else
19	if useNativeMatlabNetcdf;
20	varid=netcdf.inqVarID(ncid,varname);
21	netcdf.putAtt(ncid,varid,attrname,attrvalue);
22	else;%try to use old mex stuff
23	if strcmp(attrname,'_FillValue'),
24	attrname='FillValue_';
25	end
26	if ischar(attrvalue)
27	attrvalue(find(double(attrvalue)==10))=[];
28	eval(['ncid{''' varname '''}.' attrname '=''' attrvalue ''';']);
29	else
30	eval(['ncid{''' varname '''}.' attrname '=' num2str(attrvalue) ';']);
31	end
32	end;
33	end
34
35
36