def lookup_carbon_intensity_kg(
    iso_code: str,
    default_intensity_grams: float = 89.0,
) -> float:
    """Find the carbon intensity in kg/km of a flight given the departure country

    Data sourced from Graver et al. ICCT Report: CO2 Emissions From Commercial
    Aviation: 2013, 2018, and 2019. (Table 4)

    See: https://theicct.org/publication/co2-emissions-from-commercial-aviation-2013-2018-and-2019

    Args:
        iso_code: Flight departure country
        default_intensity_grams: Default value departure countries not explicitly
            enumerated in paper

    Returns: CO2 intensity in kg/km flown of a flight originating from departure country
    """
    intensities = {
        "AE": 89,
        "AI": 87,
        "AS": 95,
        "AU": 90,
        "BM": 87,
        "CC": 90,
        "CN": 88,
        "CX": 90,
        "DE": 91,
        "ES": 79,
        "FK": 87,
        "FR": 87,
        "GB": 87,
        "GF": 87,
        "GG": 87,
        "GI": 87,
        "GP": 87,
        "GU": 95,
        "HK": 88,
        "IM": 87,
        "IN": 85,
        "JP": 95,
        "KY": 87,
        "MP": 95,
        "MQ": 87,
        "MS": 87,
        "NC": 87,
        "NF": 90,
        "PF": 87,
        "PM": 87,
        "PR": 95,
        "RE": 87,
        "SH": 87,
        "TC": 87,
        "UM": 95,
        "US": 95,
        "VG": 87,
        "VI": 95,
        "WF": 87,
        "YT": 87,
    }
    intensity_grams = intensities.get(iso_code, default_intensity_grams)
    return intensity_grams / 1000