pycontrails 0.50.2__cp311-cp311-macosx_11_0_arm64.whl → 0.51.0__cp311-cp311-macosx_11_0_arm64.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,38 +1,42 @@
-pycontrails-0.50.2.dist-info/RECORD,,
-pycontrails-0.50.2.dist-info/LICENSE,sha256=gJ-h7SFFD1mCfR6a7HILvEtodDT6Iig8bLXdgqR6ucA,10175
-pycontrails-0.50.2.dist-info/WHEEL,sha256=sieEctgmsyAnWfDYOiunmkigyyjGmYuUaApm_YItwoI,110
-pycontrails-0.50.2.dist-info/NOTICE,sha256=gKI8DcN1WhiXB2SFRKDogcjONldGubTvBxiOYdC4CXU,1926
-pycontrails-0.50.2.dist-info/top_level.txt,sha256=Z8J1R_AiBAyCVjNw6jYLdrA68PrQqTg0t3_Yek_IZ0Q,29
-pycontrails-0.50.2.dist-info/METADATA,sha256=RSGkFvfH8frp65X3U_UldVqHY50lwg3A6w-MyHDVkGQ,8367
-pycontrails/_version.py,sha256=UlHN6EjrMB7VOwF7Q-DjBtQVs_CzfLcy0S2X4HVxI-U,413
+pycontrails-0.51.0.dist-info/RECORD,,
+pycontrails-0.51.0.dist-info/LICENSE,sha256=gJ-h7SFFD1mCfR6a7HILvEtodDT6Iig8bLXdgqR6ucA,10175
+pycontrails-0.51.0.dist-info/WHEEL,sha256=sieEctgmsyAnWfDYOiunmkigyyjGmYuUaApm_YItwoI,110
+pycontrails-0.51.0.dist-info/NOTICE,sha256=gKI8DcN1WhiXB2SFRKDogcjONldGubTvBxiOYdC4CXU,1926
+pycontrails-0.51.0.dist-info/top_level.txt,sha256=Z8J1R_AiBAyCVjNw6jYLdrA68PrQqTg0t3_Yek_IZ0Q,29
+pycontrails-0.51.0.dist-info/METADATA,sha256=ffNiyz3UxTW2Snj6wwHdKIse5R-GTWYebWi6Gx3ZJ3Y,8366
+pycontrails/_version.py,sha256=PduUz-IZe7h3dOvDzyQFqtIuB8bpKtQaiE1lK2u0IVU,413
 pycontrails/__init__.py,sha256=c_Vtz7CvdiVAL8ggluash9-8tGcLO_5Vvu-3_Ie47CE,1985
 pycontrails/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-pycontrails/core/rgi_cython.cpython-311-darwin.so,sha256=iI508Soe4v04Dhb2J-RaGpS4JsvaIjHANHXlIqJE4u4,317693
+pycontrails/core/rgi_cython.cpython-311-darwin.so,sha256=tu0YjTRoOXf5JKEf-AqDWrt55hzIdcqzEIcrh9SHKBg,310896
 pycontrails/core/vector.py,sha256=P61Rv1M040A1SRHmdWVR_RfXcWa6faM330n-9499ilI,71722
 pycontrails/core/models.py,sha256=sipUAH_r0A1TxQ0yBjaPj8WC33Piw34ABHyrmRePfhg,39003
 pycontrails/core/interpolation.py,sha256=FFYdUnTzGnkjoGs-FKsK3Y3nSK5lQI-mpq5HPALZ-y4,25495
 pycontrails/core/fleet.py,sha256=84oTb8RJ3-bPVvZn3O2ljMEZLwJ9Q-E455ZzQJ4QaQU,16075
-pycontrails/core/flight.py,sha256=zjtyP9iaS6ZdVQhtvj56lujQ-JG7r6-3hRPLmThnFf0,79999
+pycontrails/core/flight.py,sha256=tkR7p8ry5Gxd6Das-84Ny3ClXApwbcVFiqkfihy63T4,81467
 pycontrails/core/fuel.py,sha256=kJZ3P1lPm1L6rdPREM55XQ-VfJ_pt35cP4sO2Nnvmjs,4332
-pycontrails/core/polygon.py,sha256=eX0gLL6OA_3658hjT0uraoyQWn4WvKYS-uxh3YCkTZY,17685
-pycontrails/core/datalib.py,sha256=tqYR69vVeY6LsKgtOYkZU1sEJNrbrwjqRPC1eK8mQaU,23270
+pycontrails/core/polygon.py,sha256=bLddqnsaA7XVXFcEVXR3V1tVIu5qQ8sxPUjFbwt0lO8,17993
+pycontrails/core/datalib.py,sha256=m39gKxdhy9KSBhwWlrzEc5z-tTqPfcs7bPcqSxbU0k0,23971
 pycontrails/core/cache.py,sha256=rCBZiRSVoFkRwf_ay7O-eba2PaYoXdlJSz1znZPuOQk,27957
 pycontrails/core/__init__.py,sha256=4ZE7x1gMa_Q7GcgcWpPP9fQ-sTulCXWmMjsCJ0odakY,840
 pycontrails/core/flightplan.py,sha256=s7tHbjMFbHAJkSWV6hPkghuW6jDb1n5UhWAo9XbJ9z0,7349
-pycontrails/core/met.py,sha256=1TGotmJK6Ks_hmC8QGNknFuAWoNLiSFQuNAveCLcPyk,92303
+pycontrails/core/met.py,sha256=HORTWNre9rSyOgNnVnVaW0A6jDjU5W5A5s2L89KIZyM,93602
 pycontrails/core/aircraft_performance.py,sha256=4yNGfArt741HNFjxpWvEu86BTzie9LaPjC4KH3dXYts,21966
 pycontrails/core/airports.py,sha256=aeyAXVkioIRomrP79UtNrxindL4f1DJyXFaojZCuBBw,6758
 pycontrails/core/met_var.py,sha256=EhrLGdrCAp8cKb-3Whd_ttLMZn4_zLMhE-QyFinESqo,9197
 pycontrails/core/coordinates.py,sha256=vVITA90x0Jx-UQG2XMm3JAKKsIrFUcU861xV-L9czTI,5291
 pycontrails/datalib/goes.py,sha256=Myz53yENqopOOsVT7K-M16SLvoBPPnOg-14fq6DvzdQ,26359
 pycontrails/datalib/__init__.py,sha256=s5b8W6scXgespKZjPcaC-8jyLd3FqZqN0BZccLxaGuE,237
-pycontrails/datalib/ecmwf/arco_era5.py,sha256=wwBIXZxRRV98g42kNYalirsn9hmq3QBSD_4zdoIEwEg,20258
-pycontrails/datalib/ecmwf/era5.py,sha256=dGEqXtjuE0OgWxOtkMmBtegeMmo1n6nb2cUrNkXx7HM,18285
+pycontrails/datalib/ecmwf/arco_era5.py,sha256=6UCy75rPTJZokzQK6jccxLy6nOPBbAP8Z-6QbN_eMsY,18568
+pycontrails/datalib/ecmwf/era5.py,sha256=_TV0g24HVmU93fNO8ykvhWaaHSwnc2QT3cFHhhmsSjU,18178
+pycontrails/datalib/ecmwf/era5_model_level.py,sha256=R_Z30DvJU7i7_sPgYV5GDhHBPIKd2ChoePYfWP0rTYw,18881
 pycontrails/datalib/ecmwf/hres.py,sha256=yJapluaSaOrHw23iUc-OIIW2Jd7iVLlDcvMSQpd7ToU,28172
-pycontrails/datalib/ecmwf/variables.py,sha256=gihPeicRzaaxinjjL_tWzaQqg2m0ixba8xplwagnGV8,9296
-pycontrails/datalib/ecmwf/__init__.py,sha256=wQ2c-wEAQvOdAgHYAK0qj2FJkh1Jfj0J3kmcytvYtG4,1215
-pycontrails/datalib/ecmwf/common.py,sha256=0D2dQBi9mBOQnxog_A95yjkQ2brG0z0__teFc1ghDMU,3702
+pycontrails/datalib/ecmwf/variables.py,sha256=S5SV_kPrWRhwOpOYcWw3V5VRUB_HOVCTw3Sf7UClJEk,9605
+pycontrails/datalib/ecmwf/hres_model_level.py,sha256=CBrtZnX9Hbu1xBk8ALypLDUycg_7o30wJS-G3VELo2Q,19065
+pycontrails/datalib/ecmwf/__init__.py,sha256=kwgk9P4RYdLgOYcBLXX5rWz1T_yL7aO8nt2Eb1PB-eE,1455
+pycontrails/datalib/ecmwf/common.py,sha256=veOzDmAjLJcfSRQrQNl4hz8g4rusOEl8Ph9GBwoQOok,3832
+pycontrails/datalib/ecmwf/model_levels.py,sha256=x83WJtjC6OnHcUsiNgvYIrVX4lY-pkXR-YlvUo9vYis,2712
 pycontrails/datalib/ecmwf/ifs.py,sha256=s__NOQcUfbvxVPyCE7qgO8cxYLJoh7PWj7zZS-nKR7g,10594
+pycontrails/datalib/ecmwf/static/model_level_dataframe_v20240418.csv,sha256=PmvGLRzn6uuCKSwiasSuVcehvvmSaqP7cnLuN6hhCQQ,9788
 pycontrails/datalib/gfs/gfs.py,sha256=wOvOl1d-7B024uaaak9gnXg99Lbnv7yhE00ajSGRw6U,21569
 pycontrails/datalib/gfs/variables.py,sha256=KESzB1sTD3hsU8T-qZFD29oFM3l2ZcEtjAE_H7BHKIE,2861
 pycontrails/datalib/gfs/__init__.py,sha256=tWxgqmlW8Uo07J-3fBTXPrteatzTka9mSXomhWy3NVA,684
@@ -44,7 +48,7 @@ pycontrails/ext/empirical_grid.py,sha256=WSC266aKsQLzCmtrZJCpLdDBykZ9rlFE9xEXmZj
 pycontrails/ext/bada.py,sha256=j4Tj7oWSV_6UxYYa9_OjC1yTVzJMQdNRDI4aUQam_xM,1063
 pycontrails/utils/iteration.py,sha256=q_vb39VjxRr4hqTyPYko3gK4sboJOJf_Evq6m_2DL-g,319
 pycontrails/utils/__init__.py,sha256=Gt_57sBgfliFSxx9sDpuchykFDxmM11Wg9xAeSqPcnI,32
-pycontrails/utils/types.py,sha256=rMKWRDfQCGkh6QL0x3TZhReemTL9u-DmYnoqcpx8PII,4609
+pycontrails/utils/types.py,sha256=skrZuOB1yQivtg7rUi8bZwBzOD9ZINu93Cx0NaLzCkw,4773
 pycontrails/utils/temp.py,sha256=lGU0b_R8ze4yKlsOusHIIBaoNFBrmrB3vBjgHRlfcXk,1109
 pycontrails/utils/json.py,sha256=4PL7xkt-O9VGvngHh8F1nfbSbgoxjNlDZDUQRxYfcCM,5982
 pycontrails/utils/dependencies.py,sha256=ATP45xYdUbIyGFzgbOe5SbokMytvB84TcexUEFnEUZE,2559
@@ -63,9 +67,10 @@ pycontrails/models/emissions/black_carbon.py,sha256=F2SCUiV39zg2mUxbWsct6vvr_JgH
 pycontrails/models/emissions/static/edb-nvpm-v29b-engines.csv,sha256=NatpVI1D2tTDLK7uVvlanm9DhfFB44nmFA4aocUcXco,77318
 pycontrails/models/emissions/static/edb-gaseous-v29b-engines.csv,sha256=Oub-FkyR4cPlTEPo9wDPAn4i4CpGpKH2tet30MNitI0,127518
 pycontrails/models/emissions/static/default-engine-uids.csv,sha256=3blb0aqtM8YRsyT1WDo0UYTBtv1h4BwXRIC_Ll9fhnI,6217
-pycontrails/models/humidity_scaling/humidity_scaling.py,sha256=LbkxfSjxg9RdwekVPB2x7mlzchXtupZf2LJJXQ3iao4,34714
+pycontrails/models/humidity_scaling/humidity_scaling.py,sha256=9619kSgNIBjkM8vFWcwny7t4mhg00pX0aKabrDGPn7I,36658
 pycontrails/models/humidity_scaling/__init__.py,sha256=nqsab_j9BCwMbTfCn4BjXMdhItlvNKkgUJ9-lb8RyIo,1119
-pycontrails/models/humidity_scaling/quantiles/era5-quantiles.pq,sha256=tfYhbafF9Z-gGCg6VQ1YBlOaK_01e65Dc6s9b-hQ6Zo,286375
+pycontrails/models/humidity_scaling/quantiles/era5-pressure-level-quantiles.pq,sha256=tfYhbafF9Z-gGCg6VQ1YBlOaK_01e65Dc6s9b-hQ6Zo,286375
+pycontrails/models/humidity_scaling/quantiles/era5-model-level-quantiles.pq,sha256=pShCvNUo0NYtAHhT9IBRuj38X9jejdlKfv-ZoOKmtKI,35943
 pycontrails/models/cocip/radiative_forcing.py,sha256=ERuFcYMo0_1iiOricnZ8D4ext23bMnTCeZwg9vd6Vzs,44944
 pycontrails/models/cocip/wind_shear.py,sha256=p8d3iaNzxPA3MoxFEM1ZDKt0aticoD6U9cv0QmbuBzs,3860
 pycontrails/models/cocip/cocip.py,sha256=0ILMcjbgsM00rCXAeo9UkSnUsmvq1XvORhNPI-ReNcM,97591
@@ -90,5 +95,5 @@ pycontrails/physics/geo.py,sha256=lqEpTLhex2r-o-6EHSSEiPC8xqZ-NlwwCBob9-cJA_w,30
 pycontrails/physics/units.py,sha256=j-G5AC9eWIvv2MTOq9lUOoOQKFNJJuHzWLanHRji2tE,12272
 pycontrails/physics/constants.py,sha256=pHQQmccMUwuNnY4hFtm3L8G2rnUQcfJnroyQr8HAVeM,3146
 pycontrails/physics/__init__.py,sha256=_1eWbEy6evEWdfJCEkwDiSdpiDNzNWEPVqaPekHyhwU,44
-pycontrails/physics/thermo.py,sha256=x3IsnUkEIA1zrDTFkS4GbV5B-ZKe_sxLbkCwQJ5Ign4,12810
+pycontrails/physics/thermo.py,sha256=UA7jLApcbu53fo1ACtgF578Np17efn0UN1yxwus4gEQ,12789
 pycontrails/physics/jet.py,sha256=JX-o5dDjyITMNUIOcD_7UBt_4ldc9v-5gI6bmwf0wQ4,25624