def get_utm_zone(lon):
    """ Calculate UTM zone.
    Arguments:
        lon: float
            Longitude, in degrees. West: negative, East: positive.
    Returns:
        zone: int
            UTM zone number.
    """

    zone = int(1+(lon+180.0)/6.0)

    return zone