pycontrails 0.50.2__cp312-cp312-win_amd64.whl → 0.51.0__cp312-cp312-win_amd64.whl

pycontrails/datalib/ecmwf/static/model_level_dataframe_v20240418.csv

@@ -0,0 +1,139 @@
+n,a [Pa],b,ph [hPa],pf [hPa],Geopotential Altitude [m],Geometric Altitude [m],Temperature [K],Density [kg/m^3]
+0,0.0,0.0,0.0,,,,,
+1,2.000365,0.0,0.02,0.01,79301.79,80301.65,198.05,1.8e-05
+2,3.102241,0.0,0.031,0.0255,73721.58,74584.91,209.21,4.2e-05
+3,4.666084,0.0,0.0467,0.0388,71115.75,71918.79,214.42,6.3e-05
+4,6.827977,0.0,0.0683,0.0575,68618.43,69365.77,221.32,9e-05
+5,9.746966,0.0,0.0975,0.0829,66210.99,66906.53,228.06,0.000127
+6,13.605424,0.0,0.1361,0.1168,63890.03,64537.43,234.56,0.000173
+7,18.608931,0.0,0.1861,0.1611,61651.77,62254.39,240.83,0.000233
+8,24.985718,0.0,0.2499,0.218,59492.5,60053.46,246.87,0.000308
+9,32.98571,0.0,0.3299,0.2899,57408.61,57930.78,252.71,0.0004
+10,42.879242,0.0,0.4288,0.3793,55396.62,55882.68,258.34,0.000512
+11,54.955463,0.0,0.5496,0.4892,53453.2,53905.62,263.78,0.000646
+12,69.520576,0.0,0.6952,0.6224,51575.15,51996.21,269.04,0.000806
+13,86.895882,0.0,0.869,0.7821,49767.41,50159.36,270.65,0.001007
+14,107.415741,0.0,1.0742,0.9716,48048.7,48413.94,270.65,0.001251
+15,131.425507,0.0,1.3143,1.1942,46416.22,46756.98,269.02,0.001546
+16,159.279404,0.0,1.5928,1.4535,44881.17,45199.69,264.72,0.001913
+17,191.338562,0.0,1.9134,1.7531,43440.23,43738.55,260.68,0.002343
+18,227.968948,0.0,2.2797,2.0965,42085.0,42364.93,256.89,0.002843
+19,269.539581,0.0,2.6954,2.4875,40808.05,41071.2,253.31,0.003421
+20,316.420746,0.0,3.1642,2.9298,39602.76,39850.56,249.94,0.004084
+21,368.982361,0.0,3.6898,3.427,38463.25,38696.94,246.75,0.004838
+22,427.592499,0.0,4.2759,3.9829,37384.22,37604.95,243.73,0.005693
+23,492.616028,0.0,4.9262,4.601,36360.94,36569.72,240.86,0.006655
+24,564.413452,0.0,5.6441,5.2851,35389.15,35586.89,238.14,0.007731
+25,643.339905,0.0,6.4334,6.0388,34465.0,34652.52,235.55,0.008931
+26,729.744141,0.0,7.2974,6.8654,33585.02,33763.05,233.09,0.010261
+27,823.967834,0.0,8.2397,7.7686,32746.04,32915.27,230.74,0.011729
+28,926.34491,0.0,9.2634,8.7516,31945.53,32106.57,228.6,0.013337
+29,1037.201172,0.0,10.372,9.8177,31177.59,31330.96,227.83,0.015012
+30,1156.853638,0.0,11.5685,10.9703,30438.54,30584.71,227.09,0.016829
+31,1285.610352,0.0,12.8561,12.2123,29726.69,29866.09,226.38,0.018793
+32,1423.770142,0.0,14.2377,13.5469,29040.48,29173.5,225.69,0.02091
+33,1571.622925,0.0,15.7162,14.977,28378.46,28505.47,225.03,0.023186
+34,1729.448975,0.0,17.2945,16.5054,27739.29,27860.64,224.39,0.025624
+35,1897.519287,0.0,18.9752,18.1348,27121.74,27237.73,223.77,0.028232
+36,2076.095947,0.0,20.761,19.8681,26524.63,26635.56,223.17,0.031013
+37,2265.431641,0.0,22.6543,21.7076,25946.9,26053.04,222.6,0.033972
+38,2465.770508,0.0,24.6577,23.656,25387.55,25489.15,222.04,0.037115
+39,2677.348145,0.0,26.7735,25.7156,24845.63,24942.93,221.5,0.040445
+40,2900.391357,0.0,29.0039,27.8887,24320.28,24413.5,220.97,0.043967
+41,3135.119385,0.0,31.3512,30.1776,23810.67,23900.02,220.46,0.047685
+42,3381.743652,0.0,33.8174,32.5843,23316.04,23401.71,219.97,0.051604
+43,3640.468262,0.0,36.4047,35.1111,22835.68,22917.85,219.49,0.055727
+44,3911.490479,0.0,39.1149,37.7598,22368.91,22447.75,219.02,0.060059
+45,4194.930664,0.0,41.9493,40.5321,21915.16,21990.82,218.57,0.064602
+46,4490.817383,0.0,44.9082,43.4287,21473.98,21546.62,218.12,0.069359
+47,4799.149414,0.0,47.9915,46.4498,21045.0,21114.77,217.7,0.07433
+48,5119.89502,0.0,51.199,49.5952,20627.87,20694.9,217.28,0.079516
+49,5452.990723,0.0,54.5299,52.8644,20222.24,20286.66,216.87,0.084916
+50,5798.344727,0.0,57.9834,56.2567,19827.95,19889.88,216.65,0.090458
+51,6156.074219,0.0,61.5607,59.7721,19443.55,19503.09,216.65,0.09611
+52,6526.946777,0.0,65.2695,63.4151,19068.35,19125.61,216.65,0.101968
+53,6911.870605,0.0,69.1187,67.1941,18701.27,18756.34,216.65,0.108045
+54,7311.869141,0.0,73.1187,71.1187,18341.27,18394.25,216.65,0.114355
+55,7727.412109,7e-06,77.281,75.1999,17987.41,18038.35,216.65,0.120917
+56,8159.354004,2.4e-05,81.6182,79.4496,17638.78,17687.77,216.65,0.127751
+57,8608.525391,5.9e-05,86.145,83.8816,17294.53,17341.62,216.65,0.134877
+58,9076.400391,0.000112,90.8774,88.5112,16953.83,16999.08,216.65,0.142321
+59,9562.682617,0.000199,95.828,93.3527,16616.09,16659.55,216.65,0.150106
+60,10065.978516,0.00034,101.0047,98.4164,16281.1,16322.83,216.65,0.158248
+61,10584.631836,0.000562,106.4153,103.71,15948.85,15988.88,216.65,0.16676
+62,11116.662109,0.00089,112.0681,109.2417,15619.3,15657.7,216.65,0.175655
+63,11660.067383,0.001353,117.9714,115.0198,15292.44,15329.24,216.65,0.184946
+64,12211.547852,0.001992,124.1337,121.0526,14968.24,15003.5,216.65,0.194646
+65,12766.873047,0.002857,130.5637,127.3487,14646.68,14680.44,216.65,0.20477
+66,13324.668945,0.003971,137.2703,133.917,14327.75,14360.05,216.65,0.215331
+67,13881.331055,0.005378,144.2624,140.7663,14011.41,14042.3,216.65,0.226345
+68,14432.139648,0.007133,151.5493,147.9058,13697.65,13727.18,216.65,0.237825
+69,14975.615234,0.009261,159.1403,155.3448,13386.45,13414.65,216.65,0.249786
+70,15508.256836,0.011806,167.045,163.0927,13077.79,13104.7,216.65,0.262244
+71,16026.115234,0.014816,175.2731,171.1591,12771.64,12797.3,216.65,0.275215
+72,16527.322266,0.018318,183.8344,179.5537,12467.99,12492.44,216.65,0.288713
+73,17008.789063,0.022355,192.7389,188.2867,12166.81,12190.1,216.65,0.302755
+74,17467.613281,0.026964,201.9969,197.3679,11868.08,11890.24,216.65,0.317357
+75,17901.621094,0.032176,211.6186,206.8078,11571.79,11592.86,216.65,0.332536
+76,18308.433594,0.038026,221.6146,216.6166,11277.92,11297.93,216.65,0.348308
+77,18685.71875,0.044548,231.9954,226.805,10986.7,11005.69,216.74,0.364545
+78,19031.289063,0.051773,242.7719,237.3837,10696.22,10714.22,218.62,0.378253
+79,19343.511719,0.059728,253.9549,248.3634,10405.61,10422.64,220.51,0.392358
+80,19620.042969,0.068448,265.5556,259.7553,10114.89,10130.98,222.4,0.406868
+81,19859.390625,0.077958,277.5852,271.5704,9824.08,9839.26,224.29,0.42179
+82,20059.931641,0.088286,290.0548,283.82,9533.2,9547.49,226.18,0.43713
+83,20219.664063,0.099462,302.9762,296.5155,9242.26,9255.7,228.08,0.452897
+84,20337.863281,0.111505,316.3607,309.6684,8951.3,8963.9,229.97,0.469097
+85,20412.308594,0.124448,330.2202,323.2904,8660.32,8672.11,231.86,0.485737
+86,20442.078125,0.138313,344.5663,337.3932,8369.35,8380.36,233.75,0.502825
+87,20425.71875,0.153125,359.4111,351.9887,8078.41,8088.67,235.64,0.520367
+88,20361.816406,0.16891,374.7666,367.0889,7787.51,7797.04,237.53,0.53837
+89,20249.511719,0.185689,390.645,382.7058,7496.68,7505.51,239.42,0.556842
+90,20087.085938,0.203491,407.0583,398.8516,7205.93,7214.09,241.31,0.57579
+91,19874.025391,0.222333,424.019,415.5387,6915.29,6922.8,243.2,0.595219
+92,19608.572266,0.242244,441.5395,432.7792,6624.76,6631.66,245.09,0.615138
+93,19290.226563,0.263242,459.6321,450.5858,6334.38,6340.68,246.98,0.635553
+94,18917.460938,0.285354,478.3096,468.9708,6044.15,6049.89,248.86,0.656471
+95,18489.707031,0.308598,497.5845,487.947,5754.1,5759.3,250.75,0.677899
+96,18006.925781,0.332939,517.4198,507.5021,5464.6,5469.3,252.63,0.699815
+97,17471.839844,0.358254,537.7195,527.5696,5176.77,5180.98,254.5,0.722139
+98,16888.6875,0.384363,558.343,548.0312,4892.26,4896.02,256.35,0.744735
+99,16262.046875,0.411125,579.1926,568.7678,4612.58,4615.92,258.17,0.767472
+100,15596.695313,0.438391,600.1668,589.6797,4338.77,4341.73,259.95,0.790242
+101,14898.453125,0.466003,621.1624,610.6646,4071.8,4074.41,261.68,0.812937
+102,14173.324219,0.4938,642.0764,631.6194,3812.53,3814.82,263.37,0.835453
+103,13427.769531,0.521619,662.8084,652.4424,3561.7,3563.69,265.0,0.857686
+104,12668.257813,0.549301,683.262,673.0352,3319.94,3321.67,266.57,0.879541
+105,11901.339844,0.576692,703.3467,693.3043,3087.75,3089.25,268.08,0.900929
+106,11133.304688,0.603648,722.9795,713.1631,2865.54,2866.83,269.52,0.921768
+107,10370.175781,0.630036,742.0855,732.5325,2653.58,2654.69,270.9,0.941988
+108,9617.515625,0.655736,760.5996,751.3426,2452.04,2452.99,272.21,0.961527
+109,8880.453125,0.680643,778.4661,769.5329,2260.99,2261.8,273.45,0.980334
+110,8163.375,0.704669,795.6396,787.0528,2080.41,2081.09,274.63,0.998368
+111,7470.34375,0.727739,812.0847,803.8622,1910.19,1910.76,275.73,1.015598
+112,6804.421875,0.749797,827.7756,819.9302,1750.14,1750.63,276.77,1.032005
+113,6168.53125,0.770798,842.6959,835.2358,1600.04,1600.44,277.75,1.047576
+114,5564.382813,0.790717,856.8376,849.7668,1459.58,1459.91,278.66,1.06231
+115,4993.796875,0.809536,870.2004,863.519,1328.43,1328.7,279.52,1.076209
+116,4457.375,0.827256,882.791,876.4957,1206.21,1206.44,280.31,1.089286
+117,3955.960938,0.843881,894.6222,888.7066,1092.54,1092.73,281.05,1.101558
+118,3489.234375,0.859432,905.7116,900.1669,987.0,987.15,281.73,1.113047
+119,3057.265625,0.873929,916.0815,910.8965,889.17,889.29,282.37,1.123777
+120,2659.140625,0.887408,925.7571,920.9193,798.62,798.72,282.96,1.133779
+121,2294.242188,0.8999,934.7666,930.2618,714.94,715.02,283.5,1.143084
+122,1961.5,0.911448,943.1399,938.9532,637.7,637.76,284.0,1.151724
+123,1659.476563,0.922096,950.9082,947.024,566.49,566.54,284.47,1.159733
+124,1387.546875,0.931881,958.1037,954.5059,500.91,500.95,284.89,1.167147
+125,1143.25,0.94086,964.7584,961.4311,440.58,440.61,285.29,1.173999
+126,926.507813,0.949064,970.9046,967.8315,385.14,385.16,285.65,1.180323
+127,734.992188,0.95655,976.5737,973.7392,334.22,334.24,285.98,1.186154
+128,568.0625,0.963352,981.7968,979.1852,287.51,287.52,286.28,1.191523
+129,424.414063,0.969513,986.6036,984.2002,244.68,244.69,286.56,1.196462
+130,302.476563,0.975078,991.023,988.8133,205.44,205.44,286.81,1.201001
+131,202.484375,0.980072,995.0824,993.0527,169.5,169.51,287.05,1.205168
+132,122.101563,0.984542,998.8081,996.9452,136.62,136.62,287.26,1.208992
+133,62.78125,0.9885,1002.225,1000.5165,106.54,106.54,287.46,1.212498
+134,22.835938,0.991984,1005.3562,1003.7906,79.04,79.04,287.64,1.21571
+135,3.757813,0.995003,1008.2239,1006.79,53.92,53.92,287.8,1.21865
+136,0.0,0.99763,1010.8487,1009.5363,30.96,30.96,287.95,1.221341
+137,0.0,1.0,1013.25,1012.0494,10.0,10.0,288.09,1.223803

pycontrails/datalib/ecmwf/variables.py

@@ -240,5 +240,17 @@ SURFACE_VARIABLES = [
     CloudAreaFraction,
     SurfaceSolarDownwardRadiation,
 ]
+MODEL_LEVEL_VARIABLES = [
+    met_var.AirTemperature,
+    met_var.SpecificHumidity,
+    met_var.VerticalVelocity,
+    met_var.EastwardWind,
+    met_var.NorthwardWind,
+    RelativeVorticity,
+    Divergence,
+    CloudAreaFractionInLayer,
+    SpecificCloudIceWaterContent,
+    SpecificCloudLiquidWaterContent,
+]
 
 ECMWF_VARIABLES = PRESSURE_LEVEL_VARIABLES + SURFACE_VARIABLES

pycontrails/models/humidity_scaling/humidity_scaling.py

@@ -591,9 +591,15 @@ class HumidityScalingByLevel(HumidityScaling):
 
 
 @functools.cache
-def _load_quantiles() -> pd.DataFrame:
+def _load_quantiles(level_type: str) -> pd.DataFrame:
     """Load precomputed ERA5 and IAGOS quantiles.
 
+    Parameters
+    ----------
+    level_type : {"pressure", "model"}
+        Select whether to load precomputed quantiles from pressure- vs
+        model-level ERA5 data.
+
     Returns
     -------
     pd.DataFrame
@@ -602,17 +608,20 @@ def _load_quantiles() -> pd.DataFrame:
         and the interpolation methodology. The IAOGS quantiles are in the
         ``("iagos", "iagos")`` column.
     """
-    path = pathlib.Path(__file__).parent / "quantiles" / "era5-quantiles.pq"
-    return pd.read_parquet(path)
+    path = pathlib.Path(__file__).parent / "quantiles" / f"era5-{level_type}-level-quantiles.pq"
+    df = pd.read_parquet(path)
+    df.attrs["path"] = str(path)
+    return df
 
 
 def histogram_matching(
     era5_rhi: ArrayLike,
     product_type: str,
+    level_type: str,
     member: int | None,
     q_method: str | None,
 ) -> npt.NDArray[np.float64]:
-    """Map ERA5-derived RHi to it's corresponding IAGOS quantile via histogram matching.
+    """Map ERA5-derived RHi to its corresponding IAGOS quantile via histogram matching.
 
     This matching is performed on a **single** ERA5 ensemble member.
 
@@ -622,6 +631,11 @@ def histogram_matching(
         ERA5-derived RHi values for the given ensemble member.
     product_type : {"reanalysis", "ensemble_members"}
         The ERA5 product type.
+    level_type : {"pressure", "model"}
+        Select whether to perform quantile mapping based on quantiles from
+        pressure- or model-level ERA5 data. Selecting ``level_type == "model"``
+        when ``product_type == "ensemble_members"`` will produce a warning
+        and change ``product_type`` to  ``"reanalysis"``.
     member : int | None
         The ERA5 ensemble member to use. Must be in the range ``[0, 10)``.
         Only used if ``product_type == "ensemble_members"``.
@@ -634,18 +648,30 @@ def histogram_matching(
         The IAGOS quantiles corresponding to the ERA5-derived RHi values. Returned
         as a numpy array with the same shape and dtype as ``era5_rhi``.
     """
-    df = _load_quantiles()
+    if level_type not in ["pressure", "model"]:
+        msg = f"Invalid 'level_type' value '{level_type}'. " "Must be one of ['pressure', 'model']."
+        raise ValueError(msg)
+    df = _load_quantiles(level_type)
     iagos_quantiles = df[("iagos", "iagos")]
 
+    if product_type == "ensemble_members" and level_type == "model":
+        msg = (
+            "No quantiles available for model-level ensemble data. "
+            "Switching to product_type = 'reanalysis'."
+        )
+        warnings.warn(msg)
+        product_type = "reanalysis"
+
     if product_type == "ensemble_members":
         col = f"ensemble{member}", q_method or "linear-q"
     elif product_type == "reanalysis":
         col = "reanalysis", q_method or "linear-q"
     else:
-        raise ValueError(
+        msg = (
             f"Invalid 'product_type' value '{product_type}'. "
             "Must be one of ['reanalysis', 'ensemble_members']."
         )
+        raise ValueError(msg)
 
     try:
         era5_quantiles = df[col]
@@ -691,7 +717,7 @@ def histogram_matching_all_members(
 
     # Perform histogram matching on the given ensemble member
     ensemble_member_rhi = histogram_matching(
-        era5_rhi_all_members[:, member], "ensemble_members", member, q_method
+        era5_rhi_all_members[:, member], "ensemble_members", "pressure", member, q_method
     )
 
     # Perform histogram matching on all other ensemble members
@@ -702,7 +728,7 @@ def histogram_matching_all_members(
             ensemble_mean_rhi += ensemble_member_rhi
         else:
             ensemble_mean_rhi += histogram_matching(
-                era5_rhi_all_members[:, r], "ensemble_members", r, q_method
+                era5_rhi_all_members[:, r], "ensemble_members", "pressure", r, q_method
             )
 
     # Divide by the number of ensemble members to get the mean
@@ -766,6 +792,9 @@ class HistogramMatchingParams(models.ModelParams):
     #: The ERA5 product. Must be one of ``"reanalysis"`` or ``"ensemble_members"``.
     product_type: str = "reanalysis"
 
+    #: The ERA5 vertical level type. Must be one of ``"pressure"`` or ``"model"``.
+    level_type: str = "pressure"
+
     #: The ERA5 ensemble member to use. Must be in the range ``[0, 10)``.
     #: Only used if ``product_type`` is ``"ensemble_members"``.
     member: int | None = None
@@ -779,6 +808,23 @@ class HistogramMatching(HumidityScaling):
     formula = "era5_quantiles -> iagos_quantiles"
     default_params = HistogramMatchingParams
 
+    def __init__(
+        self,
+        met: MetDataset | None = None,
+        params: dict[str, Any] | None = None,
+        **params_kwargs: Any,
+    ):
+        if (params is None or "level_type" not in params) and (
+            params_kwargs is None or "level_type" not in params_kwargs
+        ):
+            msg = (
+                "The default level_type will change from 'pressure' to 'model' "
+                "in a future release. To silence this warning, "
+                "provide a 'level_type' value when instantiating HistogramMatching."
+            )
+            warnings.warn(msg, DeprecationWarning)
+        super().__init__(met, params, **params_kwargs)
+
     @overrides
     def scale(
         self,
@@ -793,6 +839,7 @@ class HistogramMatching(HumidityScaling):
         rhi_1 = histogram_matching(
             rhi,
             self.params["product_type"],
+            self.params["level_type"],
             self.params["member"],
             self.params["interpolation_q_method"],
         )

pycontrails/physics/thermo.py

@@ -66,7 +66,7 @@ def c_pm(q: ArrayScalarLike) -> ArrayScalarLike:
     Some models (including CoCiP) use a constant value here (1004 :math:`J \ kg^{-1} \ K^{-1}`)
 
     """
-    return constants.c_pd * (1.0 + (constants.c_pv / constants.c_pd + 1.0) * q)
+    return constants.c_pd * (1.0 - q) + constants.c_pv * q
 
 
 def p_vapor(q: ArrayScalarLike, p: ArrayScalarLike) -> ArrayScalarLike:

pycontrails/utils/types.py

@@ -58,7 +58,6 @@ def support_arraylike(
     - `numpy ufuncs <https://numpy.org/doc/stable/reference/ufuncs.html>`_
     """
 
-    @functools.wraps(func)
     def wrapped(arr: ArrayScalarLike) -> ArrayScalarLike:
         x = np.asarray(arr)
 
@@ -86,7 +85,9 @@ def support_arraylike(
         # Pass numpy `ret` through for anything else
         return ret
 
-    return wrapped
+    # this line produces a mypy error starting on mypy version 1.1.0,
+    # likely due to changes in https://github.com/python/mypy/pull/16942
+    return functools.update_wrapper(wrapped, func)  # type: ignore
 
 
 def apply_nan_mask_to_arraylike(arr: ArrayLike, nan_mask: np.ndarray) -> ArrayLike:

{pycontrails-0.50.2.dist-info → pycontrails-0.51.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: pycontrails
-Version: 0.50.2
+Version: 0.51.0
 Summary: Python library for modeling aviation climate impacts
 Author-email: Breakthrough Energy <py@contrails.org>
 License: Apache-2.0
@@ -36,7 +36,7 @@ Requires-Dist: xarray >=2022.3
 Provides-Extra: complete
 Requires-Dist: pycontrails[ecmwf,gcp,gfs,goes,jupyter,pyproj,vis,zarr] ; extra == 'complete'
 Provides-Extra: dev
-Requires-Dist: black[jupyter] ==24.2.0 ; extra == 'dev'
+Requires-Dist: black[jupyter] ==24.4.1 ; extra == 'dev'
 Requires-Dist: dep-license ; extra == 'dev'
 Requires-Dist: fastparquet >=0.8 ; extra == 'dev'
 Requires-Dist: ipdb >=0.13 ; extra == 'dev'
@@ -47,10 +47,10 @@ Requires-Dist: platformdirs >=3.0 ; extra == 'dev'
 Requires-Dist: pre-commit >=2.10 ; extra == 'dev'
 Requires-Dist: psutil ; extra == 'dev'
 Requires-Dist: pyarrow >=5.0 ; extra == 'dev'
-Requires-Dist: pytest >=8.1 ; extra == 'dev'
+Requires-Dist: pytest >=8.2 ; extra == 'dev'
 Requires-Dist: pytest-cov >=2.11 ; extra == 'dev'
 Requires-Dist: requests >=2.25 ; extra == 'dev'
-Requires-Dist: ruff ==0.1.15 ; extra == 'dev'
+Requires-Dist: ruff ==0.4.1 ; extra == 'dev'
 Provides-Extra: docs
 Requires-Dist: doc8 >=1.1 ; extra == 'docs'
 Requires-Dist: furo >=2023.3 ; extra == 'docs'

{pycontrails-0.50.2.dist-info → pycontrails-0.51.0.dist-info}/RECORD

@@ -1,32 +1,36 @@
 pycontrails/__init__.py,sha256=8WUs6hZAAIH1yKfwYJ8UEqNsk5voRyLmmD3Dje9DZaE,2055
-pycontrails/_version.py,sha256=Z3flk-6qKL4selVxji0mrZLN1YfKORVEfpHjsiUr8UQ,429
+pycontrails/_version.py,sha256=oEDnT0Uuy9HWKUTC_b3rQmjMj8G7tQ7o2sa9xwltcUo,429
 pycontrails/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 pycontrails/core/__init__.py,sha256=eNypVTz1kHBSKAJX3CgfKw-VKrMRRkKTutmjSlrfeUs,870
 pycontrails/core/aircraft_performance.py,sha256=ikeJmdvFRDa1RdfR-JKfhQbiiIzL0c2vzcBmobmoxMs,22511
 pycontrails/core/airports.py,sha256=nGKXN3jOtzsDCaJZVFNO3e3w-U3lqMTz5Ww5jALiRJY,6984
 pycontrails/core/cache.py,sha256=NNrFNUvmDEkeS0d1VpDNWmF225a8u2d-TGcpRgQd9Qs,28838
 pycontrails/core/coordinates.py,sha256=cb8RprpoSgRTFbAXTPNfuUHVnOxyV3zZ0Ac88P5YbBw,5465
-pycontrails/core/datalib.py,sha256=u2dNc8HxkeJd1jMPTZsOTEagTI0wQkzF1ComctmxJO8,23993
+pycontrails/core/datalib.py,sha256=j7fGMAIdvcg4caq_t2SHvHphgk_yfVM2lQVjBrdnhuA,24716
 pycontrails/core/fleet.py,sha256=WKF_s_gRXHmB9b1OW7RUkM1TzfvVD8Ab0Md-qKRwkzs,16544
-pycontrails/core/flight.py,sha256=B7yOtmQgyu-8v5je4faWVh9XOm_2Yb9UvJJaQ-fpj9c,82183
+pycontrails/core/flight.py,sha256=0fj3DPkFnTrbmizXgSRDGkx9saRadQrWt8znO7IpkL4,83693
 pycontrails/core/flightplan.py,sha256=cpMZ6VCYbfwh3vnew2XgVEHnqBx1NzeAhrTVCvlbbss,7569
 pycontrails/core/fuel.py,sha256=06YUDhvC8Rx6KbUXRB9qLTsJX2V7tLbzjwAfDH0R6l8,4472
 pycontrails/core/interpolation.py,sha256=Sp9s17i1QQx3jZwnwvt3vo6enWwlkYwTVKCE27N2Wdk,26214
-pycontrails/core/met.py,sha256=J7m0fJNJrFfcXP4d0ctf_ZK4eDa2kT9kWFHfwnEECo8,94924
+pycontrails/core/met.py,sha256=l0A1MEFrOUtYWO_ZHoHzYgqa1rmIh7fxA8PiwE26yT8,96251
 pycontrails/core/met_var.py,sha256=JzB7UhBLQyU4TuKZqemhpBHA6Dbt89BPYO2sYBLMkL4,9504
 pycontrails/core/models.py,sha256=VS-ct4xkojJIuqdPpT1ke1ZetNzv10nNx_Z_XalZyeo,40175
-pycontrails/core/polygon.py,sha256=ukcYh4TzGxz-ggipYe1m6DR50aEBgQhxWkEtjoTqC3w,18227
-pycontrails/core/rgi_cython.cp312-win_amd64.pyd,sha256=DGfxzri-C_BWh5eRvXDO4PNNSNfDzZvIPaO3C4oQ3po,264704
+pycontrails/core/polygon.py,sha256=iliIQv-EUpstP3nOEqgpqLhztW4UI_A4NeCGQuj0Ku8,18542
+pycontrails/core/rgi_cython.cp312-win_amd64.pyd,sha256=2sy8MF_JVEIUvP7_GHHzIEB45DFWxXuBnicFa_GnLLk,264704
 pycontrails/core/vector.py,sha256=MF0oWX0Ghp_G41A712VbU1GEKRNU_Pj9PtzMREPG5Z8,73916
 pycontrails/datalib/__init__.py,sha256=WnXqgv20SrHZLjFZ9qpQEwRnf0QvfZU-YvMqh7PAUwg,246
 pycontrails/datalib/goes.py,sha256=NVk6QFfHbGpOINlKlG6PR1Hn_LCmzP7k20VJQZXpVG8,27124
-pycontrails/datalib/ecmwf/__init__.py,sha256=MZBpXU2ZD4GOJX9Pr43_N-8TLF-KtV6aZQ5hADmGFCE,1262
-pycontrails/datalib/ecmwf/arco_era5.py,sha256=aWI69KhTfR42mZYssXCrB43lJelcV85GC8FPZDPxMak,20835
-pycontrails/datalib/ecmwf/common.py,sha256=XWPlWPQskuTMHBxWeKvRSuYWHRcKeiJAxfl5lScyTPc,3806
-pycontrails/datalib/ecmwf/era5.py,sha256=r4RVTfrl_3gDwf84oo1McYOn3HBwYtD6OUrloRRq6ME,18826
+pycontrails/datalib/ecmwf/__init__.py,sha256=a3prI136k2NCxCicLMWznvr52cTx6JjHpmV6ts4OGxI,1508
+pycontrails/datalib/ecmwf/arco_era5.py,sha256=kFNmci84-tx7AJDKF3Kv0BQja6L-QJ4M1RAsuErraHQ,19094
+pycontrails/datalib/ecmwf/common.py,sha256=_MQkB0qvrqMhLlGj4F4MXXi2YC_kuijeQHIX365jaYY,3940
+pycontrails/datalib/ecmwf/era5.py,sha256=D-AjF1f8nusL1NN5bfMgDWkbyZiNAfTtYWw8nlsEryA,18715
+pycontrails/datalib/ecmwf/era5_model_level.py,sha256=Ve2enkyoARQkG5CSD9o8CYw1ageJiU31aVAeoDojwW0,19362
 pycontrails/datalib/ecmwf/hres.py,sha256=agGc9Oope2yd0u_7nH4RuxOHkiU6dc_JvxAxdSIAoCY,28953
+pycontrails/datalib/ecmwf/hres_model_level.py,sha256=_A2QMEWOj8BairM30EpkuW_p0NMlWm93pxC8X6abnPk,19559
 pycontrails/datalib/ecmwf/ifs.py,sha256=PTdP-sRUYPAxDmAhxMDffdfUH4H3qMpTWahu4mn87GA,10877
-pycontrails/datalib/ecmwf/variables.py,sha256=gkXy3g7LOIxhp7RGGxu6MU26NAhYHpd2ewAmbmuljEo,9540
+pycontrails/datalib/ecmwf/model_levels.py,sha256=ZqgGmCDp_vcYNsTqPqUmYz6pPtILBEglWU7lMfynlO8,2791
+pycontrails/datalib/ecmwf/variables.py,sha256=vJIaxDM1Yg0-2hCR9BbD3NcJENk-Mh81NCnkzCh58tA,9861
+pycontrails/datalib/ecmwf/static/model_level_dataframe_v20240418.csv,sha256=9u7CVA3QnPUmNLIWUkF5b9wFunczkvx1zSudwGmtOv8,9927
 pycontrails/datalib/gfs/__init__.py,sha256=VcE2j62ITbY7F3tEtgaLrfyjHWci-4mvLtnVg3SVgtE,712
 pycontrails/datalib/gfs/gfs.py,sha256=u8IJTXiA8Gx1CTKYPzHkopVi_BBSuUqLItTBpzmC0E8,22213
 pycontrails/datalib/gfs/variables.py,sha256=u2rAsjCJTehB6MhQj9YhlZYovdwoscqAH33T5aBDzSw,2961
@@ -66,8 +70,9 @@ pycontrails/models/emissions/static/default-engine-uids.csv,sha256=6e-0Fjbka1www
 pycontrails/models/emissions/static/edb-gaseous-v29b-engines.csv,sha256=Np_ej6WYFo7PGYtvIt_IAcEVCwTTbNNImKSYQ3Lhcpg,128114
 pycontrails/models/emissions/static/edb-nvpm-v29b-engines.csv,sha256=MwLLrcATd38KPddTpHpMGBrZuA4I7he-1B5otTp4ar8,77533
 pycontrails/models/humidity_scaling/__init__.py,sha256=-xqDCJzKJx2nX6yl-gglHheQHWDhkvb8X7atbMJT2LA,1156
-pycontrails/models/humidity_scaling/humidity_scaling.py,sha256=0e2hOpVDK36ghKqBNjJo19mSlkR1SJq1qtb48NX5Coo,35692
-pycontrails/models/humidity_scaling/quantiles/era5-quantiles.pq,sha256=tfYhbafF9Z-gGCg6VQ1YBlOaK_01e65Dc6s9b-hQ6Zo,286375
+pycontrails/models/humidity_scaling/humidity_scaling.py,sha256=WEe-0rMTJFPzbsXuHVHLeUCHr552C73TlxIUhdKCOmA,37683
+pycontrails/models/humidity_scaling/quantiles/era5-model-level-quantiles.pq,sha256=pShCvNUo0NYtAHhT9IBRuj38X9jejdlKfv-ZoOKmtKI,35943
+pycontrails/models/humidity_scaling/quantiles/era5-pressure-level-quantiles.pq,sha256=tfYhbafF9Z-gGCg6VQ1YBlOaK_01e65Dc6s9b-hQ6Zo,286375
 pycontrails/models/ps_model/__init__.py,sha256=QggqLRpqUh6imcn7GFPcKFSU4s3WjgfdvO8hH_OO8NY,512
 pycontrails/models/ps_model/ps_aircraft_params.py,sha256=QP6jcK0L1X28GwuVKX8WS0kEaaB72MAS3kMFzcmaJao,12911
 pycontrails/models/ps_model/ps_grid.py,sha256=oAjo_SY9R-vyDsaTAQpHl4yP1VLjQi5ApwaiNdFT_M4,19172
@@ -78,17 +83,17 @@ pycontrails/physics/__init__.py,sha256=AScCMSMSZjKxfL6mssdSLwcja1ml7MzREThQp5PLr
 pycontrails/physics/constants.py,sha256=SWG7H7eJCvQXfUR3qS6_fYzNvEeRZga50qT2RuaHoYU,3262
 pycontrails/physics/geo.py,sha256=Dwc9dUzUlHGCDVqFdh492vJZsC5CB8J45cJnsHaMRF4,31229
 pycontrails/physics/jet.py,sha256=iQ0g5cocj6j0JZPSxpTF3LLypSZGx3QQc6UyXsFfJGc,26461
-pycontrails/physics/thermo.py,sha256=mwtsee1nzlIP6OUdPPMrlrLTHoWqhdm6Qxdb9xS50R0,13261
+pycontrails/physics/thermo.py,sha256=5S8QvjaGDPR5vCbDsAsvwmplJPhTeDnITjCEU_sGMRY,13240
 pycontrails/physics/units.py,sha256=GnCJCKz25jcEYGcQVWfuYViVha6mukf08Eo4pK_3iYo,12744
 pycontrails/utils/__init__.py,sha256=VmklFC-5I5lGFQEzuomlPk_bM6CoM9XDljfjCovG3vw,33
 pycontrails/utils/dependencies.py,sha256=SjEdbDDKfGmmYResWZndMCUySO0W0ptWAeY1aA_kcx8,2625
 pycontrails/utils/iteration.py,sha256=En2YY4NiNwCNtAVO8HL6tv9byBGKs8MKSI7R8P-gZy4,332
 pycontrails/utils/json.py,sha256=xCv71CKVZNHk4MyoYC-hl7dXObXXbI7P8gcNCn3AUoU,6172
 pycontrails/utils/temp.py,sha256=5XXqQoEfWjz1OrhoOBZD5vkkCFeuq9LpZkyhc38gIeY,1159
-pycontrails/utils/types.py,sha256=gf374r1iY_DxjWITbllF9m4-6zP91uP4gLqYy7tbsGA,4772
-pycontrails-0.50.2.dist-info/LICENSE,sha256=HVr8JnZfTaA-12BfKUQZi5hdrB3awOwLWs5X_ga5QzA,10353
-pycontrails-0.50.2.dist-info/METADATA,sha256=uMuGB8pL2WBIsW3tSYNGSxUH_qlRxPHLi6NXk160ilY,8528
-pycontrails-0.50.2.dist-info/NOTICE,sha256=qYeNEp8OjDK5jSW3hTlr9LQRjZeEhXQm0zDei5UFaYs,1969
-pycontrails-0.50.2.dist-info/WHEEL,sha256=fZWyj_84lK0cA-ZNCsdwhbJl0OTrpWkxInEn424qrSs,102
-pycontrails-0.50.2.dist-info/top_level.txt,sha256=Z8J1R_AiBAyCVjNw6jYLdrA68PrQqTg0t3_Yek_IZ0Q,29
-pycontrails-0.50.2.dist-info/RECORD,,
+pycontrails/utils/types.py,sha256=gNG9cSZ3djW7jufg0h1fXM3kD24sBY6ENE6wsxY_Q6o,4937
+pycontrails-0.51.0.dist-info/LICENSE,sha256=HVr8JnZfTaA-12BfKUQZi5hdrB3awOwLWs5X_ga5QzA,10353
+pycontrails-0.51.0.dist-info/METADATA,sha256=hh94W9zQKbUU4aH_DBn9RysO5Dr-nkw14-tyDhMkbBg,8527
+pycontrails-0.51.0.dist-info/NOTICE,sha256=qYeNEp8OjDK5jSW3hTlr9LQRjZeEhXQm0zDei5UFaYs,1969
+pycontrails-0.51.0.dist-info/WHEEL,sha256=fZWyj_84lK0cA-ZNCsdwhbJl0OTrpWkxInEn424qrSs,102
+pycontrails-0.51.0.dist-info/top_level.txt,sha256=Z8J1R_AiBAyCVjNw6jYLdrA68PrQqTg0t3_Yek_IZ0Q,29
+pycontrails-0.51.0.dist-info/RECORD,,