# UranusPoleExample2.py

import spiceypy as spice

kernels_dir = '../../../../kernels/'                                       
kernels = ['naif0012.tls','uranus_ringframes_rfrench20201201_v1.tf']     
kernels_list = [kernels_dir + k for k in kernels]                          
                                                                           
spice.furnsh(kernels_list)                                                 

ETsec_ring6_ingress = spice.str2et('UTC 1986 Jan 24 20:01:51.081')
ring_pole	= [0.0, 0.0, 1.0]

rotate		= spice.pxform('URING_6','B1950',ETsec_ring6_ingress)
pole_B1950	= spice.mxv(rotate,ring_pole)
len,RA,DE	= spice.recrad(pole_B1950)
RAdeg		= RA * spice.dpr()
DEdeg		= DE * spice.dpr()
print('B1950 pole direction for Voyager 2 RSS ring 6 ingress:')
print("RA = ",'{:07.4f}'.format(RAdeg)," Dec = ",'{:07.4f}'.format(DEdeg))

rotate		= spice.pxform('URING_6','J2000',ETsec_ring6_ingress)
pole_J2000	= spice.mxv(rotate,ring_pole)
len,RA,DE	= spice.recrad(pole_J2000)
RAdeg		= RA * spice.dpr()
DEdeg		= DE * spice.dpr()
print('J2000 pole direction for Voyager 2 RSS ring 6 ingress:')
print("RA = ",'{:07.4f}'.format(RAdeg)," Dec = ",'{:07.4f}'.format(DEdeg))