Load(): Floating point exception
Hello @pabretonniere and @omula ,
As reported by @obellprat , when running the following Load() call, you get a 'Floating point exception':
data <- Load('tos', 'm04o', 'ESA', paste(1998:2008, '1101', sep = ''), output = 'lonlat', configfile = '/home/Earth/obellpra/CFU/s2dverification/inst/config/IC3.conf', nprocs = 1)
After investigating, the problem happens when using NCO's ncrename. Load() interpolates and renames one of the ESA files with these system calls:
cdo remapcon,t106grid -selname,tos /esnas/obs/esa/monthly_mean/tos/tos_199201.nc tmp.nc
ncrename -v lat,latitude -d lat,latitude -v lon,longitude -d lon,longitude -v .lev,ensemble -d .lev,ensemble tmp.nc
After that, if running
cdo sinfo tmp.nc
you get the following error:
cdf_inq_dimlen : ncid = 65536 dimid = 0 length = 0
Error (cdf_inq_dimlen) : NetCDF: Invalid dimension ID or name
HDF5-DIAG: Error detected in HDF5 (1.8.14) thread 0:
#000: H5T.c line 1763 in H5Tclose(): not a datatype
major: Invalid arguments to routine
minor: Inappropriate type
Error (cdf_close) : NetCDF: HDF error
When Load() tries to read the longitudes and latitudes from this file via the 'ncdf4' library, it gets nothing and crashes.
Is it possible it's a bug in the NCO libraries? Is it possible it happens because the data is not formatted properly? If you see the 'comments' field in the NetCDF file you can see they warn it uses signed byte expressions that may not be supported by some libraries.