CL_fr_G502J2000 — Gamma50 (Veis) to EME2000 (J2000) frame vector transformation
[pos_J2000,[vel_J2000,jacob]] = CL_fr_G502J2000(cjdutc,pos_G50,[vel_G50,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 Gamma50 (Veis) frame (3xN)
(optional) velocity in Gamma50 (Veis) 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 J2000 frame (3xN)
(optional) velocity in J2000 frame (3xN)
(optional) jacobian of the transformation (6x6xN)
// Conversion of position G50 to J2000 pos_G50 = [ [3952930.5;3127929.25;4128420.75] , [3945680.5;3125978.25;5687420.75]]; cjd = [21010 , 21011]; pos_J2000=CL_fr_G502J2000(cjd,pos_G50); // Conversion of position and velocity : G50 to J2000 + jacobian pos_G50 = [ [3952930.5;3127929.25;4128420.75] , [3945680.5;3125978.25;5687420.75]]; vel_G50 = [ [1.e3;3.e3;7.e3] , [1.e3;3.e3;7.e3]]; cjd = [21010 , 21011]; [pos_J2000,vel_J2000,jacob]=CL_fr_G502J2000(cjd,pos_G50,vel_G50);