Project

General

Profile

CDO error :

Added by Alice Guccione 3 days ago

I am new to cdo and I am trying to compute the 95 percentile with timpctl over a list of 42 years of data. Firs I compute the minfile and maxfile like this:

cdo -b F32 -timmin -setrtomiss,0,0.099999999999 [ -cat /.../{1981..2022}/????????00/PREC_ACC_NC_*.nc ] minfile_WET_1981_2022_years.nc 
cdo -b F32 -timmax -setrtomiss,0,0.099999999999 [ -cat /.../{1981..2022}/????????00/PREC_ACC_NC_*.nc ] maxfile_WET_1981_2022_years.nc

and then I compute the 95 percentile with:

cdo -b F32 -timpctl,95 -setrtomiss,0,0.099999999999 [ -cat /.../CHAPTER/{1981..2022}/????????00/PREC_ACC_NC_*.nc ] minfile_WET_1981_2022_years.nc maxfile_WET_1981_2022_years.nc perc95_WET_1981_2022_years.nc

This last line raises the error:

cdo: ../../../src/percentiles_hist.cc:218: double histGetBin(int, double, const T*) [with T = short unsigned int]: Assertion `i - 1 < nbins' failed.
/var/spool/slurmd/job64492/slurm_script: line 7: 1727709 Aborted                 (core dumped) cdo -b F32 -timpctl,95 -setrtomiss,0,0.099999999999 [ -cat /data/datasets/obs/CHAPTER/{1981..2022}/????????00/PREC_ACC_NC_*.nc ] minfile_WET_1981_2022_years.nc maxfile_WET_1981_2022_years.nc perc95_WET_1981_2022_years.nc

I don't understand why, since when I work with just the first 21 files or the second 21 files it works perfectly, for example the two lines below work perfectly:

cdo -b F32 -timpctl,95 -setrtomiss,0,0.099999999999 [ -cat /.../{1981..2001}/????????00/PREC_ACC_NC_*.nc ] minfile_WET_1981_2001_years.nc maxfile_WET_1981_2001_years.nc perc95_WET_1981_2001_years.nc
cdo -b F32 -timpctl,95 -setrtomiss,0,0.099999999999 [ -cat /.../{2002..2022}/????????00/PREC_ACC_NC_*.nc ] minfile_WET_2002_2022_years.nc maxfile_WET_2002_2022_years.nc perc95_WET_1_year_2002_2022.nc

Thanks in advance for the help!