def METHOD_NAME():
    """test the azimuth of a polygon poly"""
    data = (
        ([(0, 0, 0), (1, 0, 0), (1, 1, 1), (0, 1, 1)], 180),
        # polygon, expected,
        ([(0, 0, 0), (1, 0, 0), (1, 1, 0), (0, 1, 0)], 0),
        (
            [
                (3.571913, -9.390334, 1.487381),
                (10.905826, -6.194443, 1.487381),
                (8.998819, -1.818255, 0.0),
                (1.664906, -5.014146, 0.0),
            ],
            360 - 23.546134,
        ),
    )
    for poly, expected in data:
        result = surface.azimuth(poly)
        assert almostequal(expected, result, places=3) == True