In Duffett-Smith, this is section 15.
# - - - S i d e r e a l T i m e . g s t
def gst ( self, eLong ):
"""Convert LST to GST.
[ self is local sidereal time at longitude eLong
radians east of Greenwich ->
return the equivalent GST as a SiderealTime instance ]
"""
If we convert the longitude eLong to
hours, this function becomes a simple matter of
subtracting the longitude from the LST and normalizing to
the interval [0,24).
#-- 1 --
# [ deltaHours := eLong expressed in hours ]
deltaHours = radiansToHours ( eLong )
#-- 2 --
gstHours = ( self.hours - deltaHours ) % 24.0
#-- 3 --
return SiderealTime ( gstHours )