def ha_from_lst(lst, ra):
    """Given local sidereal time (lst) and right ascension (ra),
        both in same units, (degrees, radians, or decimal hours)
        calculate the local hour angle.

        Return value is in same units as input arguments.
    """
    hourangle = lst - ra
    return hourangle