cloudnetpy 1.85.5__py3-none-any.whl → 1.86.1__py3-none-any.whl

cloudnetpy/categorize/containers.py

@@ -92,10 +92,18 @@ class ClassData:
         rain_from_disdrometer = self._find_rain_from_disdrometer()
         ind = ~rain_from_disdrometer.mask
         is_rain[ind] = rain_from_disdrometer[ind]
-        mask = ma.getmaskarray(self.tw)
-        first_valid_ind = np.nonzero(np.all(~mask, axis=0))[0][0]
-        is_positive_temp = self.tw[:, first_valid_ind] > T0 + 5  # Filter snowfall
-        return is_rain & is_positive_temp
+        # Filter out snowfall:
+        if (
+            self.data.disdrometer is not None
+            and "synop_WaWa" in self.data.disdrometer.data
+        ):
+            wawa = self.data.disdrometer.data["synop_WaWa"].data
+            liquid_rain_only = (wawa >= 57) & (wawa <= 68)
+        else:
+            mask = ma.getmaskarray(self.tw)
+            first_valid_ind = np.nonzero(np.all(~mask, axis=0))[0][0]
+            liquid_rain_only = self.tw[:, first_valid_ind] > T0 + 5
+        return is_rain & liquid_rain_only
 
     def _find_rain_from_radar_echo(self) -> npt.NDArray:
         first_gate_with_data = np.argmin(self.z.mask.all(axis=0))

cloudnetpy/categorize/disdrometer.py

@@ -24,14 +24,17 @@ class Disdrometer(DataSource):
 
     def interpolate_to_grid(self, time_grid: npt.NDArray) -> None:
         for key, array in self.data.items():
+            method = "nearest" if key == "synop_WaWa" else "linear"
             self.data[key].data = interpolate_1d(
-                self.time, array.data, time_grid, max_time=1
+                self.time, array.data, time_grid, max_time=1, method=method
             )
 
     def _init_rainfall_rate(self) -> None:
-        keys = ("rainfall_rate", "n_particles")
+        keys = ("rainfall_rate", "n_particles", "synop_WaWa")
         for key in keys:
             if key not in self.dataset.variables:
+                if key == "synop_WaWa":
+                    continue
                 msg = f"variable {key} is missing"
                 raise DisdrometerDataError(msg)
             self.append_data(self.dataset.variables[key][:], key)

cloudnetpy/categorize/radar.py

@@ -401,7 +401,6 @@ class Radar(DataSource):
 
         if fvel.shape == (n_time, n_height):
             # Folding velocity is already expanded in radar file
-            # Not yet in current files
             return fvel
         if utils.isscalar(fvel):
             # e.g. MIRA

cloudnetpy/instruments/bowtie.py

@@ -127,4 +127,9 @@ ATTRIBUTES = RPG_ATTRIBUTES | {
         ),
         dimensions=("time", "range"),
     ),
+    "nyquist_velocity": MetaData(
+        long_name="Nyquist velocity",
+        units="m s-1",
+        dimensions=("time", "chirp_sequence"),
+    ),
 }

cloudnetpy/instruments/fd12p.py

@@ -198,7 +198,4 @@ ATTRIBUTES = {
         comment="Cumulated snow since 00:00 UTC",
         dimensions=("time",),
     ),
-    "synop_WaWa": MetaData(
-        long_name="Synop code WaWa", units="1", dimensions=("time",)
-    ),
 }

cloudnetpy/instruments/rpg.py

@@ -1,9 +1,7 @@
-"""This module contains RPG Cloud Radar related functions."""
-
 import datetime
 import logging
 import math
-from collections.abc import Callable, Sequence
+from collections.abc import Sequence
 from os import PathLike
 from uuid import UUID
 
@@ -15,7 +13,7 @@ from rpgpy import RPGFileError
 from cloudnetpy import output, utils
 from cloudnetpy.cloudnetarray import CloudnetArray
 from cloudnetpy.constants import G_TO_KG, HPA_TO_PA, KM_H_TO_M_S, MM_H_TO_M_S
-from cloudnetpy.exceptions import InconsistentDataError, ValidTimeStampError
+from cloudnetpy.exceptions import ValidTimeStampError
 from cloudnetpy.instruments import instruments
 from cloudnetpy.instruments.cloudnet_instrument import CloudnetInstrument
 from cloudnetpy.instruments.instruments import Instrument
@@ -68,9 +66,8 @@ def rpg2nc(
     l1_files = utils.get_sorted_filenames(path_to_l1_files, ".LV1")
     fmcw94_objects, valid_files = _get_fmcw94_objects(l1_files, date)
     one_day_of_data = create_one_day_data_record(fmcw94_objects)
-    if not valid_files:
-        return uuid, []
-    print_info(one_day_of_data)
+    one_day_of_data["nyquist_velocity"] = _expand_nyquist(one_day_of_data)
+    _print_info(one_day_of_data)
     fmcw = Fmcw(one_day_of_data, site_meta)
     fmcw.convert_time_to_fraction_hour()
     fmcw.mask_invalid_ldr()
@@ -93,7 +90,7 @@ def rpg2nc(
     return uuid, valid_files
 
 
-def print_info(data: dict) -> None:
+def _print_info(data: dict) -> None:
     dual_pol = data["dual_polarization"]
     if dual_pol == 0:
         mode = "single polarisation"
@@ -108,57 +105,66 @@ RpgObjects = Sequence[Fmcw94Bin] | Sequence[HatproBinCombined]
 
 
 def create_one_day_data_record(rpg_objects: RpgObjects) -> dict:
-    """Concatenates all RPG data from one day."""
+    """Concatenates all RPG FMCW / HATPRO data from one day."""
     rpg_raw_data, rpg_header = _stack_rpg_data(rpg_objects)
-    if len(rpg_objects) > 1:
-        rpg_header = _reduce_header(rpg_header)
+    if "range" in rpg_header:
+        rpg_header["range"] = rpg_objects[0].header["range"]
+    should_be_constant = [
+        "model_number",
+        "dual_polarization",
+        "antenna_separation",
+        "antenna_diameter",
+        "antenna_gain",
+        "half_power_beam_width",
+        "radar_frequency",
+    ]
+    for key in should_be_constant:
+        if key not in rpg_header:
+            continue
+        unique_values = np.unique(rpg_header[key])
+        if len(unique_values) > 1:
+            msg = f"More than one value for {key} found: {unique_values}"
+            raise ValueError(msg)
+        rpg_header[key] = unique_values[0]
+
     rpg_raw_data = _mask_invalid_data(rpg_raw_data)
     return {**rpg_header, **rpg_raw_data}
 
 
-def _stack_rpg_data(rpg_objects: RpgObjects) -> tuple[dict, dict]:
-    """Combines data from hourly RPG objects.
-
-    Notes:
-        Ignores variable names starting with an underscore.
+def _expand_nyquist(data: dict) -> npt.NDArray:
+    """Expands Nyquist velocity from time X chirp => time X range."""
+    nyquist_velocity = ma.array(data["nyquist_velocity"])
+    chirp_start_indices = ma.array(data["chirp_start_indices"])
+    n_time = chirp_start_indices.shape[0]
+    n_range = len(data["range"])
+    expanded_nyquist = np.empty((n_time, n_range))
+    for t in range(n_time):
+        starts = chirp_start_indices[t].compressed()
+        v_nyq = nyquist_velocity[t].compressed()
+        ends = np.r_[starts[1:], n_range]
+        seg_lengths = ends - starts
+        expanded_nyquist[t, :] = np.repeat(v_nyq, seg_lengths)
+    return expanded_nyquist
 
-    """
-
-    def _stack(source: dict, target: dict, fun: Callable) -> None:
-        for name, value in source.items():
-            if not name.startswith("_"):
-                target[name] = fun((target[name], value)) if name in target else value
 
+def _stack_rpg_data(rpg_objects: RpgObjects) -> tuple[dict, dict]:
     data: dict = {}
     header: dict = {}
     for rpg in rpg_objects:
-        _stack(rpg.data, data, ma.concatenate)
-        _stack(rpg.header, header, ma.vstack)
+        for src, dst in ((rpg.data, data), (rpg.header, header)):
+            for name, value in src.items():
+                if name.startswith("_"):
+                    continue
+                arr = dst.get(name)
+                fun = (
+                    ma.concatenate
+                    if any(isinstance(x, ma.MaskedArray) for x in (value, arr))
+                    else np.concatenate
+                )
+                dst[name] = fun((arr, value)) if arr is not None else value
     return data, header
 
 
-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 ("sample_duration", "calibration_interval", "noise_threshold"):
-                logging.warning(msg)
-            else:
-                raise InconsistentDataError(msg)
-        reduced_header[key] = first_profile_value
-    return reduced_header
-
-
 def _mask_invalid_data(data_in: dict) -> dict:
     """Masks zeros and other fill values from data."""
     data = data_in.copy()
@@ -193,33 +199,68 @@ def _get_fmcw94_objects(
             continue
         objects.append(obj)
         valid_files.append(file)
-    if objects:
-        objects, valid_files = _remove_files_with_bad_height(objects, valid_files)
-    if not valid_files:
-        raise ValidTimeStampError
+    if not objects:
+        msg = "No valid files found"
+        raise ValidTimeStampError(msg)
+    objects = _interpolate_to_common_height(objects)
+    objects = _pad_chirp_related_fields(objects)
+    objects = _expand_time_related_fields(objects)
     return objects, valid_files
 
 
-def _remove_files_with_bad_height(objects: list, files: list) -> tuple[list, list]:
-    lengths = [obj.data["Zh"].shape[1] for obj in objects]
-    most_common = np.bincount(lengths).argmax()
-    files = [
-        file
-        for file, obj, length in zip(files, objects, lengths, strict=True)
-        if length == most_common
-    ]
-    objects = [
-        obj
-        for obj, length in zip(objects, lengths, strict=True)
-        if length == most_common
-    ]
-    n_removed = len(lengths) - len(files)
-    if n_removed > 0:
-        logging.warning(
-            "Removed %s RPG-FMCW-94 files due to inconsistent height vector",
-            n_removed,
-        )
-    return objects, files
+def _interpolate_to_common_height(objects: list[Fmcw94Bin]) -> list[Fmcw94Bin]:
+    range_arrays = [obj.header["range"] for obj in objects]
+    if all(np.array_equal(range_arrays[0], r) for r in range_arrays[1:]):
+        return objects
+    # Use range with the highest range gate for interpolation
+    target_height = max(range_arrays, key=lambda r: r[-1])
+    for obj in objects:
+        src_range = obj.header["range"]
+        if np.array_equal(src_range, target_height):
+            continue
+        for key, arr in obj.data.items():
+            if arr.ndim == 2 and arr.shape[1] == src_range.size:
+                obj.data[key] = utils.interpolate_2D_along_y(
+                    src_range, arr, target_height
+                )
+        obj.header["range"] = target_height
+    return objects
+
+
+def _pad_chirp_related_fields(objects: list[Fmcw94Bin]) -> list[Fmcw94Bin]:
+    """Pads chirp-related header fields with masked values to have the same length."""
+    chirp_lens = [len(obj.header["chirp_start_indices"]) for obj in objects]
+    if all(chirp_lens[0] == length for length in chirp_lens[1:]):
+        return objects
+    max_chirp_len = max(chirp_lens)
+    for obj in objects:
+        n_chirps = len(obj.header["chirp_start_indices"])
+        if n_chirps == max_chirp_len:
+            continue
+        for key, arr in obj.header.items():
+            if not isinstance(arr, str) and arr.ndim == 1 and arr.size == n_chirps:
+                pad_len = max_chirp_len - n_chirps
+                masked_arr = ma.array(arr, dtype=arr.dtype)
+                pad = ma.masked_all(pad_len, dtype=arr.dtype)
+                obj.header[key] = ma.concatenate([masked_arr, pad])
+    return objects
+
+
+def _expand_time_related_fields(objects: list[Fmcw94Bin]) -> list[Fmcw94Bin]:
+    for obj in objects:
+        n_time = obj.data["time"].size
+        for key in obj.header:
+            if key in ("range", "time") or key.startswith("_"):
+                continue
+            arr = obj.header[key]
+            # Handle outliers in latitude and longitude (e.g. Galati 2024-02-11):
+            if key in ("latitude", "longitude"):
+                arr = ma.median(arr)
+            if utils.isscalar(arr):
+                obj.header[key] = np.repeat(arr, n_time)
+            else:
+                obj.header[key] = np.tile(arr, (n_time, 1))
+    return objects
 
 
 def _validate_date(obj: Fmcw94Bin, expected_date: datetime.date) -> None:
@@ -443,9 +484,11 @@ RPG_ATTRIBUTES = {
         long_name="File code",
         units="1",
         comment="Indicates the RPG software version.",
-        dimensions=None,
+        dimensions=("time",),
+    ),
+    "program_number": MetaData(
+        long_name="Program number", units="1", dimensions=("time",)
     ),
-    "program_number": MetaData(long_name="Program number", units="1", dimensions=None),
     "model_number": MetaData(
         long_name="Model number",
         units="1",
@@ -469,54 +512,51 @@ RPG_ATTRIBUTES = {
         dimensions=None,
     ),
     "sample_duration": MetaData(
-        long_name="Sample duration", units="s", dimensions=None
+        long_name="Sample duration", units="s", dimensions=("time",)
     ),
     "calibration_interval": MetaData(
-        long_name="Calibration interval in samples", units="1", dimensions=None
+        long_name="Calibration interval in samples", units="1", dimensions=("time",)
     ),
     "number_of_spectral_samples": MetaData(
         long_name="Number of spectral samples in each chirp sequence",
         units="1",
-        dimensions=("chirp_sequence",),
-    ),
-    "nyquist_velocity": MetaData(
-        long_name="Nyquist velocity", units="m s-1", dimensions=("chirp_sequence",)
+        dimensions=("time", "chirp_sequence"),
     ),
     "number_of_averaged_chirps": MetaData(
         long_name="Number of averaged chirps in sequence",
         units="1",
-        dimensions=("chirp_sequence",),
+        dimensions=("time", "chirp_sequence"),
     ),
     "chirp_start_indices": MetaData(
         long_name="Chirp sequences start indices",
         units="1",
-        dimensions=("chirp_sequence",),
+        dimensions=("time", "chirp_sequence"),
     ),
     "integration_time": MetaData(
         long_name="Integration time",
         units="s",
         comment="Effective integration time of chirp sequence",
-        dimensions=("chirp_sequence",),
+        dimensions=("time", "chirp_sequence"),
     ),
     "range_resolution": MetaData(
         long_name="Vertical resolution of range",
         units="m",
-        dimensions=("chirp_sequence",),
+        dimensions=("time", "chirp_sequence"),
     ),
     "FFT_window": MetaData(
         long_name="FFT window type",
         units="1",
         definition=DEFINITIONS["FFT_window"],
-        dimensions=None,
+        dimensions=("time",),
     ),
     "input_voltage_range": MetaData(
-        long_name="ADC input voltage range (+/-)", units="mV", dimensions=None
+        long_name="ADC input voltage range (+/-)", units="mV", dimensions=("time",)
     ),
     "noise_threshold": MetaData(
         long_name="Noise filter threshold factor",
         units="1",
         comment="Multiple of the standard deviation of Doppler spectra.",
-        dimensions=None,
+        dimensions=("time",),
     ),
     "time_ms": MetaData(long_name="Time ms", units="ms", dimensions=("time",)),
     "quality_flag": MetaData(
@@ -545,7 +585,4 @@ RPG_ATTRIBUTES = {
     "pc_temperature": MetaData(
         long_name="PC temperature", units="K", dimensions=("time",)
     ),
-    "correlation_coefficient": MetaData(
-        long_name="Correlation coefficient", units="1", dimensions=None
-    ),
 }

cloudnetpy/metadata.py

@@ -224,4 +224,7 @@ COMMON_ATTRIBUTES = {
         comment="SNR threshold used in data screening.",
         dimensions=None,
     ),
+    "synop_WaWa": MetaData(
+        long_name="Synop code WaWa", units="1", dimensions=("time",)
+    ),
 }

cloudnetpy/output.py

@@ -391,7 +391,6 @@ def _write_vars2nc(nc: netCDF4.Dataset, cloudnet_variables: dict) -> None:
         else:
             fill_value = False
         size = obj.dimensions if obj.dimensions is not None else ()
-
         nc_variable = nc.createVariable(
             obj.name,
             obj.data_type,
@@ -399,7 +398,11 @@ def _write_vars2nc(nc: netCDF4.Dataset, cloudnet_variables: dict) -> None:
             zlib=True,
             fill_value=fill_value,
         )
-        nc_variable[:] = obj.data
+        try:
+            nc_variable[:] = obj.data
+        except IndexError as err:
+            msg = f"Unable to write variable {obj.name} to file: {err}"
+            raise IndexError(msg) from err
         for attr in obj.fetch_attributes():
             setattr(nc_variable, attr, getattr(obj, attr))
 

cloudnetpy/plotting/plot_meta.py

@@ -564,7 +564,7 @@ ATTRIBUTES = {
             plot_range=(0, 5),
         ),
         "radar_rain_atten": PlotMeta(
-            plot_range=(0, 5),
+            plot_range=(0, 15),
         ),
         "radar_melting_atten": PlotMeta(
             plot_range=(0, 5),

cloudnetpy/products/classification.py

@@ -11,6 +11,7 @@ from numpy import ma
 
 from cloudnetpy import output, utils
 from cloudnetpy.categorize import atmos_utils
+from cloudnetpy.constants import M_S_TO_MM_H
 from cloudnetpy.datasource import DataSource
 from cloudnetpy.metadata import MetaData
 from cloudnetpy.products.product_tools import CategorizeBits, QualityBits
@@ -189,24 +190,36 @@ def _get_attenuation_classes(data_source: DataSource) -> AttenuationClass:
     rain_atten = _read_atten("radar_rain_atten")
     melting_atten = _read_atten("radar_melting_atten")
 
-    if (
-        "lwp" not in data_source.dataset.variables
-        or data_source.getvar("radar_frequency") < 90
-    ):
+    not_w_band = data_source.getvar("radar_frequency") < 90
+
+    if "lwp" not in data_source.dataset.variables or not_w_band:
         lwp = np.zeros(data_source.time.shape)
     else:
         lwp_data = data_source.getvar("lwp")
         lwp = lwp_data.filled(0) if isinstance(lwp_data, ma.MaskedArray) else lwp_data
 
+    if "rainfall_rate" not in data_source.dataset.variables or not_w_band:
+        rain_rate = np.zeros(data_source.time.shape)
+    else:
+        rain_data = data_source.getvar("rainfall_rate") * M_S_TO_MM_H
+        rain_rate = (
+            rain_data.filled(0) if isinstance(rain_data, ma.MaskedArray) else rain_data
+        )
+
     total_atten = liquid_atten + rain_atten + melting_atten
 
     threshold_moderate = 10  # dB
     threshold_severe = 15  # dB
     threshold_lwp = 1  # kg/m2
+    threshold_rain = 3  # mm/h
 
     small = total_atten > 0
     moderate = total_atten >= threshold_moderate
-    severe = (total_atten > threshold_severe) | (lwp[:, np.newaxis] > threshold_lwp)
+    severe = (
+        (total_atten > threshold_severe)
+        | (lwp[:, np.newaxis] > threshold_lwp)
+        | (rain_rate[:, np.newaxis] > threshold_rain)
+    )
 
     return AttenuationClass(small=small, moderate=moderate, severe=severe)
 

cloudnetpy/utils.py

@@ -387,7 +387,7 @@ def interpolate_2d_mask(
 def interpolate_2d_nearest(
     x: npt.NDArray,
     y: npt.NDArray,
-    z: npt.NDArray,
+    z: ma.MaskedArray,
     x_new: npt.NDArray,
     y_new: npt.NDArray,
 ) -> ma.MaskedArray:
@@ -419,17 +419,57 @@ def interpolate_2d_nearest(
     return fun((xx, yy)).T
 
 
+def interpolate_2D_along_y(
+    y: npt.NDArray,
+    z: npt.NDArray | ma.MaskedArray,
+    y_new: npt.NDArray,
+) -> ma.MaskedArray:
+    """Fast 1D nearest-neighbor interpolation along y for each x.
+
+    Args:
+        y: 1D numpy array of y-coordinates (length M).
+        z: 2D array of shape (N, M).
+        y_new: 1D numpy array of new y-coordinates.
+
+    Returns:
+        Masked 2D masked array interpolated along y.
+
+    Notes:
+        Only interpolates along y. Points outside range are masked.
+    """
+    idx = np.searchsorted(y, y_new, side="left")
+    idx = np.clip(idx, 0, len(y) - 1)
+    left = np.maximum(idx - 1, 0)
+    choose_right = (idx == 0) | (
+        (idx < len(y)) & (np.abs(y[idx] - y_new) < np.abs(y_new - y[left]))
+    )
+    idx[~choose_right] = left[~choose_right]
+    z_interp = ma.array(z[:, idx])
+    mask = (y_new < y.min()) | (y_new > y.max())
+    if z_interp.mask is ma.nomask:
+        z_mask = np.zeros(z_interp.shape, dtype=bool)
+    else:
+        z_mask = z_interp.mask.copy()
+    z_mask[:, mask] = True
+    return ma.MaskedArray(z_interp, mask=z_mask)
+
+
 def interpolate_1d(
-    time: npt.NDArray, y: ma.MaskedArray, time_new: npt.NDArray, max_time: float
-) -> npt.NDArray:
+    time: npt.NDArray,
+    y: ma.MaskedArray,
+    time_new: npt.NDArray,
+    max_time: float,
+    method: str = "linear",
+) -> ma.MaskedArray:
     """1D linear interpolation preserving the mask.
 
     Args:
         time: 1D array in fraction hour.
-        y: 1D masked array, data values.
+        y: 1D array, data values.
         time_new: 1D array, new time coordinates.
         max_time: Maximum allowed gap in minutes. Values outside this gap will
             be masked.
+        method: Interpolation method, 'linear' (default) or 'nearest'.
     """
     if np.max(time) > 24 or np.min(time) < 0:
         msg = "Time vector must be in fraction hours between 0 and 24"
@@ -439,7 +479,7 @@ def interpolate_1d(
             return ma.masked_all(time_new.shape)
         time = time[~y.mask]
         y = y[~y.mask]
-    fun = interp1d(time, y, fill_value=(y[0], y[-1]), bounds_error=False)
+    fun = interp1d(time, y, kind=method, fill_value=(y[0], y[-1]), bounds_error=False)
     interpolated = ma.array(fun(time_new))
     bad_idx = get_gap_ind(time, time_new, max_time / 60)
 

cloudnetpy/version.py

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

{cloudnetpy-1.85.5.dist-info → cloudnetpy-1.86.1.dist-info}/METADATA

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

{cloudnetpy-1.85.5.dist-info → cloudnetpy-1.86.1.dist-info}/RECORD

@@ -5,18 +5,18 @@ cloudnetpy/concat_lib.py,sha256=u4UOjYzLnThaq-89iwA837OcOJpfmj_3RmRwBCFUh74,1321
 cloudnetpy/constants.py,sha256=YnoSzZm35NDooJfhlulSJBc7g0eSchT3yGytRaTaJEI,845
 cloudnetpy/datasource.py,sha256=EMJ4UHD8Z-JJ9Q82S7RgU1I2q4Z0RcBzBMKUAIwnZBI,6356
 cloudnetpy/exceptions.py,sha256=ZB3aUwjVRznR0CcZ5sZHrB0yz13URDf52Ksv7G7C7EA,1817
-cloudnetpy/metadata.py,sha256=CFpXmdEkVPzvLPv2xHIR-aMMQ-TR26KfESYw-98j7sk,7213
-cloudnetpy/output.py,sha256=0bybnILsgKHWIuw2GYkqTz2iMCJDZLUN25IQ9o_v3Cg,14968
+cloudnetpy/metadata.py,sha256=qqLXSUVhIi6H_Z1zlfXeEkc9qBqA1XO-UwsCVxjIcfk,7317
+cloudnetpy/output.py,sha256=NpuxmtLk2JLTdubCupdL4rFbx30YqgeU6vwAe_v8thc,15133
 cloudnetpy/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-cloudnetpy/utils.py,sha256=x3hukn7wuOyF3ctG8EOhvLbSUHOVKLW9ez2LzTjeolU,33110
-cloudnetpy/version.py,sha256=LIie3vB4_JAVH9ThGGPAWNjAwLy9Nrli2ODRwPbsFfY,72
+cloudnetpy/utils.py,sha256=9MMYok1A7oJS2XBkIFhaiDWt_MYLo8tLzF3_880masM,34321
+cloudnetpy/version.py,sha256=QtEAE1savq9vVLdajaUQbntR7f8B0XJ7TQqre-v5kTA,72
 cloudnetpy/categorize/__init__.py,sha256=gtvzWr0IDRn2oA6yHBvinEhTGTuub-JkrOv93lBsgrE,61
 cloudnetpy/categorize/atmos_utils.py,sha256=uWc9TABVYPI0sn4H5Az9Jf6NVRaWyEKIi17f0pAJQxE,10679
 cloudnetpy/categorize/attenuation.py,sha256=Y_-fzmQTltWTqIZTulJhovC7a6ifpMcaAazDJcnMIOc,990
 cloudnetpy/categorize/categorize.py,sha256=UZph73hVCU-UIxB2837oyLFKvTk8fdgLC-krgJDg8zk,23227
 cloudnetpy/categorize/classify.py,sha256=skA9K6Bxh9mFZ_fM4d78zt09BPDzfHLttXle6mFCMbw,8553
-cloudnetpy/categorize/containers.py,sha256=PIJwgQos3CxF9BG4hBNLTaZq252FTH0kdgagT31mFmc,5517
-cloudnetpy/categorize/disdrometer.py,sha256=_Hc7EhoxIIsAn7HCrgSB6aPMelaIW6c4lLeMf0NOnVI,1117
+cloudnetpy/categorize/containers.py,sha256=oXvDazVtd9u3PwNP8ssKMfjW8s6O1bKNcTYHif6AHko,5808
+cloudnetpy/categorize/disdrometer.py,sha256=gJW6qD2Cn3w0JbH69wqtlxDFUhKb2ioR-e9f7VZwxvc,1283
 cloudnetpy/categorize/droplet.py,sha256=wnMN9rHNSMZLXNXuYEd-RAS_8eAIIo2vkE7pp3DSTKs,8725
 cloudnetpy/categorize/falling.py,sha256=ZscFGFPFz_Nlc7PwjVwFDSVOzMGOuCDVAFj7pLvYctQ,4475
 cloudnetpy/categorize/freezing.py,sha256=gigqpb4qfeQSlKXkrPUwCbMnMsxl74thJWSRW2iHJOg,3796
@@ -26,7 +26,7 @@ cloudnetpy/categorize/lidar.py,sha256=CkIbsMFHROYnhIc1e1ci61-POqeL8AOf6RCZicpNBV
 cloudnetpy/categorize/melting.py,sha256=vhc6zq3L4gp7oEPHMnQlT2m6YBE5-CS5gdTNA7gVHRg,6329
 cloudnetpy/categorize/model.py,sha256=iFakrXy3npbg4qUrpUGEBEdwBnmlWsMgogPCtfWl7sw,6805
 cloudnetpy/categorize/mwr.py,sha256=QKc6f_h0mCVjzSa4BrwTtVAPx0Tn1VJ4jqGmIBke5II,1710
-cloudnetpy/categorize/radar.py,sha256=2mTDa9BLxQeaORm-YPQ1lJyjAKew6NYzjtUvjpIvBYU,16044
+cloudnetpy/categorize/radar.py,sha256=CHcCWm2Nl-ycBqmbNz8LPSTmuNzgmmNxaJO2uQ-o09c,16005
 cloudnetpy/categorize/attenuations/__init__.py,sha256=kIyQEZ6VVO6jJOAndrt7jNU15pm0Cavh5GnDjFmIG1M,1040
 cloudnetpy/categorize/attenuations/gas_attenuation.py,sha256=emr-RCxQT0i2N8k6eBNhRsmsCBPHJzQsWJfjC4fVSTo,975
 cloudnetpy/categorize/attenuations/liquid_attenuation.py,sha256=bmqmPk_93J4njE16-VQ1bPI7oNSS8m9ACuUH7IErBs8,3069
@@ -34,14 +34,14 @@ cloudnetpy/categorize/attenuations/melting_attenuation.py,sha256=zmpF7Gek4W9cF-5
 cloudnetpy/categorize/attenuations/rain_attenuation.py,sha256=wJPyCiKWzsQDzMhqbA7mYwj9YRVcJIpXWhBnEYFy3uU,2843
 cloudnetpy/instruments/__init__.py,sha256=PEgrrQNoiOuN_ctYilmt4LV2QCLg1likPjJdWtuGlLs,528
 cloudnetpy/instruments/basta.py,sha256=N-kRgl5Vm52pXzr9umo4YsA0hn4zZCOa-0_zZTzhheY,4284
-cloudnetpy/instruments/bowtie.py,sha256=WZYB_o90I5QKOupCRjzZU4Mi54oX_3teyAPBDUogwRI,4301
+cloudnetpy/instruments/bowtie.py,sha256=1wavsYkVMC81CflUD-nTjeyjmATGKYfEJ6aNcwOOapM,4450
 cloudnetpy/instruments/ceilo.py,sha256=scb3n0gLOZ1gPs_WdWw9KS5Ut_mtCvSdWFMKKCvXhBI,8979
 cloudnetpy/instruments/ceilometer.py,sha256=c37uteeuGnlE-o-Smu49H2qQJw6qZ0tc3Bzhyr1FoSo,13063
 cloudnetpy/instruments/cl61d.py,sha256=JsCHqVzFGhZi-5xcnsB507FDpyuw83uSWK3IFO3DhdI,2238
 cloudnetpy/instruments/cloudnet_instrument.py,sha256=B1UkiB0ytnT3MWYalEegql5QIPaMLg5bJy5xI50JEco,4503
 cloudnetpy/instruments/copernicus.py,sha256=ygEViERBSJdMeP9OxfLelZRDEbSRzY8n17ruYie2wm4,6970
 cloudnetpy/instruments/da10.py,sha256=FcXEu6dfAVvYzMvauiJBSyEkYhGa8hBuwrrdo-3LwHE,1960
-cloudnetpy/instruments/fd12p.py,sha256=5TFyNO26VGpO4ts9UIJiuLo4LUwQPHO6aK2fTnOtaKY,7019
+cloudnetpy/instruments/fd12p.py,sha256=M5yJDsPi0tHWMDUYOsoqciyKdRQ-zHEq2t19KnpZ6Io,6915
 cloudnetpy/instruments/galileo.py,sha256=f_-GkRxhNaQPbI8HpOwSmoKfGqyjmD16A0ZFgwLOIig,5137
 cloudnetpy/instruments/hatpro.py,sha256=TGOqwW0TfoPEYk13MFvFzwgJGzm6MVE5AsPavcIoj3I,10248
 cloudnetpy/instruments/instruments.py,sha256=WZgH7HjzM9Ane1CSnYCSLidbST8hunUeSt2lPntq9As,4999
@@ -53,7 +53,7 @@ cloudnetpy/instruments/nc_radar.py,sha256=XFKxPLKivnhHTgjE5HFrxjWZ0oCifhDUAog051
 cloudnetpy/instruments/pollyxt.py,sha256=IFq_RJrhgJ79OVyuo48PwYQK_zZ6VZFB_S5bEirRyzs,10566
 cloudnetpy/instruments/radiometrics.py,sha256=QKfnrZlQ0sFcFjmv1ShnCMTJQv64w4akjK-JAIY4gCg,16116
 cloudnetpy/instruments/rain_e_h3.py,sha256=fjv3SgeUNx9GisYqLrBnX9AjnO17VtouyoPh12VE9uo,5465
-cloudnetpy/instruments/rpg.py,sha256=d7qlfez1Pe326d4xkQrTk2xRxvs3g68LEGehErdoBoc,19122
+cloudnetpy/instruments/rpg.py,sha256=LmvViLSevzjR2HmkP6MCHHDzafpTWoCvw7GaaEC37s0,21029
 cloudnetpy/instruments/rpg_reader.py,sha256=EcsUUyKNiZ-kEEz48zibs-uXH7CMNpedqlrfWOY4Dks,11899
 cloudnetpy/instruments/toa5.py,sha256=CfmmBMv5iMGaWHIGBK01Rw24cuXC1R1RMNTXkmsm340,1760
 cloudnetpy/instruments/vaisala.py,sha256=tu7aljkMKep0uCWz-Sd-GuBXF_Yy421a4nHy0ffpMoc,4725
@@ -103,10 +103,10 @@ cloudnetpy/model_evaluation/tests/unit/test_plotting.py,sha256=5hkhtqX-JQ8-Yy6DA
 cloudnetpy/model_evaluation/tests/unit/test_statistical_methods.py,sha256=Ra3r4V0qbqkpDuaTYvEIbaasl0nZ5gmTLR4eGC0glBQ,9724
 cloudnetpy/model_evaluation/tests/unit/test_tools.py,sha256=Ia_VrLdV2NstX5gbx_3AZTOAlrgLAy_xFZ8fHYVX0xI,3817
 cloudnetpy/plotting/__init__.py,sha256=lg9Smn4BI0dVBgnDLC3JVJ4GmwoSnO-qoSd4ApvwV6Y,107
-cloudnetpy/plotting/plot_meta.py,sha256=TpUPmBsME3VhB-Atct6EG4ntzlo7Essre83ySp40Ulc,18454
+cloudnetpy/plotting/plot_meta.py,sha256=4AyKaTgp4ngxcUd3d0d7usJK-VpVrXGOCaXmpkGueNo,18455
 cloudnetpy/plotting/plotting.py,sha256=v1dD6kmMJ-t88p4rqrnFzS94qUjpI_v9fCAp49fA2V0,42015
 cloudnetpy/products/__init__.py,sha256=cBJdJBYltz5ZTKDqnRo-0StytAZK8gE3RYxxriFA4ak,295
-cloudnetpy/products/classification.py,sha256=6WxiGGJXqlAPgJ9hVNTKm3f7iGsqA7e3L3xqVxFFs7w,16894
+cloudnetpy/products/classification.py,sha256=Vbx2A43XP8F6q82wOPA_7pyUhuoyDty199bLU_7aNd0,17382
 cloudnetpy/products/der.py,sha256=UXdAxmmwChVVWSI4QSGAXphfMnbymGRTtGdKWEvh-J4,13162
 cloudnetpy/products/drizzle.py,sha256=0h1N_WVjC2GgIkAN-4ydOwl7WJn3psxeqmPHfX8WHhQ,11935
 cloudnetpy/products/drizzle_error.py,sha256=QN98Io9UsBoEYxKBqfwoS88OGBiK5U5RYnVQjyTWHCI,6220
@@ -118,10 +118,10 @@ cloudnetpy/products/lwc.py,sha256=xsNiiG6dGKIkWaFk0xWTabc1bZ4ULf6SqcqHs7itAUk,19
 cloudnetpy/products/mie_lu_tables.nc,sha256=It4fYpqJXlqOgL8jeZ-PxGzP08PMrELIDVe55y9ob58,16637951
 cloudnetpy/products/mwr_tools.py,sha256=MMWnp68U7bv157-CPB2VeTQvaR6zl7sexbBT_kJ_pn8,6734
 cloudnetpy/products/product_tools.py,sha256=eyqIw_0KhlpmmYQE69RpGdRIAOW7JVPlEgkTBp2kdps,11302
-cloudnetpy-1.85.5.dist-info/licenses/LICENSE,sha256=wcZF72bdaoG9XugpyE95Juo7lBQOwLuTKBOhhtANZMM,1094
+cloudnetpy-1.86.1.dist-info/licenses/LICENSE,sha256=wcZF72bdaoG9XugpyE95Juo7lBQOwLuTKBOhhtANZMM,1094
 docs/source/conf.py,sha256=IKiFWw6xhUd8NrCg0q7l596Ck1d61XWeVjIFHVSG9Og,1490
-cloudnetpy-1.85.5.dist-info/METADATA,sha256=dR72eUK34l0HmYW2R-Sn-8pqkt53ZslGXuMTwyJGKEc,5836
-cloudnetpy-1.85.5.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-cloudnetpy-1.85.5.dist-info/entry_points.txt,sha256=HhY7LwCFk4qFgDlXx_Fy983ZTd831WlhtdPIzV-Y3dY,51
-cloudnetpy-1.85.5.dist-info/top_level.txt,sha256=ibSPWRr6ojS1i11rtBFz2_gkIe68mggj7aeswYfaOo0,16
-cloudnetpy-1.85.5.dist-info/RECORD,,
+cloudnetpy-1.86.1.dist-info/METADATA,sha256=jSsYS8WOmfToTJzLKqwqAUh3PCJv6BPRv6lUAryjRtM,5836
+cloudnetpy-1.86.1.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+cloudnetpy-1.86.1.dist-info/entry_points.txt,sha256=HhY7LwCFk4qFgDlXx_Fy983ZTd831WlhtdPIzV-Y3dY,51
+cloudnetpy-1.86.1.dist-info/top_level.txt,sha256=ibSPWRr6ojS1i11rtBFz2_gkIe68mggj7aeswYfaOo0,16
+cloudnetpy-1.86.1.dist-info/RECORD,,