diff options
author | Volker Hoffmann <volker@cheleb.net> | 2016-05-14 17:07:57 +0200 |
---|---|---|
committer | Volker Hoffmann <volker@cheleb.net> | 2016-05-15 17:48:10 +0200 |
commit | 636af5d113e8c4977d8e4314e5c188e4215c3d6f (patch) | |
tree | 73aecae50e1b6154e4896c6b3d948c950d86508f | |
parent | 5d2a39d986350f304023bb245016503db51d688d (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.py | 32 |
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 \ |