diff --git a/examples/symba_chambers_2013/init_cond.py b/examples/symba_chambers_2013/init_cond.py index 5cda5056e..5f6782491 100755 --- a/examples/symba_chambers_2013/init_cond.py +++ b/examples/symba_chambers_2013/init_cond.py @@ -1,7 +1,7 @@ #!/usr/bin/env python3 import swiftest import numpy as np -import os +from numpy.random import default_rng # Initialize simulation object sim = swiftest.Simulation() @@ -10,9 +10,11 @@ MU2KG = swiftest.MSun TU2S = swiftest.YR2S DU2M = swiftest.AU2M +GU = swiftest.GC / (DU2M**3 / (MU2KG * TU2S**2)) sim.param['MU2KG'] = MU2KG sim.param['TU2S'] = TU2S sim.param['DU2M'] = DU2M +sim.param['GU'] = GU # Simulation time parameters sim.param['T0'] = 0.0 @@ -35,10 +37,14 @@ sim.param['ROTATION'] = "YES" sim.param['CHK_RMAX'] = 1000.0 sim.param['CHK_EJECT'] = 1000.0 +sim.param['IN_FORM'] = 'EL' +sim.param['OUT_FORM'] = 'EL' # Add central body sim.add("Sun") -sim.add("Earth") +GMcb = sim.ds['GMass'].values[0] +sim.add("Jupiter") +sim.add("Saturn") # Add bodies described in Chambers (2013) Sec. 2.1, with the uniform spatial distribution and two bodies sizes (big and small) Nb = 14 @@ -49,7 +55,35 @@ Rb = (3 * Mb / (4 * np.pi * dens) )**(1.0 / 3.0) Rs = (3 * Ms / (4 * np.pi * dens) )**(1.0 / 3.0) +# Define the initial orbital elements of the big and small bodies +avalb = default_rng().uniform(0.3, 2.0, Nb) +avals = default_rng().uniform(0.3, 2.0, Ns) +evalb = default_rng().uniform(0.0, 0.001, Nb) +evals = default_rng().uniform(0.0, 0.001, Ns) +incvalb = default_rng().uniform(0.0, 0.0005 * 180 / np.pi, Nb) +incvals = default_rng().uniform(0.0, 0.0005 * 180 / np.pi, Ns) +capomvalb = default_rng().uniform(0.0, 360.0, Nb) +capomvals = default_rng().uniform(0.0, 360.0, Ns) +omegavalb = default_rng().uniform(0.0, 360.0, Nb) +omegavals = default_rng().uniform(0.0, 360.0, Ns) +capmvalb = default_rng().uniform(0.0, 360.0, Nb) +capmvals = default_rng().uniform(0.0, 360.0, Ns) +GMvalb = np.full(Nb, Mb * GU) +GMvals = np.full(Ns, Ms * GU) +Rvalb = np.full(Nb, Rb) +Rvals = np.full(Ns, Rs) +Rhb = avalb * (GMvalb / (3 * GMcb))**(1.0/3.0) +Rhs = avals * (GMvals / (3 * GMcb))**(1.0/3.0) +# Give the bodies unique ids +idb = np.arange(100, 100 + Nb) +ids = np.arange(100 + Nb, 100 + Nb + Ns) + +# Populate the simulation object with the two types of bodies +sim.addp(idb, avalb, evalb, incvalb, capomvalb, omegavalb, capmvalb, GMpl=GMvalb, Rpl=Rvalb, Rhill=Rhb) +sim.addp(ids, avals, evals, incvals, capomvals, omegavals, capmvals, GMpl=GMvals, Rpl=Rvals, Rhill=Rhs) + +# Save everything to a set of initial conditions files sim.save("param.in") diff --git a/examples/symba_chambers_2013/param.in b/examples/symba_chambers_2013/param.in index d096d13fe..19e98e87a 100644 --- a/examples/symba_chambers_2013/param.in +++ b/examples/symba_chambers_2013/param.in @@ -21,7 +21,7 @@ CHK_QMIN_RANGE -1.0 -1.0 MU2KG 1.988409870698051e+30 TU2S 31557600.0 DU2M 149597870700.0 -IN_FORM XV +IN_FORM EL ENC_OUT enc.dat EXTRA_FORCE NO DISCARD_OUT discard.out @@ -34,4 +34,5 @@ ROTATION YES TIDES NO ENERGY YES GR NO +GU 39.476926408897626 ENERGY_OUT energy.dat diff --git a/examples/symba_chambers_2013/pl_chambers_2013.in b/examples/symba_chambers_2013/pl_chambers_2013.in index 56d3f87e4..67d6989d2 100644 --- a/examples/symba_chambers_2013/pl_chambers_2013.in +++ b/examples/symba_chambers_2013/pl_chambers_2013.in @@ -1,7 +1,937 @@ -1 -3 0.000120026935827952456416 0.010044948626436044442 -4.25875607065040958e-05 -0.87405601066282501144 -0.50863684909036299686 2.1125685010193420466e-05 -3.058408673972023793 5.407722806790176349 -0.00026983464502404899522 -0.0 0.0 0.33069999999999999396 -4.774949887540911433 0.034051061217069346033 2301.2115366660092022 +156 +5 0.03769225108898567778 0.35527124629249239585 +0.00046732617030490929307 +5.2035240151214816606 0.048518677336151047608 1.303568872392780964 +100.51672855018520636 273.38350380059858935 317.24733983021741324 +0.0 0.0 0.27560000000000001164 +-80.96882015576837074 -2387.9972166758053216 5008.7356633679251408 +6 0.01128589982009127331 0.43766275227962878915 +0.00038925687730393611812 +9.581843040035108672 0.052193289679895830957 2.486259885154824989 +113.5952699440524043 335.70375572759667193 225.29914557362599226 +0.0 0.0 0.22000000000000000111 +441.93102714484842597 378.52937229654817222 5135.9113502360652523 +100 1.105353939449133635e-05 0.0076985326356758955327 +2.365334758890906669e-05 +1.6971843938499941995 0.0005816364466293669124 0.026187138610341612749 +199.66295678003845637 349.1838313868397563 239.5201319664316486 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +101 1.105353939449133635e-05 0.0033890547724920599881 +2.365334758890906669e-05 +0.7471359987644425882 0.0001484221862068907502 0.015912119531244738047 +7.72758829713570794 98.74578965353445881 292.48405218591960875 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +102 1.105353939449133635e-05 0.0032403861730961017119 +2.365334758890906669e-05 +0.71436117806925114593 0.0006667054143762592925 0.02623486192401118175 +42.268600860678539277 152.49698993249327827 151.45780278371265126 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +103 1.105353939449133635e-05 0.00848716742472066395 +2.365334758890906669e-05 +1.8710433251236489038 0.0009470568874946718764 0.022705939228401819663 +38.80730095910011812 143.34784343975428555 218.97185967315348876 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +104 1.105353939449133635e-05 0.004200287321695052973 +2.365334758890906669e-05 +0.9259767321153884456 0.00016780856600745607692 0.008111571611109064106 +275.52852024631755512 255.91035494922905968 283.72241775520802776 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +105 1.105353939449133635e-05 0.008200656366547455989 +2.365334758890906669e-05 +1.8078803667251066578 0.0005806820778782817395 0.0043050823832896005852 +262.76594876192444872 197.55827100651990236 359.12597183160295344 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +106 1.105353939449133635e-05 0.0031948304428970366643 +2.365334758890906669e-05 +0.70431816363998178687 0.0006166870106254423523 0.028312766571177344632 +100.38580777052690962 247.08615435134802851 28.274821850602428697 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +107 1.105353939449133635e-05 0.00645602527104512983 +2.365334758890906669e-05 +1.4232667256020528068 0.0003621838900236009598 0.012197930522699866954 +189.06372569967658137 282.1051532665543391 175.25556542499867874 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +108 1.105353939449133635e-05 0.0041177716184015826697 +2.365334758890906669e-05 +0.90778568578166862757 0.00088526942168412907277 0.0052077153285733724167 +355.61137098230727815 88.33816057346382422 35.585620637422806567 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +109 1.105353939449133635e-05 0.00862613046390682299 +2.365334758890906669e-05 +1.9016785010188288485 0.00044486325805035754379 0.0175194201364592822 +74.14924815238387623 142.82151728160823723 178.34909057749027284 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +110 1.105353939449133635e-05 0.004510164863223143653 +2.365334758890906669e-05 +0.9942909619010631417 0.00081444817122918432024 0.00707240791737578315 +68.31335729815678803 294.6246189285234891 122.92252049539324332 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +111 1.105353939449133635e-05 0.004651281910070013023 +2.365334758890906669e-05 +1.0254010007810474558 0.0004312636045460260673 0.01674019444109412641 +337.31604521335424351 171.74369728670052382 266.0027373696372024 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +112 1.105353939449133635e-05 0.0038244126458011504717 +2.365334758890906669e-05 +0.84311306651055351047 2.7875407573919153787e-05 0.017986416975313480882 +225.7229782864748131 308.33509788803951324 112.41779338822905743 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +113 1.105353939449133635e-05 0.0029643452153518679503 +2.365334758890906669e-05 +0.6535064116198798789 0.00037425733872322575592 0.0041403864217484038016 +325.7386399373389736 300.94512487157606984 333.7954171286211249 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +114 1.1053539394491335305e-06 0.0038455497484053390705 +1.0978911404624207699e-05 +1.826471255610046418 0.00087780859820767002683 0.009722642750510790757 +120.71268570369510087 352.6610169497902234 335.05498321080250435 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +115 1.1053539394491335305e-06 0.0017443709230869269861 +1.0978911404624207699e-05 +0.82850139995235116963 0.00031729137906426232822 0.01699303328109756997 +182.6703622068744437 213.3557405525360764 250.28249688543755269 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +116 1.1053539394491335305e-06 0.0020702769104938987164 +1.0978911404624207699e-05 +0.98329277101103618186 0.00058251749032120785 0.011797523288670572553 +312.4959308001525642 359.33855098907400816 162.46970124565356741 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +117 1.1053539394491335305e-06 0.00080530024674720453575 +1.0978911404624207699e-05 +0.38248309059826302025 0.0005969943570699224548 0.01380590580233940004 +48.708972362298950998 330.19391063831329802 102.38386339213280962 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +118 1.1053539394491335305e-06 0.002285280534438147877 +1.0978911404624207699e-05 +1.085410274275427911 0.00025130462779814554623 0.006036721527856849477 +45.867981468353910657 64.60926497341390018 253.24516188255682891 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +119 1.1053539394491335305e-06 0.002898627766350067431 +1.0978911404624207699e-05 +1.3767239126595514076 0.0002530591745570577502 0.006701891775872563864 +256.03674703821604908 136.96778530539529584 171.0933268407595449 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +120 1.1053539394491335305e-06 0.0014542058638660382582 +1.0978911404624207699e-05 +0.6906854374182271883 0.0006003363904860686438 0.004833197414767716714 +56.14180364632849063 273.3261312610867435 199.6250014185157795 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +121 1.1053539394491335305e-06 0.0024981790226714743672 +1.0978911404624207699e-05 +1.1865279283331482407 0.00082336466808019992576 0.018265656692325986377 +132.47659416127314103 48.61694507533226073 207.43492631988763719 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +122 1.1053539394491335305e-06 0.0028828564897507973385 +1.0978911404624207699e-05 +1.3692332324558205148 0.00040127628444019838947 0.007585985018489528435 +343.91044751539715207 196.83240087837353371 62.629566549905092643 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +123 1.1053539394491335305e-06 0.0038943677929657182489 +1.0978911404624207699e-05 +1.8496577337414494835 0.00043859595304759446763 0.011023337424484119101 +72.91657450189477174 344.12080798821472172 100.72699073828245275 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +124 1.1053539394491335305e-06 0.001675555086167977485 +1.0978911404624207699e-05 +0.7958168278400463702 0.00034926098253153968964 0.0099212607809325507185 +135.68666178212228601 296.83136368103197356 107.90664998734122548 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +125 1.1053539394491335305e-06 0.0023748127187634624592 +1.0978911404624207699e-05 +1.1279342232088616083 0.0007746490841826023793 0.017642382003092794973 +268.43319860669271293 252.98510234252850637 120.5395807841271818 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +126 1.1053539394491335305e-06 0.0012043684437486966761 +1.0978911404624207699e-05 +0.5720233744428824263 0.00034480510921434492748 0.017251867284702383787 +154.90093862136961889 356.5658098586184792 57.388355737939328094 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +127 1.1053539394491335305e-06 0.0007890385972653567634 +1.0978911404624207699e-05 +0.37475950429965448896 0.00031644523056148786714 0.027616371909429979908 +171.3578942788073789 299.01930445458822305 225.76712007928702519 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +128 1.1053539394491335305e-06 0.001777100120316001991 +1.0978911404624207699e-05 +0.84404636539789934524 0.00087953695475614056525 0.0006577653108544906282 +70.992958768969344874 130.55955257913041123 40.84173864503155471 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +129 1.1053539394491335305e-06 0.0040678503543026272943 +1.0978911404624207699e-05 +1.9320545124498424983 0.0005834355620338948882 0.0002608548084170523012 +270.88252512263818517 195.45419716249062958 140.3411838625577559 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +130 1.1053539394491335305e-06 0.0026829522072340253071 +1.0978911404624207699e-05 +1.2742872689972433164 0.00060055603766298343727 0.0064920160946518451484 +162.55438763596669105 178.06818403372230364 68.91417398625686985 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +131 1.1053539394491335305e-06 0.00071194132454311543046 +1.0978911404624207699e-05 +0.33814160524074998726 0.00091270005510230359236 0.0011893404294311912296 +250.26053020749216671 223.06302666485407826 192.24977287221634015 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +132 1.1053539394491335305e-06 0.002812359973542913556 +1.0978911404624207699e-05 +1.3357504097390577869 0.00026337894924635706691 0.027464957742126525075 +199.22420478899164209 171.93696705940942593 94.48037873747355775 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +133 1.1053539394491335305e-06 0.0028963111872931898259 +1.0978911404624207699e-05 +1.3756236369290164578 0.0006573629715563674607 0.021333593911504113161 +87.075565489431383526 153.68132922246732619 173.6496245283211124 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +134 1.1053539394491335305e-06 0.002785993277843155817 +1.0978911404624207699e-05 +1.3232273597327499637 0.00073604309841295738944 0.011001193704390446401 +339.40562860328299166 83.8138037375181284 177.19280309751334812 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +135 1.1053539394491335305e-06 0.0019038535805104330423 +1.0978911404624207699e-05 +0.9042488245366071542 0.0009361201833287777522 0.011953628791727587133 +153.92620389616467946 44.723571716565317047 318.60665806354387541 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +136 1.1053539394491335305e-06 0.0010869785583050439427 +1.0978911404624207699e-05 +0.5162682118549849486 0.00022103242271550004863 0.017892470794599633216 +219.23592171773313453 303.7464473661763691 122.639027547649121175 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +137 1.1053539394491335305e-06 0.0032939926784115201725 +1.0978911404624207699e-05 +1.5645052949330442438 0.000991170480938109201 0.019190182292000351055 +346.07086810942291777 54.79631567078042309 335.26549358631774567 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +138 1.1053539394491335305e-06 0.0010896488129816246724 +1.0978911404624207699e-05 +0.5175364683413183897 0.0002088585855651833357 0.025066384881598170048 +68.72020084550582908 232.28654212667152024 35.431071301701685172 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +139 1.1053539394491335305e-06 0.0008099901996317258652 +1.0978911404624207699e-05 +0.38471061714041626667 0.0005249334614017682433 0.023288096472191086028 +254.41357781115257808 117.993033232943290045 239.40340264257432068 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +140 1.1053539394491335305e-06 0.0026223707179148065143 +1.0978911404624207699e-05 +1.2455136589537159075 0.00084338241955415730895 0.0023232936840040861563 +269.3199834748557464 357.25116637136363806 20.298945789952838936 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +141 1.1053539394491335305e-06 0.0026159634140416132347 +1.0978911404624207699e-05 +1.2424704643219992573 0.00044638238521322806543 0.0012966745044005083872 +154.64047440242666198 236.78631050296556282 350.38453034745816694 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +142 1.1053539394491335305e-06 0.003456693657444609839 +1.0978911404624207699e-05 +1.6417812842988768196 0.0006303353332627267564 0.016812028430699579862 +144.4420533567427185 153.14985299997016455 271.69228258365166084 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +143 1.1053539394491335305e-06 0.0007301401336089313031 +1.0978911404624207699e-05 +0.34678525928756914753 0.00045501001474390856362 0.02594919197139722647 +116.332998235860372915 256.36327622207414834 161.64230688982479478 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +144 1.1053539394491335305e-06 0.0034349976066789206724 +1.0978911404624207699e-05 +1.6314765903860701801 0.000856330141764382333 0.014786839646800843326 +42.332000671474226294 138.1316800890192269 216.26130379036857221 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +145 1.1053539394491335305e-06 0.0025439979715687957412 +1.0978911404624207699e-05 +1.2082899646084366285 0.0004070939489201934802 0.004503954508818640984 +48.07928625757347163 276.96063567133671768 187.66327227579773762 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +146 1.1053539394491335305e-06 0.003859828194739848691 +1.0978911404624207699e-05 +1.8332529054419246517 0.0003598790557617636025 0.014970896975960788777 +271.0689558389924514 52.94488824830970941 1.9029339518958510169 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +147 1.1053539394491335305e-06 0.0040866242766629467682 +1.0978911404624207699e-05 +1.9409713206540264441 0.00011426903199890425854 0.011184164105792333138 +67.28912042774895497 289.72221855724853867 184.15766072550803756 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +148 1.1053539394491335305e-06 0.00316122359322290065 +1.0978911404624207699e-05 +1.5014456718372271027 0.00086962582925130050334 0.01914542568965043115 +250.74751474204416013 182.77230460121793953 167.04339861860248106 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +149 1.1053539394491335305e-06 0.00097271678791404052593 +1.0978911404624207699e-05 +0.46199876980165455187 0.00034328565960485866526 0.021012970927227814233 +297.8391700642537785 90.09431163085170624 45.254280687068330735 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +150 1.1053539394491335305e-06 0.0033482251435429574658 +1.0978911404624207699e-05 +1.5902633906967305055 0.00019727479448159425065 0.025485920919319147299 +39.29235960954184037 250.38592799994987104 142.50585121059737048 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +151 1.1053539394491335305e-06 0.0025466847865406780281 +1.0978911404624207699e-05 +1.2095660865250288385 0.00016751618714245874271 0.0016713720493250405416 +333.37472688335139992 145.95098571742994409 234.96494950660738255 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +152 1.1053539394491335305e-06 0.000719807963238783997 +1.0978911404624207699e-05 +0.3418779213453242427 0.00074674271745437542386 0.0011890436742735912591 +195.44114509251284062 283.06165279842355176 315.09159268499945483 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +153 1.1053539394491335305e-06 0.00157025465995761385 +1.0978911404624207699e-05 +0.74580364006202581884 0.0002094126871885914243 0.022074802063619617526 +153.41729982837850343 124.788130215813424684 122.51830102672674627 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +154 1.1053539394491335305e-06 0.001874785507921179697 +1.0978911404624207699e-05 +0.89044273632717141354 0.00085481426099495228573 0.021979234537361322377 +202.10509721048100573 320.04294945394877914 309.06206051991284767 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +155 1.1053539394491335305e-06 0.0019122234532341110998 +1.0978911404624207699e-05 +0.9082241552287277919 0.0004597575061389033409 0.020148347940879648976 +163.72826672182648622 305.80294187113088356 283.64822483177448476 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +156 1.1053539394491335305e-06 0.0039778936569127347484 +1.0978911404624207699e-05 +1.8893289380113937881 0.00049610582898170353487 0.0014385636633325129576 +157.23049053490680649 323.571089698826313 102.65627395667121391 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +157 1.1053539394491335305e-06 0.0019708158646589233432 +1.0978911404624207699e-05 +0.93605303855254251566 0.00019680356564931200973 0.020455589191385341319 +278.47203022847156717 255.0643652208907497 258.27522081614978333 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +158 1.1053539394491335305e-06 0.003967327173706975909 +1.0978911404624207699e-05 +1.8843103115182096019 0.0009097170816114348319 0.0018044172991811243192 +117.34832961121314554 234.13353289800707557 259.1454183769003521 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +159 1.1053539394491335305e-06 0.0013842445575880614169 +1.0978911404624207699e-05 +0.6574568164714708063 1.4540080482081997083e-05 0.009892527806868892798 +133.30970394955991765 89.751557170898067284 193.05936528552655318 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +160 1.1053539394491335305e-06 0.0035162555550013580372 +1.0978911404624207699e-05 +1.670070631969415409 0.00091422718885765909876 0.028436761433306107644 +353.37301866523387162 51.552218623340351655 332.64071506725503014 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +161 1.1053539394491335305e-06 0.0019396124736002073748 +1.0978911404624207699e-05 +0.9212327656201906301 0.00057552402003210683355 0.00950969025161889879 +283.93076060868486366 267.5628271717484381 143.52189652919906848 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +162 1.1053539394491335305e-06 0.0031661554024814644338 +1.0978911404624207699e-05 +1.503788069787651116 0.00018115407388301875608 0.02283324079244578031 +177.01583573865022458 234.49703162631655573 336.92462545405476249 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +163 1.1053539394491335305e-06 0.0033722481017907932155 +1.0978911404624207699e-05 +1.6016732658992518523 0.00053832824323583156837 0.017356963073014760462 +210.48701348074004613 252.36586170021701037 241.77353994833666206 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +164 1.1053539394491335305e-06 0.0018549151304582611073 +1.0978911404624207699e-05 +0.8810051589588920873 0.00039455774939218580828 0.014186890964832121059 +268.78131437644134394 234.86755048529201417 135.82402939345988102 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +165 1.1053539394491335305e-06 0.0011704042480178788963 +1.0978911404624207699e-05 +0.5558918376585835208 5.564504104615730423e-05 0.023585563294025682635 +48.050891171653226763 199.47056209951909977 219.93699201694707313 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +166 1.1053539394491335305e-06 0.0040699423683521412583 +1.0978911404624207699e-05 +1.9330481294299255435 0.00052600080391025812437 0.021919559654086234746 +39.096865156286462195 128.8303250981957433 331.23155076245132022 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +167 1.1053539394491335305e-06 0.0028820213749619458248 +1.0978911404624207699e-05 +1.3688365887359974504 0.00026969964684439974665 0.017275229592356205338 +330.69633084561246505 147.5243863767787218 153.22005163183507648 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +168 1.1053539394491335305e-06 0.003980615057015220282 +1.0978911404624207699e-05 +1.8906214863822126127 0.00036100510956251820067 0.01521335462518577826 +340.23428559641661195 130.00201841928748081 169.2953836810935968 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +169 1.1053539394491335305e-06 0.0012767502410615831949 +1.0978911404624207699e-05 +0.60640162485459569996 0.00020696978447336191083 0.0056857512796010872477 +181.57921819815450704 115.817629846161409546 320.0491270864600324 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +170 1.1053539394491335305e-06 0.002399637727624715802 +1.0978911404624207699e-05 +1.1397250380654750934 0.0006442965482497394836 0.002941317764456898541 +333.45417971179358574 185.18080690840253055 269.38711004629510626 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +171 1.1053539394491335305e-06 0.0035439123990565288466 +1.0978911404624207699e-05 +1.6832064471304613029 0.00072650143891821498086 0.02076627535231315505 +82.57443269249242235 19.640854808506546192 343.42020390750877823 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +172 1.1053539394491335305e-06 0.0021893049972247140684 +1.0978911404624207699e-05 +1.0398260089737605405 0.00095820846394839131327 0.0047667113246031159526 +231.34174007276843099 346.25149299156703364 359.3927866655562866 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +173 1.1053539394491335305e-06 0.003808558580359920666 +1.0978911404624207699e-05 +1.8089020367553392266 0.00080486608391464400675 0.026061626097182894607 +333.73026949841113264 290.72940981027193175 353.58605983201931622 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +174 1.1053539394491335305e-06 0.0036046989829116192528 +1.0978911404624207699e-05 +1.7120774682852628601 0.000408327401359263474 0.023141693393884565227 +190.15243963894701551 10.9019960190977514 206.20261856852459914 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +175 1.1053539394491335305e-06 0.0025289888373368178023 +1.0978911404624207699e-05 +1.201161269352922556 0.0008119307695628591708 0.02330748815658333778 +293.74129631246364625 228.8273707230899845 50.496213420126203175 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +176 1.1053539394491335305e-06 0.0035813927117425375896 +1.0978911404624207699e-05 +1.7010079887177620428 0.0008261188698801552469 0.0033817239341808955078 +270.3371253433707011 303.21090622901044753 55.41742303760840116 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +177 1.1053539394491335305e-06 0.0029864597633055024009 +1.0978911404624207699e-05 +1.4184403454864715766 0.0006100544512734605053 0.0031745002208599282408 +128.80488878502711714 46.8536787103431962 190.83918585108236243 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +178 1.1053539394491335305e-06 0.0023567486788243735888 +1.0978911404624207699e-05 +1.119354578719118809 0.00053509985437350090966 0.020556636907334471104 +93.86909946810600047 329.53011231398716063 68.53400707162742833 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +179 1.1053539394491335305e-06 0.0034140907141169699521 +1.0978911404624207699e-05 +1.6215467127855098806 0.0007155502442420363365 0.021636466572293039379 +15.484181145735732343 319.8815184503466753 126.01422430110272899 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +180 1.1053539394491335305e-06 0.0025675289816375823053 +1.0978911404624207699e-05 +1.2194661855177959531 0.00048216964694683286767 0.017262139448479693615 +335.34868702893140835 186.94351621903143723 26.64334186676326155 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +181 1.1053539394491335305e-06 0.0034248039486638383743 +1.0978911404624207699e-05 +1.6266350398738742555 0.00081883719141492166717 0.003694641256859914659 +218.92219023702784853 218.29951433993917931 259.11630406545350525 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +182 1.1053539394491335305e-06 0.0028977454974775277425 +1.0978911404624207699e-05 +1.3763048727716702402 0.00040844368595388948033 0.020281148678402180496 +289.91356665087704414 165.61033347014003425 288.23615258185753873 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +183 1.1053539394491335305e-06 0.0011783105169288593343 +1.0978911404624207699e-05 +0.5596469764163169991 0.00026931995881909574905 0.0010082118613559110359 +154.06847791408287662 39.901773134600148296 187.47166921424076236 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +184 1.1053539394491335305e-06 0.002491583136722370224 +1.0978911404624207699e-05 +1.1833951653006802296 0.00048583037529452476142 0.027591397048547783993 +149.2419352144119955 63.397457239281493457 295.88510332284619153 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +185 1.1053539394491335305e-06 0.0039004676796874942429 +1.0978911404624207699e-05 +1.8525549184065084329 0.0005185930845717148926 0.019438436095256288877 +292.5205667956682305 105.44538380282074286 123.70099114339642199 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +186 1.1053539394491335305e-06 0.0034369089438783143244 +1.0978911404624207699e-05 +1.6323843936087221174 0.00010509062853767581647 0.0032505177502269231964 +143.85173028989797217 189.31567203063778493 223.19791686246551876 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +187 1.1053539394491335305e-06 0.004128326569038383845 +1.0978911404624207699e-05 +1.9607781215797892571 0.0007065471846638792094 0.021817391680966469192 +267.38223511629621498 257.29519900708868363 77.6961655587950446 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +188 1.1053539394491335305e-06 0.0020499009381892692037 +1.0978911404624207699e-05 +0.97361505776992030015 0.00093181583258145170724 0.0006097968157406988537 +145.87863557509786006 138.78820490355309403 73.09644258927524163 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +189 1.1053539394491335305e-06 0.0013791485719345369476 +1.0978911404624207699e-05 +0.65503644177237263513 0.00014046622778450534863 0.015442926715237644333 +113.70608882665497674 243.3631092677804304 81.13225672390240106 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +190 1.1053539394491335305e-06 0.0024853665935430262675 +1.0978911404624207699e-05 +1.1804425738198291995 0.00096027971308601327814 0.023639897027866934448 +334.99323146372387328 76.84265628244520485 155.49317295086675017 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +191 1.1053539394491335305e-06 0.003654973676990583536 +1.0978911404624207699e-05 +1.7359557924853059596 0.00054611453698988137676 0.026616017411125021674 +173.99931107991119461 359.07752431447528352 78.019305599684969366 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +192 1.1053539394491335305e-06 0.003360253875484283402 +1.0978911404624207699e-05 +1.5959765226465592036 0.00089913403289640294203 0.014826170842651739679 +204.39386065798206005 17.787041065446423715 47.062099086162277217 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +193 1.1053539394491335305e-06 0.0036411616313956547123 +1.0978911404624207699e-05 +1.7293956630081683823 0.0004269520101434239978 0.011446361145646259325 +225.78458330355539374 323.21819225965253963 348.78801121716753642 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +194 1.1053539394491335305e-06 0.0026716566036685272943 +1.0978911404624207699e-05 +1.2689223415936377037 0.00043093138180629422026 0.007198866767697694423 +49.75485744273438371 283.14995135249375835 251.76161931871857291 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +195 1.1053539394491335305e-06 0.0021962918766307708317 +1.0978911404624207699e-05 +1.0431444771347480049 0.00086121939260353213164 0.004059279588003008102 +286.63158878609289104 179.04178066528515956 146.49836003344253754 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +196 1.1053539394491335305e-06 0.0029447439089289198142 +1.0978911404624207699e-05 +1.3986271031916246432 0.00018272629779208073367 0.018596268807723595473 +131.73550112256867806 122.63667230710116485 118.00875344636077102 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +197 1.1053539394491335305e-06 0.0036073925872425995972 +1.0978911404624207699e-05 +1.7133568148563382838 0.00075563849669159832353 0.0037783795914888887967 +128.81574846279312396 26.084496990571235386 338.71113698392537117 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +198 1.1053539394491335305e-06 0.0035771279460317551398 +1.0978911404624207699e-05 +1.6989824078535733509 0.0007581912267808143613 0.009405275951210755395 +94.72942154944794879 234.18898767453580945 120.281591054773343785 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +199 1.1053539394491335305e-06 0.0035858601416377149508 +1.0978911404624207699e-05 +1.7031298263804728332 0.00095813031245370230537 0.0036113709533022904374 +146.16269587061592006 115.18323006911683137 268.9501223857537866 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +200 1.1053539394491335305e-06 0.0037229358585092736563 +1.0978911404624207699e-05 +1.7682349148822829044 0.0009838787345463682714 0.0041291372775860849323 +1.7590124360204661969 28.267489335889418811 140.55607459951542637 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +201 1.1053539394491335305e-06 0.001651848468190447193 +1.0978911404624207699e-05 +0.78455720070308299086 0.00013641617523719662343 0.019427991831481188079 +130.8007266845419565 110.181352575198033605 106.00139496728822053 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +202 1.1053539394491335305e-06 0.0027063636542107174076 +1.0978911404624207699e-05 +1.2854067025640301125 0.00080079146808039282975 0.0108589657051855435765 +46.771264060493685122 32.497203081429894667 53.366093926082044163 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +203 1.1053539394491335305e-06 0.0017526730483667882186 +1.0978911404624207699e-05 +0.8324445535132765084 9.953863767625981303e-05 0.027893078159839147095 +116.997421293990342406 10.0440500665721010876 63.77086305348634454 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +204 1.1053539394491335305e-06 0.0013663552841492629148 +1.0978911404624207699e-05 +0.6489601785763914332 0.0007547835818466555732 0.028400265298084318571 +183.67148924061424964 352.2190315881923084 274.3720362883656776 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +205 1.1053539394491335305e-06 0.0016416918432609031937 +1.0978911404624207699e-05 +0.7797332393187537125 0.00076392586396015573803 0.020535102200658142402 +221.9135945165130579 194.21835875557150075 137.1052917217631375 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +206 1.1053539394491335305e-06 0.0041503320593118873817 +1.0978911404624207699e-05 +1.9712297859918352927 0.00025071053464458946825 0.014913864673346191939 +319.79234195814888153 349.28387326164579463 138.32956542869692385 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +207 1.1053539394491335305e-06 0.0017463766922070864303 +1.0978911404624207699e-05 +0.8294540542886734702 0.00021197645574354129889 0.018247549715991118885 +105.00262688872396666 232.80424007039977141 84.949104918041584256 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +208 1.1053539394491335305e-06 0.0037525194240774044083 +1.0978911404624207699e-05 +1.7822858401553345331 0.00080294190095561387456 0.023444883786946393162 +356.27577921605848132 11.694296440931752912 96.8991254541982272 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +209 1.1053539394491335305e-06 0.0033673867988775121108 +1.0978911404624207699e-05 +1.5993643554399337159 0.00045136667731674606644 0.027694134961144371149 +122.09817852172396613 331.56919146778068352 23.29494948461385917 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +210 1.1053539394491335305e-06 0.002601544338314537602 +1.0978911404624207699e-05 +1.235622021558025363 0.00052748956446136828925 0.014390768296716713662 +47.688625110281321895 218.81863989673610149 71.091061616253483635 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +211 1.1053539394491335305e-06 0.0034572364291313065697 +1.0978911404624207699e-05 +1.6420390775791540161 0.00012278039477614099018 0.013364601572277041985 +37.099447743349671214 25.711537395277808571 46.771237173257247832 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +212 1.1053539394491335305e-06 0.0014692126295102979798 +1.0978911404624207699e-05 +0.69781300769612641943 0.00020831235893638435199 0.0105463257356791569735 +45.282702057561735387 269.00333228724332457 141.05705036356778237 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +213 1.1053539394491335305e-06 0.0012506227963414040885 +1.0978911404624207699e-05 +0.5939922088058925409 0.0005111276050050357622 0.006922367726696107386 +193.68387748392143521 15.642417562054848901 17.908891439446314564 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +214 1.1053539394491335305e-06 0.0041476312180326669947 +1.0978911404624207699e-05 +1.9699470021806242404 0.00058751583022024968673 0.02496086898048604022 +247.93088919368153711 124.53408226858653052 258.30383247130868085 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +215 1.1053539394491335305e-06 0.0029875489394385165544 +1.0978911404624207699e-05 +1.4189576574520941232 0.0004909660898418238828 0.0023047035739488575846 +238.65714156827573333 164.51271027675758774 223.70978205956649276 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +216 1.1053539394491335305e-06 0.0023921676837597391677 +1.0978911404624207699e-05 +1.1361770875017920268 0.00023970389209557418876 0.026923977514779456888 +58.044978520057981086 79.883336583186263624 94.38325642624329248 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +217 1.1053539394491335305e-06 0.0026816001073412998205 +1.0978911404624207699e-05 +1.2736450795184046658 0.0007910064207286296922 0.028329584831009391177 +35.296218138003794706 187.99468430047886613 356.03440184202310093 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +218 1.1053539394491335305e-06 0.001331364345185033964 +1.0978911404624207699e-05 +0.6323409827770205549 0.00028739168058106059726 0.009019690678395121675 +338.27263574838821114 230.96954437756605216 126.72414489683163197 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +219 1.1053539394491335305e-06 0.0032121580786483427635 +1.0978911404624207699e-05 +1.5256373686388187139 0.0002481645483046961251 0.0066040479374564345708 +60.055442506513337264 346.99690826233506868 179.83160228793084912 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +220 1.1053539394491335305e-06 0.0041836951887061120445 +1.0978911404624207699e-05 +1.9870758420364929009 0.00038625231628182366255 0.014342679397268678562 +335.15357535816582413 299.7702497912664512 317.6921811746384492 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +221 1.1053539394491335305e-06 0.0007589345479948146076 +1.0978911404624207699e-05 +0.3604613715832816201 0.0008130249542857412217 0.0132530748254426049976 +199.40136765025368959 44.661857657232609142 303.65462581497462224 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +222 1.1053539394491335305e-06 0.0011646339137693681694 +1.0978911404624207699e-05 +0.55315117628902576996 0.00058320518586504638704 0.025565557231859135229 +88.51316300040842577 64.3025295238989969 307.91222823830349853 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +223 1.1053539394491335305e-06 0.004180113624911285161 +1.0978911404624207699e-05 +1.9853747527906451253 0.00094918793041903390915 0.01898810058981277582 +339.11739337814219653 326.5281938064337055 348.13403023570242567 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +224 1.1053539394491335305e-06 0.002907903028903174131 +1.0978911404624207699e-05 +1.3811292647027828551 0.00017854471163964626287 0.01757285626967412806 +324.68564254821177428 173.27665949324585881 212.81811478340242161 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +225 1.1053539394491335305e-06 0.0030424513946792162559 +1.0978911404624207699e-05 +1.4450339697923877402 0.00083102791963415834596 0.016743369452382143164 +72.37540677365963404 25.69124406916907688 36.560572732805965757 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +226 1.1053539394491335305e-06 0.0014152066797417541456 +1.0978911404624207699e-05 +0.6721624970181487857 0.0004194325219653943738 0.02466884060344504137 +312.10194954625785613 75.06433419730545609 213.15995579142392558 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +227 1.1053539394491335305e-06 0.00074807664212413927804 +1.0978911404624207699e-05 +0.35530433181877674187 0.0006196274885235493945 0.0057712733317203191233 +112.28915034524422367 308.71008190145914796 250.26404455810629202 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +228 1.1053539394491335305e-06 0.0015240975037274374223 +1.0978911404624207699e-05 +0.7238809698040012952 0.00038945125604979613423 0.0034764685260635502795 +100.18804344114703042 289.83313211239982365 254.06542835983233886 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +229 1.1053539394491335305e-06 0.0013613878149108384452 +1.0978911404624207699e-05 +0.64660084366442682757 0.0009530460362832646236 0.00080068641185022341056 +290.2519853457544059 126.69300624511200226 129.6226360067927601 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +230 1.1053539394491335305e-06 0.0007497560275515179048 +1.0978911404624207699e-05 +0.3561019678944690181 6.2453564361704131855e-05 0.008646200757580674995 +51.77354471090342969 304.1036845586841082 280.42024752572388024 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +231 1.1053539394491335305e-06 0.0020677692895160400756 +1.0978911404624207699e-05 +0.9821017585588054377 0.00092676675850360832123 0.006735246833832509404 +263.97158005173145057 104.45971633042557869 156.07773457283062157 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +232 1.1053539394491335305e-06 0.00095787667212301207604 +1.0978911404624207699e-05 +0.4549503510590604538 0.00052515340653565799334 0.019061642884057490577 +284.31553463421840888 213.15965776812532795 3.9051805628102265189 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +233 1.1053539394491335305e-06 0.0037735814895644029364 +1.0978911404624207699e-05 +1.7922894182423769749 0.00017129142152418842175 0.025294321601488489043 +70.01841692477869117 111.37830955236358932 314.1754665499188377 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +234 1.1053539394491335305e-06 0.0007847310095064422164 +1.0978911404624207699e-05 +0.37271358479856425694 0.0003231942299594923215 0.01628512097921978874 +224.60279807470539026 103.680766227009058866 197.02880532870963748 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +235 1.1053539394491335305e-06 0.0015320077326663605468 +1.0978911404624207699e-05 +0.72763798940522450565 0.00094252324011187993847 0.01692070785353958634 +73.3265875193739447 291.12276980910894508 97.10035705120415628 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +236 1.1053539394491335305e-06 0.0017944863754330241527 +1.0978911404624207699e-05 +0.85230409115664484965 3.7676763751743336793e-05 0.021560198984918479276 +216.79211114566814445 74.30371365646033155 295.38458279934900474 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +237 1.1053539394491335305e-06 0.0022077289096507573257 +1.0978911404624207699e-05 +1.0485765774655604332 0.00092835655556198343457 0.002473122352105824439 +207.85994871205122081 242.10322219988140091 164.79329124450524091 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +238 1.1053539394491335305e-06 0.0017780507509334266778 +1.0978911404624207699e-05 +0.84449787418364552494 0.00012149108796338048226 0.022688746545291040707 +105.00259291531052952 236.90787466426584729 258.88778385063903897 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +239 1.1053539394491335305e-06 0.0009633437960757620414 +1.0978911404624207699e-05 +0.4575470005380317673 0.00021936857586009873258 0.024584544078017318142 +141.46008032785664454 156.81449724386195044 299.68311915392922629 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +240 1.1053539394491335305e-06 0.0035009114125315640493 +1.0978911404624207699e-05 +1.6627828221641500761 0.00083380110219331913645 0.016595030400960369271 +68.173689361921276486 154.75502404200838669 311.70470139259805364 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +241 1.1053539394491335305e-06 0.0006470573938361051598 +1.0978911404624207699e-05 +0.30732452000176924134 0.00021338219268733039035 0.006499688048810576234 +359.85986519138003814 244.23351243313086911 98.01858239568616682 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +242 1.1053539394491335305e-06 0.00065901961137865321196 +1.0978911404624207699e-05 +0.31300606046393031345 0.00084676255106667984527 0.00028332565690039374258 +44.286246649765288907 20.245198329963866257 276.88008765639568765 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +243 1.1053539394491335305e-06 0.003806998050224602309 +1.0978911404624207699e-05 +1.8081608518475495995 0.00027745564680026148143 0.021943110426273559899 +266.24251634239237774 252.3739085379230005 283.488252705583875 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +244 1.1053539394491335305e-06 0.0029050019908176934909 +1.0978911404624207699e-05 +1.3797513959918767945 0.0008840097714271473667 0.011119175742524006062 +218.83470399188448141 307.5328062303744332 53.607438100582484708 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +245 1.1053539394491335305e-06 0.0041506206755627920196 +1.0978911404624207699e-05 +1.9713668663368231204 0.00072844746599102529454 0.014679484165514528754 +216.30736619477607974 127.43423126496449527 280.14881244016112305 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +246 1.1053539394491335305e-06 0.0007277677857464638177 +1.0978911404624207699e-05 +0.34565849576542212462 0.00059173314782810218417 0.007831481826486607831 +243.4161800124451247 113.710960693445571223 254.21623664849960278 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +247 1.1053539394491335305e-06 0.0017717047625134947208 +1.0978911404624207699e-05 +0.84148380176337678193 0.00052751195938902041473 0.014082285132383451365 +194.45249461995499018 318.31107915541781495 157.54667905036907882 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +248 1.1053539394491335305e-06 0.004169267887052278984 +1.0978911404624207699e-05 +1.980223492309989286 0.00067309606461422974606 0.017396151114924100073 +184.39865388000126245 210.62388541157199029 167.14339720153236613 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +249 1.1053539394491335305e-06 0.0019522337566574509366 +1.0978911404624207699e-05 +0.92722733394492307113 0.00058907724893334585843 0.021547757024762807104 +190.98028164050376176 93.657838091120083845 154.55544413858271469 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +250 1.1053539394491335305e-06 0.004103199565563551535 +1.0978911404624207699e-05 +1.9488438722295042815 0.00023535689304344910105 0.020747164931323787418 +235.4258950080400723 210.09500044028584398 329.76890939133602387 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +251 1.1053539394491335305e-06 0.00096993472194017628167 +1.0978911404624207699e-05 +0.46067740774293175487 6.722535556689789992e-05 0.0099933940728082100435 +339.95817378663616637 181.27958061831904502 34.322050877546928405 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +252 1.1053539394491335305e-06 0.0009661516009396302622 +1.0978911404624207699e-05 +0.45888058746597204518 0.0005025170114595179727 0.0024236270953698216615 +57.831023059303113598 16.636326610996761133 83.793613608447557795 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 +253 1.1053539394491335305e-06 0.0032298402438196435668 +1.0978911404624207699e-05 +1.5340356389864386877 0.00062884613640746395407 0.013163328547706066446 +348.2310414025072305 241.09038942507976344 294.8270165591735008 +0.4000000000000000222 0.4000000000000000222 0.4000000000000000222 +0.0 0.0 0.0 diff --git a/python/swiftest/swiftest/init_cond.py b/python/swiftest/swiftest/init_cond.py index 297b9df68..cf8f109cc 100644 --- a/python/swiftest/swiftest/init_cond.py +++ b/python/swiftest/swiftest/init_cond.py @@ -102,7 +102,6 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds): 'Pluto': np.longdouble(0.4) } - # Unit conversion factors DCONV = swiftest.AU2M / param['DU2M'] VCONV = (swiftest.AU2M / swiftest.JD2S) / (param['DU2M'] / param['TU2S']) @@ -274,8 +273,22 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds): return ds -def vec2xr(param, idvals, v1, v2, v3, v4, v5, v6, GMpl=None, Rpl=None, Rhill=None, t=0.0): - +def vec2xr(param, idvals, v1, v2, v3, v4, v5, v6, GMpl=None, Rpl=None, Rhill=None, Ip_x=None, Ip_y=None, Ip_z=None, rot_x=None, rot_y=None, rot_z=None, t=0.0): + + if param['ROTATION'] == 'YES': + if Ip_x is None: + Ip_x = np.full_like(v1, 0.4) + if Ip_y is None: + Ip_y = np.full_like(v1, 0.4) + if Ip_z is None: + Ip_z = np.full_like(v1, 0.4) + if rot_x is None: + rot_x = np.full_like(v1, 0.0) + if rot_y is None: + rot_y = np.full_like(v1, 0.0) + if rot_z is None: + rot_z = np.full_like(v1, 0.0) + dims = ['time', 'id', 'vec'] if GMpl is not None: ispl = True @@ -290,18 +303,18 @@ def vec2xr(param, idvals, v1, v2, v3, v4, v5, v6, GMpl=None, Rpl=None, Rhill=Non return None clab, plab, tlab = swiftest.io.make_swiftest_labels(param) + vec = np.vstack([v1, v2, v3, v4, v5, v6]) if ispl: + vec = np.vstack([vec, GMpl, Rpl]) if param['RHILL_PRESENT'] == 'YES': - vec = np.vstack([v1, v2, v3, v4, v5, v6, GMpl, Rpl, Rhill]).T - else: - vec = np.vstack([v1, v2, v3, v4, v5, v6, GMpl, Rpl]).T - else: - vec = np.vstack([v1, v2, v3, v4, v5, v6]).T - bodyframe = np.expand_dims(vec, axis=0) + vec = np.vstack([vec, Rhill]) + if param['ROTATION'] == 'YES': + vec = np.vstack([vec, Ip_x, Ip_y, Ip_z, rot_x, rot_y, rot_z]) + bodyframe = np.expand_dims(vec.T, axis=0) if ispl: - bodyxr = xr.DataArray(bodyframe, dims=dims, coords={'time': [t], 'id': tpid, 'vec': plab}) + bodyxr = xr.DataArray(bodyframe, dims=dims, coords={'time': [t], 'id': idvals, 'vec': plab}) else: - bodyxr = xr.DataArray(bodyframe, dims=dims, coords={'time': [t], 'id': tpid, 'vec': tlab}) + bodyxr = xr.DataArray(bodyframe, dims=dims, coords={'time': [t], 'id': idvals, 'vec': tlab}) ds = bodyxr.to_dataset(dim='vec') return ds \ No newline at end of file diff --git a/python/swiftest/swiftest/simulation_class.py b/python/swiftest/swiftest/simulation_class.py index e7b8c560b..e0ec698a2 100644 --- a/python/swiftest/swiftest/simulation_class.py +++ b/python/swiftest/swiftest/simulation_class.py @@ -4,6 +4,7 @@ from swiftest import constants from datetime import date import xarray as xr +import numpy as np class Simulation: """ @@ -52,6 +53,7 @@ def __init__(self, codename="Swiftest", param_file=""): self.codename = codename if param_file != "" : self.read_param(param_file, codename) + self.param['GU'] = constants.GC / (param['DU2M'])**3 / (param['MU2KG'] * param['TU2S']**2) return @@ -73,7 +75,7 @@ def add(self, plname, date=date.today().isoformat(), idval=None): return - def addp(self, idvals, t1, t2, t3, t4, t5, t6, Gmass=None, radius=None, Rhill=None, t=0.0): + def addp(self, idvals, t1, t2, t3, t4, t5, t6, GMpl=None, Rpl=None, Rhill=None, Ip_x=None, Ip_y=None, Ip_z=None, rot_x=None, rot_y=None, rot_z=None): """ Adds a body (test particle or massive body) to the internal DataSet given a set up 6 vectors (orbital elements or cartesian state vectors, depending on the value of self.param). Input all angles in degress @@ -90,16 +92,18 @@ def addp(self, idvals, t1, t2, t3, t4, t5, t6, Gmass=None, radius=None, Rhill=No Gmass : Optional: Array of G*mass values if these are massive bodies radius : Optional: Array radius values if these are massive bodies Rhill : Optional: Array Rhill values if these are massive bodies + Ip_x,y,z : Optional: Principal axes moments of inertia + rot_x,y,z: Optional: Rotation rate vector components Returns ------- self.ds : xarray dataset """ - dsnew = init_cond.vec2xr(self.param, idvals, t1, t2, t3, t4, t5, t6, Gmass, radius, Rhill) + t = self.param['T0'] + + dsnew = init_cond.vec2xr(self.param, idvals, t1, t2, t3, t4, t5, t6, GMpl, Rpl, Rhill, Ip_x, Ip_y, Ip_z, rot_x, rot_y, rot_z, t) if dsnew is not None: self.ds = xr.combine_by_coords([self.ds, dsnew]) return - - def read_param(self, param_file, codename="Swiftest"):