Mercurial > forge

annotate main/netcdf/inst/private/test_netcdf_high_level_interface.m @ 11988:2f62d37d5abc octave-forge

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
restructure test (3)
author abarth93
date Mon, 29 Jul 2013 19:39:36 +0000
parents
children 6f80e9b5a5cb
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
11988
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
1
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
2 fname = [tempname '-octave-netcdf.nc'];
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
3
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
4
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
5 nccreate(fname,'temp','Dimensions',{'lon',10,'lat',20});
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
6 nccreate(fname,'salt','Dimensions',{'lon',10,'lat',20});
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
7 nccreate(fname,'u','Dimensions',{'lon','lat'});
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
8 u = randn(10,20);
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
9 ncwrite(fname,'u',u);
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
10
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
11 u2 = ncread(fname,'u');
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
12 assert(isequalwithequalnans(u,u2));
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
13
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
14 u2 = ncread(fname,'u',[10 5],[inf inf],[1 1]);
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
15 assert(isequalwithequalnans(u(10:end,5:end),u2));
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
16
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
17 ncwriteatt(fname,'temp','units','degree Celsius');
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
18 assert(strcmp(ncreadatt(fname,'temp','units'),'degree Celsius'));
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
19
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
20 ncwriteatt(fname,'temp','range',[0 10]);
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
21 assert(isequal(ncreadatt(fname,'temp','range'),[0 10]));
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
22
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
23 ncwriteatt(fname,'temp','float_range',single([0 10]));
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
24 assert(isequal(ncreadatt(fname,'temp','float_range'),[0 10]));
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
25
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
26 ncwriteatt(fname,'temp','int_range',int32([0 10]));
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
27 assert(isequal(ncreadatt(fname,'temp','int_range'),[0 10]));
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
28
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
29 info = ncinfo(fname);
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
30 assert(length(info.Variables) == 3)
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
31 assert(strcmp(info.Variables(1).Name,'temp'));
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
32 assert(isequal(info.Variables(1).Size,[10 20]));
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
33 delete(fname);
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
34
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
35
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
36 nccreate(fname,'temp','Dimensions',{'lon',10,'lat',20},'Format','64bit');
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
37 delete(fname);
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
38
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
39 nccreate(fname,'temp','Dimensions',{'lon',10,'lat',20},'Format','classic');
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
40 info = ncinfo(fname);
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
41 assert(strcmp(info.Format,'classic'));
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
42
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
43 delete(fname);
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
44
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
45 % netcdf4
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
46
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
47 nccreate(fname,'temp','Dimensions',{'lon',10,'lat',20},'Format','netcdf4');
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
48 ncwriteatt(fname,'temp','uint_range',uint32([0 10]));
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
49 assert(isequal(ncreadatt(fname,'temp','uint_range'),[0 10]));
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
50
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
51 info = ncinfo(fname);
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
52 assert(strcmp(info.Format,'netcdf4'));
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
53 delete(fname)
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
54
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
55 % scalar variable
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
56 nccreate(fname,'temp','Format','netcdf4','Datatype','double');
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
57 ncwrite(fname,'temp',123);
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
58 assert(ncread(fname,'temp') == 123)
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
59 delete(fname)
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
60
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
61
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
62 % test unlimited dimension with nccreate
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
63 fname = [tempname '-octave-netcdf.nc'];
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
64 nccreate(fname,'temp','Dimensions',{'lon',10,'lat',inf});
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
65 %system(['ncdump -h ' fname])
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
66
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
67 info = ncinfo(fname);
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
68 assert(~info.Dimensions(1).Unlimited)
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
69 assert(info.Dimensions(2).Unlimited)
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
70
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
71 delete(fname)
2f62d37d5abc restructure test (3)
abarth93
parents:
diff changeset
72