pycontrails 0.58.0__cp312-cp312-macosx_10_13_x86_64.whl → 0.59.0__cp312-cp312-macosx_10_13_x86_64.whl

pycontrails/_version.py

@@ -28,7 +28,7 @@ version_tuple: VERSION_TUPLE
 commit_id: COMMIT_ID
 __commit_id__: COMMIT_ID
 
-__version__ = version = '0.58.0'
-__version_tuple__ = version_tuple = (0, 58, 0)
+__version__ = version = '0.59.0'
+__version_tuple__ = version_tuple = (0, 59, 0)
 
-__commit_id__ = commit_id = 'g142ff7355'
+__commit_id__ = commit_id = 'g3ff518b95'

pycontrails/core/interpolation.py

@@ -26,15 +26,15 @@ class PycontrailsRegularGridInterpolator(scipy.interpolate.RegularGridInterpolat
 
     This class is a thin wrapper around the
     :class:`scipy.interpolate.RegularGridInterpolator` in order to make typical
-    ``pycontrails`` linear interpolation use-cases more performant:
+    pycontrails linear interpolation use-cases more performant:
 
-    #. Avoid ``RegularGridInterpolator`` constructor validation when `method="linear"`.
+    #. Avoid ``RegularGridInterpolator`` constructor validation when ``method="linear"``.
        In :func:`interp`, parameters are carefully crafted to fit into the intended form,
        thereby making validation unnecessary.
     #. Override the :meth:`_evaluate_linear` method with a faster implementation. See
        the :meth:`_evaluate_linear` docstring for more information.
 
-    **This class should not be used directly. Instead, use the ``interp`` function.**
+    **This class should not be used directly. Instead, use the** :func:`interp` **function.**
 
     .. versionchanged:: 0.40.0
 
@@ -428,8 +428,8 @@ def interp(
         In particular, the dimensions of ``da`` must be ``longitude``, ``latitude``,
         ``level``, and ``time``. The three spatial dimensions must be monotonically
         increasing with ``float64`` dtype. The ``time`` dimension must be
-        monotonically increasing with ``datetime64`` dtype.
-        Assumed to be cheap to load into memory (:attr:`xr.DataArray.values` is
+        monotonically increasing with :class:`numpy.datetime64` dtype.
+        Assumed to be cheap to load into memory (:attr:`xarray.DataArray.values` is
         used without hesitation).
     method : str
         Passed into :class:`scipy.interpolate.RegularGridInterpolator`.
@@ -442,7 +442,7 @@ def interp(
         ``coords``.
     indices : tuple | None, optional
         Experimental. Provide intermediate artifacts computed by
-        :meth:``scipy.interpolate.RegularGridInterpolator._find_indices`
+        :meth:`scipy.interpolate.RegularGridInterpolator._find_indices`
         to avoid redundant computation. If known and provided, this can speed
         up interpolation by avoiding an unnecessary call to ``_find_indices``.
         By default, None. Must be used precisely.

pycontrails/core/met.py

@@ -354,8 +354,7 @@ class MetBase(ABC, Generic[XArrayType]):
     def hash(self) -> str:
         """Generate a unique hash for this met instance.
 
-        Note this is not as robust as it could be since `repr`
-        cuts off.
+        Note this is not as robust as it could be since :func:`repr` cuts off.
 
         Returns
         -------
@@ -686,7 +685,7 @@ class MetBase(ABC, Generic[XArrayType]):
 class MetDataset(MetBase):
     """Meteorological dataset with multiple variables.
 
-    Composition around xr.Dataset to enforce certain
+    Composition around :class:`xarray.Dataset` to enforce certain
     variables and dimensions for internal usage
 
     Parameters
@@ -698,17 +697,17 @@ class MetDataset(MetBase):
         Defaults to None.
     wrap_longitude : bool, optional
         Wrap data along the longitude dimension. If True, duplicate and shift longitude
-        values (ie, -180 -> 180) to ensure that the longitude dimension covers the entire
+        values (ie, ``-180 -> 180``) to ensure that the longitude dimension covers the entire
         interval ``[-180, 180]``. Defaults to False.
     copy : bool, optional
         Copy data on construction. Defaults to True.
     attrs : dict[str, Any], optional
-        Attributes to add to :attr:`data.attrs`. Defaults to None.
-        Generally, pycontrails :class:`Models` may use the following attributes:
+        Attributes to add to :attr:`data.attrs`. Defaults to None. Generally, pycontrails
+        :class:`pycontrails.core.models.Models` may use the following attributes:
 
-        - ``provider``: Name of the data provider (e.g. "ECMWF").
-        - ``dataset``: Name of the dataset (e.g. "ERA5").
-        - ``product``: Name of the product type (e.g. "reanalysis").
+        - ``provider``: Name of the data provider (e.g. ``"ECMWF"``).
+        - ``dataset``: Name of the dataset (e.g. ``"ERA5"``).
+        - ``product``: Name of the product type (e.g. ``"reanalysis"``).
 
     **attrs_kwargs : Any
         Keyword arguments to add to :attr:`data.attrs`. Defaults to None.
@@ -812,7 +811,7 @@ class MetDataset(MetBase):
         return MetDataArray._from_fastpath(da)
 
     def get(self, key: str, default_value: Any = None) -> Any:
-        """Shortcut to :meth:`data.get(k, v)` method.
+        """Shortcut to :meth:`xarray.Dataset.get` method.
 
         Parameters
         ----------
@@ -889,7 +888,7 @@ class MetDataset(MetBase):
 
         See Also
         --------
-        - :meth:`xarray.Dataset.update`
+        xarray.Dataset.update
         """
         other = other or {}
         other.update(kwargs)
@@ -957,8 +956,8 @@ class MetDataset(MetBase):
         Returns
         -------
         list[str]
-            List of met keys verified in MetDataset.
-            Returns an empty list if any MetVariable is missing.
+            List of met keys verified in :class:`MetDataset`.
+            Returns an empty list if any :class:`MetVariable` is missing.
 
         Raises
         ------
@@ -1141,12 +1140,12 @@ class MetDataset(MetBase):
 
     @property
     def provider_attr(self) -> str:
-        """Look up the 'provider' attribute with a custom error message.
+        """Look up the ``"provider"`` attribute with a custom error message.
 
         Returns
         -------
         str
-            Provider of the data. If not one of 'ECMWF' or 'NCEP',
+            Provider of the data. If not one of ``"ECMWF"`` or ``"NCEP"``,
             a warning is issued.
         """
         supported = ("ECMWF", "NCEP")
@@ -1155,13 +1154,13 @@ class MetDataset(MetBase):
 
     @property
     def dataset_attr(self) -> str:
-        """Look up the 'dataset' attribute with a custom error message.
+        """Look up the ``"dataset"`` attribute with a custom error message.
 
         Returns
         -------
         str
-            Dataset of the data. If not one of 'ERA5', 'HRES', 'IFS',
-            or 'GFS', a warning is issued.
+            Dataset of the data. If not one of ``"ERA5"``, ``"HRES"``, ``"IFS"``,
+            or ``"GFS"``, a warning is issued.
         """
         supported = ("ERA5", "HRES", "IFS", "GFS")
         examples = {
@@ -1173,13 +1172,13 @@ class MetDataset(MetBase):
 
     @property
     def product_attr(self) -> str:
-        """Look up the 'product' attribute with a custom error message.
+        """Look up the ``"product"`` attribute with a custom error message.
 
         Returns
         -------
         str
-            Product of the data. If not one of 'forecast', 'ensemble', or 'reanalysis',
-            a warning is issued.
+            Product of the data. If not one of ``"forecast"``, ``"ensemble"``,
+            or ``"reanalysis"``, a warning is issued.
 
         """
         supported = ("reanalysis", "forecast", "ensemble")
@@ -1208,6 +1207,7 @@ class MetDataset(MetBase):
 
             By default, this method returns a new :class:`MetDataset` instead
             of renaming in place. To retain the old behavior, set ``inplace=True``.
+            The ``inplace`` behavior is deprecated and will be removed in a future release.
 
         Parameters
         ----------
@@ -1224,6 +1224,11 @@ class MetDataset(MetBase):
         data_renamed = standardize_variables(self.data, variables)
 
         if inplace:
+            warnings.warn(
+                "The inplace behavior is deprecated and will be removed in a future release. ",
+                DeprecationWarning,
+                stacklevel=2,
+            )
             self.data = data_renamed
             return None
 

pycontrails/datalib/ecmwf/model_levels.py

@@ -139,9 +139,9 @@ def model_level_pressure(sp: xr.DataArray, model_levels: npt.ArrayLike) -> xr.Da
             [564.02437124, 560.81744834, 557.61052544, 554.40360254],
             [551.19667964, 547.98975674, 544.78283384, 541.57591094]]])
     Coordinates:
+      * model_level  (model_level) int64 16B 80 100
       * longitude    (longitude) float64 32B -180.0 -60.0 60.0 180.0
       * latitude     (latitude) float64 32B -90.0 -30.0 30.0 90.0
-      * model_level  (model_level) int64 16B 80 100
 
     See Also
     --------

pycontrails/datalib/goes.py

@@ -308,9 +308,14 @@ def gcs_goes_path(
     fs = fs or gcsfs.GCSFileSystem(token="anon")
     rpaths = fs.glob(rpath)
 
-    out = [r for r in rpaths if _extract_band_from_rpath(r) in bands]
-    if not out:
-        raise RuntimeError(f"No data found for {time} in {region} for bands {bands}")
+    out = []
+    for r in rpaths:
+        if (band := _extract_band_from_rpath(r)) in bands:
+            out.append(r)
+            bands.remove(band)
+
+    if bands:
+        raise FileNotFoundError(f"No data found for {time} in {region} for bands {bands}")
     return out
 
 
@@ -427,6 +432,8 @@ class GOES:
 
     """
 
+    __slots__ = ("bands", "bucket", "cachestore", "fs", "region")
+
     __marker = object()
 
     def __init__(

pycontrails/datalib/gruan.py

@@ -0,0 +1,343 @@
+"""Support for accessing `GRUAN <https://www.gruan.org/>`_ data over FTP."""
+
+import datetime
+import ftplib
+import functools
+import os
+import tempfile
+from concurrent import futures
+
+import xarray as xr
+
+from pycontrails.core import cache
+
+#: GRUAN FTP server address
+FTP_SERVER = "ftp.ncdc.noaa.gov"
+
+#: Base path for GRUAN data on the FTP server
+FTP_BASE_PATH = "/pub/data/gruan/processing/level2"
+
+#: All available GRUAN products and sites on the FTP server as of 2025-10
+#: This is simply the hardcoded output of :func:`available_sites` at that time to
+#: avoid a lookup that changes infrequently.
+AVAILABLE_PRODUCTS_TO_SITES = {
+    "RS-11G-GDP.1": ["SYO", "TAT", "NYA", "LIN"],
+    "RS41-EDT.1": ["LIN", "POT", "SNG"],
+    "RS92-GDP.1": ["BOU", "CAB", "LIN", "PAY", "POT", "SOD", "TAT"],
+    "RS92-GDP.2": [
+        "BAR",
+        "BEL",
+        "BOU",
+        "CAB",
+        "DAR",
+        "GRA",
+        "LAU",
+        "LIN",
+        "MAN",
+        "NAU",
+        "NYA",
+        "PAY",
+        "POT",
+        "REU",
+        "SGP",
+        "SOD",
+        "TAT",
+        "TEN",
+        "GVN",
+    ],
+    "RS92-PROFILE-BETA.2": ["BOU", "CAB", "LIN", "POT", "SOD", "TAT"],
+    "RS92-PROFILE-BETA.3": ["BOU", "CAB", "LIN", "POT", "SOD", "TAT"],
+}
+
+
+def extract_gruan_time(filename: str) -> tuple[datetime.datetime, int]:
+    """Extract launch time and revision number from a GRUAN filename.
+
+    Parameters
+    ----------
+    filename : str
+        GRUAN filename, e.g. "LIN-RS-01_2_RS92-GDP_002_20210125T132400_1-000-001.nc"
+
+    Returns
+    -------
+    tuple[datetime.datetime, int]
+        Launch time as a datetime object and revision number as an integer.
+    """
+    parts = filename.split("_")
+    if len(parts) != 6:
+        raise ValueError(f"Unexpected filename format: {filename}")
+    time_part = parts[4]
+    try:
+        time = datetime.datetime.strptime(time_part, "%Y%m%dT%H%M%S")
+    except ValueError as e:
+        raise ValueError(f"Unexpected time segment: {time_part}") from e
+
+    revision_part = parts[5].removesuffix(".nc")
+    if not revision_part[-3:].isdigit():
+        raise ValueError(f"Unexpected revision segment: {revision_part}")
+    revision = int(revision_part[-3:])
+
+    return time, revision
+
+
+def _fetch_product_tree(prod: str) -> dict[str, list[str]]:
+    result = {}
+    with ftplib.FTP(FTP_SERVER) as ftp:
+        ftp.login()
+        prod_path = f"{FTP_BASE_PATH}/{prod}"
+        versions = [v.split("/")[-1] for v in ftp.nlst(prod_path)]
+
+        for v in versions:
+            version_path = f"{prod_path}/{v}"
+            sites = [s.split("/")[-1] for s in ftp.nlst(version_path)]
+
+            key = f"{prod}.{int(v.split('-')[-1])}"
+            result[key] = sites
+    return result
+
+
+@functools.cache
+def available_sites() -> dict[str, list[str]]:
+    """Get a list of available GRUAN sites for each supported product.
+
+    The :attr:`GRUAN.AVAILABLE` is a hardcoded snapshot of this data. The data returned
+    by this function does not change frequently, so it is cached for efficiency.
+
+    Returns
+    -------
+    dict[str, list[str]]
+        Mapping of product names to lists of available site identifiers.
+    """
+    with ftplib.FTP(FTP_SERVER) as ftp:
+        ftp.login()
+        files = [p.split("/")[-1] for p in ftp.nlst(FTP_BASE_PATH)]
+        products = [p for p in files if "." not in p]  # crude filter to exclude non-directories
+
+    # Compute each product tree in separate thread to speed up retrieval
+    # The FTP server only allows up to 5 connections from the same client
+    out = {}
+    with futures.ThreadPoolExecutor(max_workers=min(len(products), 5)) as tpe:
+        result = tpe.map(_fetch_product_tree, products)
+        for r in result:
+            out.update(r)
+
+    return out
+
+
+class GRUAN:
+    """Access `GRUAN <https://www.gruan.org/>`_ data over anonymous FTP.
+
+    GRUAN is the Global Climate Observing System Reference Upper-Air Network. It provides
+    high-quality measurements of atmospheric variables from ground to stratosphere
+    through a global network of radiosonde stations.
+
+    .. versionadded:: 0.59.0
+
+    Parameters
+    ----------
+    product : str
+        GRUAN data product. See :attr:`AVAILABLE` for available products. These currently
+        include:
+        - ``RS92-GDP.2``
+        - ``RS92-GDP.1``
+        - ``RS92-PROFILE-BETA.2``
+        - ``RS92-PROFILE-BETA.3``
+        - ``RS41-EDT.1``
+        - ``RS-11G-GDP.1``
+    site : str
+        GRUAN station identifier. See :attr:`AVAILABLE` for available sites for each product.
+    cachestore : cache.CacheStore | None, optional
+        Cache store to use for downloaded files. If not provided, a disk cache store
+        will be created in the user cache directory under ``gruan/``. Set to ``None``
+        to disable caching.
+
+    Notes
+    -----
+    The FTP files have the following hierarchy::
+
+        /pub/data/gruan/processing/level2/
+            {product-root}/
+                version-{NNN}/
+                    {SITE}/
+                        {YYYY}/
+                            <filename>.nc
+
+    - {product-root} is the product name without the trailing version integer (e.g. ``RS92-GDP``)
+    - version-{NNN} zero-pads to three digits (suffix ``.2`` -> ``version-002``)
+    - {SITE} is the station code (e.g. ``LIN``)
+    - {YYYY} is launch year
+    - Filenames encode launch time and revision (parsed by :func:`extract_gruan_time`)
+
+    Discovery helpers methods:
+
+        - :attr:`AVAILABLE` or :func:`available_sites` -> products and sites
+        - :meth:`years` -> list available years for (product, site)
+        - :meth:`list_files` -> list available NetCDF files for the given year
+        - :meth:`get` -> download and open a single NetCDF file as an :class:`xarray.Dataset`
+
+    Typical workflow:
+
+        1. Inspect :attr:`AVAILABLE` (fast) or call :func:`available_sites` (live)
+        2. Instantiate ``GRUAN(product, site)``
+        3. Call ``years()``
+        4. Call ``list_files(year)``
+        5. Call ``get(filename)`` for an ``xarray.Dataset``
+
+    """
+
+    # Convenience access to available sites
+    available_sites = staticmethod(available_sites)
+    AVAILABLE = AVAILABLE_PRODUCTS_TO_SITES
+
+    __slots__ = ("_ftp", "cachestore", "product", "site")
+
+    __marker = object()
+
+    def __init__(
+        self,
+        product: str,
+        site: str,
+        cachestore: cache.CacheStore | None = __marker,  # type: ignore[assignment]
+    ) -> None:
+        known = AVAILABLE_PRODUCTS_TO_SITES
+
+        if product not in known:
+            known = available_sites()  # perhaps AVAILABLE_PRODUCTS_TO_SITES is outdated
+            if product not in known:
+                raise ValueError(f"Unknown GRUAN product: {product}. Known products: {list(known)}")
+        self.product = product
+
+        if site not in known[product]:
+            known = available_sites()  # perhaps AVAILABLE_PRODUCTS_TO_SITES is outdated
+            if site not in known[product]:
+                raise ValueError(
+                    f"Unknown GRUAN site '{site}' for product '{product}'. "
+                    f"Known sites: {known[product]}"
+                )
+        self.site = site
+
+        if cachestore is self.__marker:
+            cache_root = cache._get_user_cache_dir()
+            cache_dir = f"{cache_root}/gruan"
+            cachestore = cache.DiskCacheStore(cache_dir=cache_dir)
+        self.cachestore = cachestore
+
+        self._ftp: ftplib.FTP | None = None
+
+    def __repr__(self) -> str:
+        return f"GRUAN(product='{self.product}', site='{self.site}')"
+
+    def _connect(self) -> ftplib.FTP:
+        """Connect to the GRUAN FTP server."""
+        if self._ftp is None or self._ftp.sock is None:
+            self._ftp = ftplib.FTP(FTP_SERVER)
+            self._ftp.login()
+            return self._ftp
+
+        try:
+            self._ftp.pwd()  # check if connection is still alive
+        except (*ftplib.all_errors, ConnectionError):  # type: ignore[misc]
+            # If we encounter any error, reset the connection and retry
+            self._ftp = None
+            return self._connect()
+        return self._ftp
+
+    @property
+    def base_path_product(self) -> str:
+        """Get the base path for GRUAN data product on the FTP server."""
+        product, version = self.product.rsplit(".")
+        return f"/pub/data/gruan/processing/level2/{product}/version-{version.zfill(3)}"
+
+    @property
+    def base_path_site(self) -> str:
+        """Get the base path for GRUAN data site on the FTP server."""
+        return f"{self.base_path_product}/{self.site}"
+
+    def years(self) -> list[int]:
+        """Get a list of available years for the selected product and site."""
+        ftp = self._connect()
+        ftp.cwd(self.base_path_site)
+        years = ftp.nlst()
+        return sorted(int(year) for year in years)
+
+    def list_files(self, year: int | None = None) -> list[str]:
+        """List available files for a given year.
+
+        Parameters
+        ----------
+        year : int | None, optional
+            Year to list files for. If ``None``, list files for all available years. The later
+            may be time-consuming.
+
+        Returns
+        -------
+        list[str]
+            List of available GRUAN filenames for the specified year.
+        """
+        if year is None:
+            years = self.years()
+            return sorted(file for y in years for file in self.list_files(y))
+
+        path = f"{self.base_path_site}/{year}"
+
+        ftp = self._connect()
+        try:
+            ftp.cwd(path)
+        except ftplib.error_perm as e:
+            available = self.years()
+            if year not in available:
+                msg = f"No data available for year {year}. Available years are: {available}"
+                raise ValueError(msg) from e
+            raise
+        return sorted(ftp.nlst())
+
+    def get(self, filename: str) -> xr.Dataset:
+        """Download a GRUAN dataset by filename.
+
+        Parameters
+        ----------
+        filename : str
+            GRUAN filename to download, e.g. "LIN-RS-01_2_RS92-GDP_002_20210125T132400_1-000-001.nc"
+
+        Returns
+        -------
+        xr.Dataset
+            The GRUAN dataset retrieved from the FTP server. If caching is enabled,
+            the file is downloaded to the cache store and loaded from there on subsequent calls.
+        """
+        if self.cachestore is None:
+            return self._get_no_cache(filename)
+        return self._get_with_cache(filename)
+
+    def _get_no_cache(self, filename: str) -> xr.Dataset:
+        t, _ = extract_gruan_time(filename)
+        path = f"{self.base_path_site}/{t.year}/{filename}"
+
+        ftp = self._connect()
+
+        try:
+            # On windows, NamedTemporaryFile cannot be reopened while still open.
+            # After python 3.11 support is dropped, we can use delete_on_close=False
+            # in NamedTemporaryFile to streamline this.
+            with tempfile.NamedTemporaryFile(delete=False) as tmp:
+                ftp.retrbinary(f"RETR {path}", tmp.write)
+            return xr.load_dataset(tmp.name)
+        finally:
+            os.remove(tmp.name)
+
+    def _get_with_cache(self, filename: str) -> xr.Dataset:
+        if self.cachestore is None:
+            raise ValueError("Cachestore is not configured.")
+
+        lpath = self.cachestore.path(filename)
+        if self.cachestore.exists(lpath):
+            return xr.open_dataset(lpath)
+
+        t, _ = extract_gruan_time(filename)
+        path = f"{self.base_path_site}/{t.year}/{filename}"
+
+        ftp = self._connect()
+        with open(lpath, "wb") as f:
+            ftp.retrbinary(f"RETR {path}", f.write)
+
+        return xr.open_dataset(lpath)

pycontrails/datalib/himawari/himawari.py

@@ -372,6 +372,8 @@ class Himawari:
     This interface requires the ``s3fs`` package to download data from the
     `AWS Public Dataset <https://registry.opendata.aws/himawari/>`_.
 
+    .. versionadded:: 0.57.0
+
     Parameters
     ----------
     region : HimawariRegion | str, optional
@@ -395,6 +397,8 @@ class Himawari:
     HimawariRegion
     """
 
+    __slots__ = ("bands", "bucket", "cachestore", "fs", "region")
+
     __marker = object()
 
     def __init__(

pycontrails/models/cocip/contrail_properties.py

@@ -236,7 +236,7 @@ def initial_ice_particle_number(
         phase, [:math:`# m^{-1}`]
     """
     if min_aei is not None:
-        aei = np.clip(aei, min=min_aei)  # type: ignore[call-overload]
+        aei = np.clip(aei, min_aei, None)
     return fuel_dist * aei
 
 

pycontrails/physics/thermo.py

@@ -84,7 +84,7 @@ def p_vapor(q: ArrayScalarLike, p: ArrayScalarLike) -> ArrayScalarLike:
     ArrayScalarTypeVar
         Vapor pressure, [:math:`Pa`]
     """
-    return q * p * (constants.R_v / constants.R_d)
+    return q * p / constants.epsilon
 
 
 def water_vapor_partial_pressure_along_mixing_line(
@@ -399,7 +399,7 @@ def rh(q: ArrayScalarLike, T: ArrayScalarLike, p: ArrayScalarLike) -> ArrayScala
     ArrayScalarLike
         Relative Humidity, :math:`[0 - 1]`
     """
-    return (q * p * (constants.R_v / constants.R_d)) / e_sat_liquid(T)
+    return (q * p) / (constants.epsilon * e_sat_liquid(T))
 
 
 def rhi(q: ArrayScalarLike, T: ArrayScalarLike, p: ArrayScalarLike) -> ArrayScalarLike:
@@ -419,7 +419,7 @@ def rhi(q: ArrayScalarLike, T: ArrayScalarLike, p: ArrayScalarLike) -> ArrayScal
     ArrayScalarLike
         Relative Humidity over ice, :math:`[0 - 1]`
     """
-    return (q * p * (constants.R_v / constants.R_d)) / e_sat_ice(T)
+    return (q * p) / (constants.epsilon * e_sat_ice(T))
 
 
 # --------------

{pycontrails-0.58.0.dist-info → pycontrails-0.59.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pycontrails
-Version: 0.58.0
+Version: 0.59.0
 Summary: Python library for modeling aviation climate impacts
 Author-email: "Contrails.org" <py@contrails.org>
 License-Expression: Apache-2.0
@@ -35,7 +35,6 @@ Requires-Dist: xarray>=2022.3
 Provides-Extra: complete
 Requires-Dist: pycontrails[ecmwf,gcp,gfs,jupyter,pyproj,sat,vis,zarr]; extra == "complete"
 Provides-Extra: dev
-Requires-Dist: fastparquet>=0.8; extra == "dev"
 Requires-Dist: ipdb>=0.13; extra == "dev"
 Requires-Dist: memory_profiler; extra == "dev"
 Requires-Dist: mypy>=1.8; extra == "dev"

{pycontrails-0.58.0.dist-info → pycontrails-0.59.0.dist-info}/RECORD

@@ -1,15 +1,15 @@
-pycontrails-0.58.0.dist-info/RECORD,,
-pycontrails-0.58.0.dist-info/WHEEL,sha256=mer8vOuI-KlBtJIdzgyZALHdls6RZGRaE6RPt7dfKrk,138
-pycontrails-0.58.0.dist-info/top_level.txt,sha256=Z8J1R_AiBAyCVjNw6jYLdrA68PrQqTg0t3_Yek_IZ0Q,29
-pycontrails-0.58.0.dist-info/METADATA,sha256=-pOZxATs-SGOiS7FppQtQYehP9MYFbQV2l_VnbptyJg,9129
-pycontrails-0.58.0.dist-info/licenses/LICENSE,sha256=gJ-h7SFFD1mCfR6a7HILvEtodDT6Iig8bLXdgqR6ucA,10175
-pycontrails-0.58.0.dist-info/licenses/NOTICE,sha256=fiBPdjYibMpDzf8hqcn7TvAQ-yeK10q_Nqq24DnskYg,1962
-pycontrails/_version.py,sha256=ife3IDQF-DAoy1CwSUU90Z0xGG9P1GsYeBh-Ic3uecI,714
+pycontrails-0.59.0.dist-info/RECORD,,
+pycontrails-0.59.0.dist-info/WHEEL,sha256=mer8vOuI-KlBtJIdzgyZALHdls6RZGRaE6RPt7dfKrk,138
+pycontrails-0.59.0.dist-info/top_level.txt,sha256=Z8J1R_AiBAyCVjNw6jYLdrA68PrQqTg0t3_Yek_IZ0Q,29
+pycontrails-0.59.0.dist-info/METADATA,sha256=s3R-ww4kBEE-scE3ZTAxgvdndp--LPpmF1Bgoa9dofc,9081
+pycontrails-0.59.0.dist-info/licenses/LICENSE,sha256=gJ-h7SFFD1mCfR6a7HILvEtodDT6Iig8bLXdgqR6ucA,10175
+pycontrails-0.59.0.dist-info/licenses/NOTICE,sha256=fiBPdjYibMpDzf8hqcn7TvAQ-yeK10q_Nqq24DnskYg,1962
+pycontrails/_version.py,sha256=4JDTQzRvD4h--EgLI4W50oyfLqa3M8rQJJovGr52AZA,714
 pycontrails/__init__.py,sha256=9ypSB2fKZlKghTvSrjWo6OHm5qfASwiTIvlMew3Olu4,2037
 pycontrails/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 pycontrails/core/vector.py,sha256=hms3hea2Y86LOuZFfOqzGZToqpOVksgXikM_S5w355w,73498
 pycontrails/core/models.py,sha256=3mDTqp1V5aae9akuYwbMGIUEkESKSYTjZeyu2IiMW7s,43915
-pycontrails/core/interpolation.py,sha256=TpLDfx-Bx2tGDjWJOyU4n9Dkz6Obl05dB7Ol5diYdsk,24816
+pycontrails/core/interpolation.py,sha256=vttZiW78RnHmzLv2lmflwR9po3ghRxwT7MtAj49WvUY,24836
 pycontrails/core/fleet.py,sha256=a_vVwAjMbjkhszg7ejP3V0yly_wJ9Va_OQEATG-9UHw,16572
 pycontrails/core/flight.py,sha256=dmqO1PxADMHIcK9U8XSnuXbP59ftQgeKvdI2xGP04ig,81460
 pycontrails/core/fuel.py,sha256=kJZ3P1lPm1L6rdPREM55XQ-VfJ_pt35cP4sO2Nnvmjs,4332
@@ -17,18 +17,19 @@ pycontrails/core/polygon.py,sha256=kXYwj1Xy-mo8GEWXFAO_OJEtZbGx11DE_sZw6iyNvN4,1
 pycontrails/core/cache.py,sha256=JQIy1sQf0Vil7wlEk9ZIvVacnOD4wM3X8-UkMFjR2wQ,28177
 pycontrails/core/__init__.py,sha256=p0O09HxdeXU0X5Z3zrHMlTfXa92YumT3fJ8wJBI5ido,856
 pycontrails/core/flightplan.py,sha256=0mvA3IO19Sap-7gwpmEIV35_mg6ChvajwhurvjZZt_U,7521
-pycontrails/core/met.py,sha256=JosS0DXGEB4NBQ40WU1b4yG04h3DkT7Vs-P0jK-UkCY,104453
+pycontrails/core/met.py,sha256=C_2NoRovK6Irz20XyDDLNYRfOwfHcOFtI9hE-IeuvHA,104859
 pycontrails/core/aircraft_performance.py,sha256=CPIgIi5nUuCHiNVLAvZcWECRfakmMd-wUWd3lMA6oGM,28204
 pycontrails/core/airports.py,sha256=CzZrgJNZ7wtNv8vg9sJczMhFov7k0gmrGR4tRKCH8i8,6782
 pycontrails/core/met_var.py,sha256=g69vqbxpJeXEQU8vrrcoUR1PX3zCo2-k3au1Lv2TiIw,12027
-pycontrails/core/rgi_cython.cpython-312-darwin.so,sha256=OCktVJT2FYKkU7ZU5JMwrNQMR8cKaXRrrIGOYWgcYnk,330104
+pycontrails/core/rgi_cython.cpython-312-darwin.so,sha256=6p7q-UbhSrtIggDRiB8CphGikkElMQSrBWwCu4S3xOo,330104
 pycontrails/core/coordinates.py,sha256=0ySsHtqTon7GMbuwmmxMbI92j3ueMteJZh4xxNm5zto,5391
-pycontrails/datalib/goes.py,sha256=kZIhx1cKIcNtIxQ3IiuN45aCpwDzzTOxi0Y4fjyaOLs,27481
+pycontrails/datalib/goes.py,sha256=QFy7peNrfQcmujEcmUzK8zR1xBxli2vHxqak84PtYTc,27632
 pycontrails/datalib/landsat.py,sha256=6ylDkAjnyX7b4ZbHn4bprO8HB8ADPFyMkwWehIs8FLg,20915
 pycontrails/datalib/geo_utils.py,sha256=w6VYhJQeMpBXaBclqANv4Nn0yqPIxlQr6GTUpjArTj0,9070
 pycontrails/datalib/__init__.py,sha256=hW9NWdFPC3y_2vHMteQ7GgQdop3917MkDaf5ZhU2RBY,369
+pycontrails/datalib/gruan.py,sha256=n-MI2iCAi3FdlbuWJKSYJFeKCtiRn3e-o4pRrUOC_dU,11673
 pycontrails/datalib/sentinel.py,sha256=ed1l1avq8lBvQinY_vNSsWRcpqxUdAPY61AGyPcLawo,23532
-pycontrails/datalib/himawari/himawari.py,sha256=uJXhu1URdu3Hen9wMwkgxyPHdjd_xyIt9E7WblVy2xQ,23327
+pycontrails/datalib/himawari/himawari.py,sha256=GueKQsytmuy4ojaECnRq-PZKNfIcJLvftDn56fqwMLQ,23424
 pycontrails/datalib/himawari/__init__.py,sha256=SWupVbeuyK07IPDCgiNjN6hoLB7hlceabJ3fixhDkl0,619
 pycontrails/datalib/himawari/header_struct.py,sha256=WbPkNBNUVm8tGKU8wpj4rldY17g5MHQ_OfbWicZSokc,9975
 pycontrails/datalib/_met_utils/metsource.py,sha256=mlKcRko5ZKuYK5uwWn6AAgUSJLMQAYq1nFqskVMGgYo,23999
@@ -40,7 +41,7 @@ pycontrails/datalib/ecmwf/variables.py,sha256=lU3BNe265XVhCXvdMwZqfkWQwtsetZxVRL
 pycontrails/datalib/ecmwf/hres_model_level.py,sha256=CcxMKiFJyLvM9njmBVywAXJxyWE7atsgHXBubKJQqHM,17779
 pycontrails/datalib/ecmwf/__init__.py,sha256=wdfhplEaW2UKTItIoshTtVEjbPyfDYoprTJNxbKZuvA,2021
 pycontrails/datalib/ecmwf/common.py,sha256=axOxvdrey9YD34uk0Ocav08MxKvC2uVaiwvyQgFZMEw,3970
-pycontrails/datalib/ecmwf/model_levels.py,sha256=_kgpnogaS6MlfvTX9dB5ASTHFUlZuQ_DRb-VADwEa0k,16996
+pycontrails/datalib/ecmwf/model_levels.py,sha256=xdH2B5KvzAIk2_OIOgZXFIZ6ndR6MxxqZsvstO0cNmg,16996
 pycontrails/datalib/ecmwf/ifs.py,sha256=_1UarorPp9VlgFZc-NnZy8YnfEqBdp7GV1A-ye6JqS8,10733
 pycontrails/datalib/ecmwf/static/model_level_dataframe_v20240418.csv,sha256=PmvGLRzn6uuCKSwiasSuVcehvvmSaqP7cnLuN6hhCQQ,9788
 pycontrails/datalib/gfs/gfs.py,sha256=VqS0MRLawgzkBDpjDUYoswXByIy6XUqA9XP7lM1ueBk,22238
@@ -100,7 +101,7 @@ pycontrails/models/cocip/cocip_params.py,sha256=BWmTt6yE4m-LM7lyCtj05FK3wVvU9n7i
 pycontrails/models/cocip/wake_vortex.py,sha256=F5S8n4eBrBM-7qNcVUtX3IrXD7Kt9pWnrKj6UK-HGeA,14555
 pycontrails/models/cocip/cocip_uncertainty.py,sha256=TZ85xAbDc5zRgQKP7wb3AfHoUIvkfHycuX86dwZCqwM,12257
 pycontrails/models/cocip/radiative_heating.py,sha256=1U4SQWwogtyQ2u6J996kAHP0OfpZ3hH2_x4Cyt3Cy8U,18984
-pycontrails/models/cocip/contrail_properties.py,sha256=BRldMsxNIYlKHg5ozD5_wWtJ7OGTcQvinkQUioHqaqk,55732
+pycontrails/models/cocip/contrail_properties.py,sha256=V1ACeeYaD_a_piqarEnjzxOvYh7cP-1c-xaWv5Kgc8U,55703
 pycontrails/models/cocip/unterstrasser_wake_vortex.py,sha256=bIRS-Z4MRMdkYtth2RaDe5h1ZN0HvCE_Sw96PXQEHKQ,18931
 pycontrails/models/ps_model/__init__.py,sha256=Fuum5Rq8ya8qkvbeq2wh6NDo-42RCRnK1Y-2syYy0Ck,553
 pycontrails/models/ps_model/ps_model.py,sha256=fgFekJpGuAu73KvpfLhlAbIwR7JJGwQpLILWmrONywc,31925
@@ -116,7 +117,7 @@ pycontrails/physics/geo.py,sha256=ITK23l1A2lzjNPTFC8ZKyQH59I5Cy_TvuvM_gbALo94,36
 pycontrails/physics/units.py,sha256=p-6PzFLpVCMpvmfrhXVh3Hs-nMJw9Y1x-hvgnL9Lo9c,12281
 pycontrails/physics/constants.py,sha256=JHYL2IJY7del2BE_1QfKaEwtIwkbtyHvyxlm_JPHR90,3201
 pycontrails/physics/__init__.py,sha256=_1eWbEy6evEWdfJCEkwDiSdpiDNzNWEPVqaPekHyhwU,44
-pycontrails/physics/thermo.py,sha256=mGVduIY1c1DM6dVAdGcklGIJmscHXo3glDXmdeUtAZk,15376
+pycontrails/physics/thermo.py,sha256=bKn1rlHHIoLvkrd42Vpox6pQVAu4OjonYSmxXXFvPug,15338
 pycontrails/physics/jet.py,sha256=Je1d3vgbBEaVIAL1WZ3C-4p2f9fy9dWOjP5vFVsGGh8,30358
 pycontrails/physics/static/iata-cargo-load-factors-20250221.csv,sha256=ixsnQk1DyGxHMo0pDy4aOoQIwgOyrGfhMRPumEwPMBc,3841
 pycontrails/physics/static/iata-passenger-load-factors-20250221.csv,sha256=Q2olRIqUpbOaavvM5ikG8m1v1YQAN3KLNHeFDPvM53Q,3835