Resize variable dimensions NetCdf

Source MatLab: http://ecco2.jpl.nasa.gov/data1/matlab/netcdf_toolbox/netcdf/@ncvar/resize.m

function theResult = resize(self, newSize)

% ncvar/resize – Resize variable dimensions.
% resize(self, newSize) resizes the dimensions of self,
% an "ncvar" object. All variables related to the
% changed dimensions are similarly affected. The
% newSize must have the same number of dimensions
% as the existing variable. The new self is returned.

% Copyright (C) 1998 Dr. Charles R. Denham, ZYDECO.
% All Rights Reserved.
% Disclosure without explicit written consent from the
% copyright owner does not constitute publication.

% Version of 03-Nov-1998 08:52:22.
% Updated 12-Aug-1999 09:42:01.

if nargin < 1, help(mfilename), return, end
if nargout > 0, theResult = self; end

% Check for no-change.

if isequal(ncsize(self), newSize)
result = self;
if nargout > 0
theResult = result;
else
ncans(result)
end
return
end

theItemName = name(self);

% Check for writeability.

f = parent(self);
thePermission = permission(f);
theSrcName = name(f);

if isequal(thePermission, 'nowrite')
disp([' ## NetCDF source file must be writeable.'])
return
end

% Check request.

d = dim(self);
for i = 1:length(d)
if ~isrecdim(d

Unknown macro: {i}

) & newSize <= 0
disp([' ## Dimension "' name(d

) '" size requires positive integer.'])
return
end
end

% Create temporary file.

g = [];

i = 0;
while isempty(g)
i = i + 1;
theTmpName = ['tmp_' int2str(i) '.nc'];
if exist(theTmpName, 'file') ~= 2
g = netcdf(theTmpName, 'noclobber');
end
end

theTmpName = name(g);

% Copy affected dimensions first.

d = dim(self);
for i = 1:length(d)
if isrecdim(d

Unknown macro: {i}

)
g(name(d

)) = 0;
elseif newSize <= 0
error([' ## Dimension "' name(d

Unknown macro: {i}

) '" requires positive integer.'])
else
g(name(d

)) = newSize;
end
end

% Copy other dimensions.

d = dim(f);
for i = 1:length(d)
g(name(d

Unknown macro: {i}

)) = ncsize(d

);
end

% Copy global attributes.

a = att(f);
for i = 1:length(a)
copy(a

Unknown macro: {i}

, g)
end

% Copy variable definitions and attributes.

v = var(f);
for i = 1:length(v)
copy(v

, g, 0, 1)
end

% Copy variable data as minimal rectangular array.
% Note that the "()" operator is out-of-context
% inside this method, so we have to build our own
% calls to "ncvar/subsref" and "ncvar/subsasgn".
% It might be easier for us to use "virtual"
% variables instead, which could be transferred
% with the more intelligent "ncvar/copy" method.

v = var(f);
w = var(g);

for i = 1:length(v)
sv = ncsize(v

Unknown macro: {i}

);
sw = ncsize(w

);
if ~isempty(sw)
d = dim(w

Unknown macro: {i}

);
if isrecdim(d

Unknown macro: {1}

)
if sw(1) == 0
if isequal(name(d

), theItemName)
sw(1) = newSize;
else
sw(1) = sv(1);
end
end
end
end
theMinimalSize = min(sv, sw);
if prod(theMinimalSize) > 0
if isequal(sv, sw)
copy(v

, g, 1)
else
theIndices = cell(size(theMinimalSize));
for j = 1:length(theIndices)
theIndices

Unknown macro: {j}

= 1:theMinimalSize(j);
end
theStruct.type = '()';
theStruct.subs = theIndices;
theData = subsref(v

Unknown macro: {i}

, theStruct);
w

= subsasgn(w

Unknown macro: {i}

, theStruct, theData);
end
end
end

% Close both files.

f = close(f);
g = close(g);

% Delete old file.

delete(theSrcName)

% Rename new file to old file name.

fcopy(theTmpName, theSrcName)
delete(theTmpName)

% Open the new file.

g = netcdf(theSrcName, thePermission);

% Return the resized variable.

result = g

Unknown macro: {theItemName}

;

if nargout > 0
theResult = result;
else
ncans(result)
end

Blog

Resize variable dimensions NetCdf