Name

CL_mod_moonSunG50 — Moon and Sun coordinates in Gamma50 (Veis) frame

Calling Sequence

   [r_moon,rl,r_sun,rs] = CL_mod_moonSunG50(cjd[,outs])
   
   

Description

  • Calculates the moon and sun coordinates in the GAMMA50 frame using a simplified version of BROWN theories for the moon and NEWCOMB for the sun.

    [A,B]=CL_mod_moonSunG50(cjd,'m') -> A is r_moon and B is rl.

    [A,B]=CL_mod_moonSunG50(cjd,'s') -> A is r_sun and B is rs.

    [A,B]=CL_mod_moonSunG50(cjd) -> A is r_moon and B is rl

    If function called with more than two output arguments, outs is senseless and it is set to 'b'.

    NOTE: optional input variable outs is used to choose what values are actually computed; you can always ask less outputs arguments than those computed ( ex: [r_m,rm]=CL_mod_moonSunG50(cjd,'b') ).

  • Last update : 03/10/2007

Parameters

cjd :

CNES Julian day (+ut/24) (1/1/1950) (1xN)

outs :

(optional string) type can be 'm' for the moon,'s' for the sun (or 'b' for both).

r_moon :

direction of the moon [rx;ry;rz] (3xN)

rl :

earth-moon distance [m] (1xN)

r_sun :

direction of the sun [sx;sy;sz] (3xN)

rs :

earth-sun distance [m] (1xN)

Bibliography

1 CNES - MSLIB FORTRAN 90, Volume S (ms_pos_soleil_lune)

Authors

CNES - DCT/SB

Examples

// Sun position in J2000 frame on the 25th of october 2008 and 2009
cjd = CL_dat_cal2cjd([2008 2009],10,25)
[r_sun,rs] = CL_mod_moonSunG50(cjd,'s')
M = CL_fr_G502J2000Mat(cjd)
pos_sun = M*(r_sun.*(rs.*.ones(3,1))) //first column is 2008, second column is 2009

// Sun right ascension and declination at the same dates :
pos_sph = CL_co_car2sph(pos_sun);
alpha_sun = pos_sph(1,:)
delta_sun = pos_sph(2,:)

See 'CelestLab > Demos > ORBIT_PROPERTIES > beta_sun' for more examples