**would like to add more information**

Hi!

As I've battled this little annoyance myself, I'd like to share my results

I've created a vessel on the moon, roughly on the surface

using oapiEquToLocal(theMoon, lon, lat, range, &observerPosition);

I recieved:

**(lat,lng) -> (X, Y, Z) **

(0.000000,0.000008)->(1738007.021691,0.248463,0.000000)

However, using the method above

oapiEquToGlobal(theMoon, lon, lat, range, &gposSurfGlobalPoint); gposSurfGlobalPoint -= gposMoon;

I've gotten

Equatorial = (-1692619.117615,-20872.727713,-394047.280632)

Using Matlab, I've extracted the rotation matrix between these two vectors

r = vrrotvec([1738002.718203,0.155289,0.000000], [-1692612.355560,-20873.042428,-394057.328804])

>> m = vrrotvec2mat(r)

m =

** -0.9739 0.0120 0.2267**

-0.0120 0.9945 -0.1043

-0.2267 -0.1043 -0.9684

However, using orbiter method oapiGetRotationMatrix(theMoon, &rot);

I've gotten

**matrix: [-0.973888,-0.006244,0.226945]**

matrix: [-0.012009,0.999639,-0.024032]

matrix: [-0.226713,-0.026130,-0.973611]

there is ALMOST a unity between them

** 1.0000 -0.0000 -0.0000**

-0.0000 0.9968 -0.0803

0.0000 0.0803 0.9968

extracting the *Residue Angle * using matlab,

[r1 r2 r3] = dcm2angle(A'*B)

r1 = 0

r2 = 0

r3 = -0.0804

**-0.0804*180/pi = -4.6066 degress**

Wikipedia says the inclination of the moon is** 5.145° **

Unless I've gotten something wrong, there's more rotation between these frames? and the 4.6 degrees' suppose to be the inclination (open question)

fix: I think the last angle is the angle of rotation between a stationary (planet-attached) axes and a stationary axes with a moving planet.