Mercurial > forge

comparison main/netcdf/inst/private/test_netcdf_high_level_interface.m @ 12561:8476336c120c octave-forge

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
ncread: do not replace _FillValue by NaNs for characters
author abarth93
date Fri, 20 Feb 2015 16:30:08 +0000
parents 5b3d9192547f
children a4de85b2e6a0
comparison
equal deleted inserted replaced
12560:9f933cf29066 12561:8476336c120c
88 assert(~info.Dimensions(1).Unlimited) 88 assert(~info.Dimensions(1).Unlimited)
89 assert(info.Dimensions(2).Unlimited) 89 assert(info.Dimensions(2).Unlimited)
90 90
91 delete(fname) 91 delete(fname)
92 92
93
94 % test double with _FillValue
95
96 fname = [tempname '-octave-netcdf.nc'];
97 fv = 99999.;
98 nccreate(fname,'flag','Dimensions',{'lon',10,'lat',10},'Datatype','double',...
99 'FillValue',fv);
100
101 %system(['ncdump -h ' fname])
102 data = zeros(10,10);
103 data(1,2) = fv;
104 ncid = netcdf_open(fname,'NC_WRITE');
105 varid = netcdf_inqVarID(ncid, 'flag');
106 netcdf_putVar(ncid,varid,data);
107 data2 = ncread(fname,'flag');
108 data(data == fv) = NaN;
109 assert(isequaln(data,data2))
110 delete(fname)
111
112
113 % test char with _FillValue
114
115 fname = [tempname '-octave-netcdf.nc'];
116 fv = '*';
117 nccreate(fname,'flag','Dimensions',{'lon',10,'lat',10},'Datatype','char',...
118 'FillValue',fv);
119 data = repmat('.',[10 10]);
120 data(1,2) = fv;
121
122 ncid = netcdf_open(fname,'NC_WRITE');
123 varid = netcdf_inqVarID(ncid, 'flag');
124 netcdf_putVar(ncid,varid,data);
125 data2 = ncread(fname,'flag');
126
127 assert(isequal(data,data2))
128 delete(fname)
129