cloudnetpy 1.96.1__tar.gz → 1.96.2__tar.gz

{cloudnetpy-1.96.1/cloudnetpy.egg-info → cloudnetpy-1.96.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: cloudnetpy
-Version: 1.96.1
+Version: 1.96.2
 Summary: Python package for Cloudnet processing
 Author: Simo Tukiainen
 License: MIT License
@@ -43,7 +43,7 @@ Requires-Dist: ceilopyter>=0.2.0
 Requires-Dist: cloudnet-api-client
 Requires-Dist: doppy>=0.5.0
 Requires-Dist: matplotlib
-Requires-Dist: mwrpy>=1.3.5
+Requires-Dist: mwrpy>=1.7.0
 Requires-Dist: netCDF4
 Requires-Dist: requests
 Requires-Dist: rpgpy>=0.14.5

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/categorize/attenuations/__init__.py

@@ -31,7 +31,7 @@ def calc_two_way_attenuation(
     """Calculates two-way attenuation (dB) for given specific attenuation
     (dB km-1) and height above ground level (m).
     """
-    path_lengths = path_lengths_from_ground(height_agl) * con.M_TO_KM  # km
+    path_lengths = path_lengths_from_ground(height_agl) * con.M_TO_KM
     one_way_attenuation = specific_attenuation * path_lengths
     accumulated_attenuation = ma.cumsum(one_way_attenuation, axis=1)
-    return accumulated_attenuation * con.TWO_WAY
+    return accumulated_attenuation * 2  # two way

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/categorize/classify.py

@@ -1,5 +1,6 @@
 import numpy as np
 import numpy.typing as npt
+from atmoslib.constants import T0
 from numpy import ma
 from scipy import ndimage
 
@@ -18,7 +19,6 @@ from cloudnetpy.categorize.containers import (
     ClassificationResult,
     Observations,
 )
-from cloudnetpy.constants import T0
 from cloudnetpy.products.product_tools import CategoryBits, QualityBits
 
 

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/categorize/containers.py

@@ -4,10 +4,11 @@ from os import PathLike
 
 import numpy as np
 import numpy.typing as npt
+from atmoslib.constants import T0
 from numpy import ma
 
 from cloudnetpy import utils
-from cloudnetpy.constants import MM_H_TO_M_S, T0
+from cloudnetpy.constants import MM_H_TO_M_S
 from cloudnetpy.products.product_tools import CategoryBits
 
 from .disdrometer import Disdrometer

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/categorize/falling.py

@@ -2,11 +2,11 @@
 
 import numpy as np
 import numpy.typing as npt
+from atmoslib.constants import T0
 from numpy import ma
 
 from cloudnetpy.categorize import atmos_utils
 from cloudnetpy.categorize.containers import ClassData
-from cloudnetpy.constants import T0
 
 
 def find_falling_hydrometeors(

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/categorize/freezing.py

@@ -4,12 +4,12 @@ import logging
 
 import numpy as np
 import numpy.typing as npt
+from atmoslib.constants import T0
 from numpy import ma
 from scipy.interpolate import interp1d
 
 from cloudnetpy import utils
 from cloudnetpy.categorize.containers import ClassData
-from cloudnetpy.constants import T0
 
 
 def find_freezing_region(obs: ClassData, melting_layer: npt.NDArray) -> npt.NDArray:

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/categorize/melting.py

@@ -2,13 +2,13 @@
 
 import numpy as np
 import numpy.typing as npt
+from atmoslib.constants import T0
 from numpy import ma
 from scipy.ndimage import gaussian_filter
 
 from cloudnetpy import utils
 from cloudnetpy.categorize import droplet
 from cloudnetpy.categorize.containers import ClassData
-from cloudnetpy.constants import T0
 
 
 def find_melting_layer(obs: ClassData, *, smooth: bool = True) -> npt.NDArray:

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/categorize/radar.py

@@ -7,11 +7,10 @@ from os import PathLike
 import numpy as np
 import numpy.typing as npt
 from numpy import ma
-from scipy import constants
 
 from cloudnetpy import utils
 from cloudnetpy.categorize.attenuations import RadarAttenuation
-from cloudnetpy.constants import GHZ_TO_HZ, SEC_IN_HOUR, SPEED_OF_LIGHT
+from cloudnetpy.constants import GHZ_TO_HZ, LN_TO_DB, SEC_IN_HOUR, SPEED_OF_LIGHT
 from cloudnetpy.datasource import DataSource
 
 
@@ -265,8 +264,7 @@ class Radar(DataSource):
             """
             noise_threshold = 3
             n_pulses = _number_of_independent_pulses()
-            ln_to_log10 = 10 / np.log(10)
-            z_precision = ma.divide(ln_to_log10, np.sqrt(n_pulses)) * (
+            z_precision = ma.divide(LN_TO_DB, np.sqrt(n_pulses)) * (
                 1 + (utils.db2lin(z_power_min - z_power) / noise_threshold)
             )
 
@@ -439,8 +437,7 @@ class Radar(DataSource):
 
 
 def _prf_to_folding_velocity(prf: npt.NDArray, radar_frequency: float) -> npt.NDArray:
-    ghz_to_hz = 1e9
     if len(prf) != 1:
         msg = "Unable to determine folding velocity"
         raise RuntimeError(msg)
-    return prf[0] * constants.c / (4 * radar_frequency * ghz_to_hz)
+    return prf[0] * SPEED_OF_LIGHT / (4 * radar_frequency * GHZ_TO_HZ)

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/constants.py

@@ -1,24 +1,21 @@
 """Constants used in Cloudnet processing."""
 
+import math
 from typing import Final
 
-# Triple point of water
-T0: Final = 273.16
+T_FREEZING: Final = 273.15
+"""Freezing point of water at standard pressure (K), i.e. 0 degrees Celsius"""
 
-# Ratio of the molecular weight of water vapor to dry air
-MW_RATIO: Final = 0.62198
-
-# Specific gas constant for dry  air (J kg-1 K-1)
-RS: Final = 287.058
-
-# ice density kg m-3
 RHO_ICE: Final = 917
+"""Ice density (kg m-3)"""
+
+RHO_WATER: Final = 1000
+"""Liquid water density (kg m-3)"""
 
-# Standard atmospheric pressure at sea level Pa
-P0: Final = 1013_25
+SPEED_OF_LIGHT: Final = 299_792_458
+"""Speed of light in vacuum (m s-1)"""
 
-# other
-SPEED_OF_LIGHT: Final = 3.0e8
+# Unit conversions
 SEC_IN_MINUTE: Final = 60
 SEC_IN_HOUR: Final = 3600
 SEC_IN_DAY: Final = 86400
@@ -35,7 +32,10 @@ MM_H_TO_M_S: Final = 1 / M_S_TO_MM_H
 GHZ_TO_HZ: Final = 1e9
 HZ_TO_GHZ: Final = 1e-9
 HPA_TO_PA: Final = 100
-PA_TO_HPA: Final = 1 / HPA_TO_PA
 KM_H_TO_M_S: Final = 1000 / SEC_IN_HOUR
-TWO_WAY: Final = 2
-G: Final = 9.80665
+
+CM_TO_KG_M2: Final = 10
+"""Liquid water equivalent cm to kg m-2"""
+
+LN_TO_DB: Final = 10 / math.log(10)
+"""Natural logarithm of a power ratio to decibels"""

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/disdronator/parsivel.py

@@ -9,6 +9,7 @@ import netCDF4
 import numpy as np
 import numpy.typing as npt
 
+from cloudnetpy.constants import M_S_TO_MM_H
 from cloudnetpy.disdronator.process import DisdroL1
 from cloudnetpy.disdronator.utils import convert_to_numpy
 
@@ -434,7 +435,7 @@ def _read_parsivel2nc(filename: str | PathLike) -> ParsivelOutput:
         time = cftime.num2pydate(nc["time"][:], units=nc["time"].units)
         data = {num: nc[key][:] for key, num in PARSIVEL2NC_KEYS.items()}
         # The data logger converts mm/h to m/s, so we need to revert this.
-        data[1] *= 3600 * 1000
+        data[1] *= M_S_TO_MM_H
         # The data logger attempts to convert temperature from °C to K, but this
         # is incorrectly done only for the first value.
         data[12][0] -= 273

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/disdronator/process.py

@@ -4,6 +4,7 @@ import numpy as np
 import numpy.typing as npt
 from numpy import ma
 
+from cloudnetpy import constants
 from cloudnetpy.disdronator.utils import make_rain_mask
 
 
@@ -63,10 +64,9 @@ def process_l2(l1: DisdroL1) -> DisdroL2:
     area_mm2 = (
         l1.area_eff if l1.area_eff is not None else np.repeat(l1.area_nom, n_diameter)
     )
-    area_m2 = area_mm2 * 1e-6
-    diameter_m = l1.diameter / 1000
-    interval_h = interval / 3600
-    rho_w = 1e-6  # kg mm-3
+    area_m2 = area_mm2 * constants.MM2_TO_M2
+    diameter_m = l1.diameter * constants.MM_TO_M
+    interval_h = interval / constants.SEC_IN_HOUR
 
     if spec.ndim == 2:
         number_concentration = spec / (
@@ -81,10 +81,10 @@ def process_l2(l1: DisdroL1) -> DisdroL2:
         energy_flux = (
             np.pi
             / 12
-            * rho_w
+            * constants.RHO_WATER
             * np.sum(
                 spec
-                * l1.diameter**3
+                * diameter_m**3
                 * l1.velocity**2
                 / (area_m2 * interval_h[:, np.newaxis]),
                 axis=1,
@@ -138,10 +138,10 @@ def process_l2(l1: DisdroL1) -> DisdroL2:
         energy_flux = (
             np.pi
             / 12
-            * rho_w
+            * constants.RHO_WATER
             * np.sum(
                 spec_rain
-                * l1.diameter[:, np.newaxis] ** 3
+                * diameter_m[:, np.newaxis] ** 3
                 * l1.velocity**2
                 / (area_m2[:, np.newaxis] * interval_h[:, np.newaxis, np.newaxis]),
                 axis=(1, 2),

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/instruments/fd12p.py

@@ -9,6 +9,7 @@ import numpy as np
 from numpy import ma
 
 from cloudnetpy import output
+from cloudnetpy.constants import MM_H_TO_M_S, MM_TO_M
 from cloudnetpy.exceptions import ValidTimeStampError
 from cloudnetpy.instruments import instruments
 from cloudnetpy.instruments.cloudnet_instrument import CSVFile
@@ -122,11 +123,9 @@ class FD12P(CSVFile):
 
     def convert_units(self) -> None:
         precipitation_rate = self.data["precipitation_rate"][:]
-        self.data["precipitation_rate"].data = (
-            precipitation_rate / 3600 / 1000
-        )  # mm/h -> m/s
+        self.data["precipitation_rate"].data = precipitation_rate * MM_H_TO_M_S
         for key in ("precipitation_amount", "snowfall_amount"):
-            self.data[key].data = self.data[key][:] / 1000  # mm -> m
+            self.data[key].data = self.data[key][:] * MM_TO_M
 
     def screen_all_masked(self) -> None:
         is_valid = np.ones_like(self.data["time"][:], dtype=np.bool_)

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/instruments/hatpro.py

@@ -69,8 +69,8 @@ def hatpro2l1c(
 
     try:
         hatpro_raw = lev1_to_nc(
-            "1C01",
-            str(mwr_dir),
+            data_type="1C01",
+            path_to_files=str(mwr_dir),
             instrument_type=instrument_type,
             output_file=str(output_file),
             lidar_path=lidar_file,

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/instruments/mwr3c.py

@@ -10,6 +10,7 @@ import numpy as np
 from numpy import ma
 
 from cloudnetpy import output, utils
+from cloudnetpy.constants import CM_TO_KG_M2
 from cloudnetpy.datasource import DataSource
 from cloudnetpy.exceptions import ValidTimeStampError
 from cloudnetpy.instruments.arm_utils import concatenate_files, read_geolocation
@@ -17,7 +18,6 @@ from cloudnetpy.instruments.cloudnet_instrument import CloudnetInstrument
 from cloudnetpy.instruments.instruments import MWR3C
 from cloudnetpy.metadata import COMMON_ATTRIBUTES
 
-CM_TO_KG_M2 = 10
 MISSING_VALUE = -9999
 MAX_OFF_ZENITH = 1.0  # deg
 RAW_VARIABLES = (

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/instruments/mwrlos.py

@@ -10,6 +10,7 @@ import numpy as np
 from numpy import ma
 
 from cloudnetpy import output, utils
+from cloudnetpy.constants import CM_TO_KG_M2
 from cloudnetpy.datasource import DataSource
 from cloudnetpy.exceptions import ValidTimeStampError
 from cloudnetpy.instruments.arm_utils import concatenate_files, read_geolocation
@@ -17,7 +18,6 @@ from cloudnetpy.instruments.cloudnet_instrument import CloudnetInstrument
 from cloudnetpy.instruments.instruments import WVR1100
 from cloudnetpy.metadata import COMMON_ATTRIBUTES
 
-CM_TO_KG_M2 = 10
 RAW_VARIABLES = ("time", "liq", "vap", "qc_liq", "qc_vap", "wet_window")
 
 

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/instruments/radiometrics.py

@@ -18,6 +18,7 @@ from numpy import ma
 
 from cloudnetpy import output, utils
 from cloudnetpy.cloudnetarray import CloudnetArray
+from cloudnetpy.constants import CM_TO_KG_M2, G_TO_KG
 from cloudnetpy.exceptions import InconsistentDataError, ValidTimeStampError
 from cloudnetpy.instruments import instruments
 from cloudnetpy.metadata import MetaData
@@ -236,8 +237,8 @@ class RadiometricsMP:
             _process_superblock()
 
         self.data["time"] = np.array(times, dtype="datetime64[s]")
-        self.data["lwp"] = np.array(lwps)  # mm => kg m-2
-        self.data["iwv"] = np.array(iwvs) * 10  # cm => kg m-2
+        self.data["lwp"] = np.array(lwps)  # mm = kg m-2
+        self.data["iwv"] = np.array(iwvs) * CM_TO_KG_M2
         if irt_times:
             self.data["irt"] = _find_closest(
                 np.array(irt_times, dtype="datetime64[s]"),
@@ -259,7 +260,7 @@ class RadiometricsMP:
         if ah_times:
             self.data["absolute_humidity"] = _find_closest(
                 np.array(ah_times, dtype="datetime64[s]"),
-                np.array(ahs) / 1000,  # g m-3 => kg m-3
+                np.array(ahs) * G_TO_KG,  # g m-3 => kg m-3
                 self.data["time"],
             )
 
@@ -312,8 +313,8 @@ class RadiometricsWVR:
             ],
             dtype="datetime64[s]",
         )
-        self.data["lwp"] = np.array(self.raw_data["LiqCM"]) * 10  # cm => kg m-2
-        self.data["iwv"] = np.array(self.raw_data["VapCM"]) * 10  # cm => kg m-2
+        self.data["lwp"] = np.array(self.raw_data["LiqCM"]) * CM_TO_KG_M2
+        self.data["iwv"] = np.array(self.raw_data["VapCM"]) * CM_TO_KG_M2
         is_zenith = np.abs(np.array(self.raw_data["ELact"]) - 90.0) < 1.0
         tb23_valid = np.array(self.raw_data["TbSky23"]) > 0
         tb31_valid = np.array(self.raw_data["TbSky31"]) > 0

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/instruments/rain_e_h3.py

@@ -6,6 +6,7 @@ from uuid import UUID
 import numpy as np
 
 from cloudnetpy import output
+from cloudnetpy.constants import MM_H_TO_M_S, MM_TO_M
 from cloudnetpy.exceptions import ValidTimeStampError
 from cloudnetpy.instruments import instruments
 from cloudnetpy.instruments.cloudnet_instrument import CSVFile
@@ -165,7 +166,5 @@ class RainEH3(CSVFile):
 
     def convert_units(self) -> None:
         rainfall_rate = self.data["rainfall_rate"][:]
-        self.data["rainfall_rate"].data = rainfall_rate / 3600 / 1000  # mm/h -> m/s
-        self.data["rainfall_amount"].data = (
-            self.data["rainfall_amount"][:] / 1000
-        )  # mm -> m
+        self.data["rainfall_rate"].data = rainfall_rate * MM_H_TO_M_S
+        self.data["rainfall_amount"].data = self.data["rainfall_amount"][:] * MM_TO_M

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/instruments/weather_station.py

@@ -14,7 +14,13 @@ from numpy import ma
 
 from cloudnetpy import output
 from cloudnetpy.cloudnetarray import CloudnetArray
-from cloudnetpy.constants import HPA_TO_PA, MM_H_TO_M_S, SEC_IN_HOUR
+from cloudnetpy.constants import (
+    HPA_TO_PA,
+    MM_H_TO_M_S,
+    MM_TO_M,
+    SEC_IN_HOUR,
+    SEC_IN_MINUTE,
+)
 from cloudnetpy.exceptions import ValidTimeStampError
 from cloudnetpy.instruments import instruments
 from cloudnetpy.instruments.cloudnet_instrument import CSVFile
@@ -147,7 +153,9 @@ class WS(CSVFile):
         if "rainfall_rate" not in self.data:
             return
         rainfall_rate = self.data["rainfall_rate"][:]
-        self.data["rainfall_rate"].data = rainfall_rate / 60 / 1000  # mm/min -> m/s
+        self.data["rainfall_rate"].data = (
+            rainfall_rate * MM_TO_M / SEC_IN_MINUTE
+        )  # mm/min -> m/s
 
     def convert_pressure(self) -> None:
         if "air_pressure" not in self.data:
@@ -249,9 +257,7 @@ class PalaiseauWS(WS):
             self.data[key] = CloudnetArray(array_masked, key)
 
     def convert_rainfall_amount(self) -> None:
-        self.data["rainfall_amount"].data = (
-            self.data["rainfall_amount"][:] / 1000
-        )  # mm -> m
+        self.data["rainfall_amount"].data = self.data["rainfall_amount"][:] * MM_TO_M
 
     def _validate_header(self, header: list[str]) -> None:
         column_titles = [row for row in header if "Col." in row]
@@ -710,7 +716,7 @@ class LimassolWS(WS):
     def convert_rainfall_rate(self) -> None:
         rainfall_rate = self.data["rainfall_rate"][:]
         self.data["rainfall_rate"].data = (
-            rainfall_rate / (10 * 60) / 1000
+            rainfall_rate * MM_TO_M / (10 * SEC_IN_MINUTE)
         )  # mm/(10 min) -> m/s
 
 

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/model_evaluation/plotting/plotting.py

@@ -18,6 +18,7 @@ from mpl_toolkits.axes_grid1 import make_axes_locatable
 from numpy import ma
 
 import cloudnetpy.model_evaluation.plotting.plot_tools as p_tools
+from cloudnetpy.constants import T_FREEZING
 from cloudnetpy.model_evaluation.plotting.plot_meta import ATTRIBUTES, PlotMeta
 from cloudnetpy.model_evaluation.statistics.statistical_methods import DayStatistics
 from cloudnetpy.plotting.plotting import (
@@ -205,7 +206,7 @@ def plot_colormesh(
         data = data.T
     vmin, vmax = variable_info.plot_range
     if variable_info.clabel == "\u00b0C":
-        data = data - 273.15
+        data = data - T_FREEZING
     if variable_info.plot_scale == "logarithmic":
         data, vmin, vmax = lin2log(data, vmin, vmax)
     cmap = plt.get_cmap(variable_info.cbar, 22)

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/model_evaluation/products/advance_methods.py

@@ -5,6 +5,7 @@ import scipy.stats
 from numpy import ma
 
 import cloudnetpy.utils as cl_tools
+from cloudnetpy.constants import T_FREEZING
 from cloudnetpy.model_evaluation.model_metadata import MODEL_PREFIX
 from cloudnetpy.model_evaluation.products.model_products import ModelManager
 from cloudnetpy.model_evaluation.products.observation_products import ObservationManager
@@ -49,7 +50,7 @@ class AdvanceProductMethods:
         cf = self.getvar_from_object("cf")
         h = self.getvar_from_object("h")
         temperature = self._model_obj.getvar("temperature")
-        t_screened = self.remove_extra_levels(temperature - 273.15)
+        t_screened = self.remove_extra_levels(temperature - T_FREEZING)
         iwc, lwc = (self._model_obj.get_water_content(var) for var in ["iwc", "lwc"])
         coeffs = self.set_frequency_parameters()
         z_sen = self.fit_z_sensitivity(h)

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/model_evaluation/products/model_products.py

@@ -3,9 +3,9 @@ from os import PathLike
 
 import numpy as np
 import numpy.typing as npt
+from atmoslib.constants import RS
 from numpy import ma
 
-from cloudnetpy.constants import RS
 from cloudnetpy.datasource import DataSource
 from cloudnetpy.exceptions import ModelDataError
 from cloudnetpy.model_evaluation.model_metadata import (

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/model_evaluation/tests/unit/test_model_products.py

@@ -3,7 +3,7 @@ import numpy as np
 import pytest
 from numpy import testing
 
-from cloudnetpy.constants import RS
+from atmoslib.constants import RS
 from cloudnetpy.exceptions import ModelDataError
 from cloudnetpy.model_evaluation.products.model_products import ModelManager
 

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/plotting/plotting.py

@@ -12,6 +12,7 @@ import matplotlib.pyplot as plt
 import netCDF4
 import numpy as np
 import numpy.typing as npt
+from atmoslib.constants import G
 from matplotlib import rcParams
 from matplotlib.axes import Axes
 from matplotlib.colorbar import Colorbar
@@ -264,7 +265,7 @@ class FigureData:
             "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
+            return np.mean(self.file.variables["sfc_geopotential"][:]) / G
         pressure = ma.mean(self.file.variables["pressure"][:, 0])
         temperature = ma.mean(self.file.variables["temperature"][:, 0])
         return float(atmoslib.isa_altitude(temperature, pressure))
@@ -455,17 +456,17 @@ class Plot:
 
     def _convert_units(self) -> str:
         multiply, add = "multiply", "add"
-        k_to_c = (add, -273.15, "\u00b0C")
+        k_to_c = (add, -con.T_FREEZING, "\u00b0C")
         units_conversion = {
-            "rainfall_rate": (multiply, 3600000, "mm h$^{-1}$"),
-            "snowfall_rate": (multiply, 3600000, "mm h$^{-1}$"),
-            "precipitation_rate": (multiply, 3600000, "mm h$^{-1}$"),
+            "rainfall_rate": (multiply, con.M_S_TO_MM_H, "mm h$^{-1}$"),
+            "snowfall_rate": (multiply, con.M_S_TO_MM_H, "mm h$^{-1}$"),
+            "precipitation_rate": (multiply, con.M_S_TO_MM_H, "mm h$^{-1}$"),
             "air_pressure": (multiply, 0.01, "hPa"),
             "relative_humidity": (multiply, 100, "%"),
             "rh": (multiply, 100, "%"),
-            "rainfall_amount": (multiply, 1000, "mm"),
-            "snowfall_amount": (multiply, 1000, "mm"),
-            "precipitation_amount": (multiply, 1000, "mm"),
+            "rainfall_amount": (multiply, con.M_TO_MM, "mm"),
+            "snowfall_amount": (multiply, con.M_TO_MM, "mm"),
+            "precipitation_amount": (multiply, con.M_TO_MM, "mm"),
             "air_temperature": k_to_c,
             "temperature": k_to_c,
             "Tw": k_to_c,
@@ -476,8 +477,8 @@ class Plot:
             "transmitter_temperature": k_to_c,
             "receiver_temperature": k_to_c,
             "pc_temperature": k_to_c,
-            "r_accum_RT": (multiply, 1000, "mm"),
-            "r_accum_NRT": (multiply, 1000, "mm"),
+            "r_accum_RT": (multiply, con.M_TO_MM, "mm"),
+            "r_accum_NRT": (multiply, con.M_TO_MM, "mm"),
             "cloud_top_height_agl": (multiply, con.M_TO_KM, "Height (km AGL)"),
         }
         conversion_method, conversion, units = units_conversion.get(

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/products/cod.py

@@ -16,9 +16,6 @@ from cloudnetpy.products.iwc import IwcSource
 from cloudnetpy.products.lwc import CloudAdjustor, Lwc, LwcError, LwcSource
 from cloudnetpy.products.product_tools import IceClassification
 
-# Density of liquid water (kg m-3)
-RHO_WATER = 1000
-
 # Droplet effective radius (m) assumed when it cannot be retrieved from radar
 DEFAULT_ASSUMED_DER = 10e-6
 
@@ -145,7 +142,7 @@ class OpticalDepthSource(DataSource):
         der_rel_error_filled = ma.filled(der_rel_error, ASSUMED_DER_REL_ERROR).copy()
         der_rel_error_filled[self.is_der_assumed] = ASSUMED_DER_REL_ERROR
         # Includes pixels added at lidar-only cloud tops by the lwc retrieval
-        extinction = 3 * lwc / (2 * RHO_WATER * der_filled)
+        extinction = 3 * lwc / (2 * constants.RHO_WATER * der_filled)
         rel_error = utils.l2norm(ma.filled(lwc_rel_error, 0), der_rel_error_filled)
         self._rel_error["liquid"] = ma.masked_where(
             ma.getmaskarray(extinction), rel_error

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/products/der.py

@@ -10,7 +10,7 @@ import numpy as np
 import numpy.typing as npt
 from numpy import ma
 
-from cloudnetpy import output, utils
+from cloudnetpy import constants, output, utils
 from cloudnetpy.categorize import atmos_utils
 from cloudnetpy.datasource import DataSource
 from cloudnetpy.exceptions import InvalidSourceFileError
@@ -131,8 +131,6 @@ class DerSource(DataSource):
     def append_der(self) -> None:
         """Estimate liquid droplet effective radius using Frisch et al. 2002."""
         params = self.parameters
-        rho_l = 1000  # density of liquid water (kg m-3)
-
         var_x = params.sigma_x * params.sigma_x
 
         Z = self.getvar("Z")
@@ -181,9 +179,13 @@ class DerSource(DataSource):
                 A * A + B * B + C * C,
             )  # mm
 
+            # Scaled retrieval needs positive lwp
+            if lwp[ind_t] <= 0:
+                continue
+
             # der scaled formula (6)
             A = Z[ind_t, idx_layer] ** (1 / 6) / (2 * lwp[ind_t] ** (1 / 3))
-            B = (np.pi * rho_l / 6) ** (1 / 3)
+            B = (np.pi * constants.RHO_WATER / 6) ** (1 / 3)
             C = integral ** (1 / 3) * ma.exp(-2 * var_x)
             der_scaled[ind_t, idx_layer] = 1.0e-3 * A * B * C  # μm => mm
 

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/products/drizzle.py

@@ -8,7 +8,7 @@ import numpy.typing as npt
 from numpy import ma
 from scipy.special import gamma
 
-from cloudnetpy import output, utils
+from cloudnetpy import constants, output, utils
 from cloudnetpy.metadata import MetaData
 from cloudnetpy.products.drizzle_error import get_drizzle_error
 from cloudnetpy.products.drizzle_tools import (
@@ -122,13 +122,12 @@ class DrizzleProducts:
 
     def _calc_lwc(self) -> npt.NDArray:
         """Calculates drizzle liquid water content (kg m-3)."""
-        rho_water = 1000
         dia, mu, s = (self._params.get(key) for key in ("Do", "mu", "S"))
         dia = ma.array(dia)
         mu = ma.array(mu)
         s = ma.array(s)
         gamma_ratio = gamma(4 + mu) / gamma(3 + mu) / (3.67 + mu)
-        return rho_water / 3 * self._data.beta * s * dia * gamma_ratio
+        return constants.RHO_WATER / 3 * self._data.beta * s * dia * gamma_ratio
 
     def _calc_lwf(self, lwc_in: npt.NDArray) -> npt.NDArray:
         """Calculates drizzle liquid water flux."""

{cloudnetpy-1.96.1 → cloudnetpy-1.96.2}/cloudnetpy/products/drizzle_error.py

@@ -3,6 +3,7 @@ import numpy.typing as npt
 from numpy import ma
 
 from cloudnetpy import utils
+from cloudnetpy.constants import LN_TO_DB
 from cloudnetpy.products.drizzle_tools import DrizzleSolver, DrizzleSource
 
 MU_ERROR = 0.07
@@ -172,18 +173,15 @@ def _stack_errors(
     return error
 
 
-COR = 10 / np.log(10)
-
-
 def db2lin(x_in: npt.NDArray) -> ma.MaskedArray:
     x = ma.copy(x_in)
     threshold = 100
     x[x > threshold] = threshold
-    return ma.exp(x / COR) - 1
+    return ma.exp(x / LN_TO_DB) - 1
 
 
 def lin2db(x_in: npt.NDArray) -> ma.MaskedArray:
     x = ma.copy(x_in)
     threshold = -0.9
     x[x < threshold] = threshold
-    return ma.log(x + 1) * COR
+    return ma.log(x + 1) * LN_TO_DB