Add RHP product

Author: tobiasmarke

mwrpy/level2/get_ret_coeff.py

@@ -31,6 +31,8 @@ def get_mvr_coeff(
     c_list = get_coeff_list(site, prefix, coeff_files)
 
     coeff: dict = {}
+    if len(c_list) != 1:
+        return coeff
 
     if (str(c_list[0][-3:]).lower() == "ret") and (len(c_list) == 1):
         coeff = read_coeff_ascii(c_list[0])

mwrpy/level2/lev2_collocated.py

@@ -123,6 +123,8 @@ def generate_lev2_single(
                         "relative_humidity",
                         "relative_humidity_random_error",
                         "relative_humidity_systematic_error",
+                        "relative_humidity_quality_flag",
+                        "relative_humidity_quality_flag_status",
                     ),
                 ),
                 (
@@ -178,7 +180,7 @@ def generate_lev2_single(
 
                     copy_global(nc_lwp, nc_output, nc_lwp.ncattrs())
 
-            except IndexError:
+            except OSError:
                 logging.warning("No coefficient files for product " + prod)
 
         return nc_output
@@ -309,6 +311,7 @@ def generate_lev2_multi(
                 hum_file=abs_hum_file.name if prod not in ("2P02", "2P03") else None,
                 lwp_offset=(None, None),
                 coeff_files=coeff_files,
+                multi=True if prod == "2P04" else False,
             )
 
         with (

mwrpy/level2/lev2_meta_nc.py

@@ -302,6 +302,22 @@ def get_data_attributes(rpg_variables: dict, data_type: str, coeff: dict) -> dic
         units="1",
         dimensions=("time",),
     ),
+    "relative_humidity_quality_flag": MetaData(
+        long_name="Relative humidity quality flag",
+        units="1",
+        definition=DEFINITIONS_COM["quality_flag"],
+        comment="0 indicates data with good quality according to applied tests.\n"
+        "The list of (not) applied tests is encoded in quality_flag_status",
+        dimensions=("time",),
+    ),
+    "relative_humidity_quality_flag_status": MetaData(
+        long_name="Relative humidity quality flag status",
+        units="1",
+        definition=DEFINITIONS_COM["quality_flag_status"],
+        comment="Checks not executed in determination of quality_flag.\n"
+        "0 indicates quality check has been applied.",
+        dimensions=("time",),
+    ),
 }