tengai 0.0.1

data/ext/horizons/vector_ephemeris_parser.rl

@@ -0,0 +1,115 @@
+=begin
+%%{
+
+  machine ephemeris_parser;
+
+  action mark { mark = p }
+
+  action target_body_id {
+    target_body_id = data[mark..p].pack('c*')
+  }
+
+  action center_body_id {
+    center_body_id = data[mark..p].pack('c*')
+  }
+
+  action start_time {
+    start_time = data[mark..p].pack('c*')
+  }
+
+  action stop_time {
+    stop_time = data[mark..p].pack('c*')
+  }
+
+  action step_size {
+    step_size = data[mark..p].pack('c*')
+  }
+
+  action ephemeris_table {
+    ephemeris_table = data[mark..p - 1].pack('c*')
+  }
+
+  action ephemeris_columns {
+    ephemeris_columns = data[mark..p].pack('c*')
+  }
+
+  ws = [ \t\r\n];
+
+  adbc = ('A.D.'|'B.C.');
+  year = digit{4};
+  month = upper lower{2};
+  date = digit{2};
+  hours =  digit{2};
+  minutes = digit{2};
+  seconds = digit{2} '.' digit{4};
+  tz = ('C'|'U') 'T';
+  datetime = adbc ' ' year '-' month '-' date ' ' hours ':' minutes ':' seconds ' ' tz;
+
+  hr = '*'+ '\n';
+
+  time_unit = ('minute' [s]? | 'calendar year' [s]?);
+
+  body_name = (alnum | space)*;
+  body_id = digit*;
+
+  target_body = 'Target body name:' space body_name space
+    '(' body_id >mark @target_body_id ')'
+    ws* '{' any* '}' '\n';
+
+  center_body = 'Center body name:' space body_name space
+    '(' body_id >mark @center_body_id ')'
+    ws* '{' any* '}' '\n';
+
+  start_time = 'Start time' ' '* ':' ' ' datetime >mark %start_time space* '\n';
+  stop_time  = 'Stop  time' ' '* ':' ' ' datetime >mark %stop_time space* '\n';
+  step_size  = 'Step-size' ' '* ':' ' ' (digit+ ' '* time_unit) >mark $step_size '\n';
+
+  ephemeris_columns = ('JDCT' (' ' | alpha | ',')*)
+  >mark @ephemeris_columns :> '\n';
+
+  soe = '$$SOE' '\n';
+  eoe = '$$EOE' '\n';
+  ephemeris_table = any*;
+  ephemeris = soe any* >mark %ephemeris_table :> eoe;
+
+  main := (
+    any*
+    target_body
+    center_body
+    any*
+    start_time
+    stop_time
+    step_size
+    any*
+    hr
+    ephemeris_columns
+    hr
+    ephemeris
+    any*
+  );
+
+}%%
+=end
+
+require 'date'
+module Tengai
+  class VectorEphemerisParser
+    def self.parse(data, ephemeris_table_parser=EphemerisTableParser)
+      data = data.unpack('c*') if data.is_a? String
+      eof = data.length
+
+      %% write init;
+      %% write exec;
+
+      { target_body_id: target_body_id.to_i,
+        center_body_id: center_body_id.to_i,
+        start_time: DateTime.parse(start_time),
+        stop_time: DateTime.parse(stop_time),
+        step_size: step_size,
+        ephemeris_table: ephemeris_table_parser.parse(
+          "#{ephemeris_columns}\n#{ephemeris_table}") }
+    end
+
+    %% write data;
+  end
+end

data/lib/tengai/body.rb

@@ -0,0 +1,22 @@
+$:.unshift File.dirname(File.dirname(File.dirname(__FILE__)))
+require 'ext/horizons/body_data_sheet_parser'
+require 'virtus'
+
+module Tengai
+  class Body
+    include Virtus
+    attribute :revised_on, Date
+    attribute :name, String
+    attribute :id, Integer
+
+    def self.find(client, body, parser=BodyDataSheetParser)
+      data = client.cmd('String' => body.to_s, 'Match' => /<cr>:\s*$/)
+      data = parser.parse(data)
+      self.new(data)
+    end
+
+    def ==(other)
+      id == other.id
+    end
+  end
+end

data/lib/tengai/client.rb

@@ -0,0 +1,34 @@
+require 'net/telnet'
+
+module Tengai
+  class Client
+    HOST = 'horizons.jpl.nasa.gov'.freeze
+    PORT = '6775'.freeze
+    PROMPT = /Horizons>/.freeze
+
+    def initialize(telnet=Net::Telnet)
+      @telnet = telnet
+      connect!
+    end
+
+    def connect!
+      @connection = @telnet.new(
+        'Host'   => HOST,
+        'Port'   => PORT,
+        'Prompt' => PROMPT)
+      @connection.waitfor 'Match' => PROMPT
+    end
+
+    def disconnect
+      @connection.close
+    end
+
+    def connection
+      @connection
+    end
+
+    def cmd(*args, &blk)
+      connection.cmd(*args, &blk)
+    end
+  end
+end

data/lib/tengai/ephemeris.rb

@@ -0,0 +1,56 @@
+$:.unshift File.dirname(File.dirname(__FILE__))
+require 'tengai/requests/ephemeris_request'
+require 'tengai/vector_ephemeris_table'
+require 'virtus'
+
+module Tengai
+  class Ephemeris
+    include Virtus
+    attribute :target_body_id, Integer
+    attribute :center_body_id, Integer
+    attribute :start_time, DateTime
+    attribute :stop_time, DateTime
+    attribute :step_size, Integer
+    attribute :ephemeris_table, VectorEphemerisTable
+
+    def initialize(client, data)
+      super(data)
+      @client = client
+    end
+
+    def target_body
+      @target_body ||= Tengai::Body.find(@client, target_body_id)
+    end
+
+    def ephemeris_table=(table)
+      @ephemeris_table = VectorEphemerisTable.new_with_table(table)
+    end
+
+    # Public: fetch an ephemeris table for a given body using the client
+    #
+    # client - Tengai::Client
+    # body - a Tengai::Body instance whose ephemeris is needed
+    # options - A Hash of options for refining the ephemeris table
+    #           :start_time - The Time beginning of the ephemeris
+    #           :stop_time - The Time end of the ephemeris
+    #           :interval - The resolution of the table in minutes
+    def self.fetch(client, body, options={})
+      start_time = options[:start_time] || Date.today.to_time # defaults start at today
+      stop_time = options[:stop_time] || (Date.today + 1).to_time # and end tomorrow
+      interval = options[:interval] || 1440
+
+      # Inject dependencies
+      request = options[:request] || EphemerisRequest
+      parser = options[:parser] || VectorEphemerisParser
+
+      data = request.fetch(
+        client, body.id,
+        start_time: start_time,
+        stop_time: stop_time,
+        interval: interval)
+      data = parser.parse(data)
+
+      new(client, data)
+    end
+  end
+end

data/lib/tengai/parsers/ephemeris_table_parser.rb

@@ -0,0 +1,20 @@
+require 'csv'
+
+module Tengai
+  class EphemerisTableParser
+    def self.parse(table)
+      empty_nils = ->(value) { value || '' }
+      strip      = ->(value) { value.strip }
+
+      table = CSV.parse(
+        table,
+        headers:            true,
+        header_converters:  [empty_nils, strip, :symbol],
+        converters:         [empty_nils, strip])
+
+      table.map do |row|
+        row.to_hash.tap {|r| r.delete(:'') }
+      end
+    end
+  end
+end

data/lib/tengai/requests/ephemeris_request.rb

@@ -0,0 +1,131 @@
+module Tengai
+  class EphemerisRequest
+    SYSTEM_PROMPT = %r{Horizons>}.freeze
+    DEFAULT_PROMPT = %r{<cr>:\s*$}
+    TABLE_PROMPT = %r{Observe, Elements, Vectors  \[o,e,v,\?\] :\s*$}.freeze
+    FIRST_OBSERVER_PROMPT = %r{Coordinate center \[ <id>,coord,geo  \] :\s*$}.freeze
+    SUBSEQUENT_OBSERVER_PROMPT = %r{Use previous center  \[ cr=\(y\), n, \? \] :\s*$}.freeze
+    OBSERVER_PROMPT = Regexp.union(FIRST_OBSERVER_PROMPT, SUBSEQUENT_OBSERVER_PROMPT).freeze
+    CONFIRM_OBSERVER_PROMPT = %r{Confirm selected station    \[ y/n \] -->\s*$}.freeze
+    REFERNCE_PLANE_PROMPT = %r{Reference plane \[eclip, frame, body \] :\s*}.freeze
+    START_TIME_PROMPT = %r{Starting (C|U)T .* :\s*$}.freeze
+    END_TIME_PROMPT = %r{Ending \s* (C|U)T .* :\s*$}.freeze
+    INTERVAL_PROMPT = %r{Output interval \[ex: 10m, 1h, 1d, \? \] :\s*$}.freeze
+    ACCEPT_DEFAULT_OUTPUT_PROMPT = %r{Accept default output \[ cr=\(y\), n, \?\] :\s*$}.freeze
+    OUTPUT_REFERENCE_FRAME_PROMPT = %r{Output reference frame \[J2000, B1950\] :\s*$}.freeze
+    CORRECTIONS_PROMPT = %r{Corrections \[ 1=NONE, 2=LT, 3=LT\+S \]  :\s*$}.freeze
+    OUTPUT_UNITS = %r{Output units \[1=KM-S, 2=AU-D, 3=KM-D\] :\s*$}.freeze
+    CSV_FORMAT_PROMPT = %r{Spreadsheet CSV format    \[ YES, NO \] :\s*$}.freeze
+    LABEL_CARTESIAN_OUTPUT_PROMPT = %r{Label cartesian output    \[ YES, NO \] :\s*$}.freeze
+    SELECT_OUTPUT_TABLE_TYPE_PROMPT = %r{Select output table type  \[ 1-6, \?  \] :\s*$}.freeze
+    SELECT_QUANTITIES_PROMPT = %r{Select table quantities \[ <#,#\.\.>, \?\] :\s*$}.freeze
+    COMPLETED_PROMPT = %r{>>> Select\.\.\. \[A\]gain, \[N\]ew-case, \[F\]tp, \[K\]ermit, \[M\]ail, \[R\]edisplay, \? :\s*$}.freeze
+
+    ANY_PROMPT = Regexp.union(
+      DEFAULT_PROMPT, TABLE_PROMPT, OBSERVER_PROMPT, CONFIRM_OBSERVER_PROMPT,
+      REFERNCE_PLANE_PROMPT, START_TIME_PROMPT, END_TIME_PROMPT, INTERVAL_PROMPT,
+      ACCEPT_DEFAULT_OUTPUT_PROMPT, SELECT_QUANTITIES_PROMPT, COMPLETED_PROMPT,
+      OUTPUT_REFERENCE_FRAME_PROMPT, CORRECTIONS_PROMPT, OUTPUT_UNITS,
+      CSV_FORMAT_PROMPT, LABEL_CARTESIAN_OUTPUT_PROMPT,
+      SELECT_OUTPUT_TABLE_TYPE_PROMPT, SYSTEM_PROMPT).freeze
+
+    TIME_FORMAT = '%Y-%b-%d %H:%M'.freeze
+
+    def initialize(client, body, options={})
+      @client = client
+      @body = body.to_s
+      @start_time = options[:start_time]
+      @stop_time = options[:stop_time]
+      @interval = options[:interval] || 1440
+      @options = options
+      @state = :ready
+    end
+
+    def self.fetch(client, body, options={})
+      new(client, body, options).fetch
+    end
+
+    # Public: initiates the ephemeris request
+    #
+    # Examples:
+    #
+    #   Tengai::EphemerisRequest.new(Tengai::Client.new, 499).fetch
+    #   # => #<Tengai::EphemerisRequest @data=""B\n \r\n Working ...   \b\b-  \r\n\e[?1h\e=\r*...">
+    #
+    # Returns the the request (self)
+    def fetch
+      raise "Not ready" unless @state == :ready
+      @state = :fetching
+      send_command(@body)
+      @data
+    end
+
+    private
+    def table
+      case @options[:table]
+      when :observer then 'o'
+      when :vector then 'v'
+      when :orbital_elements then 'e'
+      else 'v'
+      end
+    end
+
+    def observer
+      @options[:observer] || '500@0'
+    end
+
+    def send_command(command, prompt=ANY_PROMPT)
+      Tengai.log "> #{command}"
+      result = @client.cmd('String' => command, 'Match' => prompt) do |data|
+        Tengai.log '< ' + data
+      end
+      receive_data(result)
+    end
+
+    def receive_data(data)
+      case data
+      when DEFAULT_PROMPT
+        send_command 'E'
+      when TABLE_PROMPT
+        send_command table
+      when FIRST_OBSERVER_PROMPT
+        send_command observer
+      when SUBSEQUENT_OBSERVER_PROMPT
+        send_command 'n'
+      when CONFIRM_OBSERVER_PROMPT
+        send_command 'y'
+      when REFERNCE_PLANE_PROMPT
+        send_command 'frame'
+      when START_TIME_PROMPT
+        send_command @start_time.strftime(TIME_FORMAT)
+      when END_TIME_PROMPT
+        send_command @stop_time.strftime(TIME_FORMAT)
+      when INTERVAL_PROMPT
+        send_command @interval.to_s + 'm'
+      when ACCEPT_DEFAULT_OUTPUT_PROMPT
+        send_command 'n'
+      when OUTPUT_REFERENCE_FRAME_PROMPT
+        send_command 'J2000'
+      when CORRECTIONS_PROMPT
+        send_command '1'
+      when OUTPUT_UNITS
+        send_command '2'
+      when CSV_FORMAT_PROMPT
+        send_command 'YES'
+      when LABEL_CARTESIAN_OUTPUT_PROMPT
+        send_command 'YES'
+      when SELECT_OUTPUT_TABLE_TYPE_PROMPT
+        send_command '03'
+      when SELECT_QUANTITIES_PROMPT
+        send_command 'B'
+      when COMPLETED_PROMPT
+        @data = data
+        send_command 'N'
+      when Client::PROMPT
+        @state = :ready
+      else
+        puts "Unexpected data: #{data}"
+      end
+    end
+  end
+end

data/lib/tengai/vector_ephemeris_table.rb

@@ -0,0 +1,36 @@
+require 'date'
+require 'virtus'
+
+class VectorEphemerisTable
+  include Enumerable
+
+  def initialize(rows)
+    @rows = rows
+  end
+
+  def self.new_with_table(rows)
+    new(rows.map(&Row.method(:new)))
+  end
+
+  def each(&block)
+    @rows.each(&block)
+  end
+
+  class Row
+    include Virtus
+    attribute :jdct, DateTime
+    attribute :x, Float
+    attribute :y, Float
+    attribute :z, Float
+    attribute :vx, Float
+    attribute :vy, Float
+    attribute :vz, Float
+    attribute :lt, Float
+    attribute :rg, Float
+    attribute :rr, Float
+
+    def jdct=(julian)
+      super DateTime.jd(julian.to_f)
+    end
+  end
+end

data/lib/tengai/version.rb

@@ -0,0 +1,5 @@
+module Tengai
+  unless defined?(Tengai::VERSION)
+    VERSION = '0.0.1'.freeze
+  end
+end

data/lib/tengai.rb

@@ -0,0 +1,27 @@
+$:.unshift File.dirname(__FILE__)
+$:.unshift File.join(File.dirname(File.dirname(__FILE__)), 'ext')
+
+require 'logger'
+
+require 'tengai/version'
+require 'tengai/client'
+require 'tengai/body'
+require 'tengai/ephemeris'
+
+require 'horizons/vector_ephemeris_parser'
+require 'tengai/parsers/ephemeris_table_parser'
+
+module Tengai
+  class << self
+    attr_accessor :debug
+    attr_reader :logger
+
+    def log(message)
+      logger.debug { message } if debug
+    end
+  end
+
+  @logger ||= ::Logger.new(STDOUT)
+
+  self.debug = false
+end

data/tengai.gemspec

@@ -0,0 +1,63 @@
+$LOAD_PATH.unshift 'lib'
+$LOAD_PATH.unshift 'ext'
+require 'tengai/version'
+
+Gem::Specification.new do |s|
+  s.name        = 'tengai'
+  s.version     = Tengai::VERSION.dup
+  s.summary     = 'Ruby wrapper for the NASA HORIZONS telnet system'
+  s.description = 'Ruby wrapper for the NASA HORIZONS telnet system'
+  s.homepage    = 'http://github.com/zacstewart/tengai'
+  s.authors     = ['Zac Stewart']
+  s.email       = ['zgstewart@gmail.com']
+
+  s.require_paths = ['lib']
+
+  s.add_dependency 'virtus'
+
+  s.add_development_dependency 'pry'
+  s.add_development_dependency 'rake'
+  s.add_development_dependency 'mocha'
+  s.add_development_dependency 'simplecov'
+  s.add_development_dependency 'forgery'
+
+  s.files = %w(
+    .gitignore
+    .ruby-gemset
+    .ruby-version
+    .travis.yml
+    Gemfile
+    README.md
+    Rakefile
+    ext/horizons/body_data_sheet_parser.rb
+    ext/horizons/body_data_sheet_parser.rl
+    ext/horizons/vector_ephemeris_parser.rb
+    ext/horizons/vector_ephemeris_parser.rl
+    lib/tengai.rb
+    lib/tengai/body.rb
+    lib/tengai/client.rb
+    lib/tengai/ephemeris.rb
+    lib/tengai/parsers/ephemeris_table_parser.rb
+    lib/tengai/requests/ephemeris_request.rb
+    lib/tengai/vector_ephemeris_table.rb
+    lib/tengai/version.rb
+    tengai.gemspec
+    test/fixtures.rb
+    test/fixtures/bodies/10.txt
+    test/fixtures/bodies/399.txt
+    test/fixtures/bodies/499.txt
+    test/fixtures/bodies/501.txt
+    test/fixtures/ephemerides/mars_observed_by_earth_from_2012_12_28_to_29.txt
+    test/fixtures/ephemerides/mars_vectors_from_solar_system_center.txt
+    test/fixtures/major_bodies.txt
+    test/integration/body_integration_test.rb
+    test/integration/ephemeris_integration_test.rb
+    test/test_helper.rb
+    test/unit/body_data_sheet_parser_test.rb
+    test/unit/body_test.rb
+    test/unit/ephemeris_table_parser_test.rb
+    test/unit/ephemeris_test.rb
+    test/unit/vector_ephemeris_parser_test.rb
+    test/unit/vector_ephemeris_table_test.rb
+  )
+end

data/test/fixtures/bodies/10.txt

@@ -0,0 +1,23 @@
+*******************************************************************************
+ Revised : Sep 12, 1996                  Sun                                 10
+ Physical: May 29, 2012
+
+ PHYSICAL PROPERTIES:
+  GM (10^11 km^3/s^2)   = 1.32712440018   Mass (10^30 kg)       =    1.9891
+  Radius (photosphere)  = 6.963(10^5) km  Angular diam at 1 AU  = 1919.3"
+  Solar Radius (IAU)    = 6.955(10^5) km  Mean density          = 1.408 g/cm^3
+  Surface gravity       =  274.0 m/s^2    Moment of inertia     = 0.059
+  Escape velocity       =  617.7 km/s     Adopted sidereal per  = 25.38 d
+  Pole (RA,DEC in deg.) =  286.13,63.87   Obliquity to ecliptic = 7 deg 15'        
+  Solar constant (1 AU) = 1367.6 W/m^2    Solar lumin.(erg/s)   =  3.846(10^33)
+  Mass-energy conv rate = 4.3(10^12 gm/s) Effective temp (K)    =  5778
+  Surf. temp (photosphr)= 6600 K (bottom) Surf. temp (photosphr)=  4400 K (top)
+  Photospheric depth    = ~400 km         Chromospheric depth   = ~2500 km
+  Sunspot cycle         = 11.4 yr         Cycle 22 sunspot min. =  1991 A.D.
+
+  Motn. rel to nrby strs= apex : RA=271 deg; DEC=+30 deg
+                          speed: 19.4 km/s = 0.0112 AU/day
+  Motn. rel to 2.73K BB = apex : l=264.7+-0.8; b=48.2+-0.5
+                          speed: 369 +-11 km/s
+*******************************************************************************
+ Select ... [E]phemeris, [F]tp, [M]ail, [R]edisplay, ?, <cr>: 

data/test/fixtures/bodies/399.txt

@@ -0,0 +1,29 @@
+*******************************************************************************
+Revised: Apr 03, 2003                  Earth                                399
+ 
+PHYSICAL PROPERTIES:
+ Mean radius, km          = 6371.01+-0.02   Mass, 10^24 kg = 5.9736
+ Equ. radius, km          = 6378.136        Mass layers:
+ Polar axis, km           = 6356.752          Atmos          = 5.1   x 10^18 kg
+ Flattening               = 1/298.257         oceans         = 1.4   x 10^21 kg
+ Density, gm cm^-3        = 5.515             crust          = 2.6   x 10^22 kg
+ J2  (GEM T2, 1990)       = 0.0010826265      mantle         = 4.043 x 10^24 kg
+ gp, m s^-2 (polar)       = 9.8321863685      outer core     = 1.835 x 10^24 kg
+ ge, m s^-2 (equatorial)  = 9.7803267715      inner core     = 9.675 x 10^22 kg
+ go, m s^-2               = 9.82022         Fluid core rad   = 3480 km
+ GM, km^3 s^-2            = 398600.440      Inner core rad   = 1215 km
+ Mean rot. rate, rad s^-1 = 7.292115*10^-5  Surface Area:
+ Sidereal period, hr      = 23.93419          land           = 1.48 x 10^8 km
+ Mean solar day, days     = 1.002738          sea            = 3.62 x 10^8 km
+ Moment of inertia        = 0.3308          Love no., k2     = 0.299
+ Mean Temperature, K      = 270             Atm. pressure    = 1.0 bar
+ Solar constant, W/m^2    = 1367.6          Vis. mag. V(1,0) = -3.86
+ Volume, 10^10 km^3       = 108.321         Geometric albedo = 0.367 
+
+DYNAMICAL CHARACTERISTICS:
+ Obliquity to orbit, deg  = 23.45           Sidereal period  = 1.0000174  yrs
+ Orbit velocity, km s^-1  = 29.7859         Sidereal period  = 365.25636  days
+ Mean daily motion, n     = 0.9856474 deg/d Escape velocity  = 11.186 km s^-2
+ Hill's sphere radius     = 234.9           Magnetic moment  = 0.61 gauss Rp^3
+*******************************************************************************
+ Select ... [E]phemeris, [F]tp, [M]ail, [R]edisplay, ?, <cr>: 

data/test/fixtures/bodies/499.txt

@@ -0,0 +1,25 @@
+*******************************************************************************
+ Revised: Sep 28, 2012                 Mars                             499 / 4
+ 
+ GEOPHYSICAL DATA (updated 2009-May-26):
+  Mean radius (km)      = 3389.9(2+-4)    Density (g cm^-3)     =  3.933(5+-4)
+  Mass (10^23 kg )      =    6.4185       Flattening, f         =  1/154.409
+  Volume (x10^10 km^3)  =   16.318        Semi-major axis       =  3397+-4
+  Sidereal rot. period  =   24.622962 hr  Rot. Rate (x10^5 s)   =  7.088218
+  Mean solar day        =    1.0274907 d  Polar gravity ms^-2   =  3.758
+  Mom. of Inertia       =    0.366        Equ. gravity  ms^-2   =  3.71
+  Core radius (km)      =  ~1700          Potential Love # k2   =  0.153 +-.017
+ 
+  Grav spectral fact u  =   14 (x10^5)    Topo. spectral fact t = 96 (x10^5)
+  Fig. offset (Rcf-Rcm) = 2.50+-0.07 km   Offset (lat./long.)   = 62d / 88d
+  GM (km^3 s^-2)        = 42828.3         Equatorial Radius, Re = 3394.0 km 
+  GM 1-sigma (km^3 s^-2)= +- 0.1          Mass ratio (Sun/Mars) = 3098708+-9
+  
+  Atmos. pressure (bar) =    0.0056       Max. angular diam.    =  17.9"
+  Mean Temperature (K)  =  210            Visual mag. V(1,0)    =  -1.52
+  Geometric albedo      =    0.150        Obliquity to orbit    =  25.19 deg
+  Mean sidereal orb per =    1.88081578 y Orbit vel.  km/s      =  24.1309
+  Mean sidereal orb per =  686.98 d       Escape vel. km/s      =   5.027
+  Hill's sphere rad. Rp =  319.8          Mag. mom (gauss Rp^3) = < 1x10^-4 
+*******************************************************************************
+ Select ... [E]phemeris, [F]tp, [M]ail, [R]edisplay, ?, <cr>: 

data/test/fixtures/bodies/501.txt

@@ -0,0 +1,13 @@
+*******************************************************************************
+ Revised: Jan 08, 2004             Io  / (Jupiter)                          501
+
+ SATELLITE PHYSICAL PROPERTIES:
+  Radius (km)             = 1821.3 +- 0.2   Density (g cm^-3)   =  3.530 +-.006
+  Mass (10^20 kg )        =  893.3 +- 1.5   Geometric Albedo    =  0.6  
+
+ SATELLITE ORBITAL DATA:
+  Semi-major axis, a (km) = 421.769 (10^3)  Orbital period      = 1.769138 d
+  Eccentricity, e         = 0.0041          Rotational period   = Synchronous
+  Inclination, i  (deg)   = 0.036
+*******************************************************************************
+ Select ... [E]phemeris, [F]tp, [M]ail, [R]edisplay, ?, <cr>: