pycontrails 0.54.1__cp313-cp313-win_amd64.whl → 0.54.3__cp313-cp313-win_amd64.whl

pycontrails/models/ps_model/ps_operational_limits.py

@@ -350,49 +350,47 @@ def max_usable_lift_coefficient(
 
 
 def minimum_mach_num(
-    air_pressure: float,
-    aircraft_mass: float,
+    air_pressure: ArrayOrFloat,
+    aircraft_mass: ArrayOrFloat,
     atyp_param: PSAircraftEngineParams,
-) -> float:
+) -> ArrayOrFloat:
     """
     Calculate minimum mach number to avoid stall.
 
     Parameters
     ----------
-    air_pressure : float
+    air_pressure : ArrayOrFloat
         Ambient pressure, [:math:`Pa`]
-    aircraft_mass : float
+    aircraft_mass : ArrayOrFloat
         Aircraft mass at each waypoint, [:math:`kg`]
     atyp_param : PSAircraftEngineParams
         Extracted aircraft and engine parameters.
 
     Returns
     -------
-    float
-        Maximum usable lift coefficient.
+    ArrayOrFloat
+        Minimum mach number to avoid stall.
     """
 
     def excess_mass(
-        mach_number: float,
-        air_pressure: float,
-        aircraft_mass: float,
+        mach_number: ArrayOrFloat,
+        air_pressure: ArrayOrFloat,
+        aircraft_mass: ArrayOrFloat,
         mach_num_des: float,
         c_l_do: float,
         wing_surface_area: float,
-    ) -> float:
+    ) -> ArrayOrFloat:
         amass_max = max_allowable_aircraft_mass(
             air_pressure,
             mach_number,
             mach_num_des,
             c_l_do,
             wing_surface_area,
-            1e10,
+            1e10,  # clipped to this value which we want to ignore
         )
-        if amass_max < 0:
-            return np.nan
         return amass_max - aircraft_mass
 
-    m = scipy.optimize.root_scalar(
+    m = scipy.optimize.newton(
         excess_mass,
         args=(
             air_pressure,
@@ -401,21 +399,22 @@ def minimum_mach_num(
             atyp_param.c_l_do,
             atyp_param.wing_surface_area,
         ),
-        x0=0.5,
-        x1=0.6,
-    ).root
+        x0=np.full_like(air_pressure, 0.4),
+        x1=np.full_like(air_pressure, 0.5),
+        tol=1e-4,
+    )
 
     return m
 
 
 def maximum_mach_num(
-    altitude_ft: float,
-    air_pressure: float,
-    aircraft_mass: float,
-    air_temperature: float,
-    theta: float,
+    altitude_ft: ArrayOrFloat,
+    air_pressure: ArrayOrFloat,
+    aircraft_mass: ArrayOrFloat,
+    air_temperature: ArrayOrFloat,
+    theta: ArrayOrFloat,
     atyp_param: PSAircraftEngineParams,
-) -> float:
+) -> ArrayOrFloat:
     r"""
     Return the maximum mach number at the current operating conditions.
 
@@ -424,23 +423,23 @@ def maximum_mach_num(
 
     Parameters
     ----------
-    altitude_ft  : float
+    altitude_ft  : ArrayOrFloat
         Altitude, [:math:`ft`]
-    air_pressure : float
+    air_pressure : ArrayOrFloat
         Ambient pressure, [:math:`Pa`]
-    aircraft_mass : float
+    aircraft_mass : ArrayOrFloat
         Aircraft mass at each waypoint, [:math:`kg`]
-    air_temperature : npt.NDArray[np.float64]
+    air_temperature : ArrayOrFloat
         Array of ambient temperature, [:math: `K`]
-    theta : float | npt.NDArray[np.float64]
+    theta : ArrayOrFloat
         Climb (positive value) or descent (negative value) angle, [:math:`\deg`]
     atyp_param : PSAircraftEngineParams
         Extracted aircraft and engine parameters.
 
     Returns
     -------
-    float
-        Maximum usable lift coefficient.
+    ArrayOrFloat
+        Maximum mach number given thrust limiations.
     """
     # Max speed ignoring thrust limits
     mach_num_op_lim = max_mach_number_by_altitude(
@@ -451,27 +450,15 @@ def maximum_mach_num(
         atyp_param.p_inf_co,
     )
 
-    # If the max mach number ignoring thrust limits is possible, return that value
-    if (
-        get_excess_thrust_available(
-            mach_num_op_lim, air_temperature, air_pressure, aircraft_mass, theta, atyp_param
-        )
-        > 0
-    ):
-        return mach_num_op_lim
-
-    # Numerically solve for the speed where drag == max thrust
-    try:
-        m_max = scipy.optimize.root_scalar(
-            get_excess_thrust_available,
-            args=(air_temperature, air_pressure, aircraft_mass, theta, atyp_param),
-            x0=mach_num_op_lim,
-            x1=mach_num_op_lim - 0.05,
-        ).root
-    except ValueError:
-        return np.nan
-
-    return m_max
+    max_mach = scipy.optimize.newton(
+        func=get_excess_thrust_available,
+        args=(air_temperature, air_pressure, aircraft_mass, theta, atyp_param),
+        x0=mach_num_op_lim,
+        x1=mach_num_op_lim - 0.01,
+        tol=1e-4,
+    ).clip(max=mach_num_op_lim)
+
+    return max_mach
 
 
 # ----------------

pycontrails/physics/geo.py

@@ -855,6 +855,155 @@ def advect_level(
     return (level * 100.0 + (dt_s * dp_dt)) / 100.0
 
 
+def advect_longitude_and_latitude_near_poles(
+    longitude: npt.NDArray[np.floating],
+    latitude: npt.NDArray[np.floating],
+    u_wind: npt.NDArray[np.floating],
+    v_wind: npt.NDArray[np.floating],
+    dt: npt.NDArray[np.timedelta64] | np.timedelta64,
+) -> tuple[npt.NDArray[np.floating], npt.NDArray[np.floating]]:
+    r"""Advect a particle near the poles.
+
+    This function calculates the longitude and latitude of a particle after time ``dt``
+    caused by advection due to wind near the poles (above 80 degrees North and South).
+
+    Automatically wrap over the antimeridian if necessary.
+
+    Parameters
+    ----------
+    longitude : npt.NDArray[np.floating]
+        Original longitude, [:math:`\deg`]
+    latitude : npt.NDArray[np.floating]
+        Original latitude, [:math:`\deg`]
+    u_wind : npt.NDArray[np.floating]
+        Wind speed in the longitudinal direction, [:math:`m s^{-1}`]
+    v_wind : npt.NDArray[np.floating]
+        Wind speed in the latitudinal direction, [:math:`m s^{-1}`]
+    dt : npt.NDArray[np.timedelta64] | np.timedelta64
+        Advection timestep
+
+    Returns
+    -------
+    tuple[npt.NDArray[np.floating], npt.NDArray[np.floating]]
+        New longitude and latitude values, [:math:`\deg`]
+
+    Notes
+    -----
+    Near the poles, the longitude and latitude is converted to a 2-D Cartesian-like coordinate
+    system to avoid numerical instabilities and singularities caused by convergence of meridians.
+
+    See Also
+    --------
+    advect_longitude
+    advect_latitude
+    advect_horizontal
+    """
+    # Determine hemisphere sign (1 for Northern Hemisphere, -1 for Southern Hemisphere)
+    hemisphere_sign = np.where(latitude > 0.0, 1.0, -1.0)
+
+    # Convert longitude and latitude to radians
+    sin_lon_rad = np.sin(units.degrees_to_radians(longitude))
+    cos_lon_rad = np.cos(units.degrees_to_radians(longitude))
+
+    # Convert longitude and latitude to 2-D Cartesian-like coordinate system, [:math:`\deg`]
+    polar_radius = 90.0 - np.abs(latitude)
+    x_cartesian = sin_lon_rad * polar_radius
+    y_cartesian = -cos_lon_rad * polar_radius * hemisphere_sign
+
+    # Convert winds from eastward and northward direction (u, v) to (X, Y), [:math:`\deg s^{-1}`]
+    x_wind = units.radians_to_degrees(
+        (u_wind * cos_lon_rad - v_wind * sin_lon_rad * hemisphere_sign) / constants.radius_earth
+    )
+    y_wind = units.radians_to_degrees(
+        (u_wind * sin_lon_rad * hemisphere_sign + v_wind * cos_lon_rad) / constants.radius_earth
+    )
+
+    # Advect contrails in 2-D Cartesian-like plane, [:math:`\deg`]
+    dtype = np.result_type(latitude, v_wind)
+    dt_s = units.dt_to_seconds(dt, dtype)
+    x_cartesian_new = x_cartesian + dt_s * x_wind
+    y_cartesian_new = y_cartesian + dt_s * y_wind
+
+    # Convert `y_cartesian_new` back to `latitude`, [:math:`\deg`]
+    dist_squared = x_cartesian_new**2 + y_cartesian_new**2
+    new_latitude = (90.0 - np.sqrt(dist_squared)) * hemisphere_sign
+
+    # Convert `x_cartesian_new` back to `longitude`, [:math:`\deg`]
+    new_lon_rad = np.arctan2(y_cartesian_new, x_cartesian_new)
+
+    new_longitude = np.where(
+        (x_wind == 0.0) & (y_wind == 0.0),
+        longitude,
+        90.0 + units.radians_to_degrees(new_lon_rad) * hemisphere_sign,
+    )
+    # new_longitude = 90.0 + units.radians_to_degrees(new_lon_rad) * hemisphere_sign
+    new_longitude = (new_longitude + 180.0) % 360.0 - 180.0  # wrap antimeridian
+    return new_longitude, new_latitude
+
+
+def advect_horizontal(
+    longitude: npt.NDArray[np.floating],
+    latitude: npt.NDArray[np.floating],
+    u_wind: npt.NDArray[np.floating],
+    v_wind: npt.NDArray[np.floating],
+    dt: npt.NDArray[np.timedelta64] | np.timedelta64,
+) -> tuple[npt.NDArray[np.floating], npt.NDArray[np.floating]]:
+    r"""Advect a particle in the horizontal plane.
+
+    This function calls :func:`advect_longitude` and :func:`advect_latitude` when
+    the position is far from the poles (<= 80.0 degrees). When the position is near
+    the poles (> 80.0 degrees), :func:`advect_longitude_and_latitude_near_poles`
+    is used instead.
+
+    Parameters
+    ----------
+    longitude : npt.NDArray[np.floating]
+        Original longitude, [:math:`\deg`]
+    latitude : npt.NDArray[np.floating]
+        Original latitude, [:math:`\deg`]
+    u_wind : npt.NDArray[np.floating]
+        Wind speed in the longitudinal direction, [:math:`m s^{-1}`]
+    v_wind : npt.NDArray[np.floating]
+        Wind speed in the latitudinal direction, [:math:`m s^{-1}`]
+    dt : npt.NDArray[np.timedelta64] | np.timedelta64
+        Advection timestep
+
+    Returns
+    -------
+    tuple[npt.NDArray[np.floating], npt.NDArray[np.floating]]
+        New longitude and latitude values, [:math:`\deg`]
+    """
+    near_poles = np.abs(latitude) > 80.0
+
+    longitude_out = np.empty_like(longitude)
+    latitude_out = np.empty_like(latitude)
+
+    # Use simple spherical advection if position is far from the poles (<= 80.0 degrees)
+    cond = ~near_poles
+    lon_cond = longitude[cond]
+    lat_cond = latitude[cond]
+    u_wind_cond = u_wind[cond]
+    v_wind_cond = v_wind[cond]
+    dt_cond = dt if isinstance(dt, np.timedelta64) else dt[cond]
+    longitude_out[cond] = advect_longitude(lon_cond, lat_cond, u_wind_cond, dt_cond)
+    latitude_out[cond] = advect_latitude(lat_cond, v_wind_cond, dt_cond)
+
+    # And use Cartesian-like advection if position is near the poles (> 80.0 degrees)
+    cond = near_poles
+    lon_cond = longitude[cond]
+    lat_cond = latitude[cond]
+    u_wind_cond = u_wind[cond]
+    v_wind_cond = v_wind[cond]
+    dt_cond = dt if isinstance(dt, np.timedelta64) else dt[cond]
+    lon_out_cond, lat_out_cond = advect_longitude_and_latitude_near_poles(
+        lon_cond, lat_cond, u_wind_cond, v_wind_cond, dt_cond
+    )
+    longitude_out[cond] = lon_out_cond
+    latitude_out[cond] = lat_out_cond
+
+    return longitude_out, latitude_out
+
+
 # ---------------
 # Grid properties
 # ---------------

pycontrails/physics/jet.py

@@ -1,22 +1,28 @@
 """Jet aircraft trajectory and performance parameters.
 
 This module includes common functions to calculate jet aircraft trajectory
-and performance parameters, including fuel quantities, mass, thrust setting
-and propulsion efficiency.
+and performance parameters, including fuel quantities, mass, thrust setting,
+propulsion efficiency and load factors.
 """
 
 from __future__ import annotations
 
+import functools
 import logging
+import pathlib
 
 import numpy as np
 import numpy.typing as npt
+import pandas as pd
 
 from pycontrails.core import flight
 from pycontrails.physics import constants, units
 from pycontrails.utils.types import ArrayOrFloat, ArrayScalarLike
 
 logger = logging.getLogger(__name__)
+_path_to_static = pathlib.Path(__file__).parent / "static"
+PLF_PATH = _path_to_static / "iata-passenger-load-factors-20241115.csv"
+CLF_PATH = _path_to_static / "iata-cargo-load-factors-20241115.csv"
 
 
 # -------------------
@@ -43,7 +49,7 @@ def acceleration(
 
     See Also
     --------
-    :func:`flight.segment_duration`
+    pycontrails.Flight.segment_duration
     """
     dv_dt = np.empty_like(true_airspeed)
     dv_dt[:-1] = np.diff(true_airspeed) / segment_duration[:-1]
@@ -71,8 +77,8 @@ def climb_descent_angle(
 
     See Also
     --------
-    :func:`flight.segment_rocd`
-    :func:`flight.segment_true_airspeed`
+    pycontrails.Flight.segment_rocd
+    pycontrails.Flight.segment_true_airspeed
     """
     rocd_ms = units.ft_to_m(rocd) / 60.0
     sin_theta = rocd_ms / true_airspeed
@@ -319,8 +325,8 @@ def reserve_fuel_requirements(
 
     See Also
     --------
-    :func:`flight.segment_phase`
-    :func:`fuel_burn`
+    pycontrails.Flight.segment_phase
+    fuel_burn
     """
     segment_phase = flight.segment_phase(rocd, altitude_ft)
 
@@ -345,6 +351,128 @@ def reserve_fuel_requirements(
 # -------------
 
 
+@functools.cache
+def _historical_regional_load_factor(path: pathlib.Path) -> pd.DataFrame:
+    """Load the historical regional load factor database.
+
+    Daily load factors are estimated from linearly interpolating the monthly statistics.
+
+    Returns
+    -------
+    pd.DataFrame
+        Historical regional load factor for each day.
+
+    Notes
+    -----
+    The monthly **passenger load factor** for each region is compiled from IATA's monthly
+    publication of the Air Passenger Market Analysis, where the static file will be continuously
+    updated. The report estimates the regional passenger load factor by dividing the revenue
+    passenger-km (RPK) by the available seat-km (ASK).
+
+    The monthly **cargo load factor** for each region is compiled from IATA's monthly publication
+    of the Air Cargo Market Analysis, where the static file will be continuously updated.
+    The report estimates the regional cargo load factor by dividing the freight tonne-km (FTK)
+    by the available freight tonne-km (AFTK).
+    """
+    df = pd.read_csv(path, index_col="Date", parse_dates=True, date_format="%d/%m/%Y")
+    return df.resample("D").interpolate()
+
+
+AIRPORT_TO_REGION = {
+    "A": "Asia Pacific",
+    "B": "Europe",
+    "C": "North America",
+    "D": "Africa",
+    "E": "Europe",
+    "F": "Africa",
+    "G": "Africa",
+    "H": "Africa",
+    "K": "North America",
+    "L": "Europe",
+    "M": "Latin America",
+    "N": "Asia Pacific",
+    "O": "Middle East",
+    "P": "Asia Pacific",
+    "R": "Asia Pacific",
+    "S": "Latin America",
+    "T": "Latin America",
+    "U": "Asia Pacific",
+    "V": "Asia Pacific",
+    "W": "Asia Pacific",
+    "Y": "Asia Pacific",
+    "Z": "Asia Pacific",
+}
+
+
+def aircraft_load_factor(
+    origin_airport_icao: str | None = None,
+    first_waypoint_time: pd.Timestamp | None = None,
+    *,
+    freighter: bool = False,
+) -> float:
+    """
+    Estimate passenger/cargo load factor based on historical data.
+
+    Accounts for regional and seasonal differences.
+
+    Parameters
+    ----------
+    origin_airport_icao : str | None
+        ICAO code of origin airport. If None is provided, then globally averaged values will be
+        assumed at `first_waypoint_time`.
+    first_waypoint_time : pd.Timestamp | None
+        First waypoint UTC time. If None is provided, then regionally or globally averaged values
+        from the trailing twelve months will be used.
+    freighter: bool
+        Historical cargo load factor will be used if true, otherwise use passenger load factor.
+
+    Returns
+    -------
+    float
+        Passenger/cargo load factor [0 - 1], unitless
+    """
+    # If origin airport is provided, use regional load factor
+    if origin_airport_icao is not None:
+        first_letter = origin_airport_icao[0]
+        region = AIRPORT_TO_REGION.get(first_letter, "Global")
+    else:
+        region = "Global"
+
+    # Use passenger or cargo database
+    if freighter:
+        lf_database = _historical_regional_load_factor(CLF_PATH)
+    else:
+        lf_database = _historical_regional_load_factor(PLF_PATH)
+
+    # If `first_waypoint_time` is None, global/regional averages for the trailing twelve months
+    # will be assumed.
+    if first_waypoint_time is None:
+        t1 = lf_database.index[-1]
+        t0 = t1 - pd.DateOffset(months=12) + pd.DateOffset(days=1)
+        return lf_database.loc[t0:t1, region].mean().item()
+
+    date = first_waypoint_time.floor("D")
+
+    # If `date` is more recent than the historical data, then use most recent load factors
+    # from trailing twelve months as seasonal values are stable except in COVID years (2020-22).
+    if date > lf_database.index[-1]:
+        if date.month == 2 and date.day == 29:  # remove any leap day
+            date = date.replace(day=28)
+
+        filt = (lf_database.index.month == date.month) & (lf_database.index.day == date.day)
+        date = lf_database.index[filt][-1]
+
+    # (2) If `date` is before the historical data, then use 2019 load factors.
+    elif date < lf_database.index[0]:
+        if date.month == 2 and date.day == 29:  # remove any leap day
+            date = date.replace(day=28)
+
+        filt = (lf_database.index.month == date.month) & (lf_database.index.day == date.day)
+        date = lf_database.index[filt][0]
+
+    return lf_database.at[date, region].item()
+
+
 def aircraft_weight(aircraft_mass: ArrayOrFloat) -> ArrayOrFloat:
     """Calculate the aircraft weight at each waypoint.
 
@@ -413,7 +541,8 @@ def initial_aircraft_mass(
 
     See Also
     --------
-    :func:`reserve_fuel_requirements`
+    reserve_fuel_requirements
+    aircraft_load_factor
     """
     tom = operating_empty_weight + load_factor * max_payload + total_fuel_burn + total_reserve_fuel
     return min(tom, max_takeoff_weight)
@@ -463,9 +592,10 @@ def update_aircraft_mass(
 
     See Also
     --------
-    :func:`fuel_burn`
-    :func:`reserve_fuel_requirements`
-    :func:`initial_aircraft_mass`
+    fuel_burn
+    reserve_fuel_requirements
+    initial_aircraft_mass
+    aircraft_load_factor
     """
     if takeoff_mass is None:
         takeoff_mass = initial_aircraft_mass(

pycontrails/physics/static/iata-cargo-load-factors-20241115.csv

@@ -0,0 +1,71 @@
+Date,Global,Africa,Asia Pacific,Europe,Latin America,Middle East,North America
+15/12/2018,0.488,0.381,0.54,0.567,0.291,0.488,0.414
+15/1/2019,0.451,0.354,0.501,0.501,0.299,0.421,0.4
+15/2/2019,0.447,0.363,0.473,0.53,0.297,0.466,0.379
+15/3/2019,0.495,0.384,0.556,0.56,0.323,0.488,0.416
+15/4/2019,0.463,0.374,0.518,0.496,0.325,0.458,0.405
+15/5/2019,0.468,0.386,0.52,0.513,0.353,0.469,0.398
+15/6/2019,0.454,0.324,0.522,0.498,0.337,0.44,0.382
+15/7/2019,0.45,0.323,0.519,0.485,0.354,0.453,0.373
+15/8/2019,0.446,0.302,0.516,0.477,0.372,0.435,0.377
+15/9/2019,0.464,0.329,0.539,0.501,0.379,0.459,0.381
+15/10/2019,0.477,0.361,0.539,0.533,0.364,0.477,0.394
+15/11/2019,0.496,0.404,0.538,0.569,0.403,0.497,0.413
+15/12/2019,0.467,0.368,0.519,0.53,0.3,0.47,0.395
+15/1/2020,0.45,0.356,0.474,0.501,0.311,0.426,0.424
+15/2/2020,0.464,0.368,0.543,0.531,0.342,0.461,0.372
+15/3/2020,0.545,0.425,0.656,0.63,0.411,0.532,0.429
+15/4/2020,0.58,0.486,0.691,0.648,0.554,0.525,0.487
+15/5/2020,0.576,0.612,0.643,0.625,0.561,0.483,0.526
+15/6/2020,0.573,0.547,0.645,0.62,0.512,0.494,0.521
+15/7/2020,0.564,0.489,0.639,0.594,0.464,0.53,0.506
+15/8/2020,0.548,0.502,0.616,0.568,0.478,0.535,0.489
+15/9/2020,0.569,0.507,0.642,0.62,0.456,0.579,0.484
+15/10/2020,0.576,0.502,0.617,0.651,0.443,0.606,0.496
+15/11/2020,0.582,0.496,0.631,0.655,0.436,0.6,0.5
+15/12/2020,0.573,0.51,0.639,0.653,0.367,0.597,0.482
+15/1/2021,0.589,0.48,0.665,0.627,0.39,0.569,0.532
+15/2/2021,0.575,0.476,0.692,0.641,0.429,0.598,0.453
+15/3/2021,0.588,0.499,0.661,0.685,0.453,0.613,0.472
+15/4/2021,0.578,0.504,0.633,0.681,0.457,0.598,0.473
+15/5/2021,0.572,0.502,0.646,0.656,0.423,0.589,0.469
+15/6/2021,0.565,0.48,0.676,0.626,0.381,0.581,0.458
+15/7/2021,0.544,0.455,0.654,0.598,0.387,0.536,0.443
+15/8/2021,0.542,0.43,0.698,0.575,0.404,0.529,0.437
+15/9/2021,0.553,0.428,0.68,0.604,0.37,0.558,0.447
+15/10/2021,0.561,0.45,0.661,0.626,0.421,0.572,0.449
+15/11/2021,0.559,0.434,0.654,0.631,0.446,0.572,0.444
+15/12/2021,0.542,0.502,0.634,0.623,0.413,0.556,0.43
+15/1/2022,0.541,0.492,0.609,0.584,0.417,0.513,0.474
+15/2/2022,0.532,0.502,0.592,0.636,0.476,0.529,0.429
+15/3/2022,0.549,0.494,0.638,0.671,0.448,0.526,0.442
+15/4/2022,0.516,0.49,0.631,0.578,0.419,0.504,0.419
+15/5/2022,0.505,0.495,0.627,0.548,0.387,0.487,0.411
+15/6/2022,0.492,0.447,0.608,0.507,0.383,0.488,0.404
+15/7/2022,0.472,0.452,0.563,0.493,0.374,0.469,0.398
+15/8/2022,0.467,0.418,0.547,0.502,0.374,0.466,0.393
+15/9/2022,0.481,0.451,0.572,0.528,0.381,0.478,0.396
+15/10/2022,0.487,0.437,0.561,0.558,0.384,0.48,0.401
+15/11/2022,0.491,0.458,0.545,0.569,0.382,0.475,0.419
+15/12/2022,0.472,0.432,0.528,0.559,0.322,0.454,0.406
+15/1/2023,0.448,0.439,0.452,0.541,0.325,0.411,0.423
+15/2/2023,0.456,0.468,0.464,0.574,0.361,0.445,0.4
+15/3/2023,0.462,0.489,0.485,0.57,0.366,0.456,0.393
+15/4/2023,0.427,0.482,0.442,0.497,0.364,0.431,0.373
+15/5/2023,0.415,0.448,0.422,0.489,0.333,0.41,0.373
+15/6/2023,0.432,0.446,0.468,0.476,0.337,0.446,0.374
+15/7/2023,0.421,0.417,0.457,0.472,0.322,0.411,0.37
+15/8/2023,0.42,0.388,0.443,0.484,0.326,0.407,0.377
+15/9/2023,0.438,0.436,0.466,0.5,0.319,0.424,0.392
+15/10/2023,0.452,0.416,0.472,0.53,0.354,0.46,0.392
+15/11/2023,0.467,0.421,0.479,0.57,0.363,0.469,0.408
+15/12/2023,0.459,0.41,0.479,0.562,0.316,0.455,0.403
+15/1/2024,0.457,0.431,0.446,0.555,0.344,0.439,0.435
+15/2/2024,0.451,0.451,0.432,0.584,0.376,0.463,0.396
+15/3/2024,0.473,0.473,0.475,0.581,0.402,0.496,0.404
+15/4/2024,0.439,0.429,0.445,0.515,0.387,0.447,0.387
+15/5/2024,0.446,0.438,0.453,0.518,0.362,0.461,0.397
+15/6/2024,0.458,0.385,0.496,0.507,0.336,0.473,0.388
+15/7/2024,0.444,0.4,0.48,0.496,0.338,0.458,0.382
+15/8/2024,0.44,0.378,0.466,0.501,0.359,0.445,0.387
+15/9/2024,0.456,0.392,0.485,0.525,0.368,0.474,0.389

pycontrails/physics/static/iata-passenger-load-factors-20241115.csv

@@ -0,0 +1,71 @@
+Date,Global,Africa,Asia Pacific,Europe,Latin America,Middle East,North America
+15/12/2018,0.804,0.724,0.81,0.81,0.818,0.736,0.825
+15/1/2019,0.796,0.709,0.81,0.796,0.825,0.76,0.795
+15/2/2019,0.806,0.704,0.826,0.815,0.813,0.726,0.808
+15/3/2019,0.817,0.72,0.812,0.837,0.815,0.739,0.85
+15/4/2019,0.828,0.733,0.817,0.851,0.822,0.803,0.839
+15/5/2019,0.815,0.676,0.802,0.837,0.832,0.732,0.851
+15/6/2019,0.844,0.706,0.821,0.875,0.832,0.767,0.887
+15/7/2019,0.857,0.735,0.831,0.89,0.853,0.812,0.888
+15/8/2019,0.857,0.755,0.839,0.889,0.833,0.821,0.875
+15/9/2019,0.819,0.721,0.801,0.866,0.819,0.75,0.828
+15/10/2019,0.82,0.697,0.815,0.855,0.819,0.734,0.841
+15/11/2019,0.818,0.708,0.813,0.833,0.822,0.732,0.828
+15/12/2019,0.823,0.724,0.816,0.828,0.825,0.78,0.859
+15/1/2020,0.803,0.702,0.799,0.816,0.826,0.785,0.812
+15/2/2020,0.759,0.668,0.678,0.813,0.812,0.725,0.811
+15/3/2020,0.606,0.609,0.589,0.67,0.681,0.599,0.557
+15/4/2020,0.366,0.111,0.538,0.32,0.55,0.284,0.15
+15/5/2020,0.507,0.071,0.62,0.427,0.623,0.255,0.381
+15/6/2020,0.576,0.162,0.638,0.555,0.666,0.357,0.524
+15/7/2020,0.579,0.296,0.657,0.609,0.631,0.396,0.476
+15/8/2020,0.585,0.39,0.65,0.635,0.639,0.372,0.477
+15/9/2020,0.601,0.378,0.692,0.586,0.706,0.365,0.525
+15/10/2020,0.602,0.482,0.687,0.552,0.721,0.387,0.558
+15/11/2020,0.58,0.474,0.664,0.523,0.74,0.372,0.518
+15/12/2020,0.575,0.549,0.616,0.578,0.73,0.44,0.516
+15/1/2021,0.541,0.544,0.566,0.576,0.685,0.422,0.484
+15/2/2021,0.554,0.516,0.591,0.563,0.683,0.398,0.527
+15/3/2021,0.623,0.53,0.669,0.593,0.708,0.422,0.624
+15/4/2021,0.633,0.476,0.678,0.563,0.723,0.404,0.668
+15/5/2021,0.658,0.53,0.678,0.593,0.768,0.389,0.728
+15/6/2021,0.696,0.587,0.657,0.658,0.784,0.459,0.806
+15/7/2021,0.731,0.614,0.675,0.725,0.793,0.513,0.841
+15/8/2021,0.7,0.64,0.545,0.746,0.774,0.56,0.786
+15/9/2021,0.676,0.56,0.605,0.719,0.773,0.524,0.727
+15/10/2021,0.706,0.558,0.629,0.741,0.809,0.577,0.769
+15/11/2021,0.713,0.616,0.597,0.752,0.822,0.616,0.786
+15/12/2021,0.723,0.647,0.625,0.745,0.816,0.663,0.793
+15/1/2022,0.645,0.623,0.576,0.682,0.782,0.591,0.663
+15/2/2022,0.698,0.648,0.629,0.721,0.795,0.648,0.745
+15/3/2022,0.747,0.657,0.642,0.739,0.808,0.718,0.839
+15/4/2022,0.778,0.68,0.67,0.795,0.809,0.713,0.858
+15/5/2022,0.794,0.696,0.696,0.807,0.807,0.762,0.86
+15/6/2022,0.824,0.743,0.729,0.86,0.817,0.772,0.891
+15/7/2022,0.835,0.753,0.764,0.87,0.831,0.812,0.882
+15/8/2022,0.818,0.757,0.74,0.862,0.824,0.796,0.856
+15/9/2022,0.816,0.743,0.747,0.847,0.823,0.795,0.855
+15/10/2022,0.82,0.726,0.755,0.848,0.833,0.791,0.864
+15/11/2022,0.808,0.748,0.77,0.838,0.82,0.775,0.832
+15/12/2022,0.811,0.769,0.772,0.836,0.785,0.8,0.842
+15/1/2023,0.777,0.742,0.774,0.762,0.813,0.791,0.784
+15/2/2023,0.778,0.756,0.792,0.752,0.811,0.798,0.771
+15/3/2023,0.807,0.739,0.792,0.805,0.812,0.794,0.837
+15/4/2023,0.813,0.708,0.784,0.838,0.814,0.76,0.856
+15/5/2023,0.818,0.699,0.773,0.848,0.811,0.799,0.863
+15/6/2023,0.842,0.689,0.804,0.877,0.825,0.794,0.887
+15/7/2023,0.852,0.746,0.816,0.877,0.867,0.821,0.897
+15/8/2023,0.846,0.764,0.822,0.876,0.851,0.83,0.858
+15/9/2023,0.826,0.731,0.8,0.86,0.839,0.816,0.83
+15/10/2023,0.831,0.707,0.821,0.856,0.848,0.806,0.836
+15/11/2023,0.818,0.704,0.814,0.837,0.844,0.777,0.827
+15/12/2023,0.821,0.732,0.812,0.851,0.827,0.782,0.829
+15/1/2024,0.799,0.731,0.808,0.782,0.85,0.799,0.799
+15/2/2024,0.806,0.744,0.844,0.761,0.827,0.808,0.795
+15/3/2024,0.82,0.721,0.835,0.809,0.831,0.775,0.837
+15/4/2024,0.824,0.734,0.824,0.838,0.822,0.792,0.83
+15/5/2024,0.834,0.729,0.818,0.852,0.834,0.808,0.858
+15/6/2024,0.85,0.771,0.829,0.877,0.842,0.795,0.876
+15/7/2024,0.86,0.75,0.834,0.882,0.862,0.84,0.889
+15/8/2024,0.862,0.779,0.86,0.879,0.84,0.823,0.871
+15/9/2024,0.836,0.765,0.831,0.865,0.834,0.814,0.824