DennisHeimbigner
commented
1 year ago Try this:
ncdump -hs
| fgrep -i filter
On 2/22/2023 11:14 PM, tayloj25 wrote:
I am using netcdf-c-4.9.1 with hdf5-1.12.2 and zlib-1.2.13. I am attempting to run quite a similar example to the pres_temp_4D_rd.c example on the unidata website. My dimensions are slightly different though. My code appears as below:
include
include
include
/ This is the name of the data file we will read. /
define FILE_NAME "pressure_data.nc"
define NDIMS 4
define NLAT 45
define NLON 57
define LAT_NAME "y"
define LON_NAME "x"
define NREC 1
define REC_NAME "time"
define LVL_NAME "z"
define NLVL 1
/ Names of things. /
define PRES_NAME "pressure"
define UNITS "units"
define DEGREES_EAST "degrees_east"
define DEGREES_NORTH "degrees_north"
/ These are used to calculate the values we expect to find. /
define SAMPLE_PRESSURE 101403.1
define START_LAT 26.0
define START_LON -84.0
/ For the units attributes. /
define UNITS "units"
define PRES_UNITS "Pa"
define LAT_UNITS "degrees_north"
define LON_UNITS "degrees_east"
define MAX_ATT_LEN 80
/* Handle errors by printing an error message and exiting with a
- non-zero status. */
define ERR(e) {printf("Error: %s\n", nc_strerror(e)); return 2;}
int main() { int ncid, pres_varid; int lat_varid, lon_varid;
/ The start and count arrays will tell the netCDF library where to read our data. / size_t start[NDIMS], count[NDIMS];
/ Program variables to hold the data we will read. We will only need enough space to hold one timestep of data; one record. / float pres_in[NLVL][NLAT][NLON]; // float temp_in[NLVL][NLAT][NLON];
/ These program variables hold the latitudes and longitudes. / float lats[NLAT], lons[NLON];
/ Loop indexes. / int lvl, lat, lon, rec, i = 0;
/ Error handling. / int retval;
/ Open the file. / if ((retval = nc_open(FILE_NAME, NC_NOWRITE, &ncid))) ERR(retval);
/* Get the varids of the latitude and longitude coordinate
- variables. */ if ((retval = nc_inq_varid(ncid, LAT_NAME, &lat_varid))) ERR(retval); if ((retval = nc_inq_varid(ncid, LON_NAME, &lon_varid))) ERR(retval);
/ Read the coordinate variable data. / if ((retval = nc_get_var_float(ncid, lat_varid, &lats[0]))) ERR(retval); if ((retval = nc_get_var_float(ncid, lon_varid, &lons[0]))) ERR(retval);
/ Check the coordinate variable data. / for (lat = 0; lat < NLAT; lat++) if (lats[lat] != START_LAT + 5.lat) return 2; for (lon = 0; lon < NLON; lon++) if (lons[lon] != START_LON + 5.lon) return 2;
/* Get the varids of the pressure and temperature netCDF
- variables. */ if ((retval = nc_inq_varid(ncid, PRES_NAME, &pres_varid))) ERR(retval);
/* Read the data. Since we know the contents of the file we know
- that the data arrays in this program are the correct size to
- hold one timestep. */ count[0] = 1; count[1] = NLVL; count[2] = NLAT; count[3] = NLON; start[1] = 0; start[2] = 0; start[3] = 0;
/ Read and check one record at a time. / for (rec = 0; rec < NREC; rec++) { start[0] = rec; if ((retval = nc_get_vara_float(ncid, pres_varid, start, count, &pres_in[0][0][0]))) ERR(retval);
|/ Check the data. / i = 0; for (lvl = 0; lvl < NLVL; lvl++) for (lat = 0; lat < NLAT; lat++) for (lon = 0; lon < NLON; lon++) { if (pres_in[lvl][lat][lon] != SAMPLE_PRESSURE + i) return 2; //if (temp_in[lvl][lat][lon] != SAMPLE_TEMP + i) // return 2; i++; } } / next record / / Close the file. / if ((retval = nc_close(ncid))) ERR(retval); |
return 0; }
I am getting an error indicating that an undefined filter was encountered. I checked and it seems to be NC_ENOFILTER error I believe. I also believe that this error stems from the. nc_get_vara_float() function. I ran the pres_temp_4D_rd.c example with no error similar to this.
I attempted to check the filters of the .nc file with: ncdump -h | grep FILTERS
Nothing came up for filters in the ncdump. I also tried: ncdump -s
I did not see a "chunking" section but did see a variable like this "ChunkSizes = 1, 1, 45, 57 ; for my pressure data variable.
Any reason that this is happening ? Can't seem to find much info on this issue.
— Reply to this email directly, view it on GitHub https://github.com/Unidata/netcdf-c/issues/2636, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAG47W2TVRGMZ562M2IH7W3WY3525ANCNFSM6AAAAAAVFGLL2E. You are receiving this because you are subscribed to this thread.Message ID: @.***>
tayloj25
I am using netcdf-c-4.9.1 with hdf5-1.12.2 and zlib-1.2.13. I am attempting to run quite a similar example to the pres_temp_4D_rd.c example on the unidata website. My dimensions are slightly different though. My code appears as below:
include
include
include
/ This is the name of the data file we will read. /
define FILE_NAME "pressure_data.nc"
define NDIMS 4
define NLAT 45
define NLON 57
define LAT_NAME "y"
define LON_NAME "x"
define NREC 1
define REC_NAME "time"
define LVL_NAME "z"
define NLVL 1
/ Names of things. /
define PRES_NAME "pressure"
define UNITS "units"
define DEGREES_EAST "degrees_east"
define DEGREES_NORTH "degrees_north"
/ These are used to calculate the values we expect to find. /
define SAMPLE_PRESSURE 101403.1
define START_LAT 26.0
define START_LON -84.0
/ For the units attributes. /
define UNITS "units"
define PRES_UNITS "Pa"
define LAT_UNITS "degrees_north"
define LON_UNITS "degrees_east"
define MAX_ATT_LEN 80
/* Handle errors by printing an error message and exiting with a
define ERR(e) {printf("Error: %s\n", nc_strerror(e)); return 2;}
int main() { int ncid, pres_varid; int lat_varid, lon_varid;
/ The start and count arrays will tell the netCDF library where to read our data. / size_t start[NDIMS], count[NDIMS];
/ Program variables to hold the data we will read. We will only need enough space to hold one timestep of data; one record. / float pres_in[NLVL][NLAT][NLON]; // float temp_in[NLVL][NLAT][NLON];
/ These program variables hold the latitudes and longitudes. / float lats[NLAT], lons[NLON];
/ Loop indexes. / int lvl, lat, lon, rec, i = 0;
/ Error handling. / int retval;
/ Open the file. / if ((retval = nc_open(FILE_NAME, NC_NOWRITE, &ncid))) ERR(retval);
/* Get the varids of the latitude and longitude coordinate
variables. */ if ((retval = nc_inq_varid(ncid, LAT_NAME, &lat_varid))) ERR(retval); if ((retval = nc_inq_varid(ncid, LON_NAME, &lon_varid))) ERR(retval);
/ Read the coordinate variable data. / if ((retval = nc_get_var_float(ncid, lat_varid, &lats[0]))) ERR(retval); if ((retval = nc_get_var_float(ncid, lon_varid, &lons[0]))) ERR(retval);
/ Check the coordinate variable data. / for (lat = 0; lat < NLAT; lat++) if (lats[lat] != START_LAT + 5.lat) return 2; for (lon = 0; lon < NLON; lon++) if (lons[lon] != START_LON + 5.lon) return 2;
/* Get the varids of the pressure and temperature netCDF
variables. */ if ((retval = nc_inq_varid(ncid, PRES_NAME, &pres_varid))) ERR(retval);
/* Read the data. Since we know the contents of the file we know
hold one timestep. */ count[0] = 1; count[1] = NLVL; count[2] = NLAT; count[3] = NLON; start[1] = 0; start[2] = 0; start[3] = 0;
/ Read and check one record at a time. / for (rec = 0; rec < NREC; rec++) { start[0] = rec; if ((retval = nc_get_vara_float(ncid, pres_varid, start, count, &pres_in[0][0][0]))) ERR(retval);
/ Check the data. / i = 0; for (lvl = 0; lvl < NLVL; lvl++) for (lat = 0; lat < NLAT; lat++) for (lon = 0; lon < NLON; lon++) { if (pres_in[lvl][lat][lon] != SAMPLE_PRESSURE + i) return 2; //if (temp_in[lvl][lat][lon] != SAMPLE_TEMP + i) // return 2; i++; }
} / next record /
/ Close the file. / if ((retval = nc_close(ncid))) ERR(retval);
return 0; }
I am getting an error indicating that an undefined filter was encountered. I checked and it seems to be NC_ENOFILTER error I believe. I also believe that this error stems from the. nc_get_vara_float() function. I ran the pres_temp_4D_rd.c example with no error similar to this.
I attempted to check the filters of the .nc file with: ncdump -h | grep FILTERS
Nothing came up for filters in the ncdump. I also tried: ncdump -s
I did not see a "chunking" section but did see a variable like this "ChunkSizes = 1, 1, 45, 57 ; for my pressure data variable.
Any reason that this is happening ? Can't seem to find much info on this issue.