MITgcm_contrib/mitgcm_tools/nc_add.m

function [] = nc_add(nc,shortName,longName,units,grid,var,varargin)
% nc_add(nc,shortName,longName,units,grid,var)
% nc_add(nc,shortName,longName,units,grid,var,record_number)
%
% e.g.
% >> nc=netcdf('test.nc','clobber');
%
% Add a dimension (use grid==shortName)
% >> nc_add(nc,'X','Longitude','Degrees East','X',xc)
% >> nc_add(nc,'Y','Latitude','Degrees North','Y',yc)
%
% Add a record dimension (use grid==shortName and no values)
% >> nc_add(nc,'T','Time','Years','T',[])
%
% Add a variable on a grid
% >> nc_add(nc,'H','Orography','m',{'Y' 'X'},H)
%
% Add a variable on a grid with recordable dimension
% >> nc_add(nc,'time','Time','Years','T',1978,1)
% >> nc_add(nc,'theta','Pot. Temp','Degrees K',{'T' 'Y' 'X'},theta,1)
% >> nc_add(nc,'time','Time','Years','T',1978,2)
% >> nc_add(nc,'theta','Pot. Temp','Degrees K',{'T' 'Y' 'X'},theta,2)
%
% >> close(nc);
%
% $Header: /u/gcmpack/MITgcm_contrib/mitgcm_tools/nc_add.m,v 1.2 2004/06/02 19:45:46 jmc Exp $

if strcmp(shortName,grid{1})
 nc(shortName) = length(var);
end

nc{shortName} = grid;
%nc{shortName}.uniquename = longName;
nc{shortName}.long_name = longName;
nc{shortName}.units = units;
nc{shortName}.missing_value = ncdouble(NaN);

nvargs=nargin-6;
rcn=0;
if nvargs==1
 rcn=varargin{1};
end

if ~isempty(var)
 if prod(size(var))==1
  nc{shortName}(rcn) = var;
 else
  if rcn ~= 0
   nc{shortName}(rcn,:) = permute(var,ndims(var):-1:1);
  else
   nc{shortName}(:) = permute(var,ndims(var):-1:1);
  end
 end
end

1	adcroft	1.1	function [] = nc_add(nc,shortName,longName,units,grid,var,varargin)
2			% nc_add(nc,shortName,longName,units,grid,var)
3			% nc_add(nc,shortName,longName,units,grid,var,record_number)
4			%
5			% e.g.
6			% >> nc=netcdf('test.nc','clobber');
7			%
8			% Add a dimension (use grid==shortName)
9			% >> nc_add(nc,'X','Longitude','Degrees East','X',xc)
10			% >> nc_add(nc,'Y','Latitude','Degrees North','Y',yc)
11			%
12			% Add a record dimension (use grid==shortName and no values)
13			% >> nc_add(nc,'T','Time','Years','T',[])
14			%
15			% Add a variable on a grid
16			% >> nc_add(nc,'H','Orography','m',{'Y' 'X'},H)
17			%
18			% Add a variable on a grid with recordable dimension
19			% >> nc_add(nc,'time','Time','Years','T',1978,1)
20			% >> nc_add(nc,'theta','Pot. Temp','Degrees K',{'T' 'Y' 'X'},theta,1)
21			% >> nc_add(nc,'time','Time','Years','T',1978,2)
22			% >> nc_add(nc,'theta','Pot. Temp','Degrees K',{'T' 'Y' 'X'},theta,2)
23			%
24			% >> close(nc);
25			%
26	adcroft	1.3	% $Header: /u/gcmpack/MITgcm_contrib/mitgcm_tools/nc_add.m,v 1.2 2004/06/02 19:45:46 jmc Exp $
27	adcroft	1.1
28	jmc	1.2	if strcmp(shortName,grid{1})
29	adcroft	1.1	nc(shortName) = length(var);
30			end
31
32			nc{shortName} = grid;
33			%nc{shortName}.uniquename = longName;
34			nc{shortName}.long_name = longName;
35			nc{shortName}.units = units;
36			nc{shortName}.missing_value = ncdouble(NaN);
37
38			nvargs=nargin-6;
39			rcn=0;
40			if nvargs==1
41			rcn=varargin{1};
42			end
43
44			if ~isempty(var)
45			if prod(size(var))==1
46			nc{shortName}(rcn) = var;
47			else
48			if rcn ~= 0
49			nc{shortName}(rcn,:) = permute(var,ndims(var):-1:1);
50			else
51			nc{shortName}(:) = permute(var,ndims(var):-1:1);
52			end
53			end
54			end
55