cloudnetpy 1.66.17__tar.gz → 1.66.19__tar.gz

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

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

{cloudnetpy-1.66.17 → cloudnetpy-1.66.19}/cloudnetpy/instruments/instruments.py

@@ -142,11 +142,18 @@ HALO = Instrument(
     model="StreamLine",
 )
 
-WINDCUBE = Instrument(
-    manufacturer="Leosphere",
+WINDCUBE_WLS200S = Instrument(
+    manufacturer="Vaisala",
+    domain="lidar",
+    category="Doppler lidar",
+    model="WindCube WLS200S",
+)
+
+WINDCUBE_WLS70 = Instrument(
+    manufacturer="Vaisala",
     domain="lidar",
     category="Doppler lidar",
-    model="WindCube",
+    model="WindCube WLS70",
 )
 
 PARSIVEL2 = Instrument(

{cloudnetpy-1.66.17 → cloudnetpy-1.66.19}/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.19}/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.19}/cloudnetpy/plotting/plotting.py

@@ -45,6 +45,8 @@ class PlotParameters:
         plot_meta: Additional metadata for the plot.
         raise_on_empty: Whether to raise an error if no data is found for a
             plotted variable.
+        minor_ticks: Whether to display minor ticks on the x-axis.
+
     """
 
     dpi: float = 120
@@ -58,6 +60,7 @@ class PlotParameters:
     footer_text: str | None = None
     plot_meta: PlotMeta | None = None
     raise_on_empty: bool = False
+    minor_ticks: bool = False
 
 
 class Dimensions:
@@ -210,6 +213,9 @@ class SubPlot:
         ]
         self.ax.set_xticks(np.arange(0, 25, resolution, dtype=int))
         self.ax.set_xticklabels(x_tick_labels, fontsize=12)
+        if self.options.minor_ticks:
+            self.ax.xaxis.set_minor_locator(AutoMinorLocator(4))
+            self.ax.tick_params(which="minor", length=2.5)
         self.ax.set_xlim(0, 24)
 
     def set_yax(
@@ -910,6 +916,8 @@ def _get_max_gap_in_minutes(figure_data: FigureData) -> float:
         "mwr-multi": 35,
         "weather-station": 12,
         "doppler-lidar-wind": 75,
+        "doppler-lidar": 75,
+        "radar": 5,
     }
     return max_allowed_gap.get(file_type, 10)
 

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

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

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

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