1/** \file
2
3This program benchmarks creating a netCDF file and reading records.
4
5Copyright 2011, UCAR/Unidata See COPYRIGHT file for copying and
6redistribution conditions.
7*/
8
9#include <config.h>
10#include <nc_tests.h>
11#include "err_macros.h"
12#include <netcdf.h>
13#include <stdio.h>
14#include <stdlib.h>
15#include <time.h>
16#include <sys/time.h> /* Extra high precision time info. */
17
18/* We will create this file. */
19#define FILE_NAME "bm_netcdf4_recs.nc"
20
21int main(int argc, char **argv)
22{
23    struct timeval start_timeend_timediff_time;
24    double sec;
25    int nitem = 10000; /* default number of objects of each type */
26    int i;
27    int ncid;
28    int data[] = {42};
29    int ggrpnumgrp;
30    char gname[16];
31    int vvarnumvarvnvleftnvars;
32
33    int  stat;  /* return status */
34
35    /* dimension ids */
36    int basetime_dim;
37    int forecast_dim;
38    int bounds_dim;
39    int latitude_dim;
40    int longitude_dim;
41
42    /* dimension lengths */
43    size_t basetime_len = NC_UNLIMITED;
44    size_t forecast_len = 32;
45    size_t bounds_len = 2;
46    size_t latitude_len = 121;
47    size_t longitude_len = 101;
48
49    /* variable ids */
50    int temperature_2m_id;
51
52    /* rank (number of dimensions) for each variable */
53#   define RANK_temperature_2m 4
54
55    /* variable shapes */
56    int temperature_2m_dims[RANK_temperature_2m];
57    static const float temperature_2m_FillValue_att[1] = {9.96921e+36} ;
58    static const float temperature_2m_missing_value_att[1] = {9.96921e+36} ;
59    static const float temperature_2m_valid_min_att[1] = {180} ;
60    static const float temperature_2m_valid_max_att[1] = {330} ;
61
62    /* enter define mode */
63    if (nc_create(FILE_NAMENC_CLOBBER, &ncid)) ERR;
64
65    /* define dimensions */
66    if (nc_def_dim(ncid, "basetime", basetime_len, &basetime_dim)) ERR;
67    if (nc_def_dim(ncid, "forecast", forecast_len, &forecast_dim)) ERR;
68    if (nc_def_dim(ncid, "bounds", bounds_len, &bounds_dim)) ERR;
69    if (nc_def_dim(ncid, "latitude", latitude_len, &latitude_dim)) ERR;
70    if (nc_def_dim(ncid, "longitude", longitude_len, &longitude_dim)) ERR;
71
72    /* define variables */
73    temperature_2m_dims[0] = basetime_dim;
74    temperature_2m_dims[1] = forecast_dim;
75    temperature_2m_dims[2] = latitude_dim;
76    temperature_2m_dims[3] = longitude_dim;
77    if (nc_def_var(ncid, "temperature_2m", NC_FLOATRANK_temperature_2m,
78    temperature_2m_dims, &temperature_2m_id)) ERR;
79
80    /* assign per-variable attributes */
81    if (nc_put_att_text(ncidtemperature_2m_id, "long_name", 36, "Air temperature 2m above the surface")) ERR;
82    if (nc_put_att_text(ncidtemperature_2m_id, "units", 1, "K")) ERR;
83    if (nc_put_att_float(ncidtemperature_2m_id, "_FillValue", NC_FLOAT, 1, temperature_2m_FillValue_att)) ERR;
84    if (nc_put_att_float(ncidtemperature_2m_id, "missing_value", NC_FLOAT, 1, temperature_2m_missing_value_att)) ERR;
85    if (nc_put_att_float(ncidtemperature_2m_id, "valid_min", NC_FLOAT, 1, temperature_2m_valid_min_att)) ERR;
86    if (nc_put_att_float(ncidtemperature_2m_id, "valid_max", NC_FLOAT, 1, temperature_2m_valid_max_att)) ERR;
87    if (nc_put_att_text(ncidtemperature_2m_id, "standard_name", 15, "air_temperature")) ERR;
88    if (nc_put_att_text(ncidtemperature_2m_id, "cell_methods", 10, "area: mean")) ERR;
89    if (nc_put_att_text(ncidtemperature_2m_id, "coordinates", 5, "level")) ERR;
90    if (nc_close(ncid)) ERR;
91
92    if (gettimeofday(&start_timeNULL)) ERR;
93
94    return(0);
95}


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