<< CL_mod_nutationMatrix Models CL_mod_polarMotionMatrix >>

CelestLab >> Models > CL_mod_planetsG50ecl

CL_mod_planetsG50ecl

Planet positions in Gamma50 ecliptic reference frame

Calling Sequence

[pars] = CL_mod_planetsG50ecl(cjd,planet_name)

Description

Parameters

cjd:

Modified (1950.0) julian day (1xN)

planet_name:

(string) Name of planet (Mercury,Venus,Earth,Mars,Jupiter,Saturn,Uranus,Neptune or Pluto)

pars:

Planet's keplerian parameters at time cjd [sma,ecc,inc,pom,gom,M] [6xN]

See also

Authors

Examples

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');
<< CL_mod_nutationMatrix Models CL_mod_polarMotionMatrix >>