cloudnetpy 1.66.18__py3-none-any.whl → 1.67.0__py3-none-any.whl

cloudnetpy/categorize/atmos_utils.py

@@ -340,3 +340,18 @@ def normalize_lwc_by_lwp(
     theoretical_lwp = ma.sum(lwc_adiabatic * path_lengths, axis=1)
     scaling_factors = lwp / theoretical_lwp
     return lwc_adiabatic * utils.transpose(scaling_factors)
+
+
+def calc_altitude(temperature: float, pressure: float) -> float:
+    """Calculate altitude (m) based on observed pressure (Pa) and temperature (K)
+    using the International Standard Atmosphere (ISA) model.
+
+    Args:
+        temperature: Observed temperature (K).
+        pressure: Observed atmospheric pressure (Pa).
+
+    Returns:
+        Altitude (m).
+    """
+    L = 0.0065  # Temperature lapse rate (K/m)
+    return (temperature / L) * (1 - (pressure / con.P0) ** (con.RS * L / con.G))

cloudnetpy/categorize/attenuations/melting_attenuation.py

@@ -40,7 +40,7 @@ def calc_melting_attenuation(
     amount[amount == 0] = ma.masked
 
     band = utils.get_wl_band(data.radar.radar_frequency)
-    error_factor = 0.2 if band == 0 else 0.1
+    error_factor = {"Ka": 0.2, "W": 0.1}[band]
 
     error = amount * error_factor
     error[~affected_region] = ma.masked
@@ -65,9 +65,10 @@ def _calc_melting_attenuation(
         Research: Atmospheres, 124, 9520–9533. https://doi.org/10.1029/2019JD030316
 
     """
-    if frequency > 34 and frequency < 37:
+    band = utils.get_wl_band(frequency)
+    if band == "Ka":
         a, b = 0.97, 0.61
-    elif frequency > 93 and frequency < 96:
+    elif band == "W":
         a, b = 2.9, 0.42
     else:
         msg = "Radar frequency not supported"

cloudnetpy/constants.py

@@ -14,6 +14,9 @@ RS: Final = 287.058
 # ice density kg m-3
 RHO_ICE: Final = 917
 
+# Standard atmospheric pressure at sea level Pa
+P0: Final = 1013_25
+
 # other
 SPEED_OF_LIGHT: Final = 3.0e8
 SEC_IN_MINUTE: Final = 60
@@ -30,3 +33,4 @@ PA_TO_HPA: Final = 1 / HPA_TO_PA
 KM_H_TO_M_S: Final = 1000 / SEC_IN_HOUR
 M_TO_KM: Final = 1e-3
 TWO_WAY: Final = 2
+G: Final = 9.80665

cloudnetpy/instruments/cloudnet_instrument.py

@@ -52,14 +52,12 @@ class CloudnetInstrument:
         return float(new_str)
 
     def add_height(self) -> None:
-        try:
-            zenith_angle = ma.median(self.data["zenith_angle"].data)
-        except RuntimeError:
+        zenith_angle = self._get_zenith_angle()
+        if zenith_angle is None:
             logging.warning("Assuming 0 deg zenith_angle")
             zenith_angle = 0
         height = utils.range_to_height(self.data["range"].data, zenith_angle)
         height += self.data["altitude"].data
-        height = np.array(height)
         self.data["height"] = CloudnetArray(height, "height")
 
     def linear_to_db(self, variables_to_log: tuple) -> None:
@@ -106,3 +104,11 @@ class CloudnetInstrument:
             return self.data["time"].data[:]
         except KeyError:
             return self.time
+
+    def _get_zenith_angle(self) -> float | None:
+        if "zenith_angle" not in self.data or self.data["zenith_angle"].data.size == 0:
+            return None
+        zenith_angle = ma.median(self.data["zenith_angle"].data)
+        if np.isnan(zenith_angle) or zenith_angle is ma.masked:
+            return None
+        return zenith_angle

cloudnetpy/instruments/instruments.py

@@ -75,6 +75,14 @@ CHM15KX = Instrument(
     wavelength=1064.0,
 )
 
+MIRA10 = Instrument(
+    manufacturer="METEK",
+    domain="radar",
+    category="cloud radar",
+    model="MIRA-10",
+    frequency=9.4,  # 9.2 - 9.6 GHz
+)
+
 MIRA35 = Instrument(
     manufacturer="METEK",
     domain="radar",
@@ -142,11 +150,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/instruments/mira.py

@@ -8,7 +8,7 @@ from tempfile import NamedTemporaryFile, TemporaryDirectory
 from numpy import ma
 
 from cloudnetpy import concat_lib, output, utils
-from cloudnetpy.instruments.instruments import MIRA35
+from cloudnetpy.instruments.instruments import MIRA10, MIRA35
 from cloudnetpy.instruments.nc_radar import NcRadar
 from cloudnetpy.metadata import MetaData
 
@@ -99,7 +99,13 @@ class Mira(NcRadar):
     def __init__(self, full_path: str, site_meta: dict):
         super().__init__(full_path, site_meta)
         self.date = self._init_mira_date()
-        self.instrument = MIRA35
+        if "model" not in site_meta or site_meta["model"] == "mira-35":
+            self.instrument = MIRA35
+        elif site_meta["model"] == "mira-10":
+            self.instrument = MIRA10
+        else:
+            msg = f"Invalid model: {site_meta['model']}"
+            raise ValueError(msg)
 
     def screen_by_date(self, expected_date: str) -> None:
         """Screens incorrect time stamps."""

cloudnetpy/model_evaluation/products/observation_products.py

@@ -112,7 +112,7 @@ class ObservationManager(DataSource):
             raise RuntimeError(msg)
         wband = utils.get_wl_band(float(self.radar_freq))
         rainrate_threshold = 8
-        if 90 < wband < 100:
+        if wband == "W":
             rainrate_threshold = 2
         return rainrate_threshold
 

cloudnetpy/plotting/plotting.py

@@ -21,7 +21,8 @@ from mpl_toolkits.axes_grid1 import make_axes_locatable
 from numpy import ma, ndarray
 from scipy.ndimage import uniform_filter
 
-from cloudnetpy.constants import T0
+from cloudnetpy import constants as con
+from cloudnetpy.categorize.atmos_utils import calc_altitude
 from cloudnetpy.exceptions import PlottingError
 from cloudnetpy.instruments.ceilometer import calc_sigma_units
 from cloudnetpy.plotting.plot_meta import ATTRIBUTES, PlotMeta
@@ -45,6 +46,9 @@ 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.
+        plot_above_ground: Whether to plot above ground instead of above mean sea level.
+
     """
 
     dpi: float = 120
@@ -58,6 +62,8 @@ class PlotParameters:
     footer_text: str | None = None
     plot_meta: PlotMeta | None = None
     raise_on_empty: bool = False
+    minor_ticks: bool = False
+    plot_above_ground: bool = False
 
 
 class Dimensions:
@@ -170,15 +176,35 @@ class FigureData:
         m2mm = 1e3
         file_type = getattr(self.file, "cloudnet_file_type", "")
         if file_type == "model":
-            return ma.mean(self.file.variables["height"][:], axis=0) * m2km
+            height = ma.mean(self.file.variables["height"][:], axis=0)  # height AGL
+            if not self.options.plot_above_ground:
+                site_alt = self._calc_ground_altitude()
+                height += site_alt
+            return height * m2km
         if "height" in self.file.variables:
-            return self.file.variables["height"][:] * m2km
+            height = self.file.variables["height"][:]  # height AMSL
+            if self.options.plot_above_ground:
+                if "altitude" not in self.file.variables:
+                    msg = "No altitude information in the file."
+                    raise ValueError(msg)
+                height -= self.file.variables["altitude"][:]
+            return height * m2km
         if "range" in self.file.variables:
             return self.file.variables["range"][:] * m2km
         if "diameter" in self.file.variables:
             return self.file.variables["diameter"][:] * m2mm
         return None
 
+    def _calc_ground_altitude(self) -> float:
+        if (
+            "sfc_geopotential" in self.file.variables
+            and "gdas1" not in self.file.source.lower()  # uncertain unit in gdas1
+        ):
+            return np.mean(self.file.variables["sfc_geopotential"][:]) / con.G
+        pressure = ma.mean(self.file.variables["pressure"][:, 0])
+        temperature = ma.mean(self.file.variables["temperature"][:, 0])
+        return calc_altitude(temperature, pressure)
+
     def is_mwrpy_product(self) -> bool:
         cloudnet_file_type = getattr(self.file, "cloudnet_file_type", "")
         return cloudnet_file_type in ("mwr-single", "mwr-multi")
@@ -198,7 +224,8 @@ class SubPlot:
         self.ax = ax
         self.variable = variable
         self.options = options
-        self.plot_meta = self._read_plot_meta(file_type)
+        self.file_type = file_type
+        self.plot_meta = self._read_plot_meta()
 
     def set_xax(self) -> None:
         resolution = 4
@@ -210,6 +237,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(
@@ -217,7 +247,11 @@ class SubPlot:
         ylabel: str | None = None,
         y_limits: tuple[float, float] | None = None,
     ) -> None:
-        label = ylabel if ylabel is not None else "Height (km)"
+        height_str = (
+            "Height (km AGL)" if self.options.plot_above_ground else "Height (km AMSL)"
+        )
+
+        label = ylabel if ylabel is not None else height_str
         self.ax.set_ylabel(label, fontsize=13)
         if y_limits is not None:
             self.ax.set_ylim(*y_limits)
@@ -282,12 +316,12 @@ class SubPlot:
                 va="bottom",
             )
 
-    def _read_plot_meta(self, file_type: str | None) -> PlotMeta:
+    def _read_plot_meta(self) -> PlotMeta:
         if self.options.plot_meta is not None:
             plot_meta = self.options.plot_meta
         else:
             fallback = ATTRIBUTES["fallback"].get(self.variable.name, PlotMeta())
-            file_attributes = ATTRIBUTES.get(file_type or "", {})
+            file_attributes = ATTRIBUTES.get(self.file_type or "", {})
             plot_meta = file_attributes.get(self.variable.name, fallback)
         if plot_meta.clabel is None:
             plot_meta = plot_meta._replace(clabel=_reformat_units(self.variable.units))
@@ -542,7 +576,7 @@ class Plot2D(Plot):
             self._plot_contour(
                 figure_data,
                 alt,
-                levels=np.array([T0]),
+                levels=np.array([con.T0]),
                 colors="gray",
                 linewidths=1.25,
                 linestyles="dashed",
@@ -910,6 +944,7 @@ 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/products/drizzle_tools.py

@@ -72,7 +72,12 @@ class DrizzleSource(DataSource):
         """Returns string corresponding the radar frequency."""
         radar_frequency = float(self.getvar("radar_frequency"))
         wl_band = utils.get_wl_band(radar_frequency)
-        return "35" if wl_band == 0 else "94"
+        if wl_band == "Ka":
+            return "35"
+        if wl_band == "W":
+            return "94"
+        msg = f"Unsupported band: {wl_band}"
+        raise ValueError(msg)
 
 
 class DrizzleClassification(ProductClassification):

cloudnetpy/products/ier.py

@@ -3,7 +3,7 @@
 import numpy as np
 from numpy import ma
 
-from cloudnetpy import constants, output, utils
+from cloudnetpy import constants, output
 from cloudnetpy.metadata import MetaData
 from cloudnetpy.products.iwc import DEFINITIONS as IWC_DEFINITION
 from cloudnetpy.products.product_tools import IceClassification, IceSource
@@ -83,7 +83,7 @@ class IerSource(IceSource):
 
 
 def _add_ier_comment(attributes: dict, ier: IerSource) -> dict:
-    freq = utils.get_frequency(ier.wl_band)
+    freq = ier.radar_frequency
     coeffs = ier.coefficients
     factor = np.round((coeffs[0] / 0.93), 3)
     attributes["ier"] = attributes["ier"]._replace(

cloudnetpy/products/iwc.py

@@ -113,7 +113,7 @@ def _add_iwc_error_comment(attributes: dict, lwp_prior, uncertainty: float) -> d
 
 
 def _add_iwc_comment(attributes: dict, iwc: IwcSource) -> dict:
-    freq = utils.get_frequency(iwc.wl_band)
+    freq = iwc.radar_frequency
     coeffs = iwc.coefficients
     factor = round((coeffs[0] / 0.93) * 1000) / 1000
     attributes["iwc"] = attributes["iwc"]._replace(

cloudnetpy/products/product_tools.py

@@ -181,7 +181,8 @@ class IceSource(DataSource):
 
     def __init__(self, categorize_file: str, product: str):
         super().__init__(categorize_file)
-        self.wl_band = utils.get_wl_band(float(self.getvar("radar_frequency")))
+        self.radar_frequency = float(self.getvar("radar_frequency"))
+        self.wl_band = utils.get_wl_band(self.radar_frequency)
         self.temperature = _get_temperature(categorize_file)
         self.product = product
         self.coefficients = self._get_coefficients()
@@ -216,13 +217,18 @@ class IceSource(DataSource):
         References:
             Hogan et.al. 2006, https://doi.org/10.1175/JAM2340.1
         """
+        msg = f"Unsupported band: {self.wl_band}"
         if self.product == "ier":
-            if self.wl_band == 0:
+            if self.wl_band == "Ka":
                 return IceCoefficients(0.878, -0.000205, -0.0015, 0.0016, -1.52)
-            return IceCoefficients(0.669, -0.000296, -0.00193, -0.000, -1.502)
-        if self.wl_band == 0:
+            if self.wl_band == "W":
+                return IceCoefficients(0.669, -0.000296, -0.00193, -0.000, -1.502)
+            raise ValueError(msg)
+        if self.wl_band == "Ka":
             return IceCoefficients(0.878, 0.000242, -0.0186, 0.0699, -1.63)
-        return IceCoefficients(0.669, 0.000580, -0.00706, 0.0923, -0.992)
+        if self.wl_band == "W":
+            return IceCoefficients(0.669, 0.000580, -0.00706, 0.0923, -0.992)
+        raise ValueError(msg)
 
     def _convert_z(self, z_variable: str = "Z") -> np.ndarray:
         """Calculates temperature weighted z, i.e. ice effective radius [m]."""

cloudnetpy/utils.py

@@ -691,22 +691,24 @@ def get_uuid() -> str:
     return str(uuid.uuid4())
 
 
-def get_wl_band(radar_frequency: float) -> int:
-    """Returns integer corresponding to radar frequency.
+def get_wl_band(radar_frequency: float) -> Literal["X", "Ka", "W"]:
+    """Returns IEEE radar band corresponding to radar frequency.
 
     Args:
         radar_frequency: Radar frequency (GHz).
 
     Returns:
-        0 = 35GHz radar, 1 = 94Ghz radar.
+        IEEE radar band as string.
 
     """
-    return 0 if (30 < radar_frequency < 40) else 1
-
-
-def get_frequency(wl_band: int) -> str:
-    """Returns radar frequency string corresponding to wl band."""
-    return "35.5" if wl_band == 0 else "94"
+    if 8 < radar_frequency < 12:
+        return "X"
+    if 27 < radar_frequency < 40:
+        return "Ka"
+    if 75 < radar_frequency < 110:
+        return "W"
+    msg = f"Unknown band: {radar_frequency} GHz"
+    raise ValueError(msg)
 
 
 def transpose(data: np.ndarray) -> np.ndarray:

cloudnetpy/version.py

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

{cloudnetpy-1.66.18.dist-info → cloudnetpy-1.67.0.dist-info}/METADATA

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

{cloudnetpy-1.66.18.dist-info → cloudnetpy-1.67.0.dist-info}/RECORD

@@ -2,16 +2,16 @@ cloudnetpy/__init__.py,sha256=X_FqY-4yg5GUj5Edo14SToLEos6JIsC3fN-v1FUgQoA,43
 cloudnetpy/cli.py,sha256=lHkeAErmAijI-Ugpd4DHRHfbZP4SXOake0LIY5Ovv_Q,20782
 cloudnetpy/cloudnetarray.py,sha256=Ol1ha4RPAmFZANL__U5CaMKX4oYMXYR6OnjoCZ9w3eo,7077
 cloudnetpy/concat_lib.py,sha256=lLjEC0ILPI1ghv8Wu9WVDiQoAkwwBKlW5jV8KoIde_8,11963
-cloudnetpy/constants.py,sha256=RDB9aqpBRztk3QVCFgsmi9fwhtLuit_0WJrt0D6sDcc,736
+cloudnetpy/constants.py,sha256=jjW1eO5BwmLjtCKXKGXtKF_BoKpvSNrRfpgW43fWT_0,824
 cloudnetpy/datasource.py,sha256=FcWS77jz56gIzwnbafDLdj-HjAyu0P_VtY7gkeVZThU,7952
 cloudnetpy/exceptions.py,sha256=ns48useL9RN3mPh7CqIiLA19VI9OmVbyRsKTkwbThF8,1760
 cloudnetpy/metadata.py,sha256=DOGt7EQLS-AVJEszrUrpXr3gHVQv655FzeCzKrOPvoU,5477
 cloudnetpy/output.py,sha256=lq4YSeMT_d-j4rlQkKm9KIZ8boupTBBBKV1eUawpmCI,15672
 cloudnetpy/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-cloudnetpy/utils.py,sha256=JksYOwf9ORiR_QpoKrTe1JJwXpPYJj-wlwaZKCHoh3o,29744
-cloudnetpy/version.py,sha256=YdfbC1Is18-wsauv4KijktdV_lghueKCQgNzKMMazKo,73
+cloudnetpy/utils.py,sha256=U0iMIKPiKLrLVAfs_u9pPuoWYW1RJHcM8dbLF9a4yIA,29796
+cloudnetpy/version.py,sha256=9-gHFlbKlJNv6v3rFfs2bpwbFn9E2SIntv6AB3MKCTc,72
 cloudnetpy/categorize/__init__.py,sha256=s-SJaysvVpVVo5kidiruWQO6p3gv2TXwY1wEHYO5D6I,44
-cloudnetpy/categorize/atmos_utils.py,sha256=9-ymI6i1xASf-XAFyO87FaTfvq6bF89N1i_27OkUp-M,10104
+cloudnetpy/categorize/atmos_utils.py,sha256=RcmbKxm2COkE7WEya0mK3yX5rzUbrewRVh3ekm01RtM,10598
 cloudnetpy/categorize/attenuation.py,sha256=Y_-fzmQTltWTqIZTulJhovC7a6ifpMcaAazDJcnMIOc,990
 cloudnetpy/categorize/categorize.py,sha256=HtttmUlnLE4Cw3QLdz-DRfRpd8sVr_Tp4A62XI67x24,20739
 cloudnetpy/categorize/classify.py,sha256=qovHgHsMku5kpl3cJxKteNBsG8GAkfI3Zo8QhJwZSFQ,8512
@@ -30,7 +30,7 @@ cloudnetpy/categorize/radar.py,sha256=0Wg5u2aLXRVhHiFiXb4fSqY_iGgSwcCMMVshM3wBoo
 cloudnetpy/categorize/attenuations/__init__.py,sha256=CWFHVWeTIe2hrZtgkJaX2HGftbuffsFc39Mzv5B0Lw0,1037
 cloudnetpy/categorize/attenuations/gas_attenuation.py,sha256=emr-RCxQT0i2N8k6eBNhRsmsCBPHJzQsWJfjC4fVSTo,975
 cloudnetpy/categorize/attenuations/liquid_attenuation.py,sha256=0p0G79BPkw1itCXHMwbvkNHtJGBocJzow3gNHAirChI,3036
-cloudnetpy/categorize/attenuations/melting_attenuation.py,sha256=-NDv8YZygqEsRHBs6RTkJ85Lu52o1Xo5ZUMI4cBoDbA,2334
+cloudnetpy/categorize/attenuations/melting_attenuation.py,sha256=9c9xoZHtGUbjFYJxkVc3UUDHLDy0UbNUZ32ITtnsj5w,2333
 cloudnetpy/categorize/attenuations/rain_attenuation.py,sha256=qazJzRyXf9vbjJhh4yiFmABI4L57j5W_6YZ-6qjRiBI,2839
 cloudnetpy/instruments/__init__.py,sha256=_jejVwi_viSZehmAOkEqTNI-0-exGgAJ_bHW1IRRwTI,398
 cloudnetpy/instruments/basta.py,sha256=_OTnySd36ktvxk_swWBzbv_H4AVGlkF_Ce3KtPGD1rE,3758
@@ -38,13 +38,13 @@ cloudnetpy/instruments/campbell_scientific.py,sha256=2WHfBKQjtRSl0AqvtPeX7G8Hdi3
 cloudnetpy/instruments/ceilo.py,sha256=xrI7iYNftKvGZf-3C_ESUNsu-QhXV43iWkDuKp3biZU,9552
 cloudnetpy/instruments/ceilometer.py,sha256=pdmLVljsuciyKpaGxWxJ_f1IrJK-UrkBC0lSeuirLlU,12095
 cloudnetpy/instruments/cl61d.py,sha256=g6DNBFju3wYhLFl32DKmC8pUup7y-EupXoUU0fuoGGA,1990
-cloudnetpy/instruments/cloudnet_instrument.py,sha256=6qQSyQXHBTdOvzqTxOXtx6RYDm7vLVcwNYTx1nV_kM8,3885
+cloudnetpy/instruments/cloudnet_instrument.py,sha256=NQZ_FMXh8iyzXYSCKQSpIdp0MZFAh7WqKE8mZRmVbF4,4164
 cloudnetpy/instruments/copernicus.py,sha256=nmgqGOjVQFngj7BNbpcuCwA-W3yksvBbqn__iq7MyDk,6469
 cloudnetpy/instruments/galileo.py,sha256=yQBedd7dmDnwuWi1MtXOsg4-RyRx0uRAXumCY4YuH9k,4686
 cloudnetpy/instruments/hatpro.py,sha256=DzCWzTJxTc5BSOgoeyM8RjYkSXX6NDi3QXgKRp0uxlI,8759
-cloudnetpy/instruments/instruments.py,sha256=jG5TYnZ8bdCZXnI303ZsaJBEdSKaIjKMbkGtnq6kQX0,3261
+cloudnetpy/instruments/instruments.py,sha256=cHP0RN-Z8Jl9yoDHdvaOflTrzuogDTGmd-nxSOk8Uq4,3568
 cloudnetpy/instruments/lufft.py,sha256=ugXF6pssHAAz1Y_hqPdpKuluAjxxHSR88xBmQuS6RlI,3705
-cloudnetpy/instruments/mira.py,sha256=zfpPLCJLRc11wB-ddkQC9JALdwz73EgVbknv8iDJtuE,9847
+cloudnetpy/instruments/mira.py,sha256=W8mAAc5tP4nxexiN4MhGEj4ZV2-Tz2AXp5Cpt4XPaN4,10119
 cloudnetpy/instruments/mrr.py,sha256=eeAzCp3CiHGauywjwvMUAFwZ4vBOZMcd3IlF8KsrLQo,5711
 cloudnetpy/instruments/nc_lidar.py,sha256=5gQG9PApnNPrHmS9_zanl8HEYIQuGRpbnzC3wfTcOyQ,1705
 cloudnetpy/instruments/nc_radar.py,sha256=PM4XQm2Ri9pLOh2mt9dphjdXr1_aPKMFWp2T_hDNo3c,6730
@@ -72,7 +72,7 @@ cloudnetpy/model_evaluation/products/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeu
 cloudnetpy/model_evaluation/products/advance_methods.py,sha256=rng3ZLR1Arv1AGUzq0Ehu-65628PC5LZVKpHSUpCIW8,8526
 cloudnetpy/model_evaluation/products/grid_methods.py,sha256=4no7mbKc9HlEXSNKPioqLmFZxUefuI-yqX0-Ej2jMzU,9067
 cloudnetpy/model_evaluation/products/model_products.py,sha256=uWi7zXQI7kR_ju0SL_BC1wozcq5DhaCcT-XZq33Q-bA,6861
-cloudnetpy/model_evaluation/products/observation_products.py,sha256=UJTb6-1Z-HOB9ltW8E7QHNEoByXeDQQUNi3GyS5lG_0,5500
+cloudnetpy/model_evaluation/products/observation_products.py,sha256=-vUJi54qbgM3kbp0Im3BcW8RcB-aCw397XIaqfS_4vk,5496
 cloudnetpy/model_evaluation/products/product_resampling.py,sha256=IuWvtwpya76URh1WmTTgtLxAo4HZxkz6GmftpZkMCGo,3640
 cloudnetpy/model_evaluation/products/tools.py,sha256=lXnQ9XIEf5zqk_haY3mSrekPyGbAwNWvd6ZOol1Ip1Q,2918
 cloudnetpy/model_evaluation/statistics/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
@@ -101,23 +101,23 @@ cloudnetpy/model_evaluation/tests/unit/test_statistical_methods.py,sha256=Ra3r4V
 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=d7CIT4ltUsjKw0HNnOXPkaLdZbkSmTuq6AbKdbODfZE,16730
-cloudnetpy/plotting/plotting.py,sha256=fyeOUS5Fd6t1V2Ud7aiDgCne4sNQMZrR1XLT-49SZb8,35524
+cloudnetpy/plotting/plotting.py,sha256=E_UTnraLxYz8jEROOpQr0pFbntHpI8nL6gMcwgsyFj8,37109
 cloudnetpy/products/__init__.py,sha256=2hRb5HG9hNrxH1if5laJkLeFeaZCd5W1q3hh4ewsX0E,273
 cloudnetpy/products/classification.py,sha256=KwAiBSgFwDqhM114NIgYiUjj8HoYc7gAlc8E1QgcSig,8207
 cloudnetpy/products/der.py,sha256=soypE7uSEP4uHUCCQVEhyXsKY6e9mzV9B_2S5GUizqk,12729
 cloudnetpy/products/drizzle.py,sha256=58C9Mo6oRXR8KpbVPghbJvHvFX9GfS3xUp058pbf0qw,10804
 cloudnetpy/products/drizzle_error.py,sha256=4GwlHRtNbk9ks7bGtXCco-wXbcDOKeAQwKmbhzut6Qk,6132
-cloudnetpy/products/drizzle_tools.py,sha256=SXD3GProk30NgRFIhvwLyeTEQ4PQPXFZ2RKgBR9WjbY,10870
-cloudnetpy/products/ier.py,sha256=yv5bTBtA1tOB-dAt71Gnou1ywvrnhZRhmWk84p3oJRA,6006
-cloudnetpy/products/iwc.py,sha256=_focccRSkgrS7wgKUY5ASfQ_PcmEZNQkvuvu3q9QA4o,9475
+cloudnetpy/products/drizzle_tools.py,sha256=HLxUQ89mFNo6IIe6Cj3ZH-TPkJdpMxKCOt4cOOmcLs0,11002
+cloudnetpy/products/ier.py,sha256=70jyYrhO-kAHuQ3CpDVcyKKuHT3c9Eby6ADHHlPjAFY,5986
+cloudnetpy/products/iwc.py,sha256=FQtYFYnMNuDnjuK1ee3zOfzPtKsoWH95CXtEz1iplZI,9462
 cloudnetpy/products/lwc.py,sha256=sl6Al2tuH3KkCBrPbWTmuz3jlD5UQJ4D6qBsn1tt2CQ,18962
 cloudnetpy/products/mie_lu_tables.nc,sha256=It4fYpqJXlqOgL8jeZ-PxGzP08PMrELIDVe55y9ob58,16637951
 cloudnetpy/products/mwr_tools.py,sha256=rd7UC67O4fsIE5SaHVZ4qWvUJTj41ZGwgQWPwZzOM14,5377
-cloudnetpy/products/product_tools.py,sha256=01Zc6xV8CSuYcIcLpchFf5POL3_c629-YMNDZJ51udA,10853
+cloudnetpy/products/product_tools.py,sha256=uu4l6reuGbPcW3TgttbaSrqIKbyYGhBVTdnC7opKvmg,11101
 docs/source/conf.py,sha256=IKiFWw6xhUd8NrCg0q7l596Ck1d61XWeVjIFHVSG9Og,1490
-cloudnetpy-1.66.18.dist-info/LICENSE,sha256=wcZF72bdaoG9XugpyE95Juo7lBQOwLuTKBOhhtANZMM,1094
-cloudnetpy-1.66.18.dist-info/METADATA,sha256=WGYHjn8d1aUCtLC6Kj848evl1l5I0_vgZNF5bmWpbOc,5794
-cloudnetpy-1.66.18.dist-info/WHEEL,sha256=PZUExdf71Ui_so67QXpySuHtCi3-J3wvF4ORK6k_S8U,91
-cloudnetpy-1.66.18.dist-info/entry_points.txt,sha256=HhY7LwCFk4qFgDlXx_Fy983ZTd831WlhtdPIzV-Y3dY,51
-cloudnetpy-1.66.18.dist-info/top_level.txt,sha256=ibSPWRr6ojS1i11rtBFz2_gkIe68mggj7aeswYfaOo0,16
-cloudnetpy-1.66.18.dist-info/RECORD,,
+cloudnetpy-1.67.0.dist-info/LICENSE,sha256=wcZF72bdaoG9XugpyE95Juo7lBQOwLuTKBOhhtANZMM,1094
+cloudnetpy-1.67.0.dist-info/METADATA,sha256=hslGwWAT2tUoak-gag54pZhHlz_1tS_huYydwN496RA,5793
+cloudnetpy-1.67.0.dist-info/WHEEL,sha256=PZUExdf71Ui_so67QXpySuHtCi3-J3wvF4ORK6k_S8U,91
+cloudnetpy-1.67.0.dist-info/entry_points.txt,sha256=HhY7LwCFk4qFgDlXx_Fy983ZTd831WlhtdPIzV-Y3dY,51
+cloudnetpy-1.67.0.dist-info/top_level.txt,sha256=ibSPWRr6ojS1i11rtBFz2_gkIe68mggj7aeswYfaOo0,16
+cloudnetpy-1.67.0.dist-info/RECORD,,