Planet positions in Gamma50 ecliptic reference frame
[pars] = CL_mod_planetsG50ecl(cjd,planet_name)
Computes the Keplerian orbit parameters of planets in the Gamma 50 ecliptic reference frame, using a simple model.
Modified (1950.0) julian day (1xN)
(string) Name of planet (Mercury,Venus,Earth,Mars,Jupiter,Saturn,Uranus,Neptune or Pluto)
Planet's keplerian parameters at time cjd [sma,ecc,inc,pom,gom,M] [6xN]
CNES - DCT/SB
cjd = CL_dat_cal2cjd(2016,1,1,0,0,0); dates = cjd : 1 : cjd+365*5; kep_Earth = CL_mod_planetsG50ecl(dates,'Earth'); kep_Mars = CL_mod_planetsG50ecl(dates,'Mars'); pos_Mars = CL_oe_kep2car(kep_Mars,mu=%CL_muSun); pos_Earth = CL_oe_kep2car(kep_Earth,mu=%CL_muSun); scf(); plot(pos_Earth(1,:)/%CL_au,pos_Earth(2,:)/%CL_au,'g'); plot(pos_Mars(1,:)/%CL_au,pos_Mars(2,:)/%CL_au,'r'); | ![]() | ![]() |