pds 4.109

NOAA-EMC · Dec 6, 2024 · e2b1579 · e2b1579
1 parent f8a6d93
commit e2b1579
Show file tree

Hide file tree

Showing 3 changed files with 10 additions and 5 deletions.
diff --git a/src/grib2.h.in b/src/grib2.h.in
@@ -291,7 +291,7 @@ g2int g2_gribend(unsigned char *cgrib);
 /* Constants to help with templates. */
 #define G2C_MAX_GDS_TEMPLATE 31        /**< Maximum number of grid templates. */
 #define G2C_MAX_GDS_TEMPLATE_MAPLEN 28 /**< Maximum grid template map length. */
-#define G2C_MAX_PDS_TEMPLATE 86        /**< Maximum number of PDS templates. */
+#define G2C_MAX_PDS_TEMPLATE 87        /**< Maximum number of PDS templates. */
 #define G2C_MAX_PDS_TEMPLATE_MAPLEN 45 /**< Maximum template map length. */
 #define G2C_MAX_DRS_TEMPLATE 10        /**< Maximum number of DRS templates. */
 #define G2C_MAX_DRS_TEMPLATE_MAPLEN 18 /**< Maximum DRS template map length. */

diff --git a/src/pdstemplates.c b/src/pdstemplates.c
@@ -387,6 +387,10 @@ static const struct pdstemplate templatespds[G2C_MAX_PDS_TEMPLATE] =
         /** PDT 4.108 (12/04/2024)
         4.108: Analysis or forecast at a horizontal level or in a horizontal layer at a point in time for generic optical products */
         {108, 20, 0, {1, 1, 1, -1, -4, -1, -4, 1, 1, 1, 2, 1, 1, -4, 1, -1, -4, 1, -1, -4}},
+        /** PDT 4.109 (12/04/2024)
+        4.109: Individual ensemble forecast, control and perturbed, at a horizontal level or in a horizontal layer at a 
+        point in time for generic optical products */
+        {109, 23, 0, {1, 1, 1, -1, -4, -1, -4, 1, 1, 1, 2, 1, 1, -4, 1, -1, -4, 1, -1, -4, 1, 1, 1}},
 };
 
 /**

diff --git a/tests/tst_pdstemplates.c b/tests/tst_pdstemplates.c
@@ -631,24 +631,24 @@ main()
     printf("ok!\n");
     printf("Testing all getpdstemplate() calls with extensions...\n");
     {
-#define NUM_TEST 86
+#define NUM_TEST 87
         int number[NUM_TEST] = {
             0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 20, 30, 31, 40, 41,
             42, 43, 254, 1000, 1001, 1002, 1100, 1101, 32, 44, 45, 46, 47, 48, 50, 52,
             51, 91, 33, 34, 53, 54, 57, 60, 61, 35, 49, 55, 58, 59, 62, 63, 67, 68, 70,
             71, 72, 73, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 92,
-            93, 94, 95, 96, 97, 98, 99, 100, 103, 104, 108};
+            93, 94, 95, 96, 97, 98, 99, 100, 103, 104, 108, 109};
         int expected_maplen[NUM_TEST] = {
             15, 18, 17, 31, 30, 22, 16, 15, 29, 36, 30, 32, 31, 45, 44, 18, 19, 5, 5,
             16, 19, 30, 33, 3, 9, 16, 15, 15, 22, 10, 21, 24, 35, 38, 26, 21, 15, 16,
             36, 18, 32, 19, 22, 7, 24, 38, 6, 29, 21, 7, 24, 35, 38, 7, 7, 19, 21, 32,
             34, 17, 20, 31, 34, 27, 30, 36, 40, 39, 38, 17, 31, 24, 20, 34, 27, 27, 30,
-            28, 31, 31, 34, 14, 17, 20, 23, 20};
+            28, 31, 31, 34, 14, 17, 20, 23, 20, 23};
         int expected_needext[NUM_TEST] = {
             0, 0, 0, 1, 1, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0,
             0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 0, 0, 1,
             0, 1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1,
-            1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0};
+            1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0};
         int expected_map[NUM_TEST][48] = {
             {1, 1, 1, 1, 1, 2, 1, 1, -4, 1, -1, -4, 1, -1, -4},                                                                                               /* 0 */
             {1, 1, 1, 1, 1, 2, 1, 1, -4, 1, -1, -4, 1, -1, -4, 1, 1, 1},                                                                                      /* 1 */
@@ -736,6 +736,7 @@ main()
             {1, 1, 1, -1, -4, -1, -4, 1, 1, 1, 2, 1, 1, -4, 1, -1, -4, 1, -1, -4},                                                                            /* 103 */
             {1, 1, 1, -1, -4, -1, -4, 1, 1, 1, 2, 1, 1, -4, 1, -1, -4, 1, -1, -4, 1, 1, 1},                                                                   /* 104 */
             {1, 1, 1, -1, -4, -1, -4, 1, 1, 1, 2, 1, 1, -4, 1, -1, -4, 1, -1, -4},                                                                            /* 108 */
+            {1, 1, 1, -1, -4, -1, -4, 1, 1, 1, 2, 1, 1, -4, 1, -1, -4, 1, -1, -4, 1, 1, 1},                                                                   /* 109 */
         };
         int t;
         int ext_t = 0;