From cfc96163b804d9236a01a46caf0823a131a2f543 Mon Sep 17 00:00:00 2001
From: Shahram Najm <masn@ecmwf.int>
Date: Wed, 30 Aug 2023 11:02:52 +0000
Subject: [PATCH] ECC-1674: Test added

---
 tests/CMakeLists.txt             |  1 +
 tests/grib_run_length_packing.cc | 75 ++++++++++++++++++++++++++++++++
 tests/grib_run_length_packing.sh | 13 ++++++
 3 files changed, 89 insertions(+)
 create mode 100644 tests/grib_run_length_packing.cc

diff --git a/tests/CMakeLists.txt b/tests/CMakeLists.txt
index 2e478eefc..1595e759e 100644
--- a/tests/CMakeLists.txt
+++ b/tests/CMakeLists.txt
@@ -32,6 +32,7 @@ list(APPEND test_c_bins
     grib_optimize_scaling
     grib_optimize_scaling_sh
     grib_set_force
+    grib_run_length_packing
     grib_partial_message
     grib_ecc-386
     grib_ecc-1467
diff --git a/tests/grib_run_length_packing.cc b/tests/grib_run_length_packing.cc
new file mode 100644
index 000000000..15a88619f
--- /dev/null
+++ b/tests/grib_run_length_packing.cc
@@ -0,0 +1,75 @@
+/*
+ * (C) Copyright 2005- ECMWF.
+ *
+ * This software is licensed under the terms of the Apache Licence Version 2.0
+ * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
+ *
+ * In applying this licence, ECMWF does not waive the privileges and immunities granted to it by
+ * virtue of its status as an intergovernmental organisation nor does it submit to any jurisdiction.
+ */
+#include <stdio.h>
+#include "eccodes.h"
+#undef NDEBUG
+#include <assert.h>
+
+int main(int argc, char *argv[])
+{
+    size_t str_len = 0;
+    long number_of_values=4, bits_per_value=8, max_level_value=3, number_of_level_values=5, decimal_scale_factor=1;
+    size_t level_values_length = number_of_level_values;
+    long level_values[] = {101L, 102L, 103L, 104L, 105L};
+    double *values = NULL;
+    size_t values_length = number_of_values;
+    values = (double*)malloc(values_length * sizeof(double));
+    int i = 0;
+    values[0] = 9999.0;
+
+    assert(argc == 2);
+
+    for (i=1; i<values_length; i++) {
+        values[i] = (double)((100.0 + i) / 10.0) ;
+    }
+    codes_handle *sample_handle = codes_grib_handle_new_from_samples(0, "GRIB2");
+    codes_handle *clone_handle = codes_handle_clone(sample_handle);
+    codes_handle_delete(sample_handle);
+    codes_set_long(clone_handle, "centre", 34);
+    codes_set_long(clone_handle, "significanceOfReferenceTime", 0);
+    codes_set_long(clone_handle, "dataDate", 20230830);
+    codes_set_long(clone_handle, "dataTime", 0);
+    codes_set_long(clone_handle, "productionStatusOfProcessedData", 1);
+    codes_set_long(clone_handle, "typeOfProcessedData", 0);
+    codes_set_long(clone_handle, "numberOfDataPoints", number_of_values);
+    codes_set_long(clone_handle, "interpretationOfNumberOfPoints", 0);
+    codes_set_long(clone_handle, "Ni", sqrt(number_of_values));
+    codes_set_long(clone_handle, "Nj", sqrt(number_of_values));
+
+    codes_set_long(clone_handle, "latitudeOfFirstGridPointInDegrees", 35.000000);
+    codes_set_long(clone_handle, "longitudeOfFirstGridPointInDegrees", 135.000000);
+    codes_set_long(clone_handle, "latitudeOfLastGridPointInDegrees", 35.000000 + (0.012500 * sqrt(number_of_values)));
+    codes_set_long(clone_handle, "longitudeOfLastGridPointInDegrees", 135.000000 + (0.008333 * sqrt(number_of_values)));
+    codes_set_long(clone_handle, "iDirectionIncrementInDegrees", 0.012500);
+    codes_set_long(clone_handle, "jDirectionIncrementInDegrees", 0.008333);
+
+    codes_set_long(clone_handle, "indicatorOfUnitOfTimeRange", 0);
+    codes_set_long(clone_handle, "stepUnits", 1);
+    codes_set_long(clone_handle, "forecastTime", 0);
+    codes_set_long(clone_handle, "stepRange", 0);
+
+    codes_set_string(clone_handle, "shortName", "t", &str_len);
+    codes_set_long(clone_handle, "numberOfValues", 1);
+    codes_set_string(clone_handle, "packingType", "grid_run_length", &str_len);
+    codes_set_long(clone_handle, "numberOfValues", number_of_values);
+    codes_set_long(clone_handle, "bitsPerValue", bits_per_value);
+    codes_set_long(clone_handle, "maxLevelValue", max_level_value);
+    codes_set_long(clone_handle, "numberOfLevelValues", number_of_level_values);
+    codes_set_long(clone_handle, "decimalScaleFactor", decimal_scale_factor);
+    codes_set_long_array(clone_handle, "levelValues", level_values, level_values_length);
+    //codes_set_long(clone_handle, "bitMapIndicator", 255);
+    //CODES_CHECK(codes_set_long(clone_handle, "bitmapPresent", 0), 0);
+    CODES_CHECK(codes_set_double_array(clone_handle, "values", values, values_length), 0);
+    codes_write_message(clone_handle, argv[1], "w");
+    codes_handle_delete(clone_handle);
+
+    free(values);
+    return 0;
+}
diff --git a/tests/grib_run_length_packing.sh b/tests/grib_run_length_packing.sh
index e9be4f4ac..33391ef7a 100755
--- a/tests/grib_run_length_packing.sh
+++ b/tests/grib_run_length_packing.sh
@@ -16,6 +16,8 @@ label="grib_run_length_packing_test"
 tempGrib=temp.$label.grib
 tempLog=temp.$label.log
 
+# Decoding
+# ------------------
 input=$data_dir/run_length_packing.grib2
 if [ ! -f "$input" ]; then
     echo " Input data file missing: $input"
@@ -30,4 +32,15 @@ stats=$(${tools_dir}/grib_get -F%.5f -p packingType,min,max,avg $input)
 grib_check_key_equals $input numberOfMissing 71493
 grib_check_key_equals $input missingValuesPresent 1
 
+# Encoding
+# -----------------
+$EXEC ${test_dir}/grib_run_length_packing $tempGrib
+${tools_dir}/grib_dump  $tempGrib
+grib_check_key_equals $tempGrib packingType grid_run_length
+
+stats=`${tools_dir}/grib_get -M -F%.1f -p min,max,avg $tempGrib`
+[ "$stats" = "10.1 10.3 10.2" ]
+
+
+# Clean up
 rm -f $tempGrib $tempLog