Path: Orbit/OrbitData
% Reads in JPL Horizons data.
http://ssd.jpl.nasa.gov/horizons.cgi#top
The format should look like
Ephemeris Type [change] : VECTORS
Target Body [change] : Earth [Geocenter] [399]
Coordinate Origin [change] : Sun (body center) [500@10]
Time Span [change] : Start=2010-09-03, Stop=2012-10-03, Step=1 d
Table Settings [change] : defaults
Display/Output [change] : download/save (plain text file)% Save the data to a text file. Call without the file input to get a
dialog box.
Units should be au and au/day.
For a demo type
ReadHorizons('Messenger.txt')
--------------------------------------------------------------------------
Form:
[r, v, jD, lT, rG, rR] = ReadHorizons( file )
--------------------------------------------------------------------------
------
Inputs
------
file (1,:) File name
-------
Outputs
-------
r (3,:) Position vector (AU)
v (3,:) Velocity vector (AU/day)
jD (1,:) Julian date (day)
lT (1,:) Light time (day)
rG (1,:) Range (AU)
rR (1,:) Range Rate (AU/day)
--------------------------------------------------------------------------
Orbit: OrbitData/Messenger Common: Database/Constant Common: Graphics/Plot3D
Back to the Orbit Module page