def determine_layers_below_above(layers, values, elevation):
    """Determine the layers below and above the current elevation

    Args:
        layers [<str>]: All the pressure layers to load for each parameter
        values [<LayerInfo>]: All the interpolated layers information
        elevation <float>: The elevation to determine the layers for

    Returns
        <str>: Below layer ID
        <str>: Above layer ID
        <int>: Index to the above layer
    """

    index = 0
    index_below = 0
    index_above = len(layers) - 1

    for layer in layers:
        if values[layer].hgt >= elevation:
            index_below = index - 1
            index_above = index
            break

        index += 1

    # Only need to check for the low height condition
    # We should never have a height above our highest elevation
    if index_below < 0:
        index_below = 0
        index_above = 1

    return (layers[index_below], layers[index_above], index_above)