1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
|
%% Copyright (C) 2013-2022 Alexander Barth
%%
%% This program is free software; you can redistribute it and/or modify
%% it under the terms of the GNU General Public License as published by
%% the Free Software Foundation; either version 3 of the License, or
%% (at your option) any later version.
%%
%% This program is distributed in the hope that it will be useful,
%% but WITHOUT ANY WARRANTY; without even the implied warranty of
%% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
%% GNU General Public License for more details.
%%
%% You should have received a copy of the GNU General Public License
%% along with this program; If not, see <http://www.gnu.org/licenses/>.
function test_netcdf_type(nctype,octtype)
import_netcdf
m = 5;
n = 10;
fname = [tempname '-octave-netcdf- ' nctype '.nc'];
mode = bitor(netcdf.getConstant('NC_CLOBBER'),...
netcdf.getConstant('NC_NETCDF4'));
ncid = netcdf.create(fname,mode);
dimids = [netcdf.defDim(ncid,'lon',m) ...
netcdf.defDim(ncid,'time',n)];
varid = netcdf.defVar(ncid,'variable',nctype,dimids);
netcdf.endDef(ncid)
if strcmp(octtype,'char')
z = char(floor(26*rand(m,n)) + 65);
testvals = 'a';
testvalv = 'this is a name';
else
z = zeros(m,n,octtype);
z(:) = randn(m,n);
testvals = zeros(1,1,octtype);
testvals(:) = rand(1,1);
testvalv = zeros(1,5,octtype);
testvalv(:) = rand(size(testvalv));
endif
netcdf.putVar(ncid,varid,z);
z2 = netcdf.getVar(ncid,varid);
assert(isequal(z,z2))
netcdf.putAtt(ncid,varid,'scalar_attribute',testvals);
val = netcdf.getAtt(ncid,varid,'scalar_attribute');
assert(isequal(val,testvals));
assert(strcmp(class(val),octtype))
[xtype,len] = netcdf.inqAtt(ncid,varid,'scalar_attribute');
assert(xtype == netcdf.getConstant(nctype))
assert(len == numel(testvals));
netcdf.putAtt(ncid,varid,'vector_attribute',testvalv);
val = netcdf.getAtt(ncid,varid,'vector_attribute');
assert(isequal(val,testvalv));
assert(strcmp(class(val),octtype))
[xtype,len] = netcdf.inqAtt(ncid,varid,'vector_attribute');
assert(xtype == netcdf.getConstant(nctype))
assert(len == numel(testvalv));
netcdf.close(ncid);
delete(fname);
endfunction
|
