Remove a double-declaired "obs post filters" from satwnd goes-16 YAML.

NOAA-EMC · Nov 11, 2023 · a1892c3 · a1892c3
1 parent f5320c8
commit a1892c3
Showing 1 changed file with 21 additions and 8 deletions.
diff --git a/ush/ioda/bufr2ioda/bufr2ioda_satwind_amv_goes.py b/ush/ioda/bufr2ioda/bufr2ioda_satwind_amv_goes.py
@@ -229,6 +229,7 @@ def bufr_to_ioda(config, logger):
 
     height = np.full_like(pressure, fill_value=pressure.fill_value, dtype=np.float32)
     stnelev = np.full_like(pressure, fill_value=pressure.fill_value, dtype=np.float32)
+    obserr = np.full_like(wspd, fill_value=wspd.fill_value, dtype=np.float32)
 
     end_time = time.time()
     running_time = end_time - start_time
@@ -275,6 +276,7 @@ def bufr_to_ioda(config, logger):
             pressure2 = pressure[mask]
             height2 = height[mask]
             stnelev2 = stnelev[mask]
+            obserr2 = obserr[mask]
 
             # Processing Center
             ogce2 = ogce[mask]
@@ -428,17 +430,28 @@ def bufr_to_ioda(config, logger):
                 .write_attr('long_name', 'Station Elevation') \
                 .write_data(stnelev2)
 
-            # Wind Speed
-            obsspace.create_var('ObsValue/windSpeed', dtype=wspd2.dtype, fillval=wspd2.fill_value) \
+            # Observation Error for U-Wind Component
+            obsspace.create_var('ObsError/windEastward', dtype=obserr2.dtype, fillval=wspd2.fill_value) \
+                .write_attr('units', 'm s-1') \
+                .write_attr('long_name', 'Observation Error') \
+                .write_data(obserr2)
+
+            # Observation Error for V-Wind Component
+            obsspace.create_var('ObsError/windNorthward', dtype=obserr2.dtype, fillval=wspd2.fill_value) \
                 .write_attr('units', 'm s-1') \
-                .write_attr('long_name', 'Wind Speed') \
-                .write_data(wspd2)
+                .write_attr('long_name', 'Observation Error') \
+                .write_data(obserr2)
+            # Wind Speed
+#            obsspace.create_var('ObsValue/windSpeed', dtype=wspd2.dtype, fillval=wspd2.fill_value) \
+#                .write_attr('units', 'm s-1') \
+#                .write_attr('long_name', 'Wind Speed') \
+#                .write_data(wspd2)
 
             # Wind Direction
-            obsspace.create_var('ObsValue/windDirection', dtype=wdir2.dtype, fillval=wdir2.fill_value) \
-                .write_attr('units', 'degrees') \
-                .write_attr('long_name', 'Wind Direction') \
-                .write_data(wdir2)
+#            obsspace.create_var('ObsValue/windDirection', dtype=wdir2.dtype, fillval=wdir2.fill_value) \
+#                .write_attr('units', 'degrees') \
+#                .write_attr('long_name', 'Wind Direction') \
+#                .write_data(wdir2)
 
             # U-Wind Component
             obsspace.create_var('ObsValue/windEastward', dtype=uob2.dtype, fillval=wspd2.fill_value) \