cloudnetpy 1.55.21__py3-none-any.whl → 1.55.22__py3-none-any.whl

cloudnetpy/categorize/atmos_utils.py

@@ -140,4 +140,4 @@ def k2c(temp: np.ndarray) -> np.ndarray:
 
 def mmh2ms(data: np.ndarray) -> np.ndarray:
     """Converts mm h-1 to m s-1"""
-    return data / 3600_000
+    return data / con.SEC_IN_HOUR * con.MM_TO_M

cloudnetpy/categorize/categorize.py

@@ -68,8 +68,8 @@ def generate_categorize(
         model_gap_ind = data["model"].interpolate_to_grid(time, height)
         radar_gap_ind = data["radar"].rebin_to_grid(time)
         lidar_gap_ind = data["lidar"].interpolate_to_grid(time, height)
-        gap_indices = radar_gap_ind + lidar_gap_ind + model_gap_ind
-        return list(set(range(len(time))) - set(gap_indices))
+        gap_indices = set(radar_gap_ind + lidar_gap_ind + model_gap_ind)
+        return [ind for ind in range(len(time)) if ind not in gap_indices]
 
     def _screen_bad_time_indices(valid_indices: list) -> None:
         n_time_full = len(time)

cloudnetpy/categorize/lidar.py

@@ -1,5 +1,6 @@
 """Lidar module, containing the :class:`Lidar` class."""
 import logging
+from typing import Literal
 
 import numpy as np
 from numpy import ma
@@ -22,46 +23,46 @@ class Lidar(DataSource):
         self.append_data(self.getvar("beta"), "beta")
         self._add_meta()
 
-    def interpolate_to_grid(self, time_new: np.ndarray, height_new: np.ndarray) -> list:
+    def interpolate_to_grid(
+        self, time_new: np.ndarray, height_new: np.ndarray
+    ) -> list[int]:
         """Interpolate beta using nearest neighbor."""
-        max_height = 100.0  # m
-        max_time = 1.0  # min
+        max_height = 100  # m
+        max_time = 1 / 60  # min -> fraction hour
 
-        # Remove completely masked profiles from the interpolation
-        beta = self.data["beta"][:]
-        indices = []
-        for ind, b in enumerate(beta):
-            if ma.all(b) is not ma.masked:
-                indices.append(ind)
         if self.height is None:
             msg = "Unable to interpolate lidar: no height information"
             raise RuntimeError(msg)
-        beta_interpolated = interpolate_2d_nearest(
+
+        # Interpolate beta to new grid but ignore profiles that are completely masked
+        beta = self.data["beta"][:]
+        indices = [ind for ind, b in enumerate(beta) if ma.all(b) is not ma.masked]
+        beta_interp = interpolate_2d_nearest(
             self.time[indices],
             self.height,
             beta[indices, :],
             time_new,
             height_new,
         )
+        # Mask data points that are too far from the original grid
+        time_gap_ind = _get_gap_ind(self.time[indices], time_new, max_time)
+        height_gap_ind = _get_gap_ind(self.height, height_new, max_height)
+        self._mask_profiles(beta_interp, time_gap_ind, "time")
+        self._mask_profiles(beta_interp, height_gap_ind, "height")
+        self.data["beta"].data = beta_interp
+        return time_gap_ind
 
-        # Filter profiles and range gates having data gap
-        max_time /= 60  # to fraction hour
-        bad_time_indices = _get_bad_indices(self.time[indices], time_new, max_time)
-        bad_height_indices = _get_bad_indices(self.height, height_new, max_height)
-        if bad_time_indices:
-            logging.warning(
-                "Unable to interpolate lidar for %s time steps",
-                len(bad_time_indices),
-            )
-        beta_interpolated[bad_time_indices, :] = ma.masked
-        if bad_height_indices:
-            logging.warning(
-                "Unable to interpolate lidar for %s altitudes",
-                len(bad_height_indices),
-            )
-        beta_interpolated[:, bad_height_indices] = ma.masked
-        self.data["beta"].data = beta_interpolated
-        return bad_time_indices
+    @staticmethod
+    def _mask_profiles(
+        data: ma.MaskedArray, ind: list[int], dim: Literal["time", "height"]
+    ) -> None:
+        prefix = f"Unable to interpolate lidar for {len(ind)}"
+        if dim == "time" and ind:
+            logging.warning("%s time steps", prefix)
+            data[ind, :] = ma.masked
+        elif dim == "height" and ind:
+            logging.warning("%s altitudes", prefix)
+            data[:, ind] = ma.masked
 
     def _add_meta(self) -> None:
         self.append_data(float(self.getvar("wavelength")), "lidar_wavelength")
@@ -69,15 +70,9 @@ class Lidar(DataSource):
         self.append_data(3.0, "beta_bias")
 
 
-def _get_bad_indices(
-    original_grid: np.ndarray,
-    new_grid: np.ndarray,
-    threshold: float,
-) -> list:
-    indices = []
-    for ind, value in enumerate(new_grid):
-        diffu = np.abs(original_grid - value)
-        distance = diffu[diffu.argmin()]
-        if distance > threshold:
-            indices.append(ind)
-    return indices
+def _get_gap_ind(grid: np.ndarray, new_grid: np.ndarray, threshold: float) -> list[int]:
+    return [
+        ind
+        for ind, value in enumerate(new_grid)
+        if np.min(np.abs(grid - value)) > threshold
+    ]

cloudnetpy/categorize/mwr.py

@@ -2,6 +2,7 @@
 import numpy as np
 
 from cloudnetpy import utils
+from cloudnetpy.constants import G_TO_KG
 from cloudnetpy.datasource import DataSource
 
 
@@ -36,8 +37,7 @@ class Mwr(DataSource):
 
     def _init_lwp_error(self) -> None:
         random_error, bias = 0.25, 20
-        g2kg = 1e-3
-        lwp_error = utils.l2norm(self.data["lwp"][:] * random_error, bias * g2kg)
+        lwp_error = utils.l2norm(self.data["lwp"][:] * random_error, bias * G_TO_KG)
         self.append_data(lwp_error, "lwp_error", units="kg m-2")
         self.data["lwp_error"].comment = (
             "This variable is a rough estimate of the one-standard-deviation\n"

cloudnetpy/categorize/radar.py

@@ -8,6 +8,7 @@ from scipy import constants
 
 from cloudnetpy import utils
 from cloudnetpy.categorize.classify import ClassificationResult
+from cloudnetpy.constants import SEC_IN_HOUR
 from cloudnetpy.datasource import DataSource
 
 
@@ -261,8 +262,7 @@ class Radar(DataSource):
             return z_error
 
         def _number_of_independent_pulses() -> float:
-            seconds_in_hour = 3600
-            dwell_time = utils.mdiff(self.time) * seconds_in_hour
+            dwell_time = utils.mdiff(self.time) * SEC_IN_HOUR
             return (
                 dwell_time
                 * self.radar_frequency

cloudnetpy/constants.py

@@ -18,3 +18,10 @@ RS: Final = 287.058
 
 # ice density kg m-3
 RHO_ICE: Final = 917
+
+# other
+SEC_IN_MINUTE: Final = 60
+SEC_IN_HOUR: Final = 3600
+SEC_IN_DAY: Final = 86400
+MM_TO_M: Final = 1e-3
+G_TO_KG: Final = 1e-3

cloudnetpy/instruments/disdrometer/common.py

@@ -6,6 +6,7 @@ from numpy import ma
 
 from cloudnetpy import utils
 from cloudnetpy.cloudnetarray import CloudnetArray
+from cloudnetpy.constants import MM_TO_M, SEC_IN_HOUR, SEC_IN_MINUTE
 from cloudnetpy.exceptions import DisdrometerDataError, ValidTimeStampError
 from cloudnetpy.instruments.cloudnet_instrument import CloudnetInstrument
 from cloudnetpy.instruments.vaisala import values_to_dict
@@ -28,13 +29,12 @@ class Disdrometer(CloudnetInstrument):
         self._file_data = self._read_file()
 
     def convert_units(self) -> None:
-        mm_to_m = 1e3
-        mmh_to_ms = 3600 * mm_to_m
+        mmh_to_ms = SEC_IN_HOUR / MM_TO_M
         c_to_k = 273.15
         self._convert_data(("rainfall_rate_1min_total",), mmh_to_ms)
         self._convert_data(("rainfall_rate",), mmh_to_ms)
         self._convert_data(("rainfall_rate_1min_solid",), mmh_to_ms)
-        self._convert_data(("diameter", "diameter_spread", "diameter_bnds"), mm_to_m)
+        self._convert_data(("diameter", "diameter_spread", "diameter_bnds"), 1e3)
         self._convert_data(("V_sensor_supply",), 10)
         self._convert_data(("I_mean_laser",), 100)
         self._convert_data(("T_sensor",), c_to_k, method="add")
@@ -150,7 +150,9 @@ class Disdrometer(CloudnetInstrument):
             if self.source == PARSIVEL:
                 raise NotImplementedError
             hour, minute, sec = timestamp.split(":")
-            seconds.append(int(hour) * 3600 + int(minute) * 60 + int(sec))
+            seconds.append(
+                int(hour) * SEC_IN_HOUR + int(minute) * SEC_IN_MINUTE + int(sec)
+            )
         return CloudnetArray(utils.seconds2hours(np.array(seconds)), "time")
 
     def _convert_data(self, keys: tuple, value: float, method: str = "divide") -> None:

cloudnetpy/instruments/disdrometer/parsivel.py

@@ -12,6 +12,7 @@ import numpy as np
 
 from cloudnetpy import output
 from cloudnetpy.cloudnetarray import CloudnetArray
+from cloudnetpy.constants import MM_TO_M, SEC_IN_HOUR
 from cloudnetpy.exceptions import DisdrometerDataError
 from cloudnetpy.instruments import instruments
 from cloudnetpy.instruments.cloudnet_instrument import CloudnetInstrument
@@ -148,12 +149,11 @@ class Parsivel(CloudnetInstrument):
         Disdrometer.store_vectors(self.data, n_values, spreads, "diameter")
 
     def convert_units(self) -> None:
-        mm_to_m = 1e3
-        mmh_to_ms = 3600 * mm_to_m
+        mmh_to_ms = SEC_IN_HOUR / MM_TO_M
         c_to_k = 273.15
         self._convert_data(("rainfall_rate",), mmh_to_ms)
         self._convert_data(("snowfall_rate",), mmh_to_ms)
-        self._convert_data(("diameter", "diameter_spread", "diameter_bnds"), mm_to_m)
+        self._convert_data(("diameter", "diameter_spread", "diameter_bnds"), 1e3)
         self._convert_data(("V_sensor_supply",), 10)
         self._convert_data(("T_sensor",), c_to_k, method="add")
 

cloudnetpy/instruments/nc_lidar.py

@@ -1,6 +1,6 @@
 """Module with a class for Lufft chm15k ceilometer."""
 import logging
-from typing import TYPE_CHECKING
+from typing import TYPE_CHECKING, Literal
 
 import numpy as np
 from numpy import ma
@@ -19,7 +19,7 @@ class NcLidar(Ceilometer):
         super().__init__()
         self.dataset: netCDF4.Dataset | None = None
 
-    def _fetch_range(self, reference: str) -> None:
+    def _fetch_range(self, reference: Literal["upper", "lower"]) -> None:
         if self.dataset is None:
             msg = "No dataset found"
             raise RuntimeError(msg)

cloudnetpy/instruments/rpg.py

@@ -11,6 +11,7 @@ from rpgpy import RPGFileError
 from cloudnetpy import output, utils
 from cloudnetpy.categorize.atmos_utils import mmh2ms
 from cloudnetpy.cloudnetarray import CloudnetArray
+from cloudnetpy.constants import G_TO_KG
 from cloudnetpy.exceptions import InconsistentDataError, ValidTimeStampError
 from cloudnetpy.instruments import instruments
 from cloudnetpy.instruments.cloudnet_instrument import CloudnetInstrument
@@ -162,6 +163,8 @@ def _mask_invalid_data(data_in: dict) -> dict:
     data = data_in.copy()
     fill_values = (-999, 1e-10)
     for name in data:
+        if np.issubdtype(data[name].dtype, np.integer) or name == "rainfall_rate":
+            continue
         data[name] = ma.masked_equal(data[name], 0)
         for value in fill_values:
             data[name][data[name] == value] = ma.masked
@@ -234,11 +237,13 @@ class Rpg(CloudnetInstrument):
 
     def convert_time_to_fraction_hour(self, data_type: str | None = None) -> None:
         """Converts time to fraction hour."""
-        key = "time"
-        fraction_hour = utils.seconds2hours(self.raw_data[key])
-        self.data[key] = CloudnetArray(
+        ms2s = 1e-3
+        total_time_sec = self.raw_data["time"] + self.raw_data.get("time_ms", 0) * ms2s
+        fraction_hour = utils.seconds2hours(total_time_sec)
+
+        self.data["time"] = CloudnetArray(
             np.array(fraction_hour),
-            key,
+            "time",
             data_type=data_type,
         )
 
@@ -307,7 +312,7 @@ class Fmcw(Rpg):
     def convert_units(self) -> None:
         """Converts units."""
         self.data["rainfall_rate"].data = mmh2ms(self.data["rainfall_rate"].data)
-        self.data["lwp"].data *= 1e-3  # g -> kg
+        self.data["lwp"].data *= G_TO_KG
 
     @staticmethod
     def _get_instrument(data: dict):

cloudnetpy/instruments/rpg_reader.py

@@ -6,6 +6,7 @@ from numpy import ma
 from numpy.lib import recfunctions as rfn
 from rpgpy import read_rpg
 
+from cloudnetpy.constants import G_TO_KG
 from cloudnetpy.exceptions import ValidTimeStampError
 
 
@@ -256,7 +257,7 @@ class HatproBinLwp(HatproBin):
             ],
             self.header["_n_samples"],
         )
-        self.data["lwp"] *= 1e-3  # Convert from g/m^2 to kg/m^2
+        self.data["lwp"] *= G_TO_KG
 
 
 class HatproBinIwv(HatproBin):

cloudnetpy/instruments/vaisala.py

@@ -5,13 +5,12 @@ import logging
 import numpy as np
 
 from cloudnetpy import utils
+from cloudnetpy.constants import SEC_IN_HOUR, SEC_IN_MINUTE
 from cloudnetpy.exceptions import ValidTimeStampError
 from cloudnetpy.instruments import instruments
 from cloudnetpy.instruments.ceilometer import Ceilometer, NoiseParam
 
 M2KM = 0.001
-SECONDS_IN_MINUTE = 60
-SECONDS_IN_HOUR = 3600
 
 
 class VaisalaCeilo(Ceilometer):
@@ -400,4 +399,4 @@ def values_to_dict(keys: tuple, values: list) -> dict:
 def time_to_fraction_hour(time: str) -> float:
     """Returns time (hh:mm:ss) as fraction hour"""
     hour, minute, sec = time.split(":")
-    return int(hour) + (int(minute) * SECONDS_IN_MINUTE + int(sec)) / SECONDS_IN_HOUR
+    return int(hour) + (int(minute) * SEC_IN_MINUTE + int(sec)) / SEC_IN_HOUR

cloudnetpy/products/iwc.py

@@ -3,6 +3,7 @@ import numpy as np
 from numpy import ma
 
 from cloudnetpy import output, utils
+from cloudnetpy.constants import G_TO_KG
 from cloudnetpy.metadata import MetaData
 from cloudnetpy.products.product_tools import IceClassification, IceSource
 
@@ -83,7 +84,7 @@ class IwcSource(IceSource):
         def _calc_error_in_uncorrected_ice() -> float:
             spec_liq_atten = 1.0 if self.wl_band == 0 else 4.5
             liq_atten_scaled = spec_liq_atten * self.coefficients.Z
-            return lwp_prior * liq_atten_scaled * 2 * 1e-3 * 10
+            return lwp_prior * G_TO_KG * liq_atten_scaled * 2 * 10
 
         lwp_prior = 250  # g m-2
         retrieval_uncertainty = 1.7  # dB

cloudnetpy/utils.py

@@ -7,7 +7,7 @@ import uuid
 import warnings
 from collections.abc import Iterator
 from datetime import timezone
-from typing import Final
+from typing import Literal
 
 import netCDF4
 import numpy as np
@@ -15,15 +15,12 @@ from numpy import ma
 from scipy import ndimage, stats
 from scipy.interpolate import RectBivariateSpline, RegularGridInterpolator, griddata
 
+from cloudnetpy.constants import SEC_IN_DAY, SEC_IN_HOUR, SEC_IN_MINUTE
 from cloudnetpy.exceptions import ValidTimeStampError
 
 Epoch = tuple[int, int, int]
 Date = tuple[str, str, str]
 
-SECONDS_PER_MINUTE: Final = 60
-SECONDS_PER_HOUR: Final = 3600
-SECONDS_PER_DAY: Final = 86400
-
 
 def seconds2hours(time_in_seconds: np.ndarray) -> np.ndarray:
     """Converts seconds since some epoch to fraction hour.
@@ -41,8 +38,8 @@ def seconds2hours(time_in_seconds: np.ndarray) -> np.ndarray:
         Excludes leap seconds.
 
     """
-    seconds_since_midnight = np.mod(time_in_seconds, SECONDS_PER_DAY)
-    fraction_hour = seconds_since_midnight / SECONDS_PER_HOUR
+    seconds_since_midnight = np.mod(time_in_seconds, SEC_IN_DAY)
+    fraction_hour = seconds_since_midnight / SEC_IN_HOUR
     if fraction_hour[-1] == 0:
         fraction_hour[-1] = 24
     return fraction_hour
@@ -60,10 +57,10 @@ def seconds2time(time_in_seconds: float) -> list:
         list: [hours, minutes, seconds] formatted as '05' etc.
 
     """
-    seconds_since_midnight = np.mod(time_in_seconds, SECONDS_PER_DAY)
-    hours = seconds_since_midnight // SECONDS_PER_HOUR
-    minutes = seconds_since_midnight % SECONDS_PER_HOUR // SECONDS_PER_MINUTE
-    seconds = seconds_since_midnight % SECONDS_PER_MINUTE
+    seconds_since_midnight = np.mod(time_in_seconds, SEC_IN_DAY)
+    hours = seconds_since_midnight // SEC_IN_HOUR
+    minutes = seconds_since_midnight % SEC_IN_HOUR // SEC_IN_MINUTE
+    seconds = seconds_since_midnight % SEC_IN_MINUTE
     time = [hours, minutes, seconds]
     return [str(t).zfill(2) for t in time]
 
@@ -97,7 +94,7 @@ def datetime2decimal_hours(data: np.ndarray | list) -> np.ndarray:
     output = []
     for timestamp in data:
         t = timestamp.time()
-        decimal_hours = t.hour + t.minute / 60 + t.second / 3600
+        decimal_hours = t.hour + t.minute / SEC_IN_MINUTE + t.second / SEC_IN_HOUR
         output.append(decimal_hours)
     return np.array(output)
 
@@ -124,7 +121,7 @@ def time_grid(time_step: int = 30) -> np.ndarray:
     if time_step < 1:
         msg = "Time resolution should be >= 1 seconds"
         raise ValueError(msg)
-    half_step = time_step / SECONDS_PER_HOUR / 2
+    half_step = time_step / SEC_IN_HOUR / 2
     return np.arange(half_step, 24 + half_step, half_step * 2)
 
 
@@ -223,7 +220,7 @@ def rebin_1d(
 
     Returns:
     -------
-        Rebinned data with shape (N,).
+        Re-binned data with shape (N,).
 
     """
     edges = binvec(x_new)
@@ -301,7 +298,7 @@ def _filter(array: np.ndarray, structure: np.ndarray) -> np.ndarray:
 
 
 def isbit(array: np.ndarray, nth_bit: int) -> np.ndarray:
-    """Tests if nth bit (0,1,2..) is set.
+    """Tests if nth bit (0,1,2,...) is set.
 
     Args:
     ----
@@ -671,7 +668,7 @@ def n_elements(array: np.ndarray, dist: float, var: str | None = None) -> int:
         array: Input array with arbitrary units or time in fraction hour. *x* should
             be evenly spaced or at least close to.
         dist: Distance to be covered. If x is fraction time, *dist* is in minutes.
-            Otherwise *x* and *dist* should have the same units.
+            Otherwise, *x* and *dist* should have the same units.
         var: If 'time', input is fraction hour and distance in minutes, else inputs
             have the same units. Default is None (same units).
 
@@ -701,11 +698,11 @@ def n_elements(array: np.ndarray, dist: float, var: str | None = None) -> int:
     """
     n = dist / mdiff(array)
     if var == "time":
-        n = n / 60
+        n = n / SEC_IN_MINUTE
     return int(np.round(n))
 
 
-def isscalar(array) -> bool:
+def isscalar(array: np.ndarray | float | list) -> bool:
     """Tests if input is scalar.
 
     By "scalar" we mean that array has a single value.
@@ -723,9 +720,7 @@ def isscalar(array) -> bool:
 
     """
     arr = ma.array(array)
-    if not hasattr(arr, "__len__") or arr.shape == () or len(arr) == 1:
-        return True
-    return False
+    return not hasattr(arr, "__len__") or arr.shape == () or len(arr) == 1
 
 
 def get_time() -> str:
@@ -1000,7 +995,7 @@ def append_data(data_in: dict, key: str, array: np.ndarray) -> dict:
     return data
 
 
-def edges2mid(data: np.ndarray, reference: str) -> np.ndarray:
+def edges2mid(data: np.ndarray, reference: Literal["upper", "lower"]) -> np.ndarray:
     """Shifts values half bin towards up or down.
 
     Args:
@@ -1013,8 +1008,6 @@ def edges2mid(data: np.ndarray, reference: str) -> np.ndarray:
         Shifted values.
 
     """
-    if reference not in ("lower", "upper"):
-        raise ValueError
     gaps = (data[1:] - data[0:-1]) / 2
     if reference == "lower":
         gaps = np.append(gaps, gaps[-1])
@@ -1037,31 +1030,27 @@ def get_file_type(filename: str) -> str:
     raise ValueError(msg)
 
 
-def get_files_with_common_range(files: list) -> list:
+def get_files_with_common_range(filenames: list) -> list:
     """Returns files with the same (most common) number of range gates."""
     n_range = []
-    for file in files:
+    for file in filenames:
         with netCDF4.Dataset(file) as nc:
             n_range.append(len(nc.variables["range"]))
     most_common = np.bincount(n_range).argmax()
-    n_removed = len([n for n in n_range if n != most_common])
-    if n_removed > 0:
+    if n_removed := len(filenames) - n_range.count(int(most_common)) > 0:
         logging.warning(
-            "Removing %s files due to inconsistent height vector",
-            n_removed,
+            "Removing %s files due to inconsistent height vector", n_removed
         )
-    ind = np.where(n_range == most_common)[0]
-    return [file for i, file in enumerate(files) if i in ind]
+    return [file for i, file in enumerate(filenames) if n_range[i] == most_common]
 
 
 def is_all_masked(array: np.ndarray) -> bool:
     """Tests if all values are masked."""
-    if ma.isMaskedArray(array) and hasattr(array, "mask"):
-        return array.mask.all()
-    return False
+    return ma.isMaskedArray(array) and hasattr(array, "mask") and array.mask.all()
 
 
 def find_masked_profiles_indices(array: ma.MaskedArray) -> list:
+    """Finds indices of masked profiles in a 2-D array."""
     non_masked_counts = np.ma.count(array, axis=1)
     masked_profiles_indices = np.where(non_masked_counts == 0)[0]
     return list(masked_profiles_indices)

cloudnetpy/version.py

@@ -1,4 +1,4 @@
 MAJOR = 1
 MINOR = 55
-PATCH = 21
+PATCH = 22
 __version__ = f"{MAJOR}.{MINOR}.{PATCH}"

{cloudnetpy-1.55.21.dist-info → cloudnetpy-1.55.22.dist-info}/METADATA

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

{cloudnetpy-1.55.21.dist-info → cloudnetpy-1.55.22.dist-info}/RECORD

@@ -1,29 +1,29 @@
 cloudnetpy/__init__.py,sha256=X_FqY-4yg5GUj5Edo14SToLEos6JIsC3fN-v1FUgQoA,43
 cloudnetpy/cloudnetarray.py,sha256=HKgWGi5fgW-N1LfESga2Y2TbvjmJxLhe_2eloVQjz-s,6947
 cloudnetpy/concat_lib.py,sha256=beXqjZRaYYZsEBRblOO9_Juwlj3RHIoW8Z2J8HCuZnk,9859
-cloudnetpy/constants.py,sha256=8cZpPTAuJi8gvacK5e36jaNOK-1irlGBmyXEz62tzhQ,469
+cloudnetpy/constants.py,sha256=OMp3pKHCZmdKyRvfO35E7vE3FrS_DHIs_GJuexJLCQk,600
 cloudnetpy/datasource.py,sha256=SCixsM_TWiAbQeGHHdVhPX1kzyY8dSlei2PhH7WqlVg,7980
 cloudnetpy/exceptions.py,sha256=cPh7CDgD2zHeyK1-D1jdqOKs0rPFfpyu9fj8eliGis8,1370
 cloudnetpy/metadata.py,sha256=Bcu1a9UyUq61jomuZ0_6hYIOzf61e5qCXeiwLm46ikw,5040
 cloudnetpy/output.py,sha256=YlCGSFY-xj1PHv_u1qAFH0Kr3KL7GC3jyGGRwhImjSc,14670
 cloudnetpy/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-cloudnetpy/utils.py,sha256=L-aSQsBQG73kTWvMeneBl99LlxTw4qGI_7aqm5e69jI,28158
-cloudnetpy/version.py,sha256=dBQBvia28s9CqNGuF4W6zJTlNGx7LmXKCYKJuWw6ero,73
+cloudnetpy/utils.py,sha256=LVVnsjay8j2DcgvtfC2h0OJ6CN1BC4npzV2jaXtZaGA,28071
+cloudnetpy/version.py,sha256=uQtd9awNyKlJwCZwFAIzB7JIlsT49O3XYAOfsM0mrYA,73
 cloudnetpy/categorize/__init__.py,sha256=gP5q3Vis1y9u9OWgA_idlbjfWXYN_S0IBSWdwBhL_uU,69
 cloudnetpy/categorize/atmos.py,sha256=sti9lZavQ9E59MIaHSkEEkhqvgT4OE1lgcQUznzZQG8,12681
-cloudnetpy/categorize/atmos_utils.py,sha256=Q9MeDvo0Z4P1zLL9fh5RPQeLqhY6cR5QtboceRFtrNY,3865
-cloudnetpy/categorize/categorize.py,sha256=NbfKky19n1JhZCsy6Fu4QEQTby-rE--gk3Gx3wpCZYw,17006
+cloudnetpy/categorize/atmos_utils.py,sha256=lzHIfoyi_6my3Qg-BVXqbq6DkA6AjZi5nfNwyMDi3wU,3886
+cloudnetpy/categorize/categorize.py,sha256=dVy3-D9JY6jlPyjhWI0DJDyxJbEKgt4n2SolSXlhqBU,17024
 cloudnetpy/categorize/classify.py,sha256=w4GsO9a4FkR0wmWsARRNyYZ6tu5Sg6PqNkzfckT23TI,8937
 cloudnetpy/categorize/containers.py,sha256=uXpPXqcysD_iGlC-u1QOL_kyIQN74xtyWhxNMKmOL4c,4674
 cloudnetpy/categorize/droplet.py,sha256=2b2ibrPCeA8xnsXAfqJuo1z1K6P7uRIggAa5QTz4kow,8883
 cloudnetpy/categorize/falling.py,sha256=BJYGMMJo_G8lhXo4K3Rq6qLWL819K0O-mTyFYZp5sAs,4421
 cloudnetpy/categorize/freezing.py,sha256=8ruPFVi9TI7d7RyUkNoXc2WIgRWuOeGpUke9xFaPMiA,3806
 cloudnetpy/categorize/insects.py,sha256=x_Kl5bw7hC-yO4ie9F1OtK1rREguXoLn3UkJROWTMSk,5273
-cloudnetpy/categorize/lidar.py,sha256=KCa7Zo8uWshWfiEiC4JtvhHCanoqOdRXBOlt-MvK3rU,2675
+cloudnetpy/categorize/lidar.py,sha256=9EVqQ85uR57evn1DuUa0w20YekfqHR8hAS2B4TcWl_A,2616
 cloudnetpy/categorize/melting.py,sha256=wnxuWNsMrFFTuXWpLVatExmGBiNGtonYjEWcQK8JGR4,6292
 cloudnetpy/categorize/model.py,sha256=vOr3eF2Q9BcyhDH36ru37RrkPBrCWtpi_n9LIzAuEvU,5601
-cloudnetpy/categorize/mwr.py,sha256=_RM5n6SJLlHESmLcbjpcBpW2qaTnumbtNWAh__Ejwiw,1433
-cloudnetpy/categorize/radar.py,sha256=7fBKVZvFjqlREi254O4_KJIGDvK6zDwRsEyUfL2gU0Q,13089
+cloudnetpy/categorize/mwr.py,sha256=yNvF_ZNJHjfayriovNqsohYSDvviJqEPLsN882zE71s,1457
+cloudnetpy/categorize/radar.py,sha256=cNpMDMVkjvuwYPt_y4obN90EyhTH-8W0o-qJ7DrvDZ8,13095
 cloudnetpy/instruments/__init__.py,sha256=_jejVwi_viSZehmAOkEqTNI-0-exGgAJ_bHW1IRRwTI,398
 cloudnetpy/instruments/basta.py,sha256=qnvVZMSdxvK_YuL3c33QAJYDPWh7ipQXdeGn4p7ikp0,3796
 cloudnetpy/instruments/campbell_scientific.py,sha256=2WHfBKQjtRSl0AqvtPeX7G8Hdi3Dn0WbvoAppFOMbA8,5270
@@ -38,17 +38,17 @@ cloudnetpy/instruments/instruments.py,sha256=GcUEbQFPHUhRTnp300AZ2PK0O2d2EPIGtHq
 cloudnetpy/instruments/lufft.py,sha256=nozeiMRMz7I6q_FwmlxDGhWeJlqTuNh6ru39-M4K3BI,3629
 cloudnetpy/instruments/mira.py,sha256=lq4vmu-faU8kcQum-2jB4wDVkRjIM9wUNsbkM5bpoBs,9386
 cloudnetpy/instruments/mrr.py,sha256=rqb1YZ47a5oiH2a2l5GI0f2BUXtmtKuhYUNMt-8y7oE,5878
-cloudnetpy/instruments/nc_lidar.py,sha256=L7H86wRcTToJS4AAUPabg9nnyVMeSGvSCBwl5rO6kv8,1673
+cloudnetpy/instruments/nc_lidar.py,sha256=Q4sJJwiEPthDz0Zb-laISX32jNYzlUBMafxLJiOAN5c,1704
 cloudnetpy/instruments/nc_radar.py,sha256=vO35xrzxLrLGjAxRj_adDTJ3H4qu1j3jIazl7471L1c,6905
 cloudnetpy/instruments/pollyxt.py,sha256=4X1OHSi1ttg7Y4juLfsejhs_vpomX2NhNMMy89CStvM,8594
 cloudnetpy/instruments/radiometrics.py,sha256=PXVAfT_abyq82CGQwVXzS83-gtI1XZ4YMyj8iFdqXgs,7669
-cloudnetpy/instruments/rpg.py,sha256=6DdmPu9oKHT_RcVKVAvfReTwcvGf4C0D2QqiXq6HhMM,16958
-cloudnetpy/instruments/rpg_reader.py,sha256=iWgk-RlA1hTyREzdzJW9DBaQa44iZEjWwdYqiexb8K8,11358
-cloudnetpy/instruments/vaisala.py,sha256=aB43XVPpDZUFsmkbzcBsNNokAANIHCHtJIVvE_wHKv8,14493
+cloudnetpy/instruments/rpg.py,sha256=D8TDRmypa7ahmClGy4BzHAcn8OcoHreXwd5zEwQjTJA,17185
+cloudnetpy/instruments/rpg_reader.py,sha256=SvYwW0rVuIdAcnqlHCz_Zbenxk-jnGoSYcAfm2kbDuA,11370
+cloudnetpy/instruments/vaisala.py,sha256=OPqS-wV4YvKRiyjov0fhLZzoBMZbNjATfgh1AMuNmyE,14499
 cloudnetpy/instruments/weather_station.py,sha256=OcNOvbZl2m1cEouFk5bndabTNalGVPLqwG3RJXvCiSk,6002
 cloudnetpy/instruments/disdrometer/__init__.py,sha256=lyjwttWvFvuwYxEkusoAvgRcbBmglmOp5HJOpXUqLWo,93
-cloudnetpy/instruments/disdrometer/common.py,sha256=raTlRowixhmXq76e56U0Byh_8o_odr-v7vmNUl0gaEs,15603
-cloudnetpy/instruments/disdrometer/parsivel.py,sha256=iCQK5FauTniT8f3uJtm8bpa9PCjOmwm4jVUvwpKRUiM,19695
+cloudnetpy/instruments/disdrometer/common.py,sha256=nWlVqwvlxei4wJaubBN6NoNsAOpnEqDHvKCPjrAb6Go,15701
+cloudnetpy/instruments/disdrometer/parsivel.py,sha256=I8mP4oH1CCcSnkSWJ0rqWtE0VJvlXzpvACifJJdoM5E,19730
 cloudnetpy/instruments/disdrometer/thies.py,sha256=1hJK7m7utG5QSmsBRWGu0JNsIVnZvCQZ4NZc6-Y8NaI,5095
 cloudnetpy/model_evaluation/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 cloudnetpy/model_evaluation/file_handler.py,sha256=NiBKBnV8o5DRxafKKShBCmVabDkwBH-VpcaAhA6X8j4,6486
@@ -100,15 +100,15 @@ cloudnetpy/products/drizzle.py,sha256=lxEkKcQWoPZN_LkV7TCE5AFNzqa2r1ZwO7V4UKF3MG
 cloudnetpy/products/drizzle_error.py,sha256=o6njlgyNFQ1_O8fG1d9hGCJEPQYU6yiSbMrSlNd5BcM,6153
 cloudnetpy/products/drizzle_tools.py,sha256=oYEHTUd2nEUkeScpUWw3NGAkkbeAo8DsaF7uThhfyqw,10980
 cloudnetpy/products/ier.py,sha256=BJRZtY3OfCIKNHii4OuniniVOazhoKx64fSWmV9qGMM,7826
-cloudnetpy/products/iwc.py,sha256=bJz99EmkwPUxBAsM3PtMOUr_1QztfZKnimJw5KQTrwQ,10129
+cloudnetpy/products/iwc.py,sha256=jvUEIU_PePvVeFbfbdQ8hfHLJbRiQS8d6LAnLentfho,10173
 cloudnetpy/products/lwc.py,sha256=0WOtKiYWwDlyZ6fgkVFCoOTQYNdCPPUkuBrNoE3l-zk,18968
 cloudnetpy/products/mie_lu_tables.nc,sha256=It4fYpqJXlqOgL8jeZ-PxGzP08PMrELIDVe55y9ob58,16637951
 cloudnetpy/products/mwr_multi.py,sha256=mRMzdDhTddhPLBcDiG_2kw_BSRL_01hA9-o4-cmaVaQ,3301
 cloudnetpy/products/mwr_single.py,sha256=4KyxeFg7AphEJg5P7ey8SXacyNAG3PGDOvnksvZj3R8,3168
 cloudnetpy/products/product_tools.py,sha256=dQoFIvM6flnfYXfmm1A71WH7_BZQrf3K_WrsY18Vrgs,10551
 docs/source/conf.py,sha256=IKiFWw6xhUd8NrCg0q7l596Ck1d61XWeVjIFHVSG9Og,1490
-cloudnetpy-1.55.21.dist-info/LICENSE,sha256=wcZF72bdaoG9XugpyE95Juo7lBQOwLuTKBOhhtANZMM,1094
-cloudnetpy-1.55.21.dist-info/METADATA,sha256=Z4CTwNSIDWRrpKdvtY_zOMGTkrnFbgBAurB2p8pURMg,5869
-cloudnetpy-1.55.21.dist-info/WHEEL,sha256=Xo9-1PvkuimrydujYJAjF7pCkriuXBpUPEjma1nZyJ0,92
-cloudnetpy-1.55.21.dist-info/top_level.txt,sha256=ibSPWRr6ojS1i11rtBFz2_gkIe68mggj7aeswYfaOo0,16
-cloudnetpy-1.55.21.dist-info/RECORD,,
+cloudnetpy-1.55.22.dist-info/LICENSE,sha256=wcZF72bdaoG9XugpyE95Juo7lBQOwLuTKBOhhtANZMM,1094
+cloudnetpy-1.55.22.dist-info/METADATA,sha256=75usdsxUS7PzvXmT4ViGMbI0v4J2IyWsHp1OmsFOLbc,5869
+cloudnetpy-1.55.22.dist-info/WHEEL,sha256=Xo9-1PvkuimrydujYJAjF7pCkriuXBpUPEjma1nZyJ0,92
+cloudnetpy-1.55.22.dist-info/top_level.txt,sha256=ibSPWRr6ojS1i11rtBFz2_gkIe68mggj7aeswYfaOo0,16
+cloudnetpy-1.55.22.dist-info/RECORD,,