cloudnetpy 1.86.3__py3-none-any.whl → 1.86.5__py3-none-any.whl

cloudnetpy/instruments/copernicus.py

@@ -12,7 +12,7 @@ import numpy as np
 from cloudnetpy import concat_lib, output, utils
 from cloudnetpy.instruments.instruments import COPERNICUS
 from cloudnetpy.instruments.nc_radar import ChilboltonRadar
-from cloudnetpy.metadata import MetaData
+from cloudnetpy.metadata import COMMON_ATTRIBUTES, MetaData
 
 
 def copernicus2nc(
@@ -105,7 +105,7 @@ def copernicus2nc(
             copernicus.add_site_geolocation()
             valid_indices = copernicus.add_zenith_and_azimuth_angles(
                 elevation_threshold=1.1,
-                elevation_diff_threshold=0.1,
+                elevation_diff_threshold=0.5,
                 azimuth_diff_threshold=0.1,
             )
             copernicus.screen_time_indices(valid_indices)
@@ -131,10 +131,12 @@ class Copernicus(ChilboltonRadar):
         self.instrument = COPERNICUS
 
     def calibrate_reflectivity(self) -> None:
-        default_offset = -146.8  # TODO: check this value
-        calibration_factor = self.site_meta.get("calibration_offset", default_offset)
-        self.data["Zh"].data[:] += calibration_factor
-        self.append_data(np.array(calibration_factor), "calibration_offset")
+        default_offset = (
+            -141.5
+        )  # Taken from 2016-01 raw files where this was already applied
+        zh_offset = self.site_meta.get("Zh_offset", default_offset)
+        self.data["Zh"].data[:] += zh_offset
+        self.append_data(np.array(zh_offset, dtype=np.float32), "Zh_offset")
 
     def mask_corrupted_values(self) -> None:
         """Experimental masking of corrupted Copernicus data.
@@ -169,12 +171,13 @@ class Copernicus(ChilboltonRadar):
 
 
 ATTRIBUTES = {
-    "calibration_offset": MetaData(
+    "Zh_offset": MetaData(
         long_name="Radar reflectivity calibration offset",
         units="dBZ",
         comment="Calibration offset applied.",
         dimensions=None,
     ),
+    "Zh": COMMON_ATTRIBUTES["Zh"]._replace(ancillary_variables="Zh_offset"),
     "range_offset": MetaData(
         long_name="Radar range offset",
         units="m",

cloudnetpy/instruments/galileo.py

@@ -92,6 +92,7 @@ def galileo2nc(
                 galileo.check_date(date)
             galileo.sort_timestamps()
             galileo.remove_duplicate_timestamps()
+            galileo.screen_negative_altitudes()
             snr_limit = site_meta.get("snr_limit", 3)
             galileo.screen_by_snr(snr_limit=snr_limit)
             galileo.mask_clutter()
@@ -128,7 +129,6 @@ class Galileo(ChilboltonRadar):
         self.instrument = GALILEO
 
     def mask_clutter(self) -> None:
-        """Masks clutter."""
         # Only strong Z values are valid
         n_low_gates = 15
         ind = np.where(self.data["Zh"][:, :n_low_gates] < -15) and np.where(
@@ -136,6 +136,20 @@ class Galileo(ChilboltonRadar):
         )
         self.data["v"].mask_indices(ind)
 
+    def screen_negative_altitudes(self) -> None:
+        range_var = self.data["range"][:]
+        valid_idx = np.where(range_var > 0)[0]
+        if valid_idx.size == 0:
+            msg = "No valid altitudes found."
+            raise ValueError(msg)
+        if valid_idx.size == range_var.shape[0]:
+            return
+        for arr in self.data.values():
+            if utils.isscalar(arr.data):
+                continue
+            if arr[:].shape[-1] == range_var.shape[0]:
+                arr.data = arr[..., valid_idx]
+
 
 ATTRIBUTES = {
     "antenna_diameter": MetaData(

cloudnetpy/version.py

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

{cloudnetpy-1.86.3.dist-info → cloudnetpy-1.86.5.dist-info}/METADATA

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

{cloudnetpy-1.86.3.dist-info → cloudnetpy-1.86.5.dist-info}/RECORD

@@ -9,7 +9,7 @@ cloudnetpy/metadata.py,sha256=qqLXSUVhIi6H_Z1zlfXeEkc9qBqA1XO-UwsCVxjIcfk,7317
 cloudnetpy/output.py,sha256=NpuxmtLk2JLTdubCupdL4rFbx30YqgeU6vwAe_v8thc,15133
 cloudnetpy/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 cloudnetpy/utils.py,sha256=9MMYok1A7oJS2XBkIFhaiDWt_MYLo8tLzF3_880masM,34321
-cloudnetpy/version.py,sha256=fFMhfzPgAQTnH4EcNXqfMkShM7r1U4B56xzlJp8XWrQ,72
+cloudnetpy/version.py,sha256=Va6HjzCd42qubMOwHV5hz8UCWeFMp6EzRf_6g1bNCNA,72
 cloudnetpy/categorize/__init__.py,sha256=gtvzWr0IDRn2oA6yHBvinEhTGTuub-JkrOv93lBsgrE,61
 cloudnetpy/categorize/atmos_utils.py,sha256=uWc9TABVYPI0sn4H5Az9Jf6NVRaWyEKIi17f0pAJQxE,10679
 cloudnetpy/categorize/attenuation.py,sha256=Y_-fzmQTltWTqIZTulJhovC7a6ifpMcaAazDJcnMIOc,990
@@ -39,10 +39,10 @@ cloudnetpy/instruments/ceilo.py,sha256=scb3n0gLOZ1gPs_WdWw9KS5Ut_mtCvSdWFMKKCvXh
 cloudnetpy/instruments/ceilometer.py,sha256=c37uteeuGnlE-o-Smu49H2qQJw6qZ0tc3Bzhyr1FoSo,13063
 cloudnetpy/instruments/cl61d.py,sha256=JsCHqVzFGhZi-5xcnsB507FDpyuw83uSWK3IFO3DhdI,2238
 cloudnetpy/instruments/cloudnet_instrument.py,sha256=B1UkiB0ytnT3MWYalEegql5QIPaMLg5bJy5xI50JEco,4503
-cloudnetpy/instruments/copernicus.py,sha256=ygEViERBSJdMeP9OxfLelZRDEbSRzY8n17ruYie2wm4,6970
+cloudnetpy/instruments/copernicus.py,sha256=AuGvPjn6hofzSKlh2M7knb5zBgzY3xfDvu3WF6KGlp8,7091
 cloudnetpy/instruments/da10.py,sha256=FcXEu6dfAVvYzMvauiJBSyEkYhGa8hBuwrrdo-3LwHE,1960
 cloudnetpy/instruments/fd12p.py,sha256=M5yJDsPi0tHWMDUYOsoqciyKdRQ-zHEq2t19KnpZ6Io,6915
-cloudnetpy/instruments/galileo.py,sha256=f_-GkRxhNaQPbI8HpOwSmoKfGqyjmD16A0ZFgwLOIig,5137
+cloudnetpy/instruments/galileo.py,sha256=pJ5WaB-Nwzs5zgj8c2oaH7HieYxq4HF0aUilYTSRhYE,5682
 cloudnetpy/instruments/hatpro.py,sha256=TGOqwW0TfoPEYk13MFvFzwgJGzm6MVE5AsPavcIoj3I,10248
 cloudnetpy/instruments/instruments.py,sha256=WZgH7HjzM9Ane1CSnYCSLidbST8hunUeSt2lPntq9As,4999
 cloudnetpy/instruments/lufft.py,sha256=G6KeJOeltLUlGCHHEk8ns2K7WJ9ImAr25rSB2JltawE,4286
@@ -118,10 +118,10 @@ cloudnetpy/products/lwc.py,sha256=xsNiiG6dGKIkWaFk0xWTabc1bZ4ULf6SqcqHs7itAUk,19
 cloudnetpy/products/mie_lu_tables.nc,sha256=It4fYpqJXlqOgL8jeZ-PxGzP08PMrELIDVe55y9ob58,16637951
 cloudnetpy/products/mwr_tools.py,sha256=MMWnp68U7bv157-CPB2VeTQvaR6zl7sexbBT_kJ_pn8,6734
 cloudnetpy/products/product_tools.py,sha256=eyqIw_0KhlpmmYQE69RpGdRIAOW7JVPlEgkTBp2kdps,11302
-cloudnetpy-1.86.3.dist-info/licenses/LICENSE,sha256=wcZF72bdaoG9XugpyE95Juo7lBQOwLuTKBOhhtANZMM,1094
+cloudnetpy-1.86.5.dist-info/licenses/LICENSE,sha256=wcZF72bdaoG9XugpyE95Juo7lBQOwLuTKBOhhtANZMM,1094
 docs/source/conf.py,sha256=IKiFWw6xhUd8NrCg0q7l596Ck1d61XWeVjIFHVSG9Og,1490
-cloudnetpy-1.86.3.dist-info/METADATA,sha256=QQLkcSX05Op0Lks5K_uNbwteImRdke5u9CUjLkAR63I,5836
-cloudnetpy-1.86.3.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-cloudnetpy-1.86.3.dist-info/entry_points.txt,sha256=HhY7LwCFk4qFgDlXx_Fy983ZTd831WlhtdPIzV-Y3dY,51
-cloudnetpy-1.86.3.dist-info/top_level.txt,sha256=ibSPWRr6ojS1i11rtBFz2_gkIe68mggj7aeswYfaOo0,16
-cloudnetpy-1.86.3.dist-info/RECORD,,
+cloudnetpy-1.86.5.dist-info/METADATA,sha256=tF3Zehy97hcHHsGatoLodQPTYVWUrMkQdjGzdOmEyDU,5836
+cloudnetpy-1.86.5.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+cloudnetpy-1.86.5.dist-info/entry_points.txt,sha256=HhY7LwCFk4qFgDlXx_Fy983ZTd831WlhtdPIzV-Y3dY,51
+cloudnetpy-1.86.5.dist-info/top_level.txt,sha256=ibSPWRr6ojS1i11rtBFz2_gkIe68mggj7aeswYfaOo0,16
+cloudnetpy-1.86.5.dist-info/RECORD,,