Feature #2880 point2grid qc

data/config/Point2GridConfig_default

@@ -22,11 +22,12 @@ obs_window = {
 
 ////////////////////////////////////////////////////////////////////////////////
 
-
 //
-// Observation message type
+// Point observation filtering options
 //
-message_type = [];
+message_type    = [];
+obs_quality_inc = [];
+obs_quality_exc = [];
 
 ////////////////////////////////////////////////////////////////////////////////
 
@@ -72,10 +73,6 @@ var_name_map = [
 
 ////////////////////////////////////////////////////////////////////////////////
 
-quality_mark_thresh = 2;
-
-////////////////////////////////////////////////////////////////////////////////
-
 tmp_dir = "/tmp";
 version = "V12.0.0";
 

internal/test_unit/config/Point2GridConfig_WINDS

@@ -0,0 +1,2 @@
+obs_quality_inc = [ "2" ];
+obs_quality_exc = [ "1" ];

internal/test_unit/config/Point2GridConfig_valid_time

@@ -22,11 +22,12 @@ valid_time = "20201022_173000";
 
 ////////////////////////////////////////////////////////////////////////////////
 
-
 //
-// Observation message type
+// Point observation filtering options
 //
-//message_type = [];
+message_type    = [];
+obs_quality_inc = [];
+obs_quality_exc = [];
 
 ////////////////////////////////////////////////////////////////////////////////
 
@@ -72,10 +73,6 @@ valid_time = "20201022_173000";
 
 ////////////////////////////////////////////////////////////////////////////////
 
-//quality_mark_thresh = 2;
-
-////////////////////////////////////////////////////////////////////////////////
-
 tmp_dir = "/tmp";
 version = "V12.0.0";
 

internal/test_unit/xml/unit_point2grid.xml

@@ -46,6 +46,23 @@
     </output>
   </test>
 
+  <test name="point2grid_pb2nc_WINDS">
+    <exec>&MET_BIN;/point2grid</exec>
+    <param> \
+      &OUTPUT_DIR;/pb2nc/ndas.20120409.t12z.prepbufr.tm00.nc \
+      G212 \
+      &OUTPUT_DIR;/point2grid/pb2nc_WINDS.nc \
+      -field 'name="UGRD"; level="*";' \
+      -field 'name="VGRD"; level="*";' \
+      -name UWIND,VWIND \
+    -config &CONFIG_DIR;/Point2GridConfig_WINDS \
+      -v 3 
+    </param>
+    <output>
+      <grid_nc>&OUTPUT_DIR;/point2grid/pb2nc_WINDS.nc</grid_nc>
+    </output>
+  </test>
+
   <test name="point2grid_pb2nc_TMP">
     <exec>&MET_BIN;/point2grid</exec>
     <param> \
@@ -130,7 +147,7 @@
       G212 \
       &OUTPUT_DIR;/point2grid/point2grid_GOES_16_AOD_TO_G212_compute.nc \
       -field 'name="AOD";  level="(*,*)";' \
-      -qc 0,1 -method MAX \
+      -goes_qc 0,1 -method MAX \
       -v 1
     </param>
     <output>
@@ -167,7 +184,7 @@
       &OUTPUT_DIR;/point2grid/point2grid_GOES_16_ADP.nc \
       -field 'name="AOD_Smoke";  level="(*,*)";' \
       -adp &DATA_DIR_MODEL;/goes_16/OR_ABI-L2-ADPC-M6_G16_s20192662141196_e20192662143569_c20192662144526.nc \
-      -qc 0,1 -method MAX \
+      -goes_qc 0,1 -method MAX \
       -v 1
     </param>
     <output>
@@ -185,7 +202,7 @@
       G212 \
       &OUTPUT_DIR;/point2grid/point2grid_GOES_16_AOD_TO_G212_grid_map.nc \
       -field 'name="AOD";  level="(*,*)";' \
-      -qc 0,1,2 -method MAX \
+      -goes_qc 0,1,2 -method MAX \
       -v 1
     </param>
     <output>
@@ -205,7 +222,7 @@
       G212 \
       &OUTPUT_DIR;/point2grid/point2grid_GOES_16_AOD_TO_G212.nc \
       -field 'name="AOD";  level="(*,*)";' \
-      -qc 0,1,2 -method MAX \
+      -goes_qc 0,1,2 -method MAX \
       -v 1
     </param>
     <output>
@@ -224,7 +241,7 @@
       &OUTPUT_DIR;/point2grid/point2grid_GOES_16_ADP_Enterprise_high.nc \
       -field 'name="AOD_Smoke";  level="(*,*)";' \
       -adp &DATA_DIR_MODEL;/goes_16/OR_ABI-L2-ADPC-M6_G16_s20241100001171_e20241100003544_c20241100006361.nc \
-      -qc 0,1 -method MAX \
+      -goes_qc 0,1 -method MAX \
       -v 1
     </param>
     <output>
@@ -319,6 +336,7 @@
     </output>
   </test>
 
+  <!-- Use the old -qc option to test backward compatibility -->
   <test name="point2grid_GOES_16_AOD_TO_G212_unsigned">
     <exec>&MET_BIN;/point2grid</exec>
     <env>

src/basic/vx_util/crc_array.h

@@ -22,7 +22,6 @@
 #include <vector>
 
 #include "num_array.h"
-#include "int_array.h"
 #include "is_bad_data.h"
 #include "nint.h"
 #include "vx_cal.h"

src/basic/vx_util/int_array.h

@@ -20,8 +20,6 @@
 
 #include <iostream>
 
-#include "num_array.h"
-
 #include "crc_array.h"
 
 
@@ -34,6 +32,12 @@ typedef CRC_Array<int> IntArray;
 ////////////////////////////////////////////////////////////////////////
 
 
+extern ConcatString write_css(const IntArray &);
+
+
+////////////////////////////////////////////////////////////////////////
+
+
 #endif   /*  __INT_ARRAY_H__  */
 
 

src/basic/vx_util/num_array.cc

@@ -19,6 +19,7 @@
 
 #include "num_array.h"
 
+#include "int_array.h"
 #include "is_bad_data.h"
 #include "ptile.h"
 #include "nint.h"
@@ -1382,3 +1383,22 @@ ConcatString write_css_hhmmss(const NumArray &na)
 
 
 ////////////////////////////////////////////////////////////////////////
+
+
+ConcatString write_css(const IntArray &ia)
+
+{
+
+   ConcatString css;
+
+   for ( int i=0; i<ia.n(); ++i ) {
+      css << (i == 0 ? "" : ",") << ia[i];
+   }
+
+   return css;
+
+}
+
+
+////////////////////////////////////////////////////////////////////////
+