Merge branch 'develop' of github.com:ecmwf/eccodes into develop

joobog · joobog · commit 79100b4c57d5 · 2026-03-31T12:21:24.000+02:00
diff --git a/definitions/grib2/name.def b/definitions/grib2/name.def
@@ -847,6 +847,27 @@
 	 parameterCategory = 6 ;
 	 parameterNumber = 5 ;
 	}
+#Probability of a tropical storm
+'Probability of a tropical storm' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 6 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a hurricane
+'Probability of a hurricane' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 5 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a tropical depression
+'Probability of a tropical depression' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 7 ;
+	 probabilityType = 9 ;
+	}
 #Total precipitation of at least 25 mm
 'Total precipitation of at least 25 mm' = {
 	 discipline = 0 ;
diff --git a/definitions/grib2/name.lte33.def b/definitions/grib2/name.lte33.def
@@ -901,6 +901,27 @@
 	 parameterCategory = 6 ;
 	 parameterNumber = 5 ;
 	}
+#Probability of a tropical storm
+'Probability of a tropical storm' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 6 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a hurricane
+'Probability of a hurricane' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 5 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a tropical depression
+'Probability of a tropical depression' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 7 ;
+	 probabilityType = 9 ;
+	}
 #Total precipitation of at least 25 mm
 'Total precipitation of at least 25 mm' = {
 	 discipline = 0 ;
diff --git a/definitions/grib2/paramId.def b/definitions/grib2/paramId.def
@@ -847,6 +847,27 @@
 	 parameterCategory = 6 ;
 	 parameterNumber = 5 ;
 	}
+#Probability of a tropical storm
+'131089' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 6 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a hurricane
+'131090' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 5 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a tropical depression
+'131091' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 7 ;
+	 probabilityType = 9 ;
+	}
 #Total precipitation of at least 25 mm
 '131098' = {
 	 discipline = 0 ;
diff --git a/definitions/grib2/paramId.lte33.def b/definitions/grib2/paramId.lte33.def
@@ -901,6 +901,27 @@
 	 parameterCategory = 6 ;
 	 parameterNumber = 5 ;
 	}
+#Probability of a tropical storm
+'131089' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 6 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a hurricane
+'131090' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 5 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a tropical depression
+'131091' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 7 ;
+	 probabilityType = 9 ;
+	}
 #Total precipitation of at least 25 mm
 '131098' = {
 	 discipline = 0 ;
diff --git a/definitions/grib2/productDefinitionTemplateConcept.def b/definitions/grib2/productDefinitionTemplateConcept.def
@@ -33,16 +33,16 @@
 'aerosol_optical_ensemble_instantaneous'      = {productDefinitionTemplateNumber=49;}
 # 156 Average, accumulation, extreme values or other statistically processed values at a horizontal layer in a continuous or non-continuous time interval for optical properties of aerosol
 'aerosol_optical_deterministic_statistically_processed' = {productDefinitionTemplateNumber=156;}
-# 157 Individual ensemble forecast, control and perturbed at a horizontal level or in a horizontal layer in a continuous or non-continuous time interval for optical properties of aerosol
+# 157 Individual ensemble forecast, scontrol and perturbed at a horizontal level or in a horizontal layer in a continuous or non-continuous time interval for optical properties of aerosol
 'aerosol_optical_ensemble_statistically_processed' = {productDefinitionTemplateNumber=157;}
 
 # 80 Analysis or forecast at a horizontal level or in a horizontal layer at a point in time for optical properties of aerosol with source or sink
 'aerosol_optical_srcsink_deterministic_instantaneous'  = {productDefinitionTemplateNumber=80;}
 # 81 Individual ensemble forecast, control and perturbed, at a horizontal level or in a horizontal layer at a point in time for optical properties of aerosol with source or sink
 'aerosol_optical_srcsink_ensemble_instantaneous'  = {productDefinitionTemplateNumber=81;}
-# 158 Average, accumulation, extreme values or other statistically processed values at a continuous or non-continuous time interval for optical properties of aerosol with source or sink
+# 158 Average, accumulation, extreme values or other statistically processed values at a continuous or non-continuous time interval for optical properties of aerosol with source or sink
 'aerosol_optical_srcsink_deterministic_statistically_processed'  = {productDefinitionTemplateNumber=158;}
-# 159 Individual ensemble forecast, control and perturbed at a horizontal level or in a horizontal layer in a continuous or non-continuous time interval for optical properties of aerosol with source or sink
+# 159 Individual ensemble forecast, control and perturbed at a horizontal level or in a horizontal layer in a continuous or non-continuous time interval for optical properties of aerosol with source or sink
 'aerosol_optical_srcsink_ensemble_statistically_processed'  = {productDefinitionTemplateNumber=159;}
 
 # 50 Analysis or forecast at a horizontal level or in a horizontal layer at a point in time for aerosol
diff --git a/definitions/grib2/shortName.def b/definitions/grib2/shortName.def
@@ -847,6 +847,27 @@
 	 parameterCategory = 6 ;
 	 parameterNumber = 5 ;
 	}
+#Probability of a tropical storm
+'pts' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 6 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a hurricane
+'ph' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 5 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a tropical depression
+'ptd' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 7 ;
+	 probabilityType = 9 ;
+	}
 #Total precipitation of at least 25 mm
 'tpg25' = {
 	 discipline = 0 ;
diff --git a/definitions/grib2/shortName.lte33.def b/definitions/grib2/shortName.lte33.def
@@ -901,6 +901,27 @@
 	 parameterCategory = 6 ;
 	 parameterNumber = 5 ;
 	}
+#Probability of a tropical storm
+'pts' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 6 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a hurricane
+'ph' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 5 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a tropical depression
+'ptd' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 7 ;
+	 probabilityType = 9 ;
+	}
 #Total precipitation of at least 25 mm
 'tpg25' = {
 	 discipline = 0 ;
diff --git a/definitions/grib2/units.def b/definitions/grib2/units.def
@@ -847,6 +847,27 @@
 	 parameterCategory = 6 ;
 	 parameterNumber = 5 ;
 	}
+#Probability of a tropical storm
+'%' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 6 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a hurricane
+'%' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 5 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a tropical depression
+'%' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 7 ;
+	 probabilityType = 9 ;
+	}
 #Total precipitation of at least 25 mm
 '%' = {
 	 discipline = 0 ;
diff --git a/definitions/grib2/units.lte33.def b/definitions/grib2/units.lte33.def
@@ -901,6 +901,27 @@
 	 parameterCategory = 6 ;
 	 parameterNumber = 5 ;
 	}
+#Probability of a tropical storm
+'%' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 6 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a hurricane
+'%' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 5 ;
+	 probabilityType = 9 ;
+	}
+#Probability of a tropical depression
+'%' = {
+	 discipline = 0 ;
+	 parameterCategory = 191 ;
+	 parameterNumber = 7 ;
+	 probabilityType = 9 ;
+	}
 #Total precipitation of at least 25 mm
 '%' = {
 	 discipline = 0 ;
diff --git a/src/eccodes/accessor/BufrDataArray.cc b/src/eccodes/accessor/BufrDataArray.cc
@@ -424,7 +424,8 @@ int BufrDataArray::decode_string_array(grib_context* c, unsigned char* data, lon
     return ret;
 }
 
-grib_darray* BufrDataArray::decode_double_array(grib_context* c, unsigned char* data, long* pos,
+// numeric (integer or double) or codetable or flagtable array
+grib_darray* BufrDataArray::decode_numeric_array(grib_context* c, unsigned char* data, long* pos,
                                                 bufr_descriptor* bd, int canBeMissing, int* err)
 {
     grib_darray* ret = NULL;
@@ -589,7 +590,7 @@ static int descriptor_get_min_max(bufr_descriptor* bd, long width, long referenc
     return GRIB_SUCCESS;
 }
 
-int BufrDataArray::encode_double_array(grib_context* c, grib_buffer* buff, long* pos, bufr_descriptor* bd,
+int BufrDataArray::encode_numeric_array(grib_context* c, grib_buffer* buff, long* pos, bufr_descriptor* bd,
                                        grib_darray* dvalues)
 {
     int err = 0;
@@ -807,7 +808,7 @@ int BufrDataArray::encode_double_array(grib_context* c, grib_buffer* buff, long*
     return err;
 }
 
-int BufrDataArray::encode_double_value(grib_context* c, grib_buffer* buff, long* pos, bufr_descriptor* bd, double value)
+int BufrDataArray::encode_numeric_value(grib_context* c, grib_buffer* buff, long* pos, bufr_descriptor* bd, double value)
 {
     size_t lval;
     double maxAllowed, minAllowed;
@@ -888,7 +889,8 @@ char* BufrDataArray::decode_string_value(grib_context* c, unsigned char* data, l
     return sval;
 }
 
-double BufrDataArray::decode_double_value(grib_context* c, unsigned char* data, long* pos,
+// numeric (integer or double) or codetable or flagtable
+double BufrDataArray::decode_numeric_value(grib_context* c, unsigned char* data, long* pos,
                                                             bufr_descriptor* bd, int canBeMissing,
                                                             int* err)
 {
@@ -919,7 +921,6 @@ double BufrDataArray::decode_double_value(grib_context* c, unsigned char* data,
     return dval;
 }
 
-
 int decode_element(grib_context* c, BufrDataArray* self, int subsetIndex,
                    grib_buffer* b, unsigned char* data, long* pos, int i, bufr_descriptor* descriptor, long elementIndex,
                    grib_darray* dval, grib_sarray* sval)
@@ -987,12 +988,12 @@ int decode_element(grib_context* c, BufrDataArray* self, int subsetIndex,
             return GRIB_DECODING_ERROR;
         }
         if (self->compressedData_) {
-            dar = self->decode_double_array(c, data, pos, bd, self->canBeMissing_[i], &err);
+            dar = self->decode_numeric_array(c, data, pos, bd, self->canBeMissing_[i], &err);
             grib_vdarray_push(self->numericValues_, dar);
         }
         else {
             /* Uncompressed */
-            cdval = self->decode_double_value(c, data, pos, bd, self->canBeMissing_[i], &err);
+            cdval = self->decode_numeric_value(c, data, pos, bd, self->canBeMissing_[i], &err);
             grib_context_log(c, GRIB_LOG_DEBUG, "BUFR data decoding: \t %s = %g",
                              bd->shortName, cdval);
             grib_darray_push(dval, cdval);
@@ -1001,7 +1002,6 @@ int decode_element(grib_context* c, BufrDataArray* self, int subsetIndex,
     return err;
 }
 
-
 int decode_replication(grib_context* c, BufrDataArray* self, int subsetIndex, grib_buffer* buff,
                        unsigned char* data, long* pos, int i, long elementIndex, grib_darray* dval, long* numberOfRepetitions)
 {
@@ -1067,7 +1067,6 @@ int decode_replication(grib_context* c, BufrDataArray* self, int subsetIndex, gr
     return ret;
 }
 
-
 int BufrDataArray::encode_new_bitmap(grib_context* c, grib_buffer* buff, long* pos, int idx)
 {
     grib_darray* doubleValues = NULL;
@@ -1081,11 +1080,11 @@ int BufrDataArray::encode_new_bitmap(grib_context* c, grib_buffer* buff, long* p
     if (compressedData_) {
         doubleValues = grib_darray_new(1, 1);
         grib_darray_push(doubleValues, cdval);
-        err = encode_double_array(c, buff, pos, expanded_->v[idx], doubleValues);
+        err = encode_numeric_array(c, buff, pos, expanded_->v[idx], doubleValues);
         grib_darray_delete(doubleValues);
     }
     else {
-        err = encode_double_value(c, buff, pos, expanded_->v[idx], cdval);
+        err = encode_numeric_value(c, buff, pos, expanded_->v[idx], cdval);
     }
     return err;
 }
@@ -1179,17 +1178,16 @@ int encode_new_element(grib_context* c, BufrDataArray* self, int subsetIndex,
         if (self->compressedData_) {
             grib_darray* doubleValues = grib_darray_new(1, 1);
             grib_darray_push(doubleValues, cdval);
-            err = self->encode_double_array(c, buff, pos, bd, doubleValues);
+            err = self->encode_numeric_array(c, buff, pos, bd, doubleValues);
             grib_darray_delete(doubleValues);
         }
         else {
-            err = self->encode_double_value(c, buff, pos, bd, cdval);
+            err = self->encode_numeric_value(c, buff, pos, bd, cdval);
         }
     }
     return err;
 }
 
-
 int encode_new_replication(grib_context* c, BufrDataArray* self, int subsetIndex,
                            grib_buffer* buff, unsigned char* data, long* pos, int i, long elementIndex, grib_darray* dval, long* numberOfRepetitions)
 {
@@ -1253,7 +1251,6 @@ int encode_new_replication(grib_context* c, BufrDataArray* self, int subsetIndex
     return err;
 }
 
-
 int encode_element(grib_context* c, BufrDataArray* self, int subsetIndex,
                    grib_buffer* buff, unsigned char* data, long* pos, int i, bufr_descriptor* descriptor,
                    long elementIndex, grib_darray* dval, grib_sarray* sval)
@@ -1300,7 +1297,7 @@ int encode_element(grib_context* c, BufrDataArray* self, int subsetIndex,
     else {
         /* numeric or codetable or flagtable */
         if (self->compressedData_) {
-            err = self->encode_double_array(c, buff, pos, bd, self->numericValues_->v[elementIndex]);
+            err = self->encode_numeric_array(c, buff, pos, bd, self->numericValues_->v[elementIndex]);
             if (err) {
                 grib_darray* varr = self->numericValues_->v[elementIndex];
                 grib_context_log(c, GRIB_LOG_ERROR, "Encoding key '%s' ( code=%6.6ld width=%ld scale=%ld reference=%ld )",
@@ -1320,7 +1317,7 @@ int encode_element(grib_context* c, BufrDataArray* self, int subsetIndex,
                 grib_context_log(c, GRIB_LOG_ERROR, "Invalid subset index %d (number of subsets=%ld)", subsetIndex, self->numberOfSubsets_);
                 return GRIB_INVALID_ARGUMENT;
             }
-            err = self->encode_double_value(c, buff, pos, bd, self->numericValues_->v[subsetIndex]->v[elementIndex]);
+            err = self->encode_numeric_value(c, buff, pos, bd, self->numericValues_->v[subsetIndex]->v[elementIndex]);
             if (err) {
                 grib_context_log(c, GRIB_LOG_ERROR, "Cannot encode %s=%g (subset=%d)", /*subsetIndex starts from 0*/
                                  bd->shortName, self->numericValues_->v[subsetIndex]->v[elementIndex], subsetIndex + 1);
diff --git a/src/eccodes/accessor/BufrDataArray.h b/src/eccodes/accessor/BufrDataArray.h
diff --git a/src/eccodes/grib_bits_any_endian.cc b/src/eccodes/grib_bits_any_endian.cc
diff --git a/tests/grib_dump_debug.sh b/tests/grib_dump_debug.sh
diff --git a/tests/unit_tests.cc b/tests/unit_tests.cc

Original file line number	Diff line number	Diff line change
`@@ -424,7 +424,8 @@ int BufrDataArray::decode_string_array(grib_context* c, unsigned char* data, lon`
`424`	`424`	`return ret;`
`425`	`425`	`}`
`426`	`426`
`427`		`-grib_darray* BufrDataArray::decode_double_array(grib_context* c, unsigned char* data, long* pos,`
	`427`	`+// numeric (integer or double) or codetable or flagtable array`
	`428`	`+grib_darray* BufrDataArray::decode_numeric_array(grib_context* c, unsigned char* data, long* pos,`
`428`	`429`	`bufr_descriptor* bd, int canBeMissing, int* err)`
`429`	`430`	`{`
`430`	`431`	`grib_darray* ret = NULL;`
`@@ -589,7 +590,7 @@ static int descriptor_get_min_max(bufr_descriptor* bd, long width, long referenc`
`589`	`590`	`return GRIB_SUCCESS;`
`590`	`591`	`}`
`591`	`592`
`592`		`-int BufrDataArray::encode_double_array(grib_context* c, grib_buffer* buff, long* pos, bufr_descriptor* bd,`
	`593`	`+int BufrDataArray::encode_numeric_array(grib_context* c, grib_buffer* buff, long* pos, bufr_descriptor* bd,`
`593`	`594`	`grib_darray* dvalues)`
`594`	`595`	`{`
`595`	`596`	`int err = 0;`
`@@ -807,7 +808,7 @@ int BufrDataArray::encode_double_array(grib_context* c, grib_buffer* buff, long*`
`807`	`808`	`return err;`
`808`	`809`	`}`
`809`	`810`
`810`		`-int BufrDataArray::encode_double_value(grib_context* c, grib_buffer* buff, long* pos, bufr_descriptor* bd, double value)`
	`811`	`+int BufrDataArray::encode_numeric_value(grib_context* c, grib_buffer* buff, long* pos, bufr_descriptor* bd, double value)`
`811`	`812`	`{`
`812`	`813`	`size_t lval;`
`813`	`814`	`double maxAllowed, minAllowed;`
`@@ -888,7 +889,8 @@ char* BufrDataArray::decode_string_value(grib_context* c, unsigned char* data, l`
`888`	`889`	`return sval;`
`889`	`890`	`}`
`890`	`891`
`891`		`-double BufrDataArray::decode_double_value(grib_context* c, unsigned char* data, long* pos,`
	`892`	`+// numeric (integer or double) or codetable or flagtable`
	`893`	`+double BufrDataArray::decode_numeric_value(grib_context* c, unsigned char* data, long* pos,`
`892`	`894`	`bufr_descriptor* bd, int canBeMissing,`
`893`	`895`	`int* err)`
`894`	`896`	`{`
`@@ -919,7 +921,6 @@ double BufrDataArray::decode_double_value(grib_context* c, unsigned char* data,`
`919`	`921`	`return dval;`
`920`	`922`	`}`
`921`	`923`
`922`		`-`
`923`	`924`	`int decode_element(grib_context* c, BufrDataArray* self, int subsetIndex,`
`924`	`925`	`grib_buffer* b, unsigned char* data, long* pos, int i, bufr_descriptor* descriptor, long elementIndex,`
`925`	`926`	`grib_darray* dval, grib_sarray* sval)`
`@@ -987,12 +988,12 @@ int decode_element(grib_context* c, BufrDataArray* self, int subsetIndex,`
`987`	`988`	`return GRIB_DECODING_ERROR;`
`988`	`989`	`}`
`989`	`990`	`if (self->compressedData_) {`
`990`		`- dar = self->decode_double_array(c, data, pos, bd, self->canBeMissing_[i], &err);`
	`991`	`+ dar = self->decode_numeric_array(c, data, pos, bd, self->canBeMissing_[i], &err);`
`991`	`992`	`grib_vdarray_push(self->numericValues_, dar);`
`992`	`993`	`}`
`993`	`994`	`else {`
`994`	`995`	`/* Uncompressed */`
`995`		`- cdval = self->decode_double_value(c, data, pos, bd, self->canBeMissing_[i], &err);`
	`996`	`+ cdval = self->decode_numeric_value(c, data, pos, bd, self->canBeMissing_[i], &err);`
`996`	`997`	`grib_context_log(c, GRIB_LOG_DEBUG, "BUFR data decoding: \t %s = %g",`
`997`	`998`	`bd->shortName, cdval);`
`998`	`999`	`grib_darray_push(dval, cdval);`
`@@ -1001,7 +1002,6 @@ int decode_element(grib_context* c, BufrDataArray* self, int subsetIndex,`
`1001`	`1002`	`return err;`
`1002`	`1003`	`}`
`1003`	`1004`
`1004`		`-`
`1005`	`1005`	`int decode_replication(grib_context* c, BufrDataArray* self, int subsetIndex, grib_buffer* buff,`
`1006`	`1006`	`unsigned char* data, long* pos, int i, long elementIndex, grib_darray* dval, long* numberOfRepetitions)`
`1007`	`1007`	`{`
`@@ -1067,7 +1067,6 @@ int decode_replication(grib_context* c, BufrDataArray* self, int subsetIndex, gr`
`1067`	`1067`	`return ret;`
`1068`	`1068`	`}`
`1069`	`1069`
`1070`		`-`
`1071`	`1070`	`int BufrDataArray::encode_new_bitmap(grib_context* c, grib_buffer* buff, long* pos, int idx)`
`1072`	`1071`	`{`
`1073`	`1072`	`grib_darray* doubleValues = NULL;`
`@@ -1081,11 +1080,11 @@ int BufrDataArray::encode_new_bitmap(grib_context* c, grib_buffer* buff, long* p`
`1081`	`1080`	`if (compressedData_) {`
`1082`	`1081`	`doubleValues = grib_darray_new(1, 1);`
`1083`	`1082`	`grib_darray_push(doubleValues, cdval);`
`1084`		`- err = encode_double_array(c, buff, pos, expanded_->v[idx], doubleValues);`
	`1083`	`+ err = encode_numeric_array(c, buff, pos, expanded_->v[idx], doubleValues);`
`1085`	`1084`	`grib_darray_delete(doubleValues);`
`1086`	`1085`	`}`
`1087`	`1086`	`else {`
`1088`		`- err = encode_double_value(c, buff, pos, expanded_->v[idx], cdval);`
	`1087`	`+ err = encode_numeric_value(c, buff, pos, expanded_->v[idx], cdval);`
`1089`	`1088`	`}`
`1090`	`1089`	`return err;`
`1091`	`1090`	`}`
`@@ -1179,17 +1178,16 @@ int encode_new_element(grib_context* c, BufrDataArray* self, int subsetIndex,`
`1179`	`1178`	`if (self->compressedData_) {`
`1180`	`1179`	`grib_darray* doubleValues = grib_darray_new(1, 1);`
`1181`	`1180`	`grib_darray_push(doubleValues, cdval);`
`1182`		`- err = self->encode_double_array(c, buff, pos, bd, doubleValues);`
	`1181`	`+ err = self->encode_numeric_array(c, buff, pos, bd, doubleValues);`
`1183`	`1182`	`grib_darray_delete(doubleValues);`
`1184`	`1183`	`}`
`1185`	`1184`	`else {`
`1186`		`- err = self->encode_double_value(c, buff, pos, bd, cdval);`
	`1185`	`+ err = self->encode_numeric_value(c, buff, pos, bd, cdval);`
`1187`	`1186`	`}`
`1188`	`1187`	`}`
`1189`	`1188`	`return err;`
`1190`	`1189`	`}`
`1191`	`1190`
`1192`		`-`
`1193`	`1191`	`int encode_new_replication(grib_context* c, BufrDataArray* self, int subsetIndex,`
`1194`	`1192`	`grib_buffer* buff, unsigned char* data, long* pos, int i, long elementIndex, grib_darray* dval, long* numberOfRepetitions)`
`1195`	`1193`	`{`
`@@ -1253,7 +1251,6 @@ int encode_new_replication(grib_context* c, BufrDataArray* self, int subsetIndex`
`1253`	`1251`	`return err;`
`1254`	`1252`	`}`
`1255`	`1253`
`1256`		`-`
`1257`	`1254`	`int encode_element(grib_context* c, BufrDataArray* self, int subsetIndex,`
`1258`	`1255`	`grib_buffer* buff, unsigned char* data, long* pos, int i, bufr_descriptor* descriptor,`
`1259`	`1256`	`long elementIndex, grib_darray* dval, grib_sarray* sval)`
`@@ -1300,7 +1297,7 @@ int encode_element(grib_context* c, BufrDataArray* self, int subsetIndex,`
`1300`	`1297`	`else {`
`1301`	`1298`	`/* numeric or codetable or flagtable */`
`1302`	`1299`	`if (self->compressedData_) {`
`1303`		`- err = self->encode_double_array(c, buff, pos, bd, self->numericValues_->v[elementIndex]);`
	`1300`	`+ err = self->encode_numeric_array(c, buff, pos, bd, self->numericValues_->v[elementIndex]);`
`1304`	`1301`	`if (err) {`
`1305`	`1302`	`grib_darray* varr = self->numericValues_->v[elementIndex];`
`1306`	`1303`	`grib_context_log(c, GRIB_LOG_ERROR, "Encoding key '%s' ( code=%6.6ld width=%ld scale=%ld reference=%ld )",`
`@@ -1320,7 +1317,7 @@ int encode_element(grib_context* c, BufrDataArray* self, int subsetIndex,`
`1320`	`1317`	`grib_context_log(c, GRIB_LOG_ERROR, "Invalid subset index %d (number of subsets=%ld)", subsetIndex, self->numberOfSubsets_);`
`1321`	`1318`	`return GRIB_INVALID_ARGUMENT;`
`1322`	`1319`	`}`
`1323`		`- err = self->encode_double_value(c, buff, pos, bd, self->numericValues_->v[subsetIndex]->v[elementIndex]);`
	`1320`	`+ err = self->encode_numeric_value(c, buff, pos, bd, self->numericValues_->v[subsetIndex]->v[elementIndex]);`
`1324`	`1321`	`if (err) {`
`1325`	`1322`	`grib_context_log(c, GRIB_LOG_ERROR, "Cannot encode %s=%g (subset=%d)", /subsetIndex starts from 0/`
`1326`	`1323`	`bd->shortName, self->numericValues_->v[subsetIndex]->v[elementIndex], subsetIndex + 1);`