pcpf_lla.py

# body 301 is the moon (note: bodvrd allows 'MOON' instead of 301)
mu   = spice.bodvcd(301, 'GM', 1)[1][0]    # gravity
req  = spice.bodvcd(301, 'RADII', 3)[1][0] # equatorial radius
rpol = spice.bodvcd(301, 'RADII', 3)[1][2] # polar radius
f    = (req - rpol) / req                  # first flattening

# Rectangular to geodetic coordinates (radians is default)
lon, lat, alt = spice.recgeo(r_lclf, req, f)

# Geodetic to rectangular coordinates
r_landing_lclf = spice.georec(lon_landing, lat_landing, alt_landing, req, f)