CL_fr_J20002G50 — EME2000 (J2000) to Gamma50 (Veis) frame vector transformation
[pos_G50,[vel_G50,jacob]] = CL_fr_J20002G50(cjdutc,pos_J2000,[vel_J2000,ut1_utc,xp,yp,dPsi,dEps,conv])
Jacobian of the transformation is also optionally computed.
Gamma50 (Veis) (Veis frame at date) is defined as:
- origin at the center of the Earth
- plan OXY is the true equatorial plan of the date
- OX axis goes through the Veis equinox of the date
- OZ axis goes through the true celestial pole of the date
- OY completes the orthonormal frame
CNES Julian Date (1950) in TUC time (1xN)
position in EME2000 frame (3xN)
(optional) velocity in EME2000 frame (3xN)
(optional) ut1-utc [seconds] (default is 0) (1xN)
(optional) x polar coordinate [radians] (default is 0) (1xN)
(optional) y polar coordinate [radians] (default is 0) (1xN)
(optional) Nutation corrections [radians] (default is 0) (1xN)
(optional) Nutation corrections [radians] (default is 0) (1xN)
(optional) Convention IERS. Only iers 1996 (Lieske/Wahr) is implemented (default is "iers_1996")
position in Gamma50 (Veis) frame (3xN)
(optional) velocity in Gamma50 (Veis) frame (3xN)
(optional) jacobian of the transformation (6x6xN)
// Conversion of position J2000 to G50 pos_J2000 = [ [3952930.5;3127929.25;4128420.75] , [3945680.5;3125978.25;5687420.75]]; cjd = [21010 , 21011]; pos_G50=CL_fr_J20002G50(cjd,pos_J2000); // Conversion of position and velocity : J2000 to G50 + jacobian pos_J2000 = [ [3952930.5;3127929.25;4128420.75] , [3945680.5;3125978.25;5687420.75]]; vel_J2000 = [ [1.e3;3.e3;7.e3] , [1.e3;3.e3;7.e3]]; cjd = [21010 , 21011]; [pos_G50,vel_G50,jacob]=CL_fr_J20002G50(cjd,pos_J2000,vel_J2000);