CL_fr_ter2J2000Mat — Terrestrial to EME2000 (J2000) matrix
M = CL_fr_ter2J2000Mat(cjdutc[,ut1_utc[,xp,yp[,dPsi,dEps[,conv]]]])
To convert a vector vect_J2000 to vect_ter, you can use the (hyper)matrix M' : vect_ter = M' * vect_J2000
For big hypermatrixes, it is faster to use the functions CL_fr_ter2J2000 and CL_fr_J20002ter
CNES Julian Date (1950) in TUC time (1xN)
(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")
terrestrial to EME2000 transfer (hyper)matrix (3x3xN)
1 IERS Conventions (1996), Dennis D. McCarthy
2 Explanatory Supplement to the Astronomical Almanac, Seidelman (1992)
// Conversion J2000 to terrestrial pos_J2000 = [ [3952930.5;3127929.25;4128420.75] , [3945680.5;3125978.25;5687420.75]]; cjd = [21010 , 21011]; M=CL_fr_ter2J2000Mat(cjd); pos_ter = M' * (pos_J2000); // Conversion terrestrial to J2000 pos_ter = CL_co_sph2car([ [0;%pi/4;%CL_eqRad+700.e3] , [0;-%pi/4;%CL_eqRad+800.e3]]); cjd = [21010 , 21011]; M=CL_fr_ter2J2000Mat(cjd); pos_J2000 = M * (pos_ter);