@@ -0,0 +1,498 @@
+from skyfield.api import load%0Afrom skyfield.positionlib import Geocentric%0Afrom skyfield.toposlib import Topos%0A%0Adef ts():%0A    yield load.timescale()%0A%0Adef test_beneath(ts):%0A    t = ts.utc(2018, 1, 19, 14, 37, 55)%0A    def f(xyz): return str(Topos.beneath(Geocentric(xyz, None, t)))%0A    assert f(%5B1, 0, 0%5D) == 'Topos 00deg 00%5C' 00.0%22 N 21deg 46%5C' 59.4%22 E'%0A    assert f(%5B0, 0, 1%5D) == 'Topos 90deg 00%5C' 00.0%22 N 21deg 46%5C' 59.4%22 E'%0A    assert f(%5B1, 0, 1%5D) == 'Topos 45deg 00%5C' 00.0%22 N 21deg 46%5C' 59.4%22 E'%0A

Add basic test of Topos.beneath()