core
/
eccodes
mirror of https://github.com/ecmwf/eccodes.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

ECC-1158: GRIB encoding: Implement data quality checks (limits) for spectral fields

This commit is contained in:
Shahram Najm 2020-10-19 13:46:11 +01:00
parent ae2e728855
commit 45975c115b
1 changed files with 9 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/grib_accessor_class_data_complex_packing.c
View File

@ -672,6 +672,15 @@ static int pack_double(grib_accessor* a, const double* val, size_t* len)
return GRIB_INTERNAL_ERROR; return GRIB_INTERNAL_ERROR;
} }
/* Data Quality checks */
if (a->context->grib_data_quality_checks) {
/* First value is the field's average */
double min_val = val[0];
double max_val = min_val;
if ((ret = grib_util_grib_data_quality_check(gh, min_val, max_val)) != GRIB_SUCCESS)
return ret;
}
if (pen_j == sub_j) { if (pen_j == sub_j) {
double* values; double* values;
d = grib_power(decimal_scale_factor, 10); d = grib_power(decimal_scale_factor, 10);
@ -819,7 +828,6 @@ static int pack_double(grib_accessor* a, const double* val, size_t* len)
} }
else { else {
/* _test(val[i]*d,0); */ /* _test(val[i]*d,0); */
grib_encode_unsigned_long(hres, encode_float(val[i++]), &hpos, 8 * bytes); grib_encode_unsigned_long(hres, encode_float(val[i++]), &hpos, 8 * bytes);
/* _test(val[i]*d,0); */ /* _test(val[i]*d,0); */
grib_encode_unsigned_long(hres, encode_float(val[i++]), &hpos, 8 * bytes); grib_encode_unsigned_long(hres, encode_float(val[i++]), &hpos, 8 * bytes);