def get_polar_motion(time):
 (xp, yp, status) = iers.IERS_Auto.open().pm_xy(time, return_status=True)
 wmsg = None
 if np.any((status == iers.TIME_BEFORE_IERS_RANGE)):
  wmsg = u'Tried to get polar motions for times before IERS data is valid. Defaulting to polar motion from the 50-yr mean for those. This may affect precision at the 10s of arcsec level'
  xp.ravel()[(status.ravel() == iers.TIME_BEFORE_IERS_RANGE)] = _DEFAULT_PM[0]
  yp.ravel()[(status.ravel() == iers.TIME_BEFORE_IERS_RANGE)] = _DEFAULT_PM[1]
  warnings.warn(wmsg, AstropyWarning)
 if np.any((status == iers.TIME_BEYOND_IERS_RANGE)):
  wmsg = u'Tried to get polar motions for times after IERS data is valid. Defaulting to polar motion from the 50-yr mean for those. This may affect precision at the 10s of arcsec level'
  xp.ravel()[(status.ravel() == iers.TIME_BEYOND_IERS_RANGE)] = _DEFAULT_PM[0]
  yp.ravel()[(status.ravel() == iers.TIME_BEYOND_IERS_RANGE)] = _DEFAULT_PM[1]
  warnings.warn(wmsg, AstropyWarning)
 return (xp.to(u.radian).value, yp.to(u.radian).value)

def get_polar_motion ( time ) : ( xp , yp , status ) = iers . IERS_Auto . open ( ) . pm_xy ( time , return_status = True ) wmsg = None if np . any ( ( status == iers . TIME_BEFORE_IERS_RANGE ) ) : wmsg = u'Tried to get polar motions for times before IERS data is valid. Defaulting to polar motion from the 50-yr mean for those. This may affect precision at the 10s of arcsec level' xp . ravel ( ) [ ( status . ravel ( ) == iers . TIME_BEFORE_IERS_RANGE ) ] = _DEFAULT_PM [ 0 ] yp . ravel ( ) [ ( status . ravel ( ) == iers . TIME_BEFORE_IERS_RANGE ) ] = _DEFAULT_PM [ 1 ] warnings . warn ( wmsg , AstropyWarning ) if np . any ( ( status == iers . TIME_BEYOND_IERS_RANGE ) ) : wmsg = u'Tried to get polar motions for times after IERS data is valid. Defaulting to polar motion from the 50-yr mean for those. This may affect precision at the 10s of arcsec level' xp . ravel ( ) [ ( status . ravel ( ) == iers . TIME_BEYOND_IERS_RANGE ) ] = _DEFAULT_PM [ 0 ] yp . ravel ( ) [ ( status . ravel ( ) == iers . TIME_BEYOND_IERS_RANGE ) ] = _DEFAULT_PM [ 1 ] warnings . warn ( wmsg , AstropyWarning ) return ( xp . to ( u . radian ) . value , yp . to ( u . radian ) . value )
gets the two polar motion components in radians for use with apio13 .