cloudnetpy 1.85.0__tar.gz → 1.85.1__tar.gz

{cloudnetpy-1.85.0/cloudnetpy.egg-info → cloudnetpy-1.85.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: cloudnetpy
-Version: 1.85.0
+Version: 1.85.1
 Summary: Python package for Cloudnet processing
 Author: Simo Tukiainen
 License: MIT License

{cloudnetpy-1.85.0 → cloudnetpy-1.85.1}/cloudnetpy/instruments/rpg.py

@@ -83,6 +83,10 @@ def rpg2nc(
     valid_ind = fmcw.add_zenith_angle()
     fmcw.screen_time_indices(valid_ind)
     fmcw.add_height()
+    if len(np.unique(fmcw.data["time"][:].astype("f4"))) != len(fmcw.data["time"][:]):
+        msg = "Convert time to f8 to keep values unique in netCDF"
+        logging.info(msg)
+        fmcw.data["time"].data_type = "f8"
     attributes = output.add_time_attribute(RPG_ATTRIBUTES, fmcw.date)
     output.update_attributes(fmcw.data, attributes)
     output.save_level1b(fmcw, output_file, uuid)
@@ -137,19 +141,17 @@ def _reduce_header(header: dict) -> dict:
     """Removes duplicate header data. Otherwise, we would need n_files dimension."""
     reduced_header = {}
     for key, data in header.items():
+        # Handle outliers in latitude and longitude (e.g. Galati 2024-02-11):
+        if key in ("latitude", "longitude"):
+            reduced_header[key] = ma.median(data)
+            continue
         first_profile_value = data[0]
         is_identical_value = bool(
             np.isclose(data, first_profile_value, rtol=1e-2).all(),
         )
         if is_identical_value is False:
             msg = f"Inconsistent header: {key}: {data}"
-            if key in (
-                "latitude",
-                "longitude",
-                "sample_duration",
-                "calibration_interval",
-                "noise_threshold",
-            ):
+            if key in ("sample_duration", "calibration_interval", "noise_threshold"):
                 logging.warning(msg)
             else:
                 raise InconsistentDataError(msg)

{cloudnetpy-1.85.0 → cloudnetpy-1.85.1}/cloudnetpy/products/epsilon.py

@@ -138,12 +138,12 @@ def _horizontal_wind_from_doppler_lidar_file(
         t = np.broadcast_to(time[:, None], mask.shape)[~mask]
         h = np.broadcast_to(height[None, :], mask.shape)[~mask]
 
-        if len(np.unique(t)) < 2 or len(np.unique(h)) < 2:
+        if len(np.unique(t)) < 3 or len(np.unique(h)) < 3:
             msg = "Not enough unique values for interpolation"
             raise ValidTimeStampError(msg)
 
-        interp_linear = LinearNDInterpolator(list(zip(t, h, strict=False)), V[~mask])
-        interp_nearest = NearestNDInterpolator(list(zip(t, h, strict=False)), V[~mask])
+        interp_linear = LinearNDInterpolator(list(zip(t, h, strict=True)), V[~mask])
+        interp_nearest = NearestNDInterpolator(list(zip(t, h, strict=True)), V[~mask])
         T, H = np.meshgrid(time, height, indexing="ij")
         V_linear = interp_linear(T, H)
         V_nearest = interp_nearest(T, H)

{cloudnetpy-1.85.0 → cloudnetpy-1.85.1}/cloudnetpy/version.py

@@ -1,4 +1,4 @@
 MAJOR = 1
 MINOR = 85
-PATCH = 0
+PATCH = 1
 __version__ = f"{MAJOR}.{MINOR}.{PATCH}"

{cloudnetpy-1.85.0 → cloudnetpy-1.85.1/cloudnetpy.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: cloudnetpy
-Version: 1.85.0
+Version: 1.85.1
 Summary: Python package for Cloudnet processing
 Author: Simo Tukiainen
 License: MIT License