ncdf4a13/cxx/nctst.cpp File Reference

#include "config.h"
#include <iostream>
#include <string.h>
#include "netcdfcpp.h"

Include dependency graph for nctst.cpp:

Go to the source code of this file.

Namespaces
namespace	std
Classes
class	DumpableNcFile
Functions
void	gen (const char *path, NcFile::FileFormat format)
static char *	cdl_name (const char *path)
void	dumpatts (NcVar &var)
void	dump (const char *path)
int	main (void)

---

Function Documentation

static char* cdl_name (  const char *  path  )  [static]

Definition at line 105 of file nctst.cpp. References MAX_NC_NAME. Referenced by dump(). { const char* cp = path + strlen(path); while (*(cp-1) != '/' && cp != path) // assumes UNIX path separator cp--; static char np[MAX_NC_NAME]; strncpy(&np[0], cp, MAX_NC_NAME); char* ep = np + strlen(np); while (*ep != '.' && ep != np) ep--; if (*ep == '.') *ep = '\0'; return np; }

void dump (  const char *  path  )

Definition at line 211 of file nctst.cpp. References cdl_name(), DumpableNcFile::dumpdata(), DumpableNcFile::dumpdims(), DumpableNcFile::dumpgatts(), DumpableNcFile::dumpvars(), and NcFile::num_atts(). { DumpableNcFile nc(path); // default is open in read-only mode cout << "netcdf " << cdl_name(path) << " {" << endl << "dimensions:" << endl ; nc.dumpdims(); cout << "variables:" << endl; nc.dumpvars(); if (nc.num_atts() > 0) cout << "// global attributes" << endl ; nc.dumpgatts(); cout << "data:" << endl; nc.dumpdata(); cout << "}" << endl; }

void dumpatts (  NcVar &  var  )

Definition at line 149 of file nctst.cpp. References NcVar::get_att(), NcAtt::name(), NcVar::name(), and NcAtt::values(). Referenced by DumpableNcFile::dumpvars(). { NcToken vname = var.name(); NcAtt* ap; for(int n = 0; ap = var.get_att(n); n++) { cout << "\t\t" << vname << ":" << ap->name() << " = " ; NcValues* vals = ap->values(); cout << *vals << " ;" << endl ; delete ap; delete vals; } }

void gen (  const char *  path, NcFile::FileFormat  format )

Definition at line 17 of file nctst.cpp. References NcFile::add_att(), NcVar::add_att(), NcFile::add_dim(), NcFile::add_var(), NcFile::is_valid(), ncChar, ncFloat, ncInt, ncLong, NcVar::put(), NcFile::Replace, and NcVar::set_cur(). Referenced by main(). { NcFile nc(path, NcFile::Replace, NULL, 0, format); // Create, leave in define mode // Check if the file was opened successfully if (! nc.is_valid()) { cerr << "can't create netCDF file " << path << "\n"; return; } // Create dimensions const int NLATS = 4; const int NLONS = 3; const int NFRTIMES = 2; const int TIMESTRINGLEN = 20; NcDim* latd = nc.add_dim("lat", NLATS); NcDim* lond = nc.add_dim("lon", NLONS); NcDim* frtimed = nc.add_dim("frtime"); // unlimited dimension NcDim* timelend = nc.add_dim("timelen", TIMESTRINGLEN); // Create variables and their attributes NcVar* P = nc.add_var("P", ncFloat, frtimed, latd, lond); P->add_att("long_name", "pressure at maximum wind"); P->add_att("units", "hectopascals"); static float range[] = {0., 1500.}; P->add_att("valid_range", 2, range); P->add_att("_FillValue", -9999.0f); NcVar* lat = nc.add_var("lat", ncFloat, latd); lat->add_att("long_name", "latitude"); lat->add_att("units", "degrees_north"); NcVar* lon = nc.add_var("lon", ncFloat, lond); lon->add_att("long_name", "longitude"); lon->add_att("units", "degrees_east"); NcVar* frtime = nc.add_var("frtime", ncLong, frtimed); frtime->add_att("long_name", "forecast time"); frtime->add_att("units", "hours"); NcVar* reftime = nc.add_var("reftime",ncChar,timelend); reftime->add_att("long_name", "reference time"); reftime->add_att("units", "text_time"); NcVar* scalar = nc.add_var("scalarv", ncInt); scalar->add_att("scalar_att", 1); // Global attributes nc.add_att("history", "created by Unidata LDM from NPS broadcast"); nc.add_att("title", "NMC Global Product Set: Pressure at Maximum Wind"); // Start writing data, implictly leaves define mode static float lats[NLATS] = {-90, -87.5, -85, -82.5}; lat->put(lats, NLATS); static float lons[NLONS] = {-180, -175, -170}; lon->put(lons, NLONS); static int frtimes[NFRTIMES] = {12, 18}; frtime->put(frtimes, NFRTIMES); static const char* s = "1992-3-21 12:00" ; reftime->put(s, strlen(s)); static float P_data[2][4][3] = { {{950, 951, 952}, {953, 954, 955}, {956, 957, 958}, {959, 960, 961}}, {{962, 963, 964}, {965, 966, 967}, {968, 969, 970}, {971, 972, 973}} }; // We could write all P data at once with P->put(&P_data[0][0][0], P->edges()), // but instead we write one record at a time, to show use of setcur(). long rec = 0; // start at zero-th const long nrecs = 1; // # records to write P->put(&P_data[0][0][0], nrecs, NLATS, NLONS); // write zero-th record P->set_cur(++rec); // set to next record P->put(&P_data[1][0][0], nrecs, NLATS, NLONS); // write next record // close of nc takes place in destructor }

int main (  void   )

Definition at line 238 of file nctst.cpp. References NcFile::Classic, gen(), and NcFile::Offset64Bits. { cout << "Creating netCDF classic file with C++ API..."; gen("nctst_classic.nc", NcFile::Classic); // create a netCDF file cout << "ok\n"; cout << "Creating netCDF 64-bit offset file with C++ API..."; gen("nctst_64bit_offset.nc", NcFile::Offset64Bits); // create a netCDF file cout << "ok\n"; #ifdef USE_NETCDF4 // Uncomment the following line to get lots of feedback from the library. // nc_set_log_level(5); // cout << "Creating netCDF4 file with C++ API..."; // gen("nctst_netcdf4.nc", NcFile::Netcdf4); // create a netCDF file // cout << "ok\n"; // nc_set_log_level(0); // cout << "Creating netCDF4 classic file with C++ API..."; // gen("nctst_netcdf4_classic.nc", NcFile::Netcdf4Classic); // create a netCDF file // cout << "ok\n"; #endif return 0; }

---