CL_mod_moonSunG50 — Moon and Sun coordinates in Gamma50 (Veis) frame
[r_moon,rl,r_sun,rs] = CL_mod_moonSunG50(cjd[,outs])
[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') ).
CNES Julian day (+ut/24) (1/1/1950) (1xN)
(optional string) type can be 'm' for the moon,'s' for the sun (or 'b' for both).
direction of the moon [rx;ry;rz] (3xN)
earth-moon distance [m] (1xN)
direction of the sun [sx;sy;sz] (3xN)
earth-sun distance [m] (1xN)
// 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