1/** \file
2Write a file demonstrating some of the features of netCDF-4.
3
4We create two shared dimensions, "x" and "y", in a parent group, and
5some netCDF variables in different subgroups. The variables will
6include a compound and an enum type, as well as some of the new atomic
7types, like the unsigned 64-bit integer.
8
9This is part of the netCDF package. Copyright 2006-2011 University
10Corporation for Atmospheric Research/Unidata. See COPYRIGHT file for
11conditions of use. Full documentation of the netCDF can be found at
12http://www.unidata.ucar.edu/software/netcdf/docs.
13*/
14
15#include <stdlib.h>
16#include <stdio.h>
17#include <netcdf.h>
18
19/* This is the name of the data file we will create. */
20#define FILE_NAME "simple_nc4.nc"
21
22/* We are writing 2D data, a 6 x 12 grid. */
23#define NDIMS 2
24#define NX 6
25#define NY 12
26
27/* Handle errors by printing an error message and exiting with a
28 * non-zero status. */
29#define ERRCODE 2
30#define ERR(e) {printf("Error: %s\n", nc_strerror(e)); exit(ERRCODE);}
31
32int
33main()
34{
35   /* When we create netCDF variables, groups, dimensions, or types,
36    * we get back an ID for each one. */
37   int ncidx_dimidy_dimidvarid1varid2grp1idgrp2idtypeid;
38   int dimids[NDIMS];
39
40   /* This is the data array we will write. It will be filled with a
41    * progression of numbers for this example. */
42   unsigned long long data_out[NX][NY];
43
44   /* Loop indexes, and error handling. */
45   int xyretval;
46
47   /* The following struct is written as a compound type. */
48   struct s1
49   {
50         int i1;
51         int i2;
52   };
53   struct s1 compound_data[NX][NY];
54
55   /* Create some pretend data. */
56   for (x = 0; x < NXx++)
57      for (y = 0; y < NYy++)
58      {
59         data_out[x][y] = x * NY + y;
60         compound_data[x][y].i1 = 42;
61         compound_data[x][y].i2 = -42;
62      }
63
64   /* Create the file. The NC_NETCDF4 flag tells netCDF to
65    * create a netCDF-4/HDF5 file.*/
66   if ((retval = nc_create(FILE_NAMENC_NETCDF4|NC_CLOBBER, &ncid)))
67      ERR(retval);
68
69   /* Define the dimensions in the root group. Dimensions are visible
70    * in all subgroups. */
71   if ((retval = nc_def_dim(ncid, "x", NX, &x_dimid)))
72      ERR(retval);
73   if ((retval = nc_def_dim(ncid, "y", NY, &y_dimid)))
74      ERR(retval);
75
76   /* The dimids passes the IDs of the dimensions of the variable. */
77   dimids[0] = x_dimid;
78   dimids[1] = y_dimid;
79
80   /* Define two groups, "grp1" and "grp2." */
81   if ((retval = nc_def_grp(ncid, "grp1", &grp1id)))
82      ERR (retval);
83   if ((retval = nc_def_grp(ncid, "grp2", &grp2id)))
84      ERR (retval);
85
86   /* Define an unsigned 64bit integer variable in grp1, using dimensions
87    * in the root group. */
88   if ((retval = nc_def_var(grp1id, "data", NC_UINT64NDIMS,
89                            dimids, &varid1)))
90      ERR(retval);
91
92   /* Write unsigned long long data to the file. For netCDF-4 files,
93    * nc_enddef will be called automatically. */
94   if ((retval = nc_put_var_ulonglong(grp1idvarid1, &data_out[0][0])))
95      ERR(retval);
96
97   /* Create a compound type. This will cause nc_reddef to be called. */
98   if (nc_def_compound(grp2id, sizeof(struct s1), "sample_compound_type",
99                       &typeid))
100      ERR(retval);
101   if (nc_insert_compound(grp2idtypeid, "i1",
102                          offsetof(struct s1i1), NC_INT))
103      ERR(retval);
104   if (nc_insert_compound(grp2idtypeid, "i2",
105                          offsetof(struct s1i2), NC_INT))
106      ERR(retval);
107
108   /* Define a compound type variable in grp2, using dimensions
109    * in the root group. */
110   if ((retval = nc_def_var(grp2id, "data", typeidNDIMS,
111                            dimids, &varid2)))
112      ERR(retval);
113
114   /* Write the array of struct to the file. This will cause nc_endef
115    * to be called. */
116   if ((retval = nc_put_var(grp2idvarid2, &compound_data[0][0])))
117      ERR(retval);
118
119   /* Close the file. */
120   if ((retval = nc_close(ncid)))
121      ERR(retval);
122
123   printf("*** SUCCESS writing example file simple_nc4.nc!\n");
124   return 0;
125}


HyperKWIC - Version 7.20DA executed at 11:37 on 27 Oct 2017 | Polyhedron Solutions - INTERNAL USE | COMMERCIAL (Any O/S) SN 4AKIed