diff options
| author | Volker Hoffmann <volker@cheleb.net> | 2016-05-14 09:11:11 +0200 |
|---|---|---|
| committer | Volker Hoffmann <volker@cheleb.net> | 2016-05-15 17:48:02 +0200 |
| commit | e8b57b52668921a4b7576ef88acb66c85419777a (patch) | |
| tree | 487ec6412053f3fcbbab4dd6db772d2e8bb2bd8a | |
| parent | a065bb844bb6a3045d7deb152b23588f89cc250a (diff) | |
feat: AST-1 support epoch selection for planet ephemeris from Horizons
| -rw-r--r-- | Helpers/ic_helpers.py | 11 |
1 files changed, 8 insertions, 3 deletions
diff --git a/Helpers/ic_helpers.py b/Helpers/ic_helpers.py index 2b891b4..9775c7d 100644 --- a/Helpers/ic_helpers.py +++ b/Helpers/ic_helpers.py @@ -9,17 +9,21 @@ import physics_helpers as ph import vector_helpers as vh -def Solar2(): +def Solar2(epoch='2014-01-01'): """ NASA/JPL Horizon Query for Solar System. Origin = Sun, Body Center, Ecliptic (J2000.0) Epoch = 01-01-2014, 00:00 (Central Time) Earth = Earth (No Moon) + @param: epoch - Ephemeris Epoch (Date) [String] @return: pnames - List of Planet Names [List of Strings] @return: plist - List of IC Lines for Planets [List of Strings] """ + if not epoch in [ '2014-01-01' ]: + raise Exception('Unsupport Epoch/Date for Horizons Ephemeris.') + mercury = "0 1 1.6515006786989092e-07 1.6310392545626536e-05 +1.1972692892595370e-01 -4.3426096069639630e-01 -4.6466929699041107e-02 +2.1482009903780699e-02 +8.9119212513564267e-03 -1.2427988748894680e-03 0 0 0" venus = "0 2 2.4468352521240763e-06 4.0455121182840900e-05 -4.9881994854843653e-02 +7.1760434138711926e-01 +1.2711854910906791e-02 -2.0246486048184902e-02 -1.5114700566799760e-03 +1.1477577544001651e-03 0 0 0" earth = "0 3 3.0023628776833749e-06 4.2587504470568303e-05 -1.7558926311807671e-01 +9.6752446474730014e-01 -2.9820874849228900e-05 -1.7207909676081830e-02 -3.1365065792022962e-03 +1.1244539739591660e-07 0 0 0" @@ -62,11 +66,12 @@ def Solar2(): return plist, pnames -def Solar2_Kepler(): +def Solar2_Kepler(epoch='2014-01-01'): """ Return Solar System ICs as Keplerian Elements. Usess Data from Solar2() Function. + @param: epoch - Ephemeris Epoch (Date) [String] @return: a - Semi-Major Axis (AU) [Numpy Float Array] @return: e - Eccentricity [Numpy Float Array] @return: i - Inclination (Radians) [Numpy Float Array] @@ -75,7 +80,7 @@ def Solar2_Kepler(): """ # Load NASA/JPL Horizons Data in Genga IC Format - plist, pnames = Solar2() + plist, pnames = Solar2(epoch) # Allocate Arrays mass = np.zeros(len(plist)) |