Name

CL_fr_G502ter — Gamma50 (Veis) to terrestrial frame vector transformation

Calling Sequence

   [pos_ter,[vel_ter,jacob]] = CL_fr_G502ter(cjd,pos_G50,[vel_G50,ut1_utc])
   
   

Description

  • Given the positions pos_G50 and optionnaly vel_G50 in the Gamma50 (Veis) frame, this function computes positions pos_ter, and optionnaly vel_ter in the terrestrial frame.

    Jacobian of the transformation is also optionnaly 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

    Transformation from Gamma50 to terrestrial frame is a single rotation of the Veis sideral time(as calculated in CL_mod_sidTimeG50) around the Z axis.

  • Last update : 27/6/2008

Parameters

cjd:

CNES Julian dates (dates of the Gamma50 (Veis) frame) (1xN)

pos_G50:

position in Gamma50 (Veis) frame [X;Y;Z] (3xN)

vel_G50:

(optionnal) vector in Gamma50 (Veis) frame [X;Y;Z] (3xN)

ut1_utc :

(optional) ut1-utc [seconds] (default is 0) (1xN)

pos_ter:

position in terrestrial frame [X;Y;Z] (3xN)

vel_ter:

(optional) velocity in terrestrial frame [X;Y;Z] (3xN)

jacob:

(optional) jacobian of the transformation [d(x,y,z,vx,vy,vz)_ter/d(x,y,z,vx,vy,vz)_G50] (6x6xN)

Authors

CNES - DCT/SB

Bibliography

1 CNES - MSLIB FORTRAN 90, Volume R (mr_veis_TerVrai)

See also

CL_fr_ter2J2000Mat, CL_fr_G502terMat, CL_fr_ter2G50, CL_mod_sidTimeG50

Examples

// Conversion of position G50 to terrestrial
pos_G50 = [ [3952930.5;3127929.25;4128420.75] , [3945680.5;3125978.25;5687420.75]];
cjd = [21010 , 21011];
pos_ter=CL_fr_G502ter(cjd,pos_G50);

// Conversion of position and velocity : G50 to terrestrial  + 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_ter,vel_ter,jacob]=CL_fr_G502ter(cjd,pos_G50,vel_G50);