Name

CL_fr_J20002G50 — EME2000 (J2000) to Gamma50 (Veis) frame vector transformation

Calling Sequence

   [pos_G50,[vel_G50,jacob]] = CL_fr_J20002G50(cjdutc,pos_J2000,[vel_J2000,ut1_utc,xp,yp,dPsi,dEps,conv])
   
   

Description

  • Computes the frame transfer of pos_J2000 and optionally vel_J2000 (EME2000 frame) to pos_G50 and optionally vel_G50 (Gamma50 (Veis) frame) at given dates (cjdutc).

    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

  • Last update : 27/6/2008

Parameters

cjdutc:

CNES Julian Date (1950) in TUC time (1xN)

pos_J2000:

position in EME2000 frame (3xN)

vel_J2000:

(optional) velocity in EME2000 frame (3xN)

ut1_utc :

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

xp :

(optional) x polar coordinate [radians] (default is 0) (1xN)

yp :

(optional) y polar coordinate [radians] (default is 0) (1xN)

dPsi :

(optional) Nutation corrections [radians] (default is 0) (1xN)

dEps :

(optional) Nutation corrections [radians] (default is 0) (1xN)

conv :

(optional) Convention IERS. Only iers 1996 (Lieske/Wahr) is implemented (default is "iers_1996")

pos_G50:

position in Gamma50 (Veis) frame (3xN)

vel_G50:

(optional) velocity in Gamma50 (Veis) frame (3xN)

jacob:

(optional) jacobian of the transformation (6x6xN)

Authors

CNES - DCT/SB

Bibliography

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

See also

CL_fr_G502J2000Mat, CL_fr_G502J2000

Examples

// 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);