aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorVolker Hoffmann <volker@cheleb.net>2016-05-14 17:07:57 +0200
committerVolker Hoffmann <volker@cheleb.net>2016-05-15 17:48:10 +0200
commit636af5d113e8c4977d8e4314e5c188e4215c3d6f (patch)
tree73aecae50e1b6154e4896c6b3d948c950d86508f
parent5d2a39d986350f304023bb245016503db51d688d (diff)
feat: AST-1 add Solar System ICs for 19-03-2015 to IC generator scripts (called Solar3)
-rw-r--r--IC_Helpers/genic_solar.py32
1 files changed, 30 insertions, 2 deletions
diff --git a/IC_Helpers/genic_solar.py b/IC_Helpers/genic_solar.py
index b28a7f1..b561191 100644
--- a/IC_Helpers/genic_solar.py
+++ b/IC_Helpers/genic_solar.py
@@ -35,6 +35,8 @@ group3.add_argument('--solar1', action='store_true', \
help="Base Solar System: Joachim")
group3.add_argument('--solar2', action='store_true', \
help="Base Solar System: NASA Horizron Query 2014-01-01.")
+group3.add_argument('--solar3', action='store_true', \
+ help="Base Solar System: NASA Horizron Query 2015-03-19.")
args = parser.parse_args()
# Joachim's ICs
@@ -54,10 +56,10 @@ if args.solar1:
# NASA Horizon Query
# Origin = Sun, Body Center
-# Epoch = 01-01-2014, 00:00 (Central Time)
+# Epoch = 01-01-2014, 00:00 (Coordinate Time)
# Earth = Earth (No Moon)
if args.solar2:
- plist, _ = ih.Solar2()
+ plist, _ = ih.Solar2(epoch='2014-01-01')
# Spaghetti code ahead
for iplanet, planet in enumerate(plist):
if iplanet == 0:
@@ -79,6 +81,32 @@ if args.solar2:
elif iplanet == 8:
pluto = planet
+# NASA Horizon Query
+# Origin = Sun, Body Center
+# Epoch = 19-03-2015, 00:00 (Barycentric Dynamical Time)
+# Earth = Earth (No Moon)
+if args.solar3:
+ plist, _ = ih.Solar2(epoch='2015-03-19')
+ # Spaghetti code ahead
+ for iplanet, planet in enumerate(plist):
+ if iplanet == 0:
+ mercury = planet
+ elif iplanet == 1:
+ venus = planet
+ elif iplanet == 2:
+ earth = planet
+ elif iplanet == 3:
+ mars = planet
+ elif iplanet == 4:
+ jupiter = planet
+ elif iplanet == 5:
+ saturn = planet
+ elif iplanet == 6:
+ uranus = planet
+ elif iplanet == 7:
+ neptune = planet
+ elif iplanet == 8:
+ pluto = planet
# Adjust Orbital Parameters
if args.cjs or args.ejs or \