aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
-rw-r--r--Helpers/ic_helpers.py35
1 files changed, 24 insertions, 11 deletions
diff --git a/Helpers/ic_helpers.py b/Helpers/ic_helpers.py
index adb91b7..43da479 100644
--- a/Helpers/ic_helpers.py
+++ b/Helpers/ic_helpers.py
@@ -14,6 +14,7 @@ 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 (Coordinate Time)
+ 19-03-2015, 00:00 (Barycentric Dynamical Time)
Earth = Earth (No Moon)
Format of IC (Initial Conditions) Lines:
@@ -25,18 +26,30 @@ def Solar2(epoch='2014-01-01'):
@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.')
+ if epoch == '2014-01-01':
+ 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 -9.947038156059579e-01 3.727217768375445E-02 3.850441538224456e-07 -9.268872371777320e-04 -1.726435348064618e-02 1.185854807637874e-06 0 0 0"
+ mars = "0 4 3.2125081695239055e-07 2.2660750299046701e-05 -1.5124387605410869e+00 +6.9681542097585980e-01 +5.1724065828569588e-02 -5.3309348873500006e-03 -1.1514537719341871e-02 -1.1040581507850680e-04 0 0 0"
+ jupiter = "0 5 9.5420039213714751e-04 4.6732616936774455e-04 -1.3307968217788140e+00 +5.0187260363223007e+00 +8.9354534335869748e-03 -7.3915532677911646e-03 -1.5771896615651140e-03 +1.7200784821425221e-04 0 0 0"
+ saturn = "0 6 2.8570710371524815e-04 3.8925687652332965e-04 -6.8851925827531151e+00 -7.0754782545976083e+00 +3.9713562554594872e-01 +3.6895716420435351e-03 -3.9070640598097067e-03 -7.8752021502553231e-05 0 0 0"
+ uranus = "0 7 4.3642853551857632e-05 1.6953449825499186e-04 +1.9645011082335920e+01 +3.9225997140108890e+00 -2.3986827787041001e-01 -8.0447801422034405e-04 +3.6714036954299551e-03 +2.4141556323116709e-05 0 0 0"
+ neptune = "0 8 5.1480569101603742e-05 1.6458790379443301e-04 +2.7063448807961450e+01 -1.2891858367498200e+01 -3.5810433631098909e-01 +1.3227724210187141e-03 +2.8500466203526960e-03 -8.9055141344433750e-05 0 0 0"
+ pluto = "0 9 6.5808657181639943e-09 7.7608056333903304e-06 +6.2567672814280284e+00 -3.1924879744181990e+01 +1.6063792641261310e+00 +3.1299238024423749e-03 -3.0015467471074471e-05 -8.9263215128274780e-04 0 0 0"
+
+ elif epoch == '2015-03-19':
+ mercury = "0 1 1.6515006786989092e-07 1.6310392545626536e-05 1.649966657839669e-01 -4.118561803126358e-01 -4.878999658120659e-02 2.048165231929986e-02 1.188092332449576e-02 -9.083224126400074e-04 0 0 0"
+ venus = "0 2 2.4468352521240763e-06 4.0455121182840900e-05 9.994563666017854e-02 7.131386043668728e-01 4.006403058090475e-03 -2.009951439735731e-02 2.705575613213286e-03 1.197005633552020e-03 0 0 0"
+ earth = "0 4 3.2125081695239055e-07 2.2660750299046701e-05 -9.946768738757796e-01 3.726128225503331e-02 1.801395888880666e-06 -9.240970641241845e-04 -1.725719489296372e-02 5.769142449771647e-07 0 0 0"
+ mars = "0 4 3.2125081695239055e-07 2.2660750299046701e-05 1.172066422620533e+00 8.400045843626092e-01 -1.116674189452623e-02 -7.615999183958316e-03 1.256925711308859e-02 4.502998204574119e-04 0 0 0"
+ jupiter = "0 5 9.5420039213714751e-04 4.6732616936774455e-04 -4.128843952022470e+00 3.392602891717091e+00 7.829934259091746e-02 -4.884365592005598e-03 -5.480051462380191e-03 1.320584982806127e-04 0 0 0"
+ saturn = "0 6 2.8570710371524815e-04 3.8925687652332965e-04 -5.066596590180945e+00 -8.576809170722843e+00 3.507980350698500e-01 4.495034108942656e-03 -2.857852798709087e-03 -1.292462726288397e-04 0 0 0"
+ uranus = "0 7 4.3642853551857632e-05 1.6953449825499186e-04 1.921964250612925e+01 5.528472084687846e+00 -2.283607905970542e-01 -1.119073386234832e-03 3.590779724777411e-03 2.787740493514162e-05 0 0 0"
+ neptune = "0 8 5.1480569101603742e-05 1.6458790379443301e-04 2.706344880796145e+01 -1.289185836749820e+01 -3.581043363109891e-01 1.322772421018714e-03 2.850046620352696e-03 -8.905514134443375e-05 0 0 0"
+ pluto = "0 9 6.5808657181639943e-09 7.7608056333903304e-06 7.637551631046986e+00 -3.191332944493440e+01 1.204987361866688e+00 3.118667965867572e-03 9.491586589905959e-05 -9.048694473505379e-04 0 0 0"
- 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"
- mars = "0 4 3.2125081695239055e-07 2.2660750299046701e-05 -1.5124387605410869e+00 +6.9681542097585980e-01 +5.1724065828569588e-02 -5.3309348873500006e-03 -1.1514537719341871e-02 -1.1040581507850680e-04 0 0 0"
- jupiter = "0 5 9.5420039213714751e-04 4.6732616936774455e-04 -1.3307968217788140e+00 +5.0187260363223007e+00 +8.9354534335869748e-03 -7.3915532677911646e-03 -1.5771896615651140e-03 +1.7200784821425221e-04 0 0 0"
- saturn = "0 6 2.8570710371524815e-04 3.8925687652332965e-04 -6.8851925827531151e+00 -7.0754782545976083e+00 +3.9713562554594872e-01 +3.6895716420435351e-03 -3.9070640598097067e-03 -7.8752021502553231e-05 0 0 0"
- uranus = "0 7 4.3642853551857632e-05 1.6953449825499186e-04 +1.9645011082335920e+01 +3.9225997140108890e+00 -2.3986827787041001e-01 -8.0447801422034405e-04 +3.6714036954299551e-03 +2.4141556323116709e-05 0 0 0"
- neptune = "0 8 5.1480569101603742e-05 1.6458790379443301e-04 +2.7063448807961450e+01 -1.2891858367498200e+01 -3.5810433631098909e-01 +1.3227724210187141e-03 +2.8500466203526960e-03 -8.9055141344433750e-05 0 0 0"
- pluto = "0 9 6.5808657181639943e-09 7.7608056333903304e-06 +6.2567672814280284e+00 -3.1924879744181990e+01 +1.6063792641261310e+00 +3.1299238024423749e-03 -3.0015467471074471e-05 -8.9263215128274780e-04 0 0 0"
+ else:
+ raise Exception('Unsupport Epoch/Date for Horizons Ephemeris.')
plist = [ mercury, venus, earth, mars, \
jupiter, saturn, uranus, neptune, \