code: def theta_GMST1982(jd_ut1):
      
      t = (jd_ut1 - T0) / 36525.0
      g = 67310.54841 + (8640184.812866 + (0.093104 + (-6.2e-6) * t) * t) * t
      dg = 8640184.812866 + (0.093104 * 2.0 + (-6.2e-6 * 3.0) * t) * t
      theta = (jd_ut1 % 1.0 + g * _second % 1.0) * tau
      theta_dot = (1.0 + dg * _second / 36525.0) * tau
      return theta, theta_dot
docs: Return the angle of Greenwich Mean Standard Time 1982 given the JD.
  
      This angle defines the difference between the idiosyncratic True
      Equator Mean Equinox (TEME) frame of reference used by SGP4 and the
      more standard Pseudo Earth Fixed (PEF) frame of reference.
  
      From AIAA 2006-6753 Appendix C.