GRIB message validity checker: Further checks

ecmwf · Dec 13, 2024 · c6caca4 · c6caca4
1 parent c7f6a16
commit c6caca4
Show file tree

Hide file tree

Showing 2 changed files with 53 additions and 20 deletions.
diff --git a/src/accessor/grib_accessor_class_message_is_valid.cc b/src/accessor/grib_accessor_class_message_is_valid.cc
@@ -20,8 +20,30 @@ void grib_accessor_message_is_valid_t::init(const long l, grib_arguments* c)
     length_ = 0;
 }
 
+static int check_field_values(grib_handle* h)
+{
+    //printf("DEBUG  %s \n", __func__);
+    int ret = GRIB_SUCCESS;
+    double* values = NULL;
+    size_t size = 0;
+    grib_context* c = h->context;
+
+    if ((ret = grib_get_size(h, "values", &size)) != GRIB_SUCCESS)
+        return ret;
+    values = (double*)grib_context_malloc_clear(c, size * sizeof(double));
+    if (!values)
+        return GRIB_OUT_OF_MEMORY;
+
+    if ((ret = grib_get_double_array(h, "values", values, &size)) != GRIB_SUCCESS) {
+        return ret;
+    }
+    grib_context_free(c, values);
+    return GRIB_SUCCESS;
+}
+
 static int check_grid_pl_array(grib_handle* h)
 {
+    // printf("DEBUG  %s \n", __func__);
     int ret = GRIB_SUCCESS;
     long Ni = 0,plpresent  = 0;
     long* pl = NULL; // pl array
@@ -31,7 +53,7 @@ static int check_grid_pl_array(grib_handle* h)
     // is there a PL array?
     ret = grib_get_long(h, "PLPresent", &plpresent);
     if (ret != GRIB_SUCCESS || plpresent == 0)
-        return GRIB_SUCCESS;
+        return GRIB_SUCCESS; // No PL array. So nothing to do
 
     if ((ret = grib_get_size(h, "pl", &plsize)) != GRIB_SUCCESS)
         return ret;
@@ -64,19 +86,31 @@ static int check_grid_pl_array(grib_handle* h)
     return GRIB_SUCCESS;
 }
 
+typedef int (*proj_func)(grib_handle*);
+static proj_func check_functions[] = {
+    check_field_values,
+    check_grid_pl_array
+};
+
 int grib_accessor_message_is_valid_t::unpack_long(long* val, size_t* len)
 {
+    int ret = 0;
     grib_handle* h = grib_handle_of_accessor(this);
     *len = 1;
-    *val = 0;
-
-    // TODO: Move all checks to dedicated classes, e.g.,
-    //  GribDataIntegrityChecker
-    //  BufrDataIntegrityChecker
-    // etc.
-    int ret = check_grid_pl_array(h); // Just prototyping
-    if (ret == GRIB_SUCCESS) {
-        *val = 1;
+    *val = 1; // Assume message is valid
+
+    if (h->product_kind != PRODUCT_GRIB) {
+        grib_context_log(h->context, GRIB_LOG_ERROR, "Validity checks only implemented for GRIB messages");
+        *val = 0;
+        return GRIB_NOT_IMPLEMENTED;
+    }
+
+    const size_t num_check_functions = sizeof(check_functions) / sizeof(check_functions[0]);
+    for (size_t i = 0; i < num_check_functions; ++i) {
+        proj_func cfunc = check_functions[i];
+        if ((ret = cfunc(h)) != GRIB_SUCCESS) {
+            *val = 0; // check failed
+        }
     }
 
     return GRIB_SUCCESS;

diff --git a/tests/grib_check_message_validity.sh b/tests/grib_check_message_validity.sh
@@ -17,6 +17,8 @@ tempFilt=temp.$label.filt
 sample=$ECCODES_SAMPLES_PATH/reduced_gg_pl_32_grib2.tmpl
 grib_check_key_equals $sample   isMessageValid 1
 
+# Check reduced Gaussian grid Ni
+# ------------------------------
 cat >$tempFilt<<EOF
    set Ni = 0;
    assert ( isMessageValid == 0 );
@@ -27,7 +29,8 @@ ${tools_dir}/grib_filter -o $tempGrib $tempFilt $sample
 grib_check_key_equals $tempGrib isMessageValid 0
 grib_check_key_equals $sample   isMessageValid 1
 
-
+# Check reduced Gaussian grid pl
+# ------------------------------
 cat >$tempFilt<<EOF
    meta pl_elem4 element(pl, 4);
    set pl_elem4 = 0;
@@ -38,15 +41,11 @@ ${tools_dir}/grib_filter -o $tempGrib $tempFilt $sample
 
 grib_check_key_equals $tempGrib isMessageValid 0
 
-# # Set wrong Nj. Should fail
-# input=$samples_dir/reduced_gg_pl_32_grib2.tmpl
-# ${tools_dir}/grib_set -s Nj=1 $input $tempGrib
-# set +e
-# ${tools_dir}/grib_check_gaussian_grid -v $tempGrib 2> $tempText
-# status=$?
-# set -e
-# [ $status -ne 0 ]
-# grep -q "Nj is 1 but should be 2\*N" $tempText
+
+# Check data values
+# ------------------------------
+${tools_dir}/grib_set -s bitsPerValue=25 $data_dir/sample.grib2 $tempGrib
+grib_check_key_equals $tempGrib isMessageValid 0
 
 
 # Clean up