1 #include <stdio.h>
2 #include <stdlib.h>
3 #include <netcdf.h>
4
5
6 static size_t var_chunksizes[4] = {4, 4, 4, 4} ;
7 static unsigned int var_0_filterparams[1] = {9U} ;
8
9 void
check_err(const int stat,const int line,const char * file)10 check_err(const int stat, const int line, const char *file) {
11 if (stat != NC_NOERR) {
12 (void)fprintf(stderr,"line %d of %s: %s\n", line, file, nc_strerror(stat));
13 fflush(stderr);
14 exit(1);
15 }
16 }
17
18 int
main()19 main() {/* create bzip2.nc */
20
21 int stat; /* return status */
22 int ncid; /* netCDF id */
23
24 /* group ids */
25 int bzip2_grp;
26
27 /* dimension ids */
28 int dim0_dim;
29 int dim1_dim;
30 int dim2_dim;
31 int dim3_dim;
32
33 /* dimension lengths */
34 size_t dim0_len = 4;
35 size_t dim1_len = 4;
36 size_t dim2_len = 4;
37 size_t dim3_len = 4;
38
39 /* variable ids */
40 int var_id;
41
42 /* rank (number of dimensions) for each variable */
43 # define RANK_var 4
44
45 /* variable shapes */
46 int var_dims[RANK_var];
47
48 /* enter define mode */
49 stat = nc_create("bzip2.nc", NC_CLOBBER|NC_NETCDF4, &ncid);
50 check_err(stat,__LINE__,__FILE__);
51 bzip2_grp = ncid;
52
53 /* define dimensions */
54 stat = nc_def_dim(bzip2_grp, "dim0", dim0_len, &dim0_dim);
55 check_err(stat,__LINE__,__FILE__);
56 stat = nc_def_dim(bzip2_grp, "dim1", dim1_len, &dim1_dim);
57 check_err(stat,__LINE__,__FILE__);
58 stat = nc_def_dim(bzip2_grp, "dim2", dim2_len, &dim2_dim);
59 check_err(stat,__LINE__,__FILE__);
60 stat = nc_def_dim(bzip2_grp, "dim3", dim3_len, &dim3_dim);
61 check_err(stat,__LINE__,__FILE__);
62
63 /* define variables */
64
65 var_dims[0] = dim0_dim;
66 var_dims[1] = dim1_dim;
67 var_dims[2] = dim2_dim;
68 var_dims[3] = dim3_dim;
69 stat = nc_def_var(bzip2_grp, "var", NC_FLOAT, RANK_var, var_dims, &var_id);
70 check_err(stat,__LINE__,__FILE__);
71 stat = nc_def_var_chunking(bzip2_grp, var_id, NC_CHUNKED, var_chunksizes);
72 check_err(stat,__LINE__,__FILE__);
73 stat = nc_def_var_fill(bzip2_grp, var_id, NC_NOFILL, NULL);
74 check_err(stat,__LINE__,__FILE__);
75 stat = nc_def_var_filter(bzip2_grp, var_id, 307, 1, var_0_filterparams);
76 check_err(stat,__LINE__,__FILE__);
77
78 /* leave define mode */
79 stat = nc_enddef (bzip2_grp);
80 check_err(stat,__LINE__,__FILE__);
81
82 /* assign variable data */
83
84 {
85 float var_data[256] = {((float)0), ((float)1), ((float)2), ((float)3), ((float)4), ((float)5), ((float)6), ((float)7), ((float)8), ((float)9), ((float)10), ((float)11), ((float)12), ((float)13), ((float)14), ((float)15), ((float)16), ((float)17), ((float)18), ((float)19), ((float)20), ((float)21), ((float)22), ((float)23), ((float)24), ((float)25), ((float)26), ((float)27), ((float)28), ((float)29), ((float)30), ((float)31), ((float)32), ((float)33), ((float)34), ((float)35), ((float)36), ((float)37), ((float)38), ((float)39), ((float)40), ((float)41), ((float)42), ((float)43), ((float)44), ((float)45), ((float)46), ((float)47), ((float)48), ((float)49), ((float)50), ((float)51), ((float)52), ((float)53), ((float)54), ((float)55), ((float)56), ((float)57), ((float)58), ((float)59), ((float)60), ((float)61), ((float)62), ((float)63), ((float)64), ((float)65), ((float)66), ((float)67), ((float)68), ((float)69), ((float)70), ((float)71), ((float)72), ((float)73), ((float)74), ((float)75), ((float)76), ((float)77), ((float)78), ((float)79), ((float)80), ((float)81), ((float)82), ((float)83), ((float)84), ((float)85), ((float)86), ((float)87), ((float)88), ((float)89), ((float)90), ((float)91), ((float)92), ((float)93), ((float)94), ((float)95), ((float)96), ((float)97), ((float)98), ((float)99), ((float)100), ((float)101), ((float)102), ((float)103), ((float)104), ((float)105), ((float)106), ((float)107), ((float)108), ((float)109), ((float)110), ((float)111), ((float)112), ((float)113), ((float)114), ((float)115), ((float)116), ((float)117), ((float)118), ((float)119), ((float)120), ((float)121), ((float)122), ((float)123), ((float)124), ((float)125), ((float)126), ((float)127), ((float)128), ((float)129), ((float)130), ((float)131), ((float)132), ((float)133), ((float)134), ((float)135), ((float)136), ((float)137), ((float)138), ((float)139), ((float)140), ((float)141), ((float)142), ((float)143), ((float)144), ((float)145), ((float)146), ((float)147), ((float)148), ((float)149), ((float)150), ((float)151), ((float)152), ((float)153), ((float)154), ((float)155), ((float)156), ((float)157), ((float)158), ((float)159), ((float)160), ((float)161), ((float)162), ((float)163), ((float)164), ((float)165), ((float)166), ((float)167), ((float)168), ((float)169), ((float)170), ((float)171), ((float)172), ((float)173), ((float)174), ((float)175), ((float)176), ((float)177), ((float)178), ((float)179), ((float)180), ((float)181), ((float)182), ((float)183), ((float)184), ((float)185), ((float)186), ((float)187), ((float)188), ((float)189), ((float)190), ((float)191), ((float)192), ((float)193), ((float)194), ((float)195), ((float)196), ((float)197), ((float)198), ((float)199), ((float)200), ((float)201), ((float)202), ((float)203), ((float)204), ((float)205), ((float)206), ((float)207), ((float)208), ((float)209), ((float)210), ((float)211), ((float)212), ((float)213), ((float)214), ((float)215), ((float)216), ((float)217), ((float)218), ((float)219), ((float)220), ((float)221), ((float)222), ((float)223), ((float)224), ((float)225), ((float)226), ((float)227), ((float)228), ((float)229), ((float)230), ((float)231), ((float)232), ((float)233), ((float)234), ((float)235), ((float)236), ((float)237), ((float)238), ((float)239), ((float)240), ((float)241), ((float)242), ((float)243), ((float)244), ((float)245), ((float)246), ((float)247), ((float)248), ((float)249), ((float)250), ((float)251), ((float)252), ((float)253), ((float)254), ((float)255)} ;
86 size_t var_startset[4] = {0, 0, 0, 0} ;
87 size_t var_countset[4] = {4, 4, 4, 4};
88 stat = nc_put_vara(bzip2_grp, var_id, var_startset, var_countset, var_data);
89 check_err(stat,__LINE__,__FILE__);
90 }
91
92
93 stat = nc_close(bzip2_grp);
94 check_err(stat,__LINE__,__FILE__);
95 return 0;
96 }
97