def get_ice_thickness(fjord):
    """
    Get the approximate ice thickness (in m) at the grounding line for a given fjord
    """

    # Note: these values are approximated manually by inspecting the ice thickness layer in BedMachine v3 ~10 km in from the terminus
    # and choosing a value at the upper end of gridded values.
    thickness = {"JI": 1500, "KB": 800}

    try:
        return thickness.pop(fjord)
    except KeyError:
        print("The current fjord does not have an ice thickness entry - using a default value!")
        return 1500