Moon and Sun positions in Gamma50 (Veis) reference frame - DEPRECATED
[u_moon,r_moon,u_sun,r_sun] = CL_mod_moonSunG50(cjd [,body="b"]) [u_moon,r_moon] = CL_mod_moonSunG50(cjd, body="m") [u_sun,r_sun] = CL_mod_moonSunG50(cjd, body="s")
This function is deprecated.
Replacement functions: CL_eph_moon and CL_eph_sun
Computes the Moon and/or Sun positions relative to Gamma50 (Veis) reference frame using a simplified version of BROWN theories for the moon and NEWCOMB for the Sun.
The Moon or Sun position is defined by a unit vector u_body and a distance r_body from the Earth centre to the centre of the body.
Either the Moon position, the Sun position or both can be computed.
Note: The time scale is TT and not UTC as wrongly indicated in previous versions.
Modified (1950.0) julian day (Time scale: TT) (1xN)
(string, optional) Type can be 'm' or 'moon' for the moon,'s' or 'sun' for the Sun (or 'b' or 'both' for both).
Unit vector from Earth to body [rx;ry;rz] (3xN)
Distance from Earth to body [m] (1xN)
1 CNES - MSLIB FORTRAN 90, Volume S (ms_pos_soleil_lune)
CNES - DCT/SB
// Sun position in J2000 frame on the 25th of october 2008 and 2009 cjd = CL_dat_cal2cjd([2008 2009],[10 10],[25 25]); [r_sun,rs] = CL_mod_moonSunG50(cjd,'s'); pos_sun_G50 = CL_dMult(r_sun,rs); M = CL_fr_G502J2000Mat(cjd); pos_sun_J2000 = M*pos_sun_G50; // Sun right ascension and declination at the same dates : pos_sph = CL_co_car2sph(pos_sun_J2000); alpha_sun = pos_sph(1,:); delta_sun = pos_sph(2,:); |  |  |