cloudnetpy 1.66.17__tar.gz → 1.66.18__tar.gz

{cloudnetpy-1.66.17/cloudnetpy.egg-info → cloudnetpy-1.66.18}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: cloudnetpy
-Version: 1.66.17
+Version: 1.66.18
 Summary: Python package for Cloudnet processing
 Author: Simo Tukiainen
 License: MIT License

{cloudnetpy-1.66.17 → cloudnetpy-1.66.18}/cloudnetpy/instruments/mira.py

@@ -72,6 +72,7 @@ def mira2nc(
             mira.screen_by_snr()
             mira.screen_invalid_ldr()
             mira.mask_invalid_data()
+            mira.mask_bad_angles()
             mira.add_time_and_range()
             mira.add_site_geolocation()
             mira.add_radar_specific_variables()
@@ -128,6 +129,18 @@ class Mira(NcRadar):
             )
             self.data["ldr"].data[:] = ma.masked
 
+    def mask_bad_angles(self) -> None:
+        """Masks clearly bad elevation and azimuth angles."""
+        limits = {
+            "elevation": (0, 180),
+            "azimuth_angle": (-360, 360),
+        }
+        for key, (lower, upper) in limits.items():
+            if (array := self.data[key].data) is not None:
+                margin = (upper - lower) * 0.05
+                array[array < (lower - margin)] = ma.masked
+                array[array > (upper + margin)] = ma.masked
+
 
 def _parse_input_files(input_files: str | list[str], temp_dir: str) -> tuple:
     if isinstance(input_files, list) or os.path.isdir(input_files):
@@ -219,7 +232,6 @@ def _get_keymap(filetype: str) -> dict:
                 ("SNRh2l", "SNR"),
                 ("elv", "elevation"),
                 ("azi", "azimuth_angle"),
-                ("aziv", "azimuth_velocity"),
                 ("nfft", "nfft"),
                 ("nave", "nave"),
                 ("prf", "prf"),
@@ -234,7 +246,6 @@ def _get_keymap(filetype: str) -> dict:
             "SNRg": "SNR",
             "elv": "elevation",
             "azi": "azimuth_angle",
-            "aziv": "azimuth_velocity",
             "nfft": "nfft",
             "nave": "nave",
             "prf": "prf",

{cloudnetpy-1.66.17 → cloudnetpy-1.66.18}/cloudnetpy/instruments/nc_radar.py

@@ -94,41 +94,26 @@ class NcRadar(DataSource, CloudnetInstrument):
         """Adds non-varying instrument zenith and azimuth angles and returns valid
         time indices.
         """
-        if "azimuth_velocity" in self.data:
-            azimuth = self.data["azimuth_velocity"].data
-            azimuth_reference = azimuth[0] if np.all(azimuth == azimuth[0]) else 0
-            azimuth_tolerance = 1e-6
-        else:
-            azimuth = self.data["azimuth_angle"].data
-            azimuth_reference = ma.median(azimuth)
-            azimuth_tolerance = 0.1
-
         elevation = self.data["elevation"].data
+        azimuth = self.data["azimuth_angle"].data
 
-        # Elevation is sometimes around -1000 indicating missing value
-        # Assume the instrument is pointing vertically in these cases
-        elevation[elevation < 0] = 90
-        zenith = 90 - elevation
+        elevation_diff = ma.diff(elevation, prepend=elevation[1])
+        azimuth_diff = ma.diff(azimuth, prepend=azimuth[1])
 
-        is_valid_zenith = np.abs(zenith) < 10
-        if not np.any(is_valid_zenith):
-            msg = "No valid zenith angles"
-            raise ValidTimeStampError(msg)
+        is_stable = np.abs(elevation - 90) < 1
+        is_stable &= np.abs(elevation_diff) < 1e-6
+        is_stable &= np.abs(azimuth_diff) < 1e-3
 
-        valid_zenith = zenith[is_valid_zenith]
-        is_stable_zenith = np.isclose(zenith, ma.median(valid_zenith), atol=0.1)
-        is_stable_azimuth = np.isclose(
-            azimuth,
-            azimuth_reference,
-            atol=azimuth_tolerance,
-        )
-        is_stable_profile = is_stable_zenith & is_stable_azimuth
+        # If scanning unit is broken, data are missing
+        # (assume it's vertically pointing)
+        missing_info = elevation.mask & azimuth.mask
+        is_stable[missing_info] = True
 
-        if ma.isMaskedArray(is_stable_profile):
-            is_stable_profile[is_stable_profile.mask] = False
-        n_removed = np.count_nonzero(~is_stable_profile)
+        if ma.isMaskedArray(is_stable):
+            is_stable[is_stable.mask] = False
+        n_removed = np.count_nonzero(~is_stable)
 
-        if n_removed >= len(zenith) - 1:
+        if n_removed >= len(elevation) - 1:
             msg = "Less than two profiles with valid zenith / azimuth angles"
             raise ValidTimeStampError(msg)
 
@@ -137,11 +122,10 @@ class NcRadar(DataSource, CloudnetInstrument):
                 "Filtering %s profiles due to varying zenith / azimuth angle",
                 n_removed,
             )
+        zenith = 90 - elevation
         self.append_data(zenith, "zenith_angle")
-        for key in ("elevation", "azimuth_velocity"):
-            if key in self.data:
-                del self.data[key]
-        return list(is_stable_profile)
+        del self.data["elevation"]
+        return list(is_stable)
 
     def add_radar_specific_variables(self) -> None:
         """Adds radar specific variables."""

{cloudnetpy-1.66.17 → cloudnetpy-1.66.18}/cloudnetpy/plotting/plotting.py

@@ -910,6 +910,7 @@ def _get_max_gap_in_minutes(figure_data: FigureData) -> float:
         "mwr-multi": 35,
         "weather-station": 12,
         "doppler-lidar-wind": 75,
+        "radar": 5,
     }
     return max_allowed_gap.get(file_type, 10)
 

{cloudnetpy-1.66.17 → cloudnetpy-1.66.18}/cloudnetpy/version.py

@@ -1,4 +1,4 @@
 MAJOR = 1
 MINOR = 66
-PATCH = 17
+PATCH = 18
 __version__ = f"{MAJOR}.{MINOR}.{PATCH}"

{cloudnetpy-1.66.17 → cloudnetpy-1.66.18/cloudnetpy.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: cloudnetpy
-Version: 1.66.17
+Version: 1.66.18
 Summary: Python package for Cloudnet processing
 Author: Simo Tukiainen
 License: MIT License