cloudnetpy 1.65.7__py3-none-any.whl → 1.66.0__py3-none-any.whl

{cloudnetpy-1.65.7.dist-info → cloudnetpy-1.66.0.dist-info}/RECORD

@@ -1,30 +1,36 @@
 cloudnetpy/__init__.py,sha256=X_FqY-4yg5GUj5Edo14SToLEos6JIsC3fN-v1FUgQoA,43
 cloudnetpy/cloudnetarray.py,sha256=Ol1ha4RPAmFZANL__U5CaMKX4oYMXYR6OnjoCZ9w3eo,7077
 cloudnetpy/concat_lib.py,sha256=8Ek059RMLAfbbXCkX90cgnhw_8ZpcDrxw1yPvwtuitU,9846
-cloudnetpy/constants.py,sha256=Kn95HXVMsKFo9Bu-UOYSBJuPEOjBR4JLaEXH0-MBpwk,830
+cloudnetpy/constants.py,sha256=RDB9aqpBRztk3QVCFgsmi9fwhtLuit_0WJrt0D6sDcc,736
 cloudnetpy/datasource.py,sha256=j7N4g59HPvOBWle-W9bOUF0BfRLgvR4zwOi_B50cI7Q,7921
 cloudnetpy/exceptions.py,sha256=ns48useL9RN3mPh7CqIiLA19VI9OmVbyRsKTkwbThF8,1760
 cloudnetpy/metadata.py,sha256=v_VDo2vbdTxB0zIsfP69IcrwSKiRlLpsGdq6JPI4CoA,5306
-cloudnetpy/output.py,sha256=YkCaxVkG_Mt2hng_IVnhygHteV4UMKzKALkeFZwFJL8,14822
+cloudnetpy/output.py,sha256=YrWRBEZg0QNZRVnd9ziAziH-eJSh7O5JuWiH4ZxM0_s,15584
 cloudnetpy/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-cloudnetpy/utils.py,sha256=JV0Fawnme1HoZgoiidV3eIzsn6vx0AEjBNmI1CcrBsA,28517
-cloudnetpy/version.py,sha256=B0fKgZkDpneBNddUtrKe13FSc_Y17mGFfki5_swqhMU,72
-cloudnetpy/categorize/__init__.py,sha256=gP5q3Vis1y9u9OWgA_idlbjfWXYN_S0IBSWdwBhL_uU,69
-cloudnetpy/categorize/atmos.py,sha256=vavMC_WQQyGH14eL4vAzKLKDDZt8tBrMYimztYHOjH8,12639
-cloudnetpy/categorize/atmos_utils.py,sha256=ns5ydiEN34Ng6mJiOBpxKBVDU2NXj6W3Q5IUynmNRYI,3586
-cloudnetpy/categorize/categorize.py,sha256=VcxSs22icYDG3R5x9jTX7uuSjPU8tyc5bIz0Zj5F9z0,18303
-cloudnetpy/categorize/classify.py,sha256=a-0bVCtynGfORnDGTsPuzqkuDeOOR_OMz5ai9NsMuic,9870
-cloudnetpy/categorize/containers.py,sha256=g3SQHoqlY1uJ8b-MG-BbW3oWz5IyacA8kJBeIPy_4EA,4859
+cloudnetpy/utils.py,sha256=1gne8Qqa-e4CNNJMynWr4LVirGK0Be-KcskVZMD72ac,28621
+cloudnetpy/version.py,sha256=M2-atDdtYzbMUo8WvBwyYGrpfWspuimBvKUAoMUh9xk,72
+cloudnetpy/categorize/__init__.py,sha256=s-SJaysvVpVVo5kidiruWQO6p3gv2TXwY1wEHYO5D6I,44
+cloudnetpy/categorize/atmos_utils.py,sha256=9-ymI6i1xASf-XAFyO87FaTfvq6bF89N1i_27OkUp-M,10104
+cloudnetpy/categorize/attenuation.py,sha256=Y_-fzmQTltWTqIZTulJhovC7a6ifpMcaAazDJcnMIOc,990
+cloudnetpy/categorize/categorize.py,sha256=UQdgLHTNc-J4-CRK3P1mrGNUMal2cKE6nwE5JYlFMcc,20644
+cloudnetpy/categorize/classify.py,sha256=qovHgHsMku5kpl3cJxKteNBsG8GAkfI3Zo8QhJwZSFQ,8512
+cloudnetpy/categorize/containers.py,sha256=4AUx95Br-aczyUB1dCmTfLn0MRKsk2GLkdCTdnqIfAs,5186
 cloudnetpy/categorize/disdrometer.py,sha256=keU3pFvKtk840A0oLwAyNDuqOCswBPJEKf2bV0YWyA8,2004
-cloudnetpy/categorize/droplet.py,sha256=s-q1nCZ40rw5-fabfj61mH2lmdcNRT-TASvwWGZeptY,8688
-cloudnetpy/categorize/falling.py,sha256=Wz49mbw0pubnbhoekYTqGT0S9UNaE88jWOjSPEsCIaI,4386
-cloudnetpy/categorize/freezing.py,sha256=c4k5AIgfBpvw64EaVVVYPi2Fx4SpHk1cyfceE1ydD28,3755
+cloudnetpy/categorize/droplet.py,sha256=t49KEsH5ZM68JQ4NvAf9kGgQ-evic1T4de2-jgJ2f4M,8683
+cloudnetpy/categorize/falling.py,sha256=lok0HMi1ewf9pS70mq62nRKL6wJzMyWbYmv1cdwrwnA,4404
+cloudnetpy/categorize/freezing.py,sha256=eSFD37R7vBrg7mgfSanrwhBjnFyWNBpjw2AtvRmSh48,3753
 cloudnetpy/categorize/insects.py,sha256=9J5agmktit8Or66GGNue-bThiaG9rB2SuPNZBXI7FCE,5243
+cloudnetpy/categorize/itu.py,sha256=fnyPE4FreZhwXHJbOPTlHKeG9hTXuJZIKt_VBUQqiJo,10234
 cloudnetpy/categorize/lidar.py,sha256=YQrM_LOz8NQrrD9l9HyujV1GSGwkQ8LMqXN13bEJRW4,2605
-cloudnetpy/categorize/melting.py,sha256=mYdOKxfTC2InB8NdOPwr_7NpbouQMm-9f2Q1kfTqIJE,6262
-cloudnetpy/categorize/model.py,sha256=rLyx9SKz56BmQ809PwwhCnbusshUuYPoqhLJbnZvrt8,5842
+cloudnetpy/categorize/melting.py,sha256=Oi80SrB4pNBhZlFhb6ayXkT2KMtXXV2rYu6jzEGn7tg,6211
+cloudnetpy/categorize/model.py,sha256=stNfirRtMLPPH-9hydgS33aXeXeE50Cv_ki3BTjFApo,6630
 cloudnetpy/categorize/mwr.py,sha256=F7cquERWL6mBkgboqeaCIPf9gOlKI-NWUQIBdQXGT_I,1635
-cloudnetpy/categorize/radar.py,sha256=NCBjyC-OvkcprDnaDD0ql3WoKtSWMHcMhspTh5Zlmi8,13730
+cloudnetpy/categorize/radar.py,sha256=0Wg5u2aLXRVhHiFiXb4fSqY_iGgSwcCMMVshM3wBoo0,14149
+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/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
 cloudnetpy/instruments/campbell_scientific.py,sha256=2WHfBKQjtRSl0AqvtPeX7G8Hdi3Dn0WbvoAppFOMbA8,5270
@@ -53,7 +59,7 @@ cloudnetpy/instruments/disdrometer/common.py,sha256=g52iK2aNp3Z88kovUmGVpC54NZom
 cloudnetpy/instruments/disdrometer/parsivel.py,sha256=HJZrEysQkx9MiIVPDV25CYHpXi_SjgZlgO-otoaKK34,25640
 cloudnetpy/instruments/disdrometer/thies.py,sha256=JYFREICrf6nas3bJeEnJUtWzDc_4mEmt3TxhqO--bfQ,10632
 cloudnetpy/model_evaluation/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-cloudnetpy/model_evaluation/file_handler.py,sha256=0Pg7V6eo1G6KYeE10B6lR8hIKIJoseO-pBaorcW42rY,6446
+cloudnetpy/model_evaluation/file_handler.py,sha256=jnpTxwxLSZdt08By0cjFDmq1KWC4fhJaN_MfjgA8ijQ,6440
 cloudnetpy/model_evaluation/metadata.py,sha256=PAsnOqcUoV3CJwplgWiVCF9Zt4io8tqj7CIgS4fEL-8,8412
 cloudnetpy/model_evaluation/model_metadata.py,sha256=CxpY6RPm7GOTBBmPhcNVVpm9ateUmHSUwGtFXTLq3To,1436
 cloudnetpy/model_evaluation/utils.py,sha256=Z9VqYVdtY9yTr2JeVfBn4nccIVWCN5Fd-BCyB_qYI-A,154
@@ -65,7 +71,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=vrbHT008T4RFXM2pKm7dWPLKb1smD4rUNOM6IhUcU_w,5500
+cloudnetpy/model_evaluation/products/observation_products.py,sha256=UJTb6-1Z-HOB9ltW8E7QHNEoByXeDQQUNi3GyS5lG_0,5500
 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
@@ -85,31 +91,31 @@ cloudnetpy/model_evaluation/tests/e2e/process_lwc/tests.py,sha256=ANBA0LVao3Xrm-
 cloudnetpy/model_evaluation/tests/unit/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 cloudnetpy/model_evaluation/tests/unit/conftest.py,sha256=WL_FgrDeoUYGp4PKjb37HLu79D9uu33PGQL40_ctqS0,7446
 cloudnetpy/model_evaluation/tests/unit/test_advance_methods.py,sha256=IkoAVtsWVFrPpFqQOLAPHKb9qgV-KjGGVEtWMudeiSo,10079
-cloudnetpy/model_evaluation/tests/unit/test_grid_methods.py,sha256=qkNPfHI25EE0-BXk73TijpeM7YwswX7e41-764o5lqE,26254
+cloudnetpy/model_evaluation/tests/unit/test_grid_methods.py,sha256=AEUXN5HBhKliPsSNGDCUtWOtIx6Y8iqkywb_-RfXYU0,26277
 cloudnetpy/model_evaluation/tests/unit/test_model_products.py,sha256=FRbYLshSHH2E527uJPwvUIyZKTsPFSZrwDsPsNrFSSU,3475
-cloudnetpy/model_evaluation/tests/unit/test_observation_products.py,sha256=P-W5QwRHMtem6p5SyyH7p9TvHGro3XW1baQcIwh6nFg,4892
+cloudnetpy/model_evaluation/tests/unit/test_observation_products.py,sha256=DN3yVqq8vFYca_9POjcrJ8XViMrJks_jM-aQznfN8QQ,4936
 cloudnetpy/model_evaluation/tests/unit/test_plot_tools.py,sha256=POdypGWjV2NA4DCU7w8Unk_IdPfOpUb1qBDhfA3x1Bw,9222
 cloudnetpy/model_evaluation/tests/unit/test_plotting.py,sha256=h9V8JKmrO4v9bOvv-UjRa06sZJQPhDNVHGBSImDdtkI,3277
 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=JHrr-4A9fhqdi_tQFe6mR4Fdry3hkI-lmmVu5Ny2vco,15979
-cloudnetpy/plotting/plotting.py,sha256=JIgia-Hujsa2ot4JMoBLzEdH7YWzmlGuKSLvvQ9tl0U,34253
+cloudnetpy/plotting/plot_meta.py,sha256=ZvaKU3eXy1KFxQomnsEu3mCYpwwBYKAYk7oAwOzAGSg,16143
+cloudnetpy/plotting/plotting.py,sha256=ML7v5iUlIBfeXX-F55Cf0GiHfDErblUWhfAXOCZZb3k,35141
 cloudnetpy/products/__init__.py,sha256=2hRb5HG9hNrxH1if5laJkLeFeaZCd5W1q3hh4ewsX0E,273
-cloudnetpy/products/classification.py,sha256=bNG8W1CMgGoUBpXopQjYAW3F-uEJGyojXb4A5jmErHo,7921
-cloudnetpy/products/der.py,sha256=1LDBbnbUg8feMUTGWJq3bObBhEcZ_Ee17MB1x0GwRdo,12669
+cloudnetpy/products/classification.py,sha256=AKb9GCatvhS5KR0c9LfN96nUvzi02175ZCQlvMH1Dws,8077
+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=LR2AtbFQRGFrJ2LGyiLxOfbnlznVLydXvb8RFDR0_4E,10848
-cloudnetpy/products/ier.py,sha256=ge1f_aYick20Nlznq8zbBl5umWlTP-UwMivy4Y05Sck,7839
-cloudnetpy/products/iwc.py,sha256=Q8dXV3JF3JUQgwkmQFOKakm21Tnf8oCWsH0CSqIEKl4,10209
-cloudnetpy/products/lwc.py,sha256=wSd3GDqELz4yyWBMiKDR-QhRK8scPheqsBcS1qzxYOI,18997
+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/lwc.py,sha256=sl6Al2tuH3KkCBrPbWTmuz3jlD5UQJ4D6qBsn1tt2CQ,18962
 cloudnetpy/products/mie_lu_tables.nc,sha256=It4fYpqJXlqOgL8jeZ-PxGzP08PMrELIDVe55y9ob58,16637951
-cloudnetpy/products/mwr_tools.py,sha256=tN_sPDS3BdpFJfa5a2mnc3eCMoi7syjVJMaTt962hmo,5004
-cloudnetpy/products/product_tools.py,sha256=VNw2diJj30POz68-3qNVkJP7r9AUspT_d1Fp0BbeIx8,10414
+cloudnetpy/products/mwr_tools.py,sha256=rd7UC67O4fsIE5SaHVZ4qWvUJTj41ZGwgQWPwZzOM14,5377
+cloudnetpy/products/product_tools.py,sha256=01Zc6xV8CSuYcIcLpchFf5POL3_c629-YMNDZJ51udA,10853
 docs/source/conf.py,sha256=IKiFWw6xhUd8NrCg0q7l596Ck1d61XWeVjIFHVSG9Og,1490
-cloudnetpy-1.65.7.dist-info/LICENSE,sha256=wcZF72bdaoG9XugpyE95Juo7lBQOwLuTKBOhhtANZMM,1094
-cloudnetpy-1.65.7.dist-info/METADATA,sha256=KHoz5LFKsFR3kBEdosbY1L_jC15RvnJCjnotAaTQF0w,5784
-cloudnetpy-1.65.7.dist-info/WHEEL,sha256=cVxcB9AmuTcXqmwrtPhNK88dr7IR_b6qagTj0UvIEbY,91
-cloudnetpy-1.65.7.dist-info/top_level.txt,sha256=ibSPWRr6ojS1i11rtBFz2_gkIe68mggj7aeswYfaOo0,16
-cloudnetpy-1.65.7.dist-info/RECORD,,
+cloudnetpy-1.66.0.dist-info/LICENSE,sha256=wcZF72bdaoG9XugpyE95Juo7lBQOwLuTKBOhhtANZMM,1094
+cloudnetpy-1.66.0.dist-info/METADATA,sha256=EwNyFbZwEjCO1ykTbgOgruQsjgKlbJNAPU7O9Myg8-w,5784
+cloudnetpy-1.66.0.dist-info/WHEEL,sha256=GV9aMThwP_4oNCtvEC2ec3qUYutgWeAzklro_0m4WJQ,91
+cloudnetpy-1.66.0.dist-info/top_level.txt,sha256=ibSPWRr6ojS1i11rtBFz2_gkIe68mggj7aeswYfaOo0,16
+cloudnetpy-1.66.0.dist-info/RECORD,,

{cloudnetpy-1.65.7.dist-info → cloudnetpy-1.66.0.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (74.1.2)
+Generator: setuptools (75.1.0)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

cloudnetpy/categorize/atmos.py

@@ -1,376 +0,0 @@
-"""This module contains functions to calculate
-various atmospheric parameters.
-"""
-
-from typing import Final
-
-import numpy as np
-import scipy.constants
-from numpy import ma
-
-from cloudnetpy import constants as con
-from cloudnetpy import utils
-from cloudnetpy.categorize.atmos_utils import (
-    calc_psychrometric_constant,
-    calc_saturation_vapor_pressure,
-)
-from cloudnetpy.categorize.containers import ClassificationResult
-from cloudnetpy.categorize.model import Model
-
-M_TO_KM: Final = 0.001
-TWO_WAY: Final = 2
-
-
-def calc_lwc_change_rate(temperature: np.ndarray, pressure: np.ndarray) -> np.ndarray:
-    """Returns rate of change of condensable water (LWC).
-
-    Calculates the theoretical adiabatic rate of increase of LWC
-    with height, given the cloud base temperature and pressure.
-
-    Args:
-        temperature: Temperature of cloud base (K).
-        pressure: Pressure of cloud base (Pa).
-
-    Returns:
-        dlwc/dz (kg m-3 m-1)
-
-    References:
-        Brenguier, 1991, https://bit.ly/2QCSJtb
-
-    """
-    svp = calc_saturation_vapor_pressure(temperature)
-    svp_mixing_ratio = calc_mixing_ratio(svp, pressure)
-    air_density = calc_air_density(pressure, temperature, svp_mixing_ratio)
-    kelvin_per_kg = calc_psychrometric_constant(temperature)
-    pressure_difference = pressure - svp
-    f1 = kelvin_per_kg - 1
-    f2 = 1 / (
-        kelvin_per_kg
-        + (con.LATENT_HEAT * svp_mixing_ratio * air_density / pressure_difference)
-    )
-    f3 = con.MW_RATIO * svp * pressure_difference**-2
-    dqs_dp = f1 * f2 * f3
-    return dqs_dp * air_density**2 * -scipy.constants.g
-
-
-def calc_mixing_ratio(svp: np.ndarray, pressure: np.ndarray) -> np.ndarray:
-    """Calculates mixing ratio from saturation vapor pressure and pressure.
-
-    Args:
-        svp: Saturation vapor pressure (Pa).
-        pressure: Atmospheric pressure (Pa).
-
-    Returns:
-        Mixing ratio (kg kg-1).
-
-    """
-    return con.MW_RATIO * svp / (pressure - svp)
-
-
-def calc_air_density(
-    pressure: np.ndarray,
-    temperature: np.ndarray,
-    svp_mixing_ratio: np.ndarray,
-) -> np.ndarray:
-    """Calculates air density (kg m-3).
-
-    Args:
-        pressure: Pressure (Pa).
-        temperature: Temperature (K).
-        svp_mixing_ratio: Saturation vapor pressure mixing ratio (kg/kg).
-
-    Returns:
-        Air density (kg m-3).
-
-    """
-    return pressure / (con.RS * temperature * (0.6 * svp_mixing_ratio + 1))
-
-
-def get_attenuations(data: dict, classification: ClassificationResult) -> dict:
-    """Calculates attenuations due to atmospheric gases and liquid water.
-
-    Args:
-        data: Containing :class:`Model` and :class:`Mwr` instances.
-        classification: A :class:`ClassificationResult` instance.
-
-    Returns:
-        Dictionary containing `radar_gas_atten`, `radar_liquid_atten`,
-            `liquid_atten_err`, `liquid_corrected` and `liquid_uncorrected` fields.
-
-    """
-    gas = GasAttenuation(data, classification)
-    liquid = LiquidAttenuation(data, classification)
-    return {
-        "radar_gas_atten": gas.atten,
-        "radar_liquid_atten": liquid.atten,
-        "liquid_atten_err": liquid.atten_err,
-        "liquid_corrected": liquid.corrected,
-        "liquid_uncorrected": liquid.uncorrected,
-    }
-
-
-class Attenuation:
-    """Base class for gas and liquid attenuations.
-
-    Args:
-        model: The :class:`Model` instance.
-        classification: The :class:`ClassificationResult` instance.
-
-    Attributes:
-        classification (ClassificationResult): The :class:`ClassificationResult`
-        instance.
-
-    """
-
-    def __init__(self, model: Model, classification: ClassificationResult):
-        self._dheight = utils.mdiff(model.height)
-        self._model = model.data_dense
-        self._liquid_in_pixel = utils.isbit(classification.category_bits, 0)
-        self.classification = classification
-
-
-class GasAttenuation(Attenuation):
-    """Radar gas attenuation class. Child of Attenuation.
-
-    Args:
-        data: Containing :class:`Model` instance.
-        classification: The :class:`ClassificationResult` instance.
-
-    Attributes:
-        atten (ndarray): Gas attenuation (dB).
-
-    """
-
-    def __init__(self, data: dict, classification: ClassificationResult):
-        super().__init__(data["model"], classification)
-        self.atten = self._calc_gas_atten()
-
-    def _calc_gas_atten(self) -> np.ndarray:
-        specific_atten = ma.copy(self._model["specific_gas_atten"])
-        specific_atten_corrected = self._fix_atten_in_liquid(specific_atten)
-        return self._specific_to_gas_atten(specific_atten_corrected)
-
-    def _fix_atten_in_liquid(self, atten: np.ndarray) -> np.ndarray:
-        saturated_atten = self._model["specific_saturated_gas_atten"]
-        atten[self._liquid_in_pixel] = saturated_atten[self._liquid_in_pixel]
-        return atten
-
-    def _specific_to_gas_atten(self, specific_atten: np.ndarray) -> np.ndarray:
-        layer1_atten = self._model["gas_atten"][:, 0]
-        atten_cumsum = ma.cumsum(specific_atten, axis=1)
-        atten = TWO_WAY * atten_cumsum * self._dheight * M_TO_KM
-        atten += utils.transpose(layer1_atten)
-        atten = np.insert(atten, 0, layer1_atten, axis=1)[:, :-1]
-        return ma.array(atten, mask=atten_cumsum.mask)
-
-
-class LiquidAttenuation(Attenuation):
-    """Radar liquid attenuation class. Child of Attenuation.
-
-    Args:
-        data: Containing :class:`Model` and :class:`Mwr` instances.
-        classification: The :class:`ClassificationResult` instance.
-
-    Attributes:
-        atten (ndarray): Radar liquid attenuation (dB).
-        atten_err (ndarray): Error of radar liquid attenuation (dB).
-        uncorrected (ndarray): Boolean array denoting uncorrected pixels.
-        corrected (ndarray): Boolean array denoting corrected pixels.
-
-    """
-
-    def __init__(self, data: dict, classification: ClassificationResult):
-        super().__init__(data["model"], classification)
-        if data["mwr"] is not None:
-            self._mwr = data["mwr"].data
-        else:
-            n_time = data["radar"].data["Z"][:].shape[0]
-            self._mwr = {
-                "lwp": ma.masked_all(n_time),
-                "lwp_error": ma.masked_all(n_time),
-            }
-        self._lwc_dz_err = self._get_lwc_change_rate_error()
-        self.atten = self._get_liquid_atten()
-        self.atten_err = self._get_liquid_atten_err()
-        self.uncorrected = self._find_pixels_hard_to_correct()
-        self.corrected = self._find_corrected_pixels()
-        self._mask_uncorrected_attenuation()
-
-    def _get_lwc_change_rate_error(self) -> np.ndarray:
-        atmosphere = (self._model["temperature"], self._model["pressure"])
-        return fill_clouds_with_lwc_dz(atmosphere, self._liquid_in_pixel)
-
-    def _get_liquid_atten(self) -> ma.MaskedArray:
-        """Finds radar liquid attenuation."""
-        lwp = ma.copy(self._mwr["lwp"][:])
-        lwp[lwp < 0] = 0
-        lwc = calc_adiabatic_lwc(self._lwc_dz_err, self._dheight)
-        lwc_scaled = distribute_lwp_to_liquid_clouds(lwc, lwp)
-        return self._calc_attenuation(lwc_scaled)
-
-    def _get_liquid_atten_err(self) -> ma.MaskedArray:
-        """Finds radar liquid attenuation error."""
-        lwc_err_scaled = distribute_lwp_to_liquid_clouds(
-            self._lwc_dz_err,
-            self._mwr["lwp_error"][:],
-        )
-        return self._calc_attenuation(lwc_err_scaled)
-
-    def _calc_attenuation(self, lwc_scaled: np.ndarray) -> ma.MaskedArray:
-        """Calculates liquid attenuation (dB)."""
-        liquid_attenuation = ma.zeros(lwc_scaled.shape)
-        spec_liq = self._model["specific_liquid_atten"]
-        lwp_cumsum = ma.cumsum(lwc_scaled[:, :-1] * spec_liq[:, :-1], axis=1)
-        liquid_attenuation[:, 1:] = TWO_WAY * lwp_cumsum * M_TO_KM
-        return liquid_attenuation
-
-    def _find_pixels_hard_to_correct(self) -> np.ndarray:
-        melting_layer = utils.isbit(self.classification.category_bits, 3)
-        hard_to_correct = np.cumsum(melting_layer, axis=1) >= 1
-        hard_to_correct[self.classification.is_rain == 1, :] = True
-        attenuated = self._find_attenuated_part_of_atmosphere()
-        hard_to_correct[attenuated & self.atten.mask] = True
-        return hard_to_correct
-
-    def _find_corrected_pixels(self) -> np.ndarray:
-        proper_values = ma.array(self.atten > 0)
-        filled = False
-        return proper_values.filled(filled) & ~self.uncorrected
-
-    def _mask_uncorrected_attenuation(self) -> None:
-        self.atten[self.uncorrected] = ma.masked
-
-    def _find_attenuated_part_of_atmosphere(self) -> np.ndarray:
-        return np.cumsum(self._lwc_dz_err, axis=1) > 0
-
-
-def fill_clouds_with_lwc_dz(atmosphere: tuple, is_liquid: np.ndarray) -> np.ndarray:
-    """Fills liquid clouds with lwc change rate at the cloud bases.
-
-    Args:
-        atmosphere: 2-element tuple containing temperature (K) and pressure (Pa).
-        is_liquid: Boolean array indicating presence of liquid clouds.
-
-    Returns:
-        Liquid water content change rate (kg/m3/m), so that for each cloud the base
-        value is filled for the whole cloud.
-
-    """
-    lwc_dz = get_lwc_change_rate_at_bases(atmosphere, is_liquid)
-    lwc_dz_filled = ma.zeros(lwc_dz.shape)
-    lwc_dz_filled[is_liquid] = utils.ffill(lwc_dz[is_liquid])
-    return lwc_dz_filled
-
-
-def get_lwc_change_rate_at_bases(
-    atmosphere: tuple,
-    is_liquid: np.ndarray,
-) -> np.ndarray:
-    """Finds LWC change rate in liquid cloud bases.
-
-    Args:
-        atmosphere: 2-element tuple containing temperature (K) and pressure (Pa).
-        is_liquid: Boolean array indicating presence of liquid clouds.
-
-    Returns:
-        Liquid water content change rate at cloud bases (kg/m3/m).
-
-    """
-    liquid_bases = find_cloud_bases(is_liquid)
-    lwc_dz = ma.zeros(liquid_bases.shape)
-    lwc_dz[liquid_bases] = calc_lwc_change_rate(
-        atmosphere[0][liquid_bases],
-        atmosphere[1][liquid_bases],
-    )
-    return lwc_dz
-
-
-def find_cloud_bases(array: np.ndarray) -> np.ndarray:
-    """Finds bases of clouds.
-
-    Args:
-        array: 2D boolean array denoting clouds or some other similar field.
-
-    Returns:
-        Boolean array indicating bases of the individual clouds.
-
-    """
-    zeros = np.zeros(array.shape[0])
-    array_padded = np.insert(array, 0, zeros, axis=1).astype(int)
-    return np.diff(array_padded, axis=1) == 1
-
-
-def find_cloud_tops(array: np.ndarray) -> np.ndarray:
-    """Finds tops of clouds.
-
-    Args:
-        array: 2D boolean array denoting clouds or some other similar field.
-
-    Returns:
-        Boolean array indicating tops of the individual clouds.
-
-    """
-    array_flipped = np.fliplr(array)
-    bases_of_flipped = find_cloud_bases(array_flipped)
-    return np.fliplr(bases_of_flipped)
-
-
-def find_lowest_cloud_bases(
-    cloud_mask: np.ndarray,
-    height: np.ndarray,
-) -> ma.MaskedArray:
-    """Finds altitudes of cloud bases."""
-    cloud_heights = cloud_mask * height
-    return _find_lowest_heights(cloud_heights)
-
-
-def find_highest_cloud_tops(
-    cloud_mask: np.ndarray,
-    height: np.ndarray,
-) -> ma.MaskedArray:
-    """Finds altitudes of cloud tops."""
-    cloud_heights = cloud_mask * height
-    cloud_heights_flipped = np.fliplr(cloud_heights)
-    return _find_lowest_heights(cloud_heights_flipped)
-
-
-def _find_lowest_heights(cloud_heights: np.ndarray) -> ma.MaskedArray:
-    inds = (cloud_heights != 0).argmax(axis=1)
-    heights = np.array([cloud_heights[i, ind] for i, ind in enumerate(inds)])
-    return ma.masked_equal(heights, 0.0)
-
-
-def calc_adiabatic_lwc(lwc_change_rate: np.ndarray, dheight: float) -> np.ndarray:
-    """Calculates adiabatic liquid water content (kg/m3).
-
-    Args:
-        lwc_change_rate: Liquid water content change rate (kg/m3/m) calculated at the
-            base of each cloud and filled to that cloud.
-        dheight: Median difference of the height vector (m).
-
-    Returns:
-        Liquid water content (kg/m3).
-
-    """
-    is_liquid = lwc_change_rate != 0
-    ind_from_base = utils.cumsumr(is_liquid, axis=1)
-    return ind_from_base * dheight * lwc_change_rate
-
-
-def distribute_lwp_to_liquid_clouds(lwc: np.ndarray, lwp: np.ndarray) -> np.ndarray:
-    """Finds LWC that would produce measured LWP.
-
-    Calculates LWP-weighted, normalized LWC. This is the measured
-    LWP distributed to liquid cloud pixels according to their
-    theoretical proportion, i.e., sum(scaled LWC) = measured LWP.
-
-    Args:
-        lwc: 2D liquid water content (kg/m3).
-        lwp: 1D liquid water path (kg/m2).
-
-    Returns:
-        2D LWP-weighted, normalized LWC (kg/m2).
-
-    """
-    lwc_sum = ma.sum(lwc, axis=1)
-    return (lwc.T / lwc_sum * lwp).T