diff --git a/.gitignore b/.gitignore index 99c34ac64..8f2237e27 100644 --- a/.gitignore +++ b/.gitignore @@ -9,7 +9,7 @@ !CHANGELOG !README.md !README.swifter -!example/*.in +!*.in !example/cleanup !example/swifter_symba_omp !Makefile diff --git a/examples/rmvs_gr_test/cb.swiftest.in b/examples/rmvs_gr_test/cb.swiftest.in new file mode 100644 index 000000000..2386b53c8 Binary files /dev/null and b/examples/rmvs_gr_test/cb.swiftest.in differ diff --git a/examples/rmvs_gr_test/config.swiftest.in b/examples/rmvs_gr_test/config.swiftest.in new file mode 100644 index 000000000..6e623a5dd --- /dev/null +++ b/examples/rmvs_gr_test/config.swiftest.in @@ -0,0 +1,29 @@ +! Swiftest input file generated using init_cond.py +T0 0 +TSTOP 1.0 +DT 0.0006844626967830253 +CB_IN cb.swiftest.in +PL_IN pl.swiftest.in +TP_IN tp.swiftest.in +IN_TYPE REAL8 +ISTEP_OUT 1 +ISTEP_DUMP 1 +BIN_OUT bin.swiftest.dat +OUT_TYPE REAL8 +OUT_FORM EL +OUT_STAT REPLACE +CHK_CLOSE yes +CHK_RMIN 0.004650467260962157 +CHK_RMAX 1000.0 +CHK_EJECT 1000.0 +CHK_QMIN 0.004650467260962157 +CHK_QMIN_COORD HELIO +CHK_QMIN_RANGE 0.004650467260962157 1000.0 +ENC_OUT enc.swiftest.dat +EXTRA_FORCE no +BIG_DISCARD no +ROTATION no +GR no +MU2KG 1.988409870698051e+30 +DU2M 149597870700.0 +TU2S 31557600.0 diff --git a/examples/rmvs_gr_test/param.swifter.in b/examples/rmvs_gr_test/param.swifter.in new file mode 100644 index 000000000..5834d2dcc --- /dev/null +++ b/examples/rmvs_gr_test/param.swifter.in @@ -0,0 +1,26 @@ +! Swifter input file generated using init_cond.py +T0 0 +TSTOP 1.0 +DT 0.0006844626967830253 +PL_IN pl.swifter.in +TP_IN tp.swifter.in +IN_TYPE ASCII +ISTEP_OUT 1 +ISTEP_DUMP 1 +BIN_OUT bin.swifter.dat +OUT_TYPE REAL8 +OUT_FORM EL +OUT_STAT NEW +J2 4.7535806948127355e-12 +J4 -2.2473967953572827e-18 +CHK_CLOSE yes +CHK_RMIN 0.004650467260962157 +CHK_RMAX 1000.0 +CHK_EJECT 1000.0 +CHK_QMIN 0.004650467260962157 +CHK_QMIN_COORD HELIO +CHK_QMIN_RANGE 0.004650467260962157 1000.0 +ENC_OUT enc.swifter.dat +EXTRA_FORCE no +BIG_DISCARD no +RHILL_PRESENT yes diff --git a/examples/rmvs_gr_test/param.tpcollider.in b/examples/rmvs_gr_test/param.tpcollider.in new file mode 100644 index 000000000..a6066b36c --- /dev/null +++ b/examples/rmvs_gr_test/param.tpcollider.in @@ -0,0 +1,26 @@ +! Swifter input file generated using init_cond.py +T0 0.0 +TSTOP 146100.0 +DT 0.1 +PL_IN pl.tpcollider.in +TP_IN tp.tpcollider.in +IN_TYPE ASCII +ISTEP_OUT 14610 +ISTEP_DUMP 14610 +BIN_OUT bin.tpcollider.dat +OUT_TYPE REAL8 +OUT_FORM XV +OUT_STAT NEW +J2 4.7535806948127355e-12 +J4 -2.2473967953572827e-18 +CHK_CLOSE yes +CHK_RMIN 0.004650467260962157 +CHK_RMAX 10000.0 +CHK_EJECT 10000.0 +CHK_QMIN 0.004650467260962157 +CHK_QMIN_COORD HELIO +CHK_QMIN_RANGE 0.004650467260962157 10000.0 +ENC_OUT enc.tpcollider.dat +EXTRA_FORCE no +BIG_DISCARD no +RHILL_PRESENT yes diff --git a/examples/rmvs_gr_test/pl.swifter.in b/examples/rmvs_gr_test/pl.swifter.in new file mode 100644 index 000000000..d0d4e7ff9 --- /dev/null +++ b/examples/rmvs_gr_test/pl.swifter.in @@ -0,0 +1,40 @@ +10 ! Planet input file generated using init_cond.py using JPL Horizons data for the major planets (and Pluto) for epoch 2021-01-28 +1 39.47692640889762629 +0.0 0.0 0.0 +0.0 0.0 0.0 +2 6.553709809565313959502e-06 0.0014751229680863789154 +1.6306381826061645943e-05 +0.1030256860922895 0.2897796047098886 0.01422904600374035 +-11.74004209950937 3.8343124110162736 1.3902496665973592 +3 9.6633133995815387361564e-05 0.006759127649782299051 +4.0453784346544178454e-05 +0.06110218027254217 -0.7245466901305982 -0.01346904300924688 +7.311995449678243 0.5941125721336201 -0.4137913843379075 +4 0.00012002693582795246295385 0.0100447565675466429165 +4.25875607065040958e-05 +-0.6061796342297583 0.7761214554702035 -3.4750047790977e-05 +-5.054824314301841 -3.891667468503358 0.00019720338148272726 +5 1.2739802010675942316241e-05 0.0072464490746299085006 +2.265740805092889601e-05 +0.2751944175855944 1.51937688993241 0.02508924593104206 +-4.835983593209577 1.344855094041679 0.14681413000004515 +6 0.037692251088985682938581 0.3552852357486060849 +0.00046732617030490929307 +3.200135438345358 -3.953498213518368 -0.05517737289975112 +2.111393749129838 1.8660266890185446 -0.05498941067210089 +7 0.011285899820091272946487 0.43763064566943408597 +0.00038925687730393611812 +5.607382165725712 -8.258649105608766 -0.07958445228024298 +1.5748468603228847 1.1414574661825514 -0.08250331331320372 +8 0.0017236589478267728883093 0.4690969274244374022 +0.00016953449859497231466 +15.28225422201768 12.53905314208462 -0.1514143582550325 +-0.9198472198098231 1.0454390993472462 0.01574538863031621 +9 0.0020336100526728304385693 0.7807192056765467829 +0.000164587904124493665 +29.47483071169769 -5.147686530859088 -0.5733441819169969 +0.19191677740340274 1.1385110364087574 -0.027844325148353527 +10 2.9242167710294538257026e-07 0.05383468172776979939 +7.943294877391593783e-06 +14.14000920780611 -31.14141812522779 -0.7565722591093476 +1.073396108697069 0.23003123192799815 -0.33424529561177047 diff --git a/examples/rmvs_gr_test/pl.swiftest.in b/examples/rmvs_gr_test/pl.swiftest.in new file mode 100644 index 000000000..7bdc4a619 Binary files /dev/null and b/examples/rmvs_gr_test/pl.swiftest.in differ diff --git a/examples/rmvs_gr_test/pl.tpcollider.in b/examples/rmvs_gr_test/pl.tpcollider.in new file mode 100644 index 000000000..95f1a048c --- /dev/null +++ b/examples/rmvs_gr_test/pl.tpcollider.in @@ -0,0 +1,36 @@ +9 ! Planet input file generated using init_cond.py using JPL Horizons data for the major planets (and Pluto) for epoch 2421-01-31 +1 0.0002959122081920778 +0.0 0.0 0.0 +0.0 0.0 0.0 +3 7.243452483873648e-10 0.006759130149202631 +4.0453784346544176e-05 +0.7019627011222136 -0.1852496375751218 -0.04300948405789237 +-0.005060711738168192 -0.01946939433513342 1.612614183356907e-06 +4 8.99701138216602e-10 0.010044559796067993 +4.25875607065041e-05 +-0.5697193205332233 0.8029811715067949 -0.0007045655351608364 +0.01430660935370767 0.01002546009563165 -1.100857588556753e-05 +5 9.549535102761467e-11 0.00724618889821577 +2.2657408050928896e-05 +0.7993560103471902 -1.148987408482768 -0.04313430212910231 +-0.01198004310618487 -0.009212920234676632 8.927002293279674e-05 +6 2.825345908631355e-07 0.3551533878488701 +0.0004673261703049093 +-4.634097296304171 -2.863163670218678 0.1152887094960275 +-0.003870142070986203 0.006066584891983585 5.910227945707925e-05 +7 8.459715183006417e-08 0.4376906786688663 +0.00038925687730393614 +-8.865265590469349 2.837147222179521 0.3099061219991513 +0.001998614714254602 0.005352670386838502 -0.0001703654134002045 +8 1.2920249163736676e-08 0.46813372268314535 +0.00016953449859497232 +15.47151038891844 -12.61881964551067 -0.2448817967119467 +-0.002462937202468351 -0.002853377332777701 2.133002929086309e-05 +9 1.5243589003230837e-08 0.7791996556565636 +0.00016458790412449367 +-24.87809567256067 17.00453520872124 0.2234899970988031 +0.001786881544518242 0.002576864647985829 -9.421562352086046e-05 +10 2.19194228290428e-12 0.05348153138291406 +7.943294877391593e-06 +-9.921482526550596 41.28101896209601 -1.54188968720156 +0.002154849797934256 0.001120743083473351 -0.0007433230817862677 diff --git a/examples/rmvs_gr_test/tp.swifter.in b/examples/rmvs_gr_test/tp.swifter.in new file mode 100644 index 000000000..3e12e22a6 --- /dev/null +++ b/examples/rmvs_gr_test/tp.swifter.in @@ -0,0 +1,31 @@ +10 +100 +0.3089139948287037 -0.277270436883302 -0.052337497440874486 +4.590828296619012 9.736244632162235 0.28280742360614397 +101 +0.3070354816876597 -0.2845400603406744 -0.05318033091888456 +4.545139601540877 9.883525912822861 0.28793674166112104 +102 +0.31066213352368505 -0.28113935471688967 -0.053108274441649984 +4.617453752609237 9.800810875613966 0.28446343275638736 +103 +0.3088067674199909 -0.2869523372878359 -0.05191462227061065 +4.6363938524337565 9.713660761176554 0.2789015757022737 +104 +0.307074553418099 -0.28073329893098137 -0.0534617731635447 +4.579194784145406 9.574048169595772 0.2861631995442269 +105 +0.30294954170866395 -0.281481925092296 -0.052234290208198066 +4.604562331430815 9.719598701455027 0.282459654106244 +106 +0.30692242784366397 -0.2878888119399549 -0.05247924519789774 +4.633660571617613 9.88377173149324 0.28490916253088794 +107 +0.30868976460301095 -0.27766108923594196 -0.05385210649760836 +4.581218660618022 9.756537560396604 0.2826382827782646 +108 +0.30635827269855037 -0.2829786891311232 -0.05244294704858842 +4.599148326746042 9.647983245028241 0.28233738073637155 +109 +0.313848717338811 -0.28389321799892847 -0.05290541615450685 +4.556080023696361 9.729321212583539 0.29047556209904657 diff --git a/examples/rmvs_gr_test/tp.swiftest.in b/examples/rmvs_gr_test/tp.swiftest.in new file mode 100644 index 000000000..d91e0934e Binary files /dev/null and b/examples/rmvs_gr_test/tp.swiftest.in differ diff --git a/examples/rmvs_gr_test/tp.tpcollider.in b/examples/rmvs_gr_test/tp.tpcollider.in new file mode 100644 index 000000000..0b3169633 --- /dev/null +++ b/examples/rmvs_gr_test/tp.tpcollider.in @@ -0,0 +1,4 @@ +1 +100 +0.3568867004762329 -0.128562827221147 -0.0429414817481014 +-0.004043366894290556 -0.02767355464504326 -0.001909056361553056 diff --git a/examples/rmvs_swifter_comparison/1pl_1tp_encounter/cb.swiftest.in b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/cb.swiftest.in new file mode 100644 index 000000000..2386b53c8 Binary files /dev/null and b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/cb.swiftest.in differ diff --git a/examples/rmvs_swifter_comparison/1pl_1tp_encounter/config.swiftest.in b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/config.swiftest.in new file mode 100644 index 000000000..914af3324 --- /dev/null +++ b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/config.swiftest.in @@ -0,0 +1,29 @@ +! Swiftest input file generated using init_cond.py +T0 0 +TSTOP 0.2 +DT 0.00034223134839151266 +CB_IN cb.swiftest.in +PL_IN pl.swiftest.in +TP_IN tp.swiftest.in +IN_TYPE REAL8 +ISTEP_OUT 1 +ISTEP_DUMP 1 +BIN_OUT bin.swiftest.dat +OUT_TYPE REAL8 +OUT_FORM XV +OUT_STAT REPLACE +CHK_CLOSE yes +CHK_RMIN 0.004650467260962157 +CHK_RMAX 1000.0 +CHK_EJECT 1000.0 +CHK_QMIN 0.004650467260962157 +CHK_QMIN_COORD HELIO +CHK_QMIN_RANGE 0.004650467260962157 1000.0 +ENC_OUT enc.swiftest.dat +EXTRA_FORCE no +BIG_DISCARD no +ROTATION no +GR no +MU2KG 1.988409870698051e+30 +DU2M 149597870700.0 +TU2S 31557600.0 diff --git a/examples/rmvs_swifter_comparison/1pl_1tp_encounter/param.swifter.in b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/param.swifter.in new file mode 100644 index 000000000..927ac2883 --- /dev/null +++ b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/param.swifter.in @@ -0,0 +1,26 @@ +! Swifter input file generated using init_cond.py +T0 0 +TSTOP 0.2 +DT 0.00034223134839151266 +PL_IN pl.swifter.in +TP_IN tp.swifter.in +IN_TYPE ASCII +ISTEP_OUT 1 +ISTEP_DUMP 1 +BIN_OUT bin.swifter.dat +OUT_TYPE REAL8 +OUT_FORM XV +OUT_STAT NEW +J2 4.7535806948127355e-12 +J4 -2.2473967953572827e-18 +CHK_CLOSE yes +CHK_RMIN 0.004650467260962157 +CHK_RMAX 1000.0 +CHK_EJECT 1000.0 +CHK_QMIN 0.004650467260962157 +CHK_QMIN_COORD HELIO +CHK_QMIN_RANGE 0.004650467260962157 1000.0 +ENC_OUT enc.swifter.dat +EXTRA_FORCE no +BIG_DISCARD no +RHILL_PRESENT yes diff --git a/examples/rmvs_swifter_comparison/1pl_1tp_encounter/pl.swifter.in b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/pl.swifter.in new file mode 100644 index 000000000..6f91ef4c9 --- /dev/null +++ b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/pl.swifter.in @@ -0,0 +1,8 @@ +2 ! Planet input file generated using init_cond.py +1 39.47692640889762629 +0.0 0.0 0.0 +0.0 0.0 0.0 +2 0.00012002693582795246295385 0.010044724833237895015 +4.25875607065041e-05 +1.0 0.0 0.0 +0.0 6.283185307179586 0.0 diff --git a/examples/rmvs_swifter_comparison/1pl_1tp_encounter/pl.swiftest.in b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/pl.swiftest.in new file mode 100644 index 000000000..d3786c3df Binary files /dev/null and b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/pl.swiftest.in differ diff --git a/examples/rmvs_swifter_comparison/1pl_1tp_encounter/tp.swifter.in b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/tp.swifter.in new file mode 100644 index 000000000..9c026369e --- /dev/null +++ b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/tp.swifter.in @@ -0,0 +1,4 @@ +1 +100 +1.01 0.0 0.0 +0.0 6.252003053624663 0.0 diff --git a/examples/rmvs_swifter_comparison/1pl_1tp_encounter/tp.swiftest.in b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/tp.swiftest.in new file mode 100644 index 000000000..e1506974a Binary files /dev/null and b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/tp.swiftest.in differ diff --git a/examples/rmvs_swifter_comparison/mars_ejecta/cb.swiftest.in b/examples/rmvs_swifter_comparison/mars_ejecta/cb.swiftest.in new file mode 100644 index 000000000..2aafc8408 --- /dev/null +++ b/examples/rmvs_swifter_comparison/mars_ejecta/cb.swiftest.in @@ -0,0 +1,4 @@ +0.00029591220828563 +0.004650467260962157 +0.0 +0.0 diff --git a/examples/rmvs_swifter_comparison/mars_ejecta/config.swiftest.in b/examples/rmvs_swifter_comparison/mars_ejecta/config.swiftest.in new file mode 100644 index 000000000..685ee4595 --- /dev/null +++ b/examples/rmvs_swifter_comparison/mars_ejecta/config.swiftest.in @@ -0,0 +1,33 @@ +! +! Parameter file for the CHO run of the 4 giant planets and Pluto. +! +!NPLMAX -1 ! not used +!NTPMAX -1 ! not used +T0 0.0e0 +TSTOP 365.25e5 ! simulation length in days +DT 2e0 ! stepsize in days +ISTEP_OUT 73500 ! output cadence +ISTEP_DUMP 73500 ! system dump cadence +CB_IN cb.swiftest.in +PL_IN pl.swiftest.in +TP_IN tp.in +IN_TYPE ASCII +BIN_OUT bin.swiftest.dat +OUT_TYPE REAL8 ! double precision real output +OUT_FORM XV ! osculating element output +OUT_STAT REPLACE +CHK_CLOSE yes ! check for planetary close encounters +CHK_RMIN 0.005e0 ! check for close solar encounters +CHK_RMAX 10000.0e0 ! discard outside of +CHK_EJECT -1.0 ! ignore this check +CHK_QMIN -1.0 ! ignore this check +!CHK_QMIN_COORD HELIO ! commented out here +!CHK_QMIN_RANGE 1.0 1000.0 ! commented out here +ENC_OUT enc.swiftest.dat +EXTRA_FORCE no ! no extra user-defined forces +BIG_DISCARD yes ! output all planets if anything discarded +ROTATION no +GR no +MU2KG 1.988409870698051e+30 +DU2M 149597870700.0 +TU2S 86400.0000 diff --git a/examples/rmvs_swifter_comparison/mars_ejecta/param.swifter.in b/examples/rmvs_swifter_comparison/mars_ejecta/param.swifter.in new file mode 100644 index 000000000..5327ebcab --- /dev/null +++ b/examples/rmvs_swifter_comparison/mars_ejecta/param.swifter.in @@ -0,0 +1,30 @@ +! +! Parameter file for the CHO run of the 4 giant planets and Pluto. +! +!NPLMAX -1 ! not used +!NTPMAX -1 ! not used +T0 0.0e0 +TSTOP 365.25e5 ! simulation length in days +DT 2e0 ! stepsize in days +ISTEP_OUT 73500 ! output cadence +ISTEP_DUMP 73500 ! system dump cadence +PL_IN pl.swifter.in +TP_IN tp.in +IN_TYPE ASCII +BIN_OUT bin.swifter.dat +OUT_TYPE REAL8 ! double precision real output +OUT_FORM XV ! osculating element output +OUT_STAT NEW +J2 0.0 ! no J2 term +J4 0.0 ! no J4 term +CHK_CLOSE yes ! check for planetary close encounters +CHK_RMIN 0.005 ! check for close solar encounters +CHK_RMAX 10000.0 ! discard outside of +CHK_EJECT -1.0 ! ignore this check +CHK_QMIN -1.0 ! ignore this check +!CHK_QMIN_COORD HELIO ! commented out here +!CHK_QMIN_RANGE 1.0 1000.0 ! commented out here +ENC_OUT enc.swifter.dat +EXTRA_FORCE no ! no extra user-defined forces +BIG_DISCARD yes ! output all planets if anything discarded +RHILL_PRESENT yes ! Hill's sphere radii in input file diff --git a/examples/rmvs_swifter_comparison/mars_ejecta/pl.swifter.in b/examples/rmvs_swifter_comparison/mars_ejecta/pl.swifter.in new file mode 100644 index 000000000..ab7f6ec08 --- /dev/null +++ b/examples/rmvs_swifter_comparison/mars_ejecta/pl.swifter.in @@ -0,0 +1,36 @@ + 9 ! Solar system. Taken from JPL Horizons on 10-31-2012 00:00 (JD 2456231.5) + 1 2.9591220828563E-04 + .0 .0 .0 + .0 .0 .0 + 2 4.912576812862670E-11 1.475E-03 ! Mercury / Mass (Msun) / Rhill (AU) + 1.63104E-05 ! Radius (AU) + 3.347381871776144E-01 -2.106110537919978E-01 -4.792146121659843E-02 + 9.457239437449661E-03 2.510612575783642E-02 1.183593414707043E-03 + 3 7.243495778974390E-10 6.758E-03 ! Venus + 4.04538E-05 + -4.641166544325086E-01 5.473160240817795E-01 3.428518329144122E-02 + -1.549753193570599E-02 -1.319081534035637E-02 7.136685719563098E-04 + 4 8.997065158526820E-10 1.004E-02 ! Earth + Moon + 4.2587571E-05 + 7.844742229036105E-01 6.083466458892074E-01 -1.996912082982283E-05 + -1.081828010368904E-02 1.352675083743791E-02 2.330162751015549E-07 + 5 9.549592181631700E-11 7.246E-03 ! Mars + 2.26601E-05 + 3.248808297463504E-01 -1.392041319192106E+00 -3.714224198825128E-02 + 1.415678377567019E-02 4.380992498628390E-03 -2.558013868476858E-04 + 6 2.825362796108150E-07 0.3553 ! Jupiter + 4.67326E-04 + 1.873343037544675E+00 4.683322528183729E+00 -6.137078942429344E-02 + -7.104028854569467E-03 3.164193073442964E-03 1.458211211387848E-04 + 7 8.459765747874750E-08 0.4356 ! Saturn + 3.89257E-04 + -8.251865479985982E+00 -5.225008639958182E+00 4.193935221241906E-01 + 2.677374423933877E-03 -4.723961471952217E-03 -2.458643348128228E-05 + 8 1.292032638982330E-08 0.4686 ! Uranus + 1.69534E-04 + 1.992234896320000E+01 2.342619564376139E+00 -2.493547557095487E-01 + -4.936420700174518E-04 3.724868172847281E-03 2.027439791557233E-05 + 9 1.524368011623570E-08 0.7757 ! Neptune + 1.64601E-04 + 2.647158936403647E+01 -1.409652133490650E+01 -3.196444184043846E-01 + 1.449117851099529E-03 2.791630281764020E-03 -9.063088654405710E-05 diff --git a/examples/rmvs_swifter_comparison/mars_ejecta/pl.swiftest.in b/examples/rmvs_swifter_comparison/mars_ejecta/pl.swiftest.in new file mode 100644 index 000000000..e1931fc7f --- /dev/null +++ b/examples/rmvs_swifter_comparison/mars_ejecta/pl.swiftest.in @@ -0,0 +1,33 @@ +8 +2 4.91257681286267e-11 +1.63104e-05 +0.3347381871776144 -0.2106110537919978 -0.04792146121659843 +0.00945723943744966 0.02510612575783642 0.001183593414707043 +3 7.24349577897439e-10 +4.04538e-05 +-0.4641166544325086 0.5473160240817795 0.03428518329144122 +-0.01549753193570599 -0.01319081534035637 0.0007136685719563098 +4 8.99706515852682e-10 +4.2587571e-05 +0.7844742229036105 0.6083466458892074 -1.996912082982283e-05 +-0.01081828010368904 0.01352675083743791 2.330162751015549e-07 +5 9.5495921816317e-11 +2.26601e-05 +0.3248808297463504 -1.392041319192106 -0.03714224198825128 +0.01415678377567019 0.00438099249862839 -0.0002558013868476858 +6 2.82536279610815e-07 +0.000467326 +1.873343037544675 4.683322528183729 -0.06137078942429344 +-0.007104028854569467 0.003164193073442964 0.0001458211211387848 +7 8.45976574787475e-08 +0.000389257 +-8.251865479985982 -5.225008639958182 0.4193935221241906 +0.002677374423933877 -0.004723961471952217 -2.458643348128228e-05 +8 1.29203263898233e-08 +0.000169534 +19.9223489632 2.342619564376139 -0.2493547557095487 +-0.0004936420700174518 0.003724868172847281 2.027439791557233e-05 +9 1.52436801162357e-08 +0.000164601 +26.47158936403647 -14.0965213349065 -0.3196444184043846 +0.001449117851099529 0.00279163028176402 -9.06308865440571e-05 diff --git a/examples/rmvs_swifter_comparison/mars_ejecta/start.in b/examples/rmvs_swifter_comparison/mars_ejecta/start.in new file mode 100644 index 000000000..d6d7c3850 --- /dev/null +++ b/examples/rmvs_swifter_comparison/mars_ejecta/start.in @@ -0,0 +1 @@ +param.in diff --git a/examples/rmvs_swifter_comparison/mars_ejecta/swiftest_rmvs_vs_swifter_rmvs.ipynb b/examples/rmvs_swifter_comparison/mars_ejecta/swiftest_rmvs_vs_swifter_rmvs.ipynb index df30cd7ee..d82127068 100644 --- a/examples/rmvs_swifter_comparison/mars_ejecta/swiftest_rmvs_vs_swifter_rmvs.ipynb +++ b/examples/rmvs_swifter_comparison/mars_ejecta/swiftest_rmvs_vs_swifter_rmvs.ipynb @@ -13,7 +13,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 2, "metadata": {}, "outputs": [ { @@ -32,9 +32,17 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 3, "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Reading Swiftest file config.swiftest.in\n" + ] + } + ], "source": [ "config_file_name = 'config.swiftest.in'\n", "config = swio.read_swiftest_config(config_file_name)\n", @@ -43,16 +51,683 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 4, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "
<xarray.Dataset>\n",
+       "Dimensions:  (id: 19, time: 5)\n",
+       "Coordinates:\n",
+       "  * id       (id) int64 0 2 3 4 5 6 ... 6000205 6000206 6000207 6000208 6000209\n",
+       "  * time     (time) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05\n",
+       "Data variables:\n",
+       "    Mass     (time, id) float64 0.0002959 4.913e-11 7.243e-10 ... nan nan nan\n",
+       "    Radius   (time, id) float64 0.00465 1.631e-05 4.045e-05 ... nan nan nan\n",
+       "    J_2      (time, id) float64 0.0 nan nan nan nan nan ... nan nan nan nan nan\n",
+       "    J_4      (time, id) float64 0.0 nan nan nan nan nan ... nan nan nan nan nan\n",
+       "    px       (time, id) float64 nan 0.3347 -0.4641 ... -0.9303 0.1796 1.422\n",
+       "    py       (time, id) float64 nan -0.2106 0.5473 ... -1.162 -1.411 -0.5606\n",
+       "    pz       (time, id) float64 nan -0.04792 0.03429 ... -0.02492 0.04703\n",
+       "    vx       (time, id) float64 nan 0.009457 -0.0155 ... 0.01503 0.007312\n",
+       "    vy       (time, id) float64 nan 0.02511 -0.01319 ... 0.001562 0.0129\n",
+       "    vz       (time, id) float64 nan 0.001184 0.0007137 ... 0.0005511 0.0005797
" + ], + "text/plain": [ + "\n", + "Dimensions: (id: 19, time: 5)\n", + "Coordinates:\n", + " * id (id) int64 0 2 3 4 5 6 ... 6000205 6000206 6000207 6000208 6000209\n", + " * time (time) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05\n", + "Data variables:\n", + " Mass (time, id) float64 0.0002959 4.913e-11 7.243e-10 ... nan nan nan\n", + " Radius (time, id) float64 0.00465 1.631e-05 4.045e-05 ... nan nan nan\n", + " J_2 (time, id) float64 0.0 nan nan nan nan nan ... nan nan nan nan nan\n", + " J_4 (time, id) float64 0.0 nan nan nan nan nan ... nan nan nan nan nan\n", + " px (time, id) float64 nan 0.3347 -0.4641 ... -0.9303 0.1796 1.422\n", + " py (time, id) float64 nan -0.2106 0.5473 ... -1.162 -1.411 -0.5606\n", + " pz (time, id) float64 nan -0.04792 0.03429 ... -0.02492 0.04703\n", + " vx (time, id) float64 nan 0.009457 -0.0155 ... 0.01503 0.007312\n", + " vy (time, id) float64 nan 0.02511 -0.01319 ... 0.001562 0.0129\n", + " vz (time, id) float64 nan 0.001184 0.0007137 ... 0.0005511 0.0005797" + ] + }, + "execution_count": 4, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "swiftestdat" ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 5, "metadata": {}, "outputs": [], "source": [ @@ -61,7 +736,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 6, "metadata": {}, "outputs": [], "source": [ @@ -70,39 +745,1415 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 7, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "[,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ,\n", + " ]" + ] + }, + "execution_count": 7, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEzCAYAAADAeS+TAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABtvElEQVR4nO3dd3hUVfrA8e+Znt5DAgkJIaEGooJSVFRQxN4Ae8Wyq2vbXcsW9WfvrljYVbFXii72joKLIEUIRHqVkN7b1HvP74+ZhCEEEkgmM4HzeZ55ZubOLecO5L5zyn2PkFKiKIqiKPtjCHYBFEVRlNCngoWiKIrSLhUsFEVRlHapYKEoiqK0SwULRVEUpV0qWCiKoijtMgW7AAciMTFRZmZmBrsYiqIoPcqKFSsqpJRJndlHjwoWmZmZLF++PNjFUBRF6VGEEDs6uw/VDKUoiqK0SwULRVEUpV0qWCiKoijt6lF9FoqiHNrcbjeFhYU4HI5gF6VHstlspKWlYTabu3zfKlgoihIyCgsLiYqKIjMzEyFEsIvTo0gpqayspLCwkH79+nX5/lUzlKIoIcPhcJCQkKACxUEQQpCQkBCwWpkKFoqihBQVKPZmd7k7tF4gvzsVLBRFUfZj7NixbS6/6qqrmDt3bsCPX2O3U1VRTmVDY8CPtT8qWCiKouzHzz//HLRja7pOfW0dUhiICQ8LWjlAdXAriqLsV2RkJA0NDUgpufnmm5k/fz79+vWjO2YZLa6rx6hrhMXEYDIE97e9qlkoiqJ0wH//+182bNjAmjVreOWVVwJe46hzuaGpCWkyExseHtBjdYQKFoqiKB2wcOFCLr74YoxGI71792b8+PEBO5YmJZW1tQgkSXGxIdHpr5qhFEVROqi7LtpFDU1Y3S6s4eFYAnCD3cFQNQtFUZQOGDduHB988AGaplFcXMwPP/wQkOPUuj3ojQ0gDMRFRwfkGAdD1SwURVE64LzzzmP+/PkMGzaMAQMGcMIJJ3T5MTy6pKy+gQhdIyY2FkOQO7X9qWChKIqyHw0NDYC3CeqFF14I6LF2OZyEOewYzWbCw4I7VLa10AlbiqIoh7EatwdPYwMGJPGxodGp7U8FC0VRlCBz6zrFTXZsbhfh4eEByRrbWSpYKIqiBJGUkl0ON+H2JgxCEB1Cndr+VLBQFEUJohqPhsNhx6xrREdHh1Sntr/QLJWiKMphwK3r7LI7iXTaMZvNhIfAndr7ooKFoihKEEgp2elwE+Z0YJCSmJiYkOvU9qeChaIoip+dO3dy0kknMXjwYIYOHcr06dMDcpwqt0aTy0WY20l4eDgWiyUgx+kq6j4LRVEUPyaTiaeffpqjjjqK+vp6RowYwSmnnMKQIUO67BguXafI6SLWZUcIQVRUVJftO1BUzUJRFMVPamoqRx11FABRUVEMHjyYXbt2ddn+vc1PLixuN0aPh+joaIxGY5ftP1BUzUJRlJB0/6e/sbaorkv3OaR3NPedNbTD62/fvp2VK1cyatSoLitDpdtDg0cjyWXHZDKFdKe2P1WzUBRFaUNDQwMXXHABzz77bJfd++DUdYqcbuLcTtD1kO/U9qdqFoqihKQDqQF0NbfbzQUXXMCll17K+eef3yX7lFKy0+7CpOuYnA7CwsKwWq1dsu/uoGoWiqIofqSUTJs2jcGDB/PnP/+5y/Zb7vLQqOnEux2IEL5Te19UsFAURfGzaNEi3n77bebPn88RRxzBEUccwRdffNGpfTo0nRKXm1ipobtcREVF9YhObX9Ba4YSQqQDbwEpgA68LKUMzIBmRVGUDjruuOOQUnbZ/qSU/O5wYQCsjiYMJhMRERFdtv/uEsw+Cw/wFynlr0KIKGCFEOJbKeXaIJZJURSlS5W5PNg1nVTpxqNpxMbH95hObX9Ba4aSUhZLKX/1va4H1gF9glUeRVGUrmbXdEpdbmIN4GlsxGazYbPZgl2sgxISfRZCiEzgSOCXNj67XgixXAixvLy8vNvLpiiKcjB0X/OTEUGk0wHQ4zq1/QU9WAghIoEPgduklHvdgSOlfFlKOVJKOTIpKan7C6goinIQylxuHJpOqkHidDqIiorCZOq5dysENVgIIcx4A8W7UsqPglkWRVGUrtKkaZS6PMSZjHgaGzAajT2yU9tf0IKF8PbwvAqsk1I+E6xyKIqidKXm5ieTEMRqLjweDzExMSE7qVFHBbP0xwKXA+OFEKt8j9ODWB5FURQcDgfHHHMMeXl5DB06lPvuu++Ati9xunFqkjSzkcaGBqxWa4/t1PYXtAY0KeX/gJ43fkxRlEOa1Wpl/vz5REZG4na7Oe644zjttNMYPXp0u9s2ejTKXR7izUb0pkakb1KjQ0HPrhcpiqJ0MSEEkZGRgDdHlNvt7tB9EZqv+clsECQKsNvtREZG9uhObX+HxlkoinLo+fJuKFnTtftMGQanPdbuapqmMWLECDZv3sxNN93UoRTlJU43Ll2SFW6hvqoKo9HYEnQOBapmoSiK0orRaGTVqlUUFhaydOlSCgoK9rt+g0ejwuUhwWJCOBx4fJMa9fRObX+qZqEoSmjqQA0g0GJjYznxxBP56quvyM3NbXOd5uYni0HQy2ykoroei8VySHRq+zt0wp6iKEoXKC8vp6amBvD2O3z33XcMGjRon+sXO924dUlfm4XG+vqWTu2emP9pf1TNQlEUxU9xcTFXXnklmqah6zpTp07lzDPPbHPdeo9GpctDksWEWdeobWoiIiICs9nczaUOPBUsFEVR/AwfPpyVK1e2u55Hl+x0uLAaBL2sZqoqKjAYDERFRXVDKbufaoZSFEU5CEVOF24pSQ+z4Ghqwu12H3Kd2v4OzbNSFEUJoFq3h2q3RrLFhA2oq6vDYrEQFhYW7KIFjAoWiqIoB8CjSwqdbmxGQS+LmfpDuFPbnwoWiqIoB2CX04VHekc/edxumg7hTm1/KlgoiqJ0UI3bQ41bo5fFjM1goLa29pDu1PangoWiKEoHuHWdQoeLMKOBZIsJu92O2+0mKirqkO3U9nfon6GiKMpB0DSNI488kjPPPBMpJYUONzrQ12YBKamrq8NsNhMeHh7sonYLFSwURVHaMH36dAYPHgxAtUejzqORYjFjMxqor69H1/VDvlPbnwoWiqIorRQWFvL5559z7bXXoktJkcNFuNFAksWE2+2msbGR8PBwLBZLsIvabdQd3IqihKTHlz7O+qr1XbrPQfGDuOuYu9pd77bbbuOJJ56grq4Ohy53Nz8BtbW1CCEOi05tf6pmoSiK4uezzz4jOTmZESNGUOfR8UhJqtWM1WjA4XDgcrmIjo7GaDQGu6jdStUsFEUJSR2pAQTCokWL+OSTT/jiiy9osNtprK/n9muu5q233qK2thaTyXTYdGr7UzULRVEUP48++ig7d+7k+7XreeL1NznppJN45513aGhoOOw6tf2pYKEoitJKpdtDg0cnwWzCIAQej4eGhgbCwsKwWq3BLl5QqGChKIrix6npFDndRJkMnH3yBD777LOWTu3o6OhgFy9oVLBQFEXxkb4pUgWQZrMghMDhcOB0OomKijrsOrX9qWChKIriU+7y0KTp9LFZsBgM6Lre0qkdERER7OIFlQoWiqIogEPTKXG5iTYZiTN5axCNjY1omnbYdmr7U8FCUZTDXnPzkwFIs5kRvk7t+vp6bDbbYdup7U8FC0VRDntlLg92TSfNZsHsyyBbV1d32Hdq+1PBQlGUw5pd0yl1uYkxG4k1e+9TdjgcOBwOIiMjMZnUvcsQ5Du4hRCvAWcCZVLK3GCWRVGUw48uJb87nBgRpFm9uZ+klAwYMICIiAisVismk4nly5cHuaTBF+yQ+QbwAvBWkMuhKEo3cTs1yn+vJyYpjPAYS1A7jktdbhyaJDPMgsngLUdjYyNSSr799lvS0tKCVrZQE9RgIaVcKITIDGYZFEUJPKlLijbVsH5JMVt+Lcft1AAIizKTlB5FYnokielR6DYdKWW3BJAmTaPM5SHObCTG1/ykaRr19fUIIbDZbAEvQ08S7JpFu4QQ1wPXA/Tt2zfIpVEU5UDUlDWxYUkJG34pob7SgdlmJHtkMpnDEqmvclBR2EDFznpWfbcTXZMcfVkcFTsbMFkMNLz4DO4tGxEGEELQFeHDOngQKX//u7f5ye7CLAS9rbvnpKirq0NKidFoZOLEiQghuOGGG7j++uu74Og9W8gHCynly8DLACNHjpRBLo6iKO1w2j1sXl7KhiUlFG+pRQhIGxzP6HOy6HdEEmbL3ndBa26dquJGiqt2YIs043FpeDw6mkdvWUcIX9AwCO9rw8EHkBKnG6cu6RdubWl+cjqd2O12IiMjWbRoEb1796asrIxTTjmFQYMGMW7cuIM82qEh5IOFoiihT9d0dq6vZsPiYrbmV6C5deJSwhlzXn8GHJNCZNz+71Mwmg0k9Y2iotFIVLy3+Sf2kXvRPDoel+4NHi4dj1tD13b/ZjSaDJgsBkwWIyaz99lo2v8gz0aPRrnLQ7zZRLTv5jspJbW1tRiNRiIjI1uGyyYnJ3PeeeexdOlSFSyCXQBFUXquyl0NrF9SwsalJTTVurBGmBgyNpWBY1JJzojqVN+DEAKT2YjJbIQIc8tybwDR8Lh3BxFnk6flc4PRt11zELEYMJoMCCHQfDffmQ2C3rbd+2xsbMTj8RAXF4fdbkfXdaKiomhsbOSbb77h3nvvPejzOFQEe+js+8CJQKIQohC4T0r5ajDLpCjK/tkbXGxc6m1mKv+9HoNBkDEsgUGjU8nITcBoDuztW0aT9+LvX1fRdbm79uF7bqpztXwuhMBkMeAwgFFIekVaMUhA7O7Utlgs2Gw2tm3bxnnnnQeAx+PhkksuYdKkSQE9p54g2KOhLg7m8RVF6RjNo7NjTSXrlxSzY00lui5J6hvFcVNzGHB0L8KiLO3vBJCaRuPPP9O4ZAlhQ4cSPno0pvj4TpfPYBBYbCYsfgOYpC69tQ+3r/bh1BAOjVgJziY75YDJbMRjsCOlJCIsCikhKyuL/Pz8TpfpUKOaoRRFaZOUkvLf61m/uIRNy0pxNLoJj7YwfEI6g0ankNAnssP7cpeWUvPhh9TMnYunqNjbWy29fQ/WIYOJHDuWiLFjkbGxXVZ+YRCYrUbMViOalOxsdCAw0N9iQfc1YbmcLjy6C6GZqa9wUo9zz34QX19Ie/0ghwMVLBRF2UNDtZONS0tYv6SE6uJGjCYD/Y5IZNDoVNIHx2EwduzCKTWNhoULqZk9h4YFC0DXiRg7hl533knkiSfi3LCBxsWLaVz0M5VvvkXlzFfxvPgizvBwjBERGCIjETZbl9xzUeR049Yl2eFWLCYjWIzIcBP2inoMGEhMikdz6wfVD3K4UMFCURTcLo1tq8pZv6SEwnVVSAkpWTGceOlAskckYw03t7+T5n0VFVHz4UfUfPghnpISjImJJFx7LbGTL8Did69UWF4eYXl5JP7hD+iNjTQtX852qxXcbtylpVBaijCaMER6A4chIgKDpWPNXf7qPBpVLg9JFhMRpt3DdpuamnC73cTGxrZ0pB9MP0hz8GgOJodqAFHBQlEOU1JKijfXeu+qXlGGy6ERGW9lxGmZDByVQmyv8I7vy+OhYcECqmfPpnHhTwBEHHccvf7xd6JOPBFh3n+wMUREEHnCCRjXrcOak4PudqM3NqI3NKA3NKDV1gIgLBaMkZEtwUO0M3OdR5cUOlxYDYIUq9+IKk2jrq4Oi8VCWFhY22XaXz+I31BeR4MbKXcP5/UfxtscTAyGzgUQu8uDzWwMaiBSwUJRDjN1FXbWLylhw5Ji6iocmKxGso9KYtDoVHrnxCIO4MLmKiykZu5caj/8CE95OabkZBL/+Adizr8AS1qfgy6jwWzGEBsLsbFIKZFOpzdoNDbiqamBqipAYAgP89Y4IiMxhIUhDHs2kRU5XbilJCfcisHvQltfX4+U8oAnNfLvB2kmpdzrfhCX3YOj0d2yjrcfxLjHPSEd6QdxazrFtQ5qmlxkJIQTE3bgNauuooKFohwGXHYPm38tY8OSEoo21YCAtIFxHHNmP7KOTN7j4tce6XZTP/8HambPpvHnn0EIIseNI3bqFCLHjUN0cUpvIQTCZsNgs2FKTETqOrrd3lLr8JSXQ3k5wmDYHTgiIqg3mqh2ayRbTYT71UBcLhdNTU1ERERgbqfG09Hytb4fREqJrsk9aiAel4azaXcA2V8/iC4llQ0uSuscSCA52kaUtfNl7QwVLBTlEKXrkl3rq1m/pJitK8vxuHVie4Uz6pwsBo5KablTuqNcO3ZQM3cuNf+dh1ZRgSklhcSbbiL2gvMxp6YG6Cz2JgwGjBERGCMioFcvpMeD3tiI5mu20urrATAYTfQJCyM2NgZpNCBMppY7tQ0GA1FRUfs8Rk1NDddeey0FBQUIIXjttdcYM2ZMx8soBEaT8N4P4teap2t7dqJ7XBpNdZ49thMmgV3XcUhJlNVISlwYVnPHg3mgqGChKIeYquJGNiwpZsMvpTTWOLGGmxg4JpVBo1Po1S/6gJpddJeLhu++o3rOHJoWLwGjkcgTTyRu6hQijjuu3T6D7iBMJowxMRhjYgBvmcuraxCNjUTZm3A31OMGDDYbrogI3FISEx2NwbDvZqBbb72VSZMmMXfu3JaaSFcwGA1YjAYstt2X3uZ+EKfDQ12jG92tYwGsCLDr1DoaMZmNRMZZ99iuu6lgoSiHAEeDm03LS1m/uJiyHfUIg6Dv0HiOm5JD5vAEbxPJAXBu20bNnLnU/ve/aNXVmHv3Jum2W4k573zMvZIDdBZdo04YKAmLJCU2DpvFhLQ70Bob0BoaaNA0jLqOYccOnOHhGCIjMUZEIsJ2D9Gtq6tj4cKFvPHGGwBYLBYsBzEKq6OkgBqXm9IGp7fJKc5KYqQVqe3ZD9LZTvLOUsFCUXooTdP5vaCS9UtK2L66Al2TJPSJ5NjJ2Qw4JoXw6AO7wOlOJ/XffEvN7Nk0LVsGJhNR48cTO2UKEceO3avzONB+mr2Rip0NB7SNRNKk6RiEYGur8no0D1FJZk48rx9GswW9sQFPaSkeShFGI4aISAyREWzesoWkpCSuvvpq8vPzGTFiBNOnTyciIqIrTw+ABoebXTUOnB6NaJuZ1Fgb1ubhvYa982IFkwoWitKDSCmp2NnA+iXFbFpWir3eTViUmWEnpjFoTAqJaftuh98X5+bN1MyZQ+28j9FqazGnp5P05z8Te965mJKSAnAWgePUfXeFt/oV7u1w1jCZwrDFxkKsb7nb3dLXoTc0otXV0rR1K7/++ivP3HsvY6ZP5/Z77uGxxx7jwQcf7LJyuj06xbV2auxuLCYDmQkRRIeFRlDYFxUsFKUHaKx1+pL3FVO5qxGDSdBvuO+u6qHxGDt4V3Uz3eGg7quvqJk9B/uvv4LZTNTJE4ibOpXwUaO6vRbRluOnDjig9avcHnbaXaRazST7jRySUlJZWYnb7SY5ec8mNGE2Y2o1RDdTGEhLSWFkVhaunTs5e8QInn7jDdylpfscottRupRUNDgpq3MC0CvaRlKkNehNTB2hgoWihCiPW2NbfgXrF5ewc20lUkKvftGccPEAskf2wnYQzROODRupmT2b2k8/Ra+rw5KRQfIddxBz3rldktAvWFy6zi6Hi3CjgSTLnpc1h8OBy+UiJiYG43465JuH6KYNHUJ6v35sF4Kcfv1Y8MYbDM7OxlNeAeXl4Bui25KSxGrt0KCB/TY59QAqWChKCJFSUrK1jvVLitm8vAyX3UNknJUjT81g0OgU4lIOvN1cb2qi7ssvvbWI/HyE2UzUqacSO2UK4ccc3ePTU0jpvUtbAn1tlj3OR9d1amtrMZlMhId3/I70559/nssuvxyXy0VWVhavv/46tujolrvKtYYG3L4husJk8tY4mlOStLp3w+VrcqrtQU1ObVHBQlFCQF2lnY2/eJP31ZbZMVkM9D8ymYGjU+gzMO6gmikc69ZRPXs2dZ9+ht7QgKV/f5LvvouYc87BFBcXgLMIjiq3Rr1Hp7fNjLVVc1xDQwO6rhMXF3dAQfGII45g+fLley03RkdjjI7GjHeIrt7Q4A0g9fVoNTUAGKxWb40jIpJqaaK0wZtPqic1ObVFBQtFCRKXw8PWleWsX1LMrg01APQZEMuISZn0PyrpoMbUaw2N1H3xOTWz5+AoKEBYrURPOpXYqVMJO+qoHl+LaM2p6xQ5XUSYDCSa9/y+PB4PDQ0NhIWFYbXuf1rXg2GwWDDEx0N8vLe/w+HYnZKkqgoqK4lAkG6xYo2OwmIy05O/fhUsFKUbSV2ya2M165eUsGVlOR6nRnRSGMec1Y+Bo1KITmw7qd1+9ykljoLfqJk9m7rPP0dvasKak0Ovf/yDmLPParlZ7VAjpWSnw/urPb1V8xNAbW0tQoiW+bQDSQiBCAvDY7ZSYginzhhDtHSTaPBgczShV5TjrCj3DdH1y6IbgCAWKCpYKEo3qCltYv3iYjb8UkJDtROLzciAo3sxaHQKKf0PLJldM62+nrrPPqN6zhyca9chbDaiTz+duKlTsOXlHXK1iNYq3B4aPTppNgvWVqOTHA4HTqeT6Ojo/XZqd5XWo5ySY8JIioxtaXKSHg9ac5NVQwNaXR0AwmzBEBnhzaQbEdHlebW6UuiWTFF6OEejm80ryli/uJjSbXUIAelD4hl7fjb98hIxWQ78IialxLF6tbcv4osvkXY71sGDSbnvXqLPPBPjfvIdHUqcmk6x002UyUB8q7vT/Tu1A3EjXWv1DjdFfqOcesfavBMs+REm055DdJv7Oxoa0Wtr0aqrATCEhbXcHGgIDw+JIczNVLBQlC6kazq/r61i/WLvXdWaRye+dwRjzu/PwGNSiIg9uGYHra6O2k8+pWb2bJwbNyLCw4k580xip07Fljv0kK9F+JNS8rvDhQFIa6P5qbGxEU3TSEhICOj34j/KyWoy0C8xgihb+6OchBAIq9XbBJWQ4L1hsMmO3ujLoltZARXlIAwYIsJb5u/o6BDdQFHBQlE6SeqSyqIG1i8pYePSUux1LmyRZoYe35tBY1JJTI88qD9yKSX2lSupmTWbuq++Qjqd2HJzSbn/fqLPOANjZOB/NYeicpeHJk0nPcyCpXVKD4+H+vp6bDZbQDq1wdfkVO+krN7b5JQSbSMxas/5Mg6EEAJjRDjGiHBITkZqmq+5qhGtsQF3SQkAlvT0oPY/qWChKO2QuqSp3kV9pYO6Srvv2UG930Pz6BiMgsxhiQwcnUJGbkKHJrdpi6e6mrpPPqF6zhxcm7dgiIgg5rxziZ0yhbChQ7v47HoWh6ZT4nITbTIS18YNbXV1dZ3u1N6wYQMXXnhhy/utW7fywAMPcNttt+3R5BQTZiY1Zu8mp84SRuOeQ3TdbvSGBgzd0KS2PypYKIc9qUua6ly+AGD3PlftHQz8hUWZiYq3kdAnkn7DE4lNCadfXiJhkQeXnVRKSdOyZdTMmUv9118jXS5secNJffghoidNCvqFIhTs2fxk3qu25nA4cDgcREVFYepER/HAgQNZtWoV4J1+tU+fPpxx1jnsqGz0NTkZO9zk1BUMZjOGELgvRgUL5ZDXZjBoFRh0j9xjm7AoM1EJYSSmRdIvL5HoBBtRCWFEJdiIircd0Mxy++OpqqJ23sfUzJmDa9s2DFFRxE6ZQuzUKdgGDuySYxwqylwe7JpORpgFc6vmJykldXV1GI1GIiMju+yY3373HX0z++G0xeN0eDrd5NSTqWCh9HidCQZJ6VFkHZEUsGDQdnl1mpYu9d4X8e134HYTdtRRpF5/PdGTTsUQduD3WhyKfnjjZcp2bAVAl2DXdIxCsN6494Va0zQ0j4bJbNrvpEbJGVmcdNX1HTp+vcPNK2+8w4QzzyPKZiI1JgzLQTYtHgpUsFBCXpcHgwQb5oMYttpZnooKav77X2rmzMX9++8YYmKIu/gi4qZMwZqT0+3l6TEkOHQdxN6pxwGkBM2jeefg7oKhpi6PRlGNg8q6RuZ//QVPPPYoGQmqGVAFCyXopC5prHW1HQgqHdRX94xg0Bap6zT+vJia2bOpnz8fPB7CR44k6eY/ETVxYo+6g7e7NdcAip0uypweMsMsxJj3vmRVV1djt9tJTk7uVF+FrkvKG5yU+0Y5FSxdwMiRI+ifkXbQ+zyUdOibFUJMk1K+6vfeCPxTSnl/wEqmHDI6FQz6RpF1ZOgGg31xl5ZR+9+PqJn7Ie7CQoyxscRffjmxUyZjzcoKdvF6jCZNo8zlIc5sbDNQOJ1O7HY7kZGRnQoUdQ43RTV2XB7dN8opjIfmfcglF1/cmeIfUjr67U4QQlwATAMSgNeBBZ09uBBiEjAdMAIzpZSPdXafSvc7qGAQbSE6wUZSRhT9j0pqCQTRCTYi40M/GLRFahqNixZRPXs2DT/8CJpG+OjRJP/5diJPPhlDAOdxPhTpUvK73YVZCHpb9/7upJTU1tZ2qlO7ucmpzrHnKKempia+/fZbXnrppc6exl6klCABXXpf64CUSH33cu97v9cSjNEWDEH8u+hQsJBSXiKEuBBYAzQBF0spF3XmwL7ayYvAKUAhsEwI8YmUcm1n9qt0PV2XNNU62w4EzX0GWvcGA+8fmQ6a5n2taUhdB11Hapr3M133LtM07x+irrX/WfO+NB2k7n3WW++7ef3mdTTcOwup+e9HeIqKMSYkkHDN1cROnowlI6NT53k4K3G6ceqSfuFWTG30VTQ2NuLxeIiLizvgvorWTU4pMTYSI3ePcgoPD6eysnKPbQ7oIi+lbz3/5bu36RAhwOC9aQ+D8B4jiDraDJUD3Ap8CAwGLhdCrJRSNnXi2McAm6WUW33H+AA4B+jyYPHq7fejW0JzhIkgNIfgmaSRFFcvImQ8urAhES3/VyUgdCexup0E3Y5Bb8Sk2TFqjZh8D2OdG7FTIqT3giqkxK1LKqWkUuotf0D4Pt/9R6W3PEvfH5x3PV9waHnG94tL95ZI+j7b67VvHe8B9l7H936vddraVzsM6RqmU3UM/YpoML5Iwycv7rmV70Lkv0zu57OW962Gae6x/b4+81u+v30e0PH2V/4O7tP/M81gQBei5VkXgphBT1JXCQ5hodKYTJRsQq/7nZpWfyZSCuyOSAwGDad9B04H++YtAEIaAIFLs9LojkBKIwkGNxFGJ8amBuxNBoQUIA3ev0vfc/N23uf2SSQIHSm8/5+kkEh0EBJp8v5/lWLv59bLWl8aZL2JOGt2h8oQCB1thvoUuElK+b3w3glzO7AM6MztpH2AnX7vC4FRrVcSQlwPXA/Qt2/fgzqQNIVRYXMd1LaBFOQfCvvlxsP28DoStEj6e3rR35OEDQtC+AKcsAI2hEgIdlG7T3NgQvc9ZMt7IXRAQwgNSj0INMD3LJrfa3sux+Nd3/d+92ce37q+5cLT9rZ77dtvO+G/bqt9iz3XBX33+yD/dlnX302000E0DpJkHd4rpgGJEYkBhAGJgTq8AwNidDA6IkEagOaH2P26ZfluYUBL0gzd6n3gd5FvvrALHYlEN7R9ce/oRf5Q0dFgcQxwrRDiJrx/Hf8DLurksdv6Sve6fkopXwZeBhg5cuRBXV+vffLOg9nssCWl5B+vfIFrzQ8YwhtZGtbAMttW+mZlkTlkMLEZfXDgodHjosntoNHjxOFx4XC7cLidODQnTo8bt9uFU3Ph8rjxaC7cmgdNc6PpbjyaB01zoesaUveg6R6k7gHd42si8iCkhtA1kBoGqWPQdYTUMEoNAxKD75efwferz4DAgPfXoe83IQbfcu/r1st86woTRswYhAkDJozCiAETJkwYhQmj79mEEaMwYhImTMKIESNmvK/DjOH0j84ixmDCLEFqvmYITYKme5sgNNnyrOu+5X7LWp416Wvyaq5ZdeM/vgGEQYBReJ9b3rP7vdFvuf+z/zoGWq3DHs8gkU4P0qmhOzR0p450STQRi0P2BwmGti4REtxoOAxNhEkLRmn1LhS6r+7rDeTNwY/mIN4S4PV9v259OP+qNMJbeOF7lr7n5sAkBEizrwblt17r7Vp/fiDrBVlHg8VMoB543vf+YmAMMLUTxy4E0v3epwFFndif0kWEEPzfNadx2avxlG0oYIpWQK3DQaGusWPLFmw2G0OHDiUvL4/09KxuyYQppcQjwaXrOKXErnlo8rio97hocttp9Lho9Dhp8riwe5w0aU4cmguHx4ldc+LUXL6H97Vbc+PSnbh0N27NhUdvfjSi6W7fw4Wuu9GlG6m7kNKDlG6QHoR0+Z49CNwAuPV+1MXdRFx4H/qFWcgMs9IvzEpWuNX32kJsGyN62j13fR+BpznA7OMz9OYApLcKXtLXX+P9fPeyvT/bY1+a33Lfds1lkG4dHK3L47evPcrjvQoLswFhM2GwWRFWI4YYE7pRp85kwGo0EGE2eoONAAzC13YPtbVVGHQDMYlxCKNhv///6uzeUU4eTSMmzERKlAWzgd1NoHs0Ser7WNa66dK/ybK5aVQD6WF3M6a+5zE6TUB8FtgCP5HTPksgO9DZIoTIl1LmtbfsgA4shAnYCEwAduFt1rpESvnbvrYZOXKkbGteXCUwqhpdnPPi/3C4NKYfIyn45AOqmxxYM3NoMJhbOhfz8vIYPnw48fHxwS5yt9CkxKlLXLqOS5fYNY35O+czY/mDaFJncOafcESMYZvdSZHTvce2cSYj/cK9QSQzzEI/X0DpF24lzmQ8LFKNN3cUi1ad1vUejWUFv5GancOACBvGNr6LxsZGamtriY2NJTw8fJ/HcHo0iv1GOfWJtRHZTbmc9uLXf7ZXENlfAGrdjxYeDyZbu4dbt24dgwcP3mOZEGKFlHJkZ06jo8HiDeA/UsolvvejgCullDd26uBCnA48i3fo7GtSyof3t74KFt1vY2k958/4mb7x4XwwbQQbfviKX/47C6fDSfLRY3FHxfN7YSHg7VPKy8tj6NCh2Gzt/6c+1BQ1FHHXwrtYVb6Kc7PP5W/H/A1hsLHD4WR7k4utdifb7U62+R67HO49fnPGmIx7BZDmoJJoNh3ygeQv639nkr2ascNyiWgjk6umaZSVlWE2m/c5V0XzKKeyeicC6BVtJSHy8MrlFOxgsQ4YCPzuW9QXWEfzQDAph3emEB2lgkVw/LC+jGlvLmPikBRmXHoUzqYGfvloFqu+/gxhMDJ04hkY+2Ty29q1VFRUYDKZGDhwIHl5efTv379bprUMFR7dw4xVM5i5ZiaZMZk8Oe5JBsa3nRDQqevssLv8AoiLbU3e14UOF/55bqOMBm/g8AsgWb6gkmTp+YHk+8o6Ll29lf8mmBgzPLfNdWpqamhqaiIpKQmzee9aQp3dTVGt98a62DALqTE2zAeZy+lf//oXM2fORAjBsGHDeP3113vMD6BgB4uM/X0updzRmUJ0lAoWwTPzp6089Pk6bh6fzV8mei9+tWUl/O+Dt1m/aAFh0TGMueBiEocMZ01BAWvWrMFutxMREcGwYcPIy8sjJSWlx1/UOuqX4l+4+6e7qXPWccfRd3DhwAsP6Nxdus5Oh2uPALLNVzP53eFqbvYHINxooJ9/jSTM20eSFW6lVw8IJDVuDycu3UC0ych/onSGDBmy1zoul4uKigoiIiKIaTUBkH+Tk81kpHcnm5x27drFcccdx9q1awkLC2Pq1KmcfvrpXHXVVQe9z+4UqGDR0ZvyuiUYKKFr2nH92Fhaz/PzN5OdHMk5R/QhJjmFM265gxGnn8OCd19j/uv/IS61D8dfciUT//IXNm/eTH5+PkuXLmXJkiUkJyeTl5fHsGHDOjU5TU8wKnUUc8+ayz8X/ZOHf3mYJcVLuH/s/cRYY9rfGLAYDPQPt9E/3ObNmeDHrUsKHa6WALLN7mRbk4t1DQ6+qqjF/2b5MINhd9NWq76SVKs5JJpn/rlpF+VuN28O74co3PtS03yntsFgIMpvjnFdl5TVOylv8DY5pcaEkRBp6ZJz8ng82O12zGbv3dy9e/fu9D57ug7VLEKFqlkEl8ujc9nMX1hVWMPsG8ZwRHpsy2dSSrb+upSF775B1a6d9Bk0hHGXXkPvAYNoamrit99+Iz8/n8LCQoQQZGVlkZeXx6BBg7AcwmkwdKnz9tq3efbXZ0kMS+SJcU9wZPKRATueR5fscrpamrW2+9VKdthduPz+3m0GQYZvlFbLyC1fU1dvq7nNDuau9mV5DVcXbOf2jF7clZW6x6/imk+34CpqRNc1PB4PJpMJg8EISDRd4vLo6BJMRoHFZOjQDa6W3hHEntW/3fWmT5/OP/7xD8LCwpg4cSLvvvtuZ0+12wS1GSpUqGARfJUNTs55cRFOj84nfzqW1Jg974zXNY2CH75l0ex3aKqtYcCoYznukiuJS/H+MquoqGD16tXk5+dTW1uLxWJhyJAh5OXlkZGR0SUppkNRQUUBdyy4g+LGYm484kam5U7DaOjevhxNSoqc7pYAsrvD3cUOuxOHvvtaYBGCDL+mLW9fifd9H6ulzfQbB6rS5eGEpetJsZr5YkQOFoOhjWDRgMvlRggwm81ICU6PjqZLDAawGI0YD6AsHQkW1dXVXHDBBcyaNYvY2FimTJnC5MmTueyyyzp1vt1FBQtUsAgVG0rqueDfP5OREM6cP4wh3LJ3a6bLYWf5p/9l+acfoXk85E08jdHnX0R4tLcZRtd1fv/9d/Lz8/ntt99wuVzExMQwfPhw8vLySExM7O7TCrgGVwMPLHmAL7d9yaiUUTxy/CMkhycHu1iAN2lfidPtCyCulv6RrU3eZ7tfIDELQV+brzYSvmdfSZrNgrmDF+/rCrbzVUUt34wcwOBI74+O1he62tpaGhsbSUhIpNqhU97gxAAkR9u6rMmptTlz5vDVV1/x6qveRNtvvfUWS5YsYcaMGV1+rEBQwQIVLELJ/PWlTHtzOZOGpvDiJUdh2McFoqG6isVz3mPN/G8w22yMOm8qR552FmbL7nkcXC4XGzZsID8/ny1btiClpE+fPuTl5ZGbm7vf8fQ9jZSSeZvn8ejSR7EZbTx83MMcn3Z8sIu1X1JKSl0eX9+If2e7dzhwk7Z73JZRQLqtdWe7haxwK+k2CxZfzXFeaTV/WLuDv/VL5dbMXi3b+1/o3G435eXlmCw2qtwmXJpObLhvlJMxcDXQX375hWuuuYZly5YRFhbGVVddxciRI7n55psDdsyupIIFKliEmlcWbuXhL9Zxy/hs/jxx//NFVxb+zsL33mDriqVEJSRx7IWXMeT4kxCtmp3q6+tZs2YN+fn5lJaWYjAYGDBgAHl5eeTk5HRqzoJQsrVmK3csvION1Ru5csiV3HrUrZiNQbpprBOklJQ3BxK/ALK9ydvM1eAXSAxAmi+Q5Nc30S/MyqdH5ezRpNV8oZNSUl5Ridvtplq3YTGb6B0bRqS1e/7977vvPmbNmoXJZOLII49k5syZWHvIRFUqWKCCRaiRUnLn3NXMWVHIcxcfydl57Y8Y2fnbaha88zqlWzeRlNGPcZddQ+bwtjt8S0pKyM/PZ/Xq1TQ2NhIWFkZubi55eXn06dMn5IeEtsfhcfDU8qeYtWEWQxOG8uS4J0mPTm9/wx5CSkmlW/M2Z/lqJc19JHUejTeH9SMnYs97F9atW8fAgYMora5HOhtolBZio6MC1uR0KFLBAhUsQpHTo3HZzF9YXVjLrFYjpPZF6jrrF//E/95/i7ryUjLzjmLcpVeTlNGvzfU1TWPr1q3k5+ezfv16PB4PCQkJLWlGYmPbP2Yo+37H99zz8z3oUufe0fdyetbpwS5SUEgpWbn6N8KT+xKhN2IwGEhOSsTcxt3cyr6pYIEKFqGqeYSUy6PzcRsjpPbF43az6uvP+OWjWTiaGhl6wgSOnXoZUQn77tx2OBysXbuW/Px8duzwjsnPzMwkLy+PwYMH95i7bFvzTxVyXvZ53H3M3YSbD52+mvZsLW/g/z5dy+WDTOT064tFukhISOgxTT+hRAULVLAIZRtK6jl/xiL6JUUw+4a2R0jti6OhgV/mzWbll58gDEZGnHEOR589GWs7HdvV1dUtw3CrqqowmUwMHjyYvLw8srKyetww3ANJFXKoaHJ5ePGHzbyycBtWk4FXz0khLSWJsLAw4uLigl28HkkFC1SwCHXfryvl2reWc1puCi9cvO8RUvtSW1bK/z54a3f6kMkXM3zCJIztdGpLKSksLCQ/P5+CggIcDgeRkZEtw3B79eq13+1DTWdThfQEUkq+Kijhwc/WUlTr4Pyj+nDXpIH8vqGAtLQ0kpOTD6ucYl1JBQtUsOgJXl64hUe+WM8tE3L48ykDDmofJVs2sfCd19i5dg1xqb05/uKryD5mTIcumB6Ph40bN5Kfn8+mTZvQdZ2UlJSWNCORkZEHVabuVmmv5J+L/sn/dv2PCX0nHFCqkFBXVu/gL7Pz+WlTBYNSonjw3FyOzoxnw4YNVFVV9ah/p1CkggUqWPQE/iOknr/4SM7qwAipfe1n28rlLHz3dSoLf6f3gMGMu+wa+gwc3P7GPo2NjRQUFJCfn09RURFCCLKzs8nLy2PgwIFtZi4NJS2pQlY8S2J44FOFdIfC6iYum/kLZfVO7jx1IJeNzsBkNOB2u3nxxRcZO3YsRx999CFXk+pOKliggkVP4T9CavYNY8jrwAipfdE1jYIfv+PnOe/SWF1FzqixHH/xlcSl9jmg/ZSVlbF69WpWr15NXV0dVqu1Zba/vn37hvTFKRRShXSFbRWNXPrKEhqcHt645hiO6ru7T2LBggX88MMPXHDBBQwbNiyIpfSaPn06r7zyClJKrrvuOm677bZgF6nDVLBABYuepKLByTkvLMKt6Xzyp+NIiencKCW3w8Hyz/7Lsk8+RPO4yTvldEZfsDt9SEfpus727dvJz89n7dq1uN1uYmNjW4bhJiQktL+TIGhwNfDA4gf4cnvopQrpiPUldVw2cylSSt6eNoohvXdnHa6urubFF19kwIABDBs2bK8LXXcrKCjgoosuYunSpVgsFiZNmsS///1vcnJyglqujlLBAhUsepr1JXVcMONnspIimX3DGMIsnf813FhTzeK577H6+68xW20cc85kjjrjnD3Sh3SUy+Vi3bp15Ofns3XrVgDS09NbZvsLC+vYEODu0hNThQCsLqzhiteWYjMZeefaUWQn79kfMWvWLDZv3sxNN91EcXFx0IPFnDlz+Prrr5k5cyYADz74IFarlTvvvDOo5eooFSxQwaInah4hdXpuKs9ffOQBj5Dal8rCnfz0/htsWf4LkQmJHDv1MoaMO8mXwvrA1dbWtqQZKS8vx2g0tsz2l52dHVIjc3pSqpCl26q45o1lxEWYee/a0aTH7zkcevPmzbzzzjuMHz+ecePG7XGh+/LLLykpKenS8qSkpHDaaaftd51169ZxzjnnsHjxYsLCwpgwYQIjR47k+eef79KyBEpQJz9SlIM1YXAv7p40iEe/XE92ciS3H+QIqdYS0tI59457KFxbwIJ3XuXrfz/Lr198zLhLryYz76gD3l9MTAzHHXccxx57LMXFxeTn57NmzRrWrl1LeHh4y2x/qampQe/fyIrN4t3T3+Wp5U/x5to3WV66PCRThSzcWM71by+nT2wY7147eq+mSI/Hw5dffkl8fDxjx44NUin3NnjwYO666y5OOeUUIiMjycvLO2RyknWGqlkoASel5K9zVvPhr4W8cMmRnDm8a2cdk7rOhiX/43/vv0ltWSkZw49k3KVXk5yZ1an9aprWMtvfhg0b0DSNpKSklmG4raf3DIZQTRXy9W8l3PzeSrKTI3lr2jEkRu7dTLho0SK+/fZbLrnkEgYM8P6IaOtXcbD9/e9/Jy0tjRtvvDHYRekQ1QyFChY9mdOjcekrv7BmVy1z/jCG4WmxXX4Mj9tN/jdfsOSjD3A0NjB03HiOvfDy/aYP6Si73d4y29/OnTsB9pjtL5hpKUItVci8lbv4y5x8hqfF8MbVxxATtncTWV1dHS+88AKZmZlccsklLctDJViUlZWRnJzM77//zsSJE1m8eHGPuaNcBQtUsOjpmkdIeXSdj2/q/AipfWlJH/LVpwgER51xDsecMxlreESX7L+ysrIlzUhNTQ1ms5ns7GwsFgtCCAwGA0KIPR6BXiaRfLPjG777/TuSIpK4Ovdq0qLS9lov0OWYtayQf35SwOh+icy8ciQR+0gp/uGHH7J27Vpuuukm4uPjW5aHSrA4/vjjqaysxGw288wzzzBhwoRgF6nDVLBABYtDwbriOib/u2tHSO1LXXkZ/5v1Nut++oGwqGhGX3AxeadMwmjqms5gKWXLbH/btm1D13WklEgp93i9v2U96e/vQO0vqDgcDsaNG8f48eP32CZUgkVPpoIFKlgcKr5bW8p1by/n9GGpvHDxkQHvMC7dupmF777G7wWriU1J5fiLryRn1LFB76gG2gweHQ00+1pW46hhxqoZ5Jflc3Svo5mWO41wU/h+t+3sMRdtquB/m8sZ1CuSSbkpCPZ/XuHh4Rx77LF73UWvgkXnqWCBChaHkv8s2MJjX67ntpNzuO3krhkhtT9SSravWsHCd1+nYucOUgcM4oTLph1Q+pCexD9VSFJ4Eo+PezwgqUKklDz65XpeXriVySPSePyC4Rg7MTxaBYvOC1Sw6Fk5nJVDxg3jsjj/qD48+90mPl9dHPDjCSHod+RILn/iOSbecAt15WV8cO8dfPL0I1QV7Qr48bubQRi4cuiVvH362xiFkau/uppXVr+Cpmtddgxdl/xzXgEvL9zKlWMyeKKTgUIJbSpYKEEhhODR84cxIiOOv8xZxZrC2m45rsFgZNj4iUx79mXGTr2U7atX8uZfb+T71/5NU21Nt5ShO+Um5jLnrDlMzJjIcyuf44Zvb6CsqazT+/VoOn+dk8+7v/zOH0/sz/+dPbTLbrhU9qTrOvVVlei63v7KAaSChRI0VpORly4fQUKElWvfWkZpnaPbjm222RhzwcVMm/4yw8afSv63X/Lqrdex5KNZuJ3dV47uEGmJ5PFxj/PA2AfIL89n8ieT+anwp4Pen9Oj8af3VvLRyl3ccepA7po0KCT6fw5FmsdDdfEuGqurcNmbglqWoAQLIcQUIcRvQghdCNGpdjSlZ0uMtDLzypHUOzxc99Zy7K6uaybpiIjYOE6+9kaufOpF0ofmsWjW27x22w0U/PAtehc22QSbEILzcs5j1pmzSAxP5Mbvb+SpZU/h1twHtB+7S+P6t1bw1W8l3HvmEG46KTtAJVY8LidVRYV4nE5ie6ViiwjuHB/BqlkUAOcDC4N0fCWEDE6NZvpFR7JmVy13zM0PynDShD7pnHvHP7nw/x4jKj6Rr/8znbfvupVtq1YcUsNbs2KzeO/097hw4IW8ufZNLv/ycnbW7ezQtg1OD1e9vpSFm8p54oLhXHNcvwCXNnhqamqYPHkygwYNYvDgwSxevJiqqipOOeUUcnJyOOWUU6iurm5Z/9FHHyU7O5uBAwfy9ddftyxfsWIFw4YNIzs7m1tuuaXl/9IzzzzDkCFDGD58OBMmTGiZTx7gzTffJDs7mwEDBvLBnDnE907DFhnJtm3bGDVqFDk5OVx44YW4XK7u+0Kg7aF73fUAfgRGdnT9ESNGSOXQNeOHzTLjrs/ks99uDGo5dF2X639eKF+5eZp8auoZcvaD/5Cl27YEtUyB8N327+SY98bIUe+Okp9v+Xy/61Y3OuXZL/xP9v/b5/LjVbsCVqa1a9cGbN8H4oorrpCvvPKKlFJKp9Mpq6ur5R133CEfffRRKaWUjz76qLzzzjullFL+9ttvcvjw4dLhcMitW7fKrKws6fF4pJRSHn300fLnn3+Wuq7LSZMmyS+++EJKKeX8+fNlY2OjlFLKGTNmyKlTp0oppaysrJSZGRly7YqlctPq1bJfv36yqqpKSinllClT5Pvvvy+llPKGG26QM2bMaLPsbX2HwHLZ2et1Z3fQqYN3IFgA1wPLgeV9+/Zt88tRDg26rsvbZ62UGXd9Jj/LLwp2caTH7ZIrPp8nX7jmIvnUhWfKL154WtaWlwa7WF1qV/0uednnl8ncN3LlPf+7Rza6Gvdap6zOIU/91wKZ8/cv5De/lQS0PKEQLGpra2VmZqbUdX2P5QMGDJBFRd7/l0VFRXLAgAFSSikfeeQR+cgjj7SsN3HiRPnzzz/LoqIiOXDgwJbl7733nrz++uv3Ot6vv/4qx44dK3Vdl6++9JK87KILZVVRodQ8Hnn99dfL9957T+q6LhMSEqTb7ZZSSvnzzz/LiRMntln+QAWLgKVSFEJ8B6S08dE/pJQfd3Q/UsqXgZfBe59FFxVPCUFCCB45bxjbKxr5y5xV9I0PZ1ha8JL1GU1mjjr9HIacMIGl8+bw65efsGHxTxx1+jmMOndKl6UPCabekb15fdLrzFg1g5lrZrKqfBVPjnuSgfEDASiqsXPZzF8ornXw2lVHc1xO5/NsddTGjQ9S37CuS/cZFTmYAQPu2e86W7duJSkpiauvvpr8/HxGjBjB9OnTKS0tJTU1FYDU1FTKyryjynbt2sXo0aNbtk9LS2PXrl2YzWbS0tL2Wt7aq6++yqRJk6gtLWHHtq1kZmYSm9IbIUTLNpWVlcTGxrZkv93XvgIpYH0WUsqTpZS5bTw6HCiUw4/NbOSly0eSEGHlureWd+sIqX2WKSKScZdezTXPvsTA0cex7JMPmXnLdfz6xcdongPrIA5FJoOJW466hVcmvkK9q55LPr+ED9Z/wPaKBqb8ZzHl9U7ennZMtwaKYPJ4PPz666/88Y9/ZOXKlURERPDYY4/tc33ZRp+WEGKfy/298847LFu2jGmXXIyjsQFLWDiW8Ig91uvovgJNJWlXQk5SlHeE1AX//pnr31rOrBvGYDMHf/Kh6MRkTvvTXzjqjHNZ+M5r/PDmK6z86jOOu/hKBowOjfQhnTEqdRRzz5rLPxf9k4d/eRiD/ROEZyrvX38SuX26v4bXXg0gUNLS0khLS2PUqFEATJ48mccee4xevXpRXFxMamoqxcXFJCcnt6zfnIkYoLCwkN69e5OWlkZhYeFey5t99913PPTQQ3z07tsYDRCTnEpWTg4//vjjHtuceOKJJCYmUlNTg8fjwWQy7bWv7hCsobPnCSEKgTHA50KIr9vbRjm8DE6N5tkLj2D1rlrumLs6pEYk9erXn8n/fIjz/3Y/JouFz559jPfvvYPq4p5/J3hCWAJ/HPQIhuqz0W2/EZfzAm7z1mAXq1ulpKSQnp7Ohg0bAPj+++8ZMmQIZ599Nm+++SbgHbF0zjnnAHD22WfzwQcf4HQ62bZtG5s2beKYY44hNTWVqKgolixZgpSSt956q2WblStXcv311/P6jBdJSkz0jniKiOTUU0/lm2++obq6murqar755htOPfVUhBCcdNJJzJ07d6/jd5vOdnp050ONhjr8vPjDJplx12dy+nfBHSG1L5rmkavnfy1fuOYiOf2KyXLtTz8Eu0idsnx7pcy97ys59tHv5deblspJcyfJvDfz5Mv5L0uP5gn48UOhg1tKKVeuXClHjBghhw0bJs855xxZVVUlKyoq5Pjx42V2drYcP368rKysbFn/oYcekllZWXLAgAEtI56klHLZsmVy6NChMisrS950000tneYnnnCCTExIkEOHDJF5w4fLs846q2WbV199Vfbv31/2799fvvbaay3Lt2zZIo8++mjZv39/OXnyZOlwONose6A6uFUiQSWkSSn5y+x8Plq5i39fehSnDUsNdpHaVFdRzufPPUnRhrXknjSR8Vdfj9kamPk6AmXR5gqufXM5qTE23rl2FL1jw2hwNfDA4gf4cvuXjEoZxaPHP0pSeFLAynCoJxKUUtJQWUFjbQ3W8AhievU66Hnj90UlElQOS0IIHjl/GEf2jeX22aso2NU9OaQOVHRiEhfe9yijzruQgh+/5d2//5mKnTva3zBEfLe2lKvfWEZGQjizbhhD79gwoI1UIZ92LlXI4UzXdWpKi2msrSE8JpbYlNQuDxSBpIKFEvJsZiMvXz6S+HAL1765nLIQGCHVFoPRyHEXXc4Ff38Ae30d7/79z6z54ZuQ6m9py6f5RfzhnRUMTonig+tHkxS15xSx/qlCEsISuPH7G3l6+dMHnCrkcKZ5PFQXFeJsbCQqMYnoxKQeNyBCBQulR/COkDqaOoeb695ajsMdunmbMocfyRVPPE/vAYP45j/P8cXzTwU9Cdy+zF62k1s/WMlRGXG8c+0oYsMt+1zXP1XIG7+9wRVfXsHO+o6lCjmcuZ1OqnbtxON2E5fSm4iY2GAX6aCoYKH0GEN6R/OvC48gvzD0Rki1FhEbxwX/eIBjL7ycDT//xNt330rp1s3BLtYeXl+0jTs/XM1xOUm8efUxRNnan27WZrLxz9H/5F8n/osd9TuY8ukUvtj6RTeUtmdyNjZSVeQdPhvfOw1rRM+9kVMFC6VHOXVoCnecOpBP84t4YX5oXXxbMxiMjD7/Qqbe9wgel4v37/krK7/6NCSC3Is/bOb+T9dy6tBevHLFiAOeC/3kjJOZe9ZccmJzuOunu7h30b00uUOz9hQsjbU1VJcUYTKbie+TjtlqbX+jEKaChdLj3Hhif847sg9Pf7uRL9cEfpa9zkobnMvljz9HxvAjmf/6S3zy9CM4GhqCUhYpJU98tZ4nv97AeUf24cVLjsJqOrhO1uZUIdcNu455m+dx0ecXsaFqQxeXuOeRUlJXUU59RTnWiAjieqdhNPX8+59VsFB6nOZZ9o7sG8ufZ+eH7Agpf+HRMZx7572ccPk0tv66jLfvvoWijeu7tQy6Lrn/07XM+HELl4zqy9NT8jAZO3cJ2FeqkFCoPXVGZ1KUf/jB+zTV1hARE8uWwiLy8vIOOEV5Tk4OOTk5LTcBArzwwgtkZ2cjhKCioqL7voxmnb1Rozsf6qY8xV9pnV2OeeQ7OfqR72RprT3Yxemwok3r5St/ukY+c/HZ8pd5c6SuaQE/pkfT5V9mr5IZd30mH/rst70yqnaFiqYK+Ydv/yBz38iVt86/VdY4ag54H6FyU97BpChvbKiXyxf9JDP6pss63w17B5OivF+/frKyslJWVVXtkaL8119/ldu2bZMZGRmyvLx8n2UP1E15qmah9FjJUTZeuXIkNU1urn97RUiPkPKXmj2Qyx6bTvbI0fz03ht89Pj9NNUFrnbk8ujc8sFK5q4o5LaTc/j76YMDMmwzISyBFye8yF9H/pUFOxcw5dMprCxb2eXHCbS6ujoWLlzItGnTALBYLMTGxvLxxx9z5ZVXAnDllVcyb948AD7++GOmTJ5MY0U5aamp5OQMoGDDBoqLi6mrq2PMmDEIIbjiiitatjnppJMIDw8HYPTo0S05pL7++mtOOeUU4uPjiYuL45RTTuGrr74C4MgjjyQzM7P7vohWen5DmnJYG9o7hn9deAR/eGcFd85dzfSLjugR49dtEZGcefvdrP7uS3548xXeuvNmzrj5r6QPHd6lx3G4NW5691e+X1/GP04fzHXjsrp0/60ZhIErh17JiF4juGPBHVz91dXcdMRNXJN7DcYDvAHtnk2FFDTYu7R8uZFhPJiTtt91DjRF+Y7t2xk2aADgHfHUNyPjgFOUn3baaYA33Xl6enq72wSDqlkoPd6kXO8IqU/yi3jxh9AeIeVPCEHeKadzyUNPYwkLZ86D/+TnOe912dzfjU4P17yxjPkbynj4vNyABwp/uYm5zDlrDhMzJvLcyue44dsbKG8q77bjd8aBpChvrKnB2diAwWjaY8TTgaQoX758OXfccQew73TnoUDVLJRDwo0n9mdTaT1PfbOR7ORIJuWGZg6ptiRnZnHZo//i+1f/zeK571G4dg2n3/xXIuMTDnqftXY3V7++lPzCWp6Zmsd5R+7/13QgNKcKGdN7DI/88giTP53MQ8c+xPFpx3do+/ZqAIHSkRTlRUVFJCYkUF9ZTnrfDKobGltGPB1IivKHH36YBQsWYPUFmbS0tDZTlIcCVbNQDglCCB67YDhHpMdy+6yeMULKn8UWxmk3/ZlJN95O8ZaNvHXXLWxbteKg9lXZ4OTil5ewZlctL15yVFACRbOemCqkvRTluqbx0gvPM3H8SUTExjH1kkuYNWvWAacov+GGG/jkk09a5sUA9pmiPCR0toe8Ox9qNJTSnj1GSNX1nBFS/ip2/i7f+OtN8qmpZ8gF77wmPb55lzuiuMYuJzz9oxz4zy/kjxvKAljKA2d32+WDix+UuW/kyvM/Pl8+uPhBOWPlDPnBug/kt9u/lStKVsjVBaulR/MEZLTWgdhXivKTTjpJ9svMlMeNGS13btvWsv6BpiifMGGCTE5Olnl5eTIvL69DKcqnT58u+/TpI41Go0xNTZXTpk1rs+wqRTkqRbnSMb8V1TL534sZ6EuMFwqz7B0ot8vJj2++wurvviI1ZyBn3noX0UnJ+91mZ1UTl878hapGF69eOZJRWQffjBVI3+34jn/n/5uypjJqnDV7fPbskGdJ6ZeCEAKTMGE0GDEZTN6H8D4bDcaW1yaDCaMwdku7vtvhoLq0CKlLYnulYvWNZgo1gUpRroKFckj6qqCYP7zzK+cc0ZtnL+wZI6TasmHxT3zz0vMIg+DUP95GztFj2lxvc1kDl838Bbtb461rjiEvPbZ7C3qQ3LqbGkcNlY5KKu2V2CptZORk4NE9aLqGR3rw6J6W95K2r1ctQaVVEGkr0BjEgbe+OxobqC0twWAyEZuSitkSuqk7AhUsVAe3ckialJvKXycO4KlvNjKgVxQ3nZQd7CIdlIFjjqdXVg6fPfs4nzz1MEdOOotxl12Dybw76d/aojouf/UXhIAPrh/N4NToIJb4wJgNZpLCk1omVFpXt47EsMQ215VSokltryDikb7A4nvd5GlC0zV0qbe5H4Mw7BFEWgca//cGYaCptob6ygrMNhuxvVIPidQdB+PwPGvlsHDTSdlsKmvgya830D8pkkm5KcEu0kGJ7ZXCxQ8+wU/vvcGKzz9m14a1nHnrncSl9mHl79Vc+dpSIq0m3rl2FFlJkcEubsA0N02ZDCastP/LXtM1NKntrpn4vW4OME7Nicfjfd2WSLsJq8uAZhW4InUcjtKgN4cFiwoWyiFLCMHjFwxnR2UTt89aRXr8GIb2jgl2sQ6K0WTmxCuuI33ocL6a8Sxv330bmWdfyd/yTSRGWXn32lGkxYVmG3qwGA1GjBixGPc9R0czXeotNRZN13BrLlyVdUiXGz3CjCdc4NE92F0ONKntM/dVoJvDgkn1WSiHvLJ6B+e8sAgBzPvTsSRH9ay5sVurqyjnvccfpvH3zexMGsad991N76SeGQRbC4U5uD1uNzUlRWhuN1GJSYRH7/nd7qs5rHXNpbmfpauawzpaa1F9FopykJKjbLxyxUim/GcxN7y9gvev65kjpJotKvbwjPlkTk1LInPXEr557G+cedtdJKZnBLtoPZ7L4aCmpAiQxKb2xhq2d23tQJvDdKnvDiL+QcUvoLTXHCaEoG9UXyItwWtm7Fn1IEU5SLl9YvjXhXms/L2Guz8M7Vn29ufDFYXc9N6vDE+P5+GH7+xx832HMkdDPdVFhQiDAUN4FJdefsVBpSj/+uuvW5avWLGCvOF5DBk4hLv/cjdRlijibHG8/Z+3mXD0BM449gymnT8Nc52ZQfGDGJIwhCWfLuGcMedwzuhzWDhvIb0iehFvi+faq65l4MCB5Obmcs011+B2d/ONjZ29UaM7H+qmPKWznvtuo8y46zP5wvxNwS7KAXtr8XaZcddn8tJXlshG5+4b9Rqqq+TsB/4un5p6hvxs+hPS2dQYxFJ2TjBSlOu6LuurKmXx5o2yovB36fG4DypFucPhkFu3bpVZWVnS4/FIKbs2Rfnnn38udV2Xuq7Liy66SM6YMaPN81EpyhWlC/xpfDZn5/Xmya838PVvJcEuToe9tGAL98wr4OTBvZh55UjCLbtbkHvCfN+hSkqduvIyGqoqsUVGEZ/ah8bGpgNOUX7RRRdhtVrp168f2dnZLF26tMtTlJ9++ukIIRBCcMwxx+yRd6o7qD4L5bAihOCJycPZUeUdITX3D2MZ0jt070uQUvKvbzfy3PzNnJXXm2em5mFuY3a75vm+0wYP5fPnnuT9e/7KCZdP44hTz+yxwznv//Q31hbVdek+h/SO5r6zhgKgaxo1pcW47HYi4+KJiItHCHHAKcp37drF6NGjW47RnFY8UCnK3W43b7/9NtOnT++ib6VjVM1COezYzEZeuXwEMWFmrn1zGeX1zmAXqU1SSh76fB3Pzd/MhSPTefbCI9oMFP5Cab7vUOZxu6kqKsTtcBCT3IvI+ISWoHogKcph32nF97Xc38GkKL/xxhsZN24cxx/fsey9XSUoNQshxJPAWYAL2AJcLaWsCUZZlMNTcrR3hNTk//zMDW8v570QGyGl6ZJ/zlvD+0t3cvWxmdxzxhAMho7VEJrn+17x+Tx+eu9N3r77Fs645U56DxgU4FJ3reYaQFdzOezUlBQDkrjU3lhajXjqSIry4uLilmyxaWlp7Ny5s2X7QKYov//++ykvL+ell17qqq+jw4JVs/gWyJVSDgc2An8LUjmUw1hunxj+NfUIfv29hr99tCZkRhK5NZ0/z17F+0t3cvP4bO49s+OBopkQgpFnnsdFDzyOEIJZ/3cXSz+ei9TbHvN/uLA31FNdtAthMBDfO32vQAHtpygHePPNN1vSjZ999tl88MEHAU9RPnPmTL7++mvef/99DIYgXLo720Pe2QdwHvBuR9ZVo6GUQJjuGyH14g/BHyHlcHvktW8ukxl3fSZn/LC5S/Zpb6iXnzz9iHxq6hly7iP3ysbami7ZbyAEajSU/4inysKd0uPZf9r3faUoHz9+vMzOzpbjx4+XlZWVLet3R4pyo9Eos7KyWra5//772yz7IZuiXAjxKTBLSvlOe+uqO7iVQJBScssHq/hsdREvXTaCiUODk0OqyeXhhrdX8NOmCh44ZyhXjMnssn1LKVvm+7ZFRgVkvu+uEIg7uL0jnsqx19cRFhVFdGIyIhi/zLtJoO7gDtg3JoT4TghR0MbjHL91/gF4gHf3s5/rhRDLhRDLy8t7xhy+Ss8ihODJycMZ3ieG22at6vIROB1R53Bz5WtLWbS5gicnD+/SQAGt5vu2hXX5fN+hStc0qouLsNfXERkXT3RSr0M6UARS0GoWQogrgT8AE6SUTR3ZRtUslEAqq3Nw9guLMBoE8246lqSo7pmzoLrRxRWvLWVdcR3TLzqSM4YHdv5wl8PO9zNnsPanH0gfMqzT8313pa6sWXjcLmqKi9E8bqKTkgmLCt0h0l2px9Us9kcIMQm4Czi7o4FCUQItOdrGzCtHUtno5A/vrMDhDvyv7rI6Bxe+vJgNpfW8fMWIgAcK8M33/ae/dMl836HKZbdTtasQXdeIS+1z2ASKQApWfewFIAr4VgixSgjxnyCVQ1H2kNsnhmemHsGKHdX8PcAjpAqrm5j60mIKq+28cfXRjB/UK2DHasvQEyZw2SPPEhEbx0eP3sfCd19H83i6tQyBYK+vo7p4FwaDkfg+aVjCwoJdpENCUIKFlDJbSpkupTzC9/hDMMqhKG05fVgqt588gI9W7uI/C7YG5BjbKhqZ+p/FVDa6eHvaKMb2b3t2uEBLSEvnkoefZvjJk1j2yYfM+r+7qCsvC0pZOktKSUNVJbVlpZhtNuL7pGEytz+XhdIxqqdHUdpwy4RszhyeyhNfr+fbtaVduu8NJfVM+c9inB6dD64fzYiMuC7d/4EyW6ycct2fOPO2u6gs3Mlbd93MpmWLg1qmAyV1nbryUhqqqwiLiiYutTcGY+jcZHkoUMFCUdoghOCpKXkM7xPDrR+sZF1x14yQWl1Yw4UvL8ZogFk3hNbMfQPHHM/lj00ntlcqnzz1MPNffwlPd6fBPgia5vGNeKonMj6B6KRkRCdnoaupqWHy5MldkqJ82LBhZGdnc8stt7Q0az7zzDMMGTKE4cOHM2HCBHbs2NGyzZtvvklOTg45OTktNwECTJs2jby8PIYPH87kyZNp6O40Lp29UaM7H+qmPKW7ldTa5aiHv5NjH/1eltc7OrWvX7ZWyqH3fiWPe/x7uaMidNOIu10uOf+Nl+VTU8+Qb911i6wqKuy2Yx/oTXlup1OW7dgmS7Zskk31dV1WjlBMUV5bW9tSvttvv72lLK2pFOWKEgS9fDmkKhud3PD2CpyegxshtXBjOVe89gvJ0VZm3zCGvgmhO1+2yWzmpCuv49w776GuvIy3776NdYsWBLtYe3HZm6gq2onUdeJ69yEsMqpL9ltXVxeSKcqjo70juqSU2O32bs8mrFKUK0o7hqXF8PSUI7jpvV/520dreHpK3gH9oX79Wwk3v7eS/smRvD3tGBIju+f+jc7qP2IUlz/+HJ8/9yRfPPckOwvyOemq6zFbu2kO8y/vhpI1bX6kaR5wu4kTApPF0vFmp5RhcNq+M8gCIZ2i/Oqrr+aLL75gyJAhPP300x075y6iahaK0gFnDE/ltpNz+OjXXby0sOMjpD5etYsb3/2VoX2i+eC60T0mUDSLTkziwvseZdR5U1nzw7e8+/c/U7FzR/sbBohEonncaG4XwmDAZLV2un+itVBOUf76669TVFTE4MGDmTVrVofPqSuomoWidNCtE3LYVNbA41+tp39SJKcM2f99Ee8v/Z2//3cNo/rFM/PKo4m09sw/N4PRyHEXXUHakGF8+cLTvPv3PzP+mhvIPfGUwDaFtKoBSF2ntrwMR0M9YVHRRCcldXmggNBOUQ5gNBq58MILefLJJ7n66qu7+vT3SdUsFKWDhBA8NTmPYX1iuO2Dlawv2fcIqZk/beVvH63hxAFJvHH1MT02UPjLHH4kVzzxPL0HDOKb/zzHF88/hcvePQkYNM1DVfEuHA31RCV0zYinfQnFFOVSSjZv9k6VK6Xk008/ZdCgbp6fpLM95N35UKOhlFBQUmuXxzz8bZsjpHRdb0l5/sd3lkunWwtSKQNH0zxy8YcfyKcvPEvOvOVaWbK1a1KpS9n2SB7/EU/2LhzxtD+hlqJc0zQ5duxYmZubK4cOHSovueSSPUZH+TtkU5QfCJVIUAkVqwtrmPrSYnJ7x/DudaOwmoxIKXnsy/W8tHArFxyVxuMXDMPUzjSoPVnhugI+f+5J7HW1XTbfd+skeE57E7UlxSAEsSm9sdi6qXO9BzukEgkqSk83PC2Wp6bksXxHNX//qABdl9zzcQEvLdzKFWMyeHLy8EM6UEDg5/tuqqulprgIg8lEfJ90FSiCrOc3pCpKkJw5vDebShuY/v0m1hbXsa64jj+c0J+7Jg3s9jHwwbLnfN9vdMl831J6czw11lRjCQ8nNjlFpe4IAYf2Tx9FCbBbJ+RwxvBU1hXX8deJAw6rQNFs93zfT3R6vm8pJbVlJTTWVBMeHUNcisrxFCpUn4WidJJb09lc1sDgVDVngqOxgW9fep6Nvywi84gRnHbTnwmP7lj+q8aaatatW0dachJRCYmEx8QedoG3K6g+C0UJUWajQQUKH1tEJGfefjcTpt3Izt9W8/adN7Pzt9Xtblexcwfv/fMvaB4PsSmpRMTGqUARYlSwUBSlSwkhOGKid75vcwfm+96+eiXv33MHmttNREwstojIbi6x0hEqWCiKEhDJmVlc9tizDD7uBBbPfY+5D/6ThqrKPdZZ/f1XfPTofUQnJXPJw09jNJuDVNo9hWKK8mY333wzkZFBCKidvVGjOx/qpjxF6ZkKfvxOPnv5+fLFay+RW1cul7qmyR/fflU+NfUMOfeRe6XDl677QFOUB0oopiiX0nuT32WXXSYjIiL2WXaVolxRlB6r9Xzfb//tNpZ/+hF5E8/gvDvvxRoeOinbQzVFuaZp3HHHHTzxxBPd+G3spu6zUBSlWzTP9/3jm6+wZv43nHjFdRx1+tn77Mh+fOnjrK9a36VlGBQ/iLuOuWu/64RqivIXXniBs88+u6UM3U0FC0VRuk3zfN8nXnkdZktopmtvTlH+/PPPM2rUKG699daApyhfsGDBfvdVVFTEnDlz9shI291UsFAUpdt1JFC0VwMIlFBMUb5y5Uo2b95MdnY2AE1NTWRnZ7dkou0Oqs9CURTFTyimKD/jjDMoKSlh+/btbN++nfDw8G4NFKBqFoqiKHt5/vnnufTSS3G5XGRlZfH666+j6zpTp07l1VdfpW/fvsyZMweAoUOHMnXqVIYMGYLJZOLFF1/E6EtR8u9//5urrroKu93Oaaed1tI3cccdd9DQ0MCUKVMA6Nu3L5988gnx8fHcc889HH300QDce++9xMfHB+Eb2JtK96EoSshoK1WFcmBUug9FURQlaFSwUBRFUdqlgoWiKIrSLhUsFEVRlHb1qNFQK1asqBBC7Gh/zTYlAhVdWZ4gU+cTug6lc4FuPJ9vv/12mKZpnkAeQ9M0k9FoDOgxulPr8ykpKTENGTJkTavVMjp7nB4VLKSUSQe7rRBieWdHA4QSdT6h61A6F+je88nPz9+em5sb0MBUUFAwODc3d10gj9GdWp+PpmmJgfj3Us1QiqIorVRUVBgnTZqU1a9fv6FZWVlDv/vuu4jS0lLj2LFjczIyMnLHjh2bU15e3jLf69/+9reUvn375mZmZuZ++OGHLTNh/fTTT+EDBgwY0rdv39yrrroqXfdNNft///d/vfr37z90wIABQ8aMGTNg48aNluZtnn/++YSMjIzcjIyM3Oeffz6hefkFF1yQ2adPn2GDBg0aMmjQoCE///xzWDd9HYAKFoqiKHu5/vrr0ydOnFi3bdu239auXbv2iCOOcNx3332pJ554Yv2OHTsKTjzxxPp77703BWDFihW2jz76KH7Dhg2/ffXVVxtvu+22vh6Pt1XoxhtvzJgxY8aO7du3F2zdutU2d+7caIARI0Y0rVq1at3GjRvXnnvuudW33357GkBpaanx8ccf77106dJ1y5cvX/f444/39g9KDz30UOH69evXrl+/fu3YsWPt3fmdHE7B4uVgF6CLqfMJXYfSucAhdj6JiYnl+/u8qqrK8Msvv0TddtttFQA2m00mJiZqX331VewNN9xQCXDDDTdUfvnll3EAc+fOjT3//POrwsLC5KBBg1wZGRnOH3/8MWLHjh3mhoYGw8knn9xoMBi49NJLK+fNmxcHcNZZZ9VHRUXpAMcdd1xDcXGxBWDevHkx48aNq+vVq5eWlJSkjRs3ru6jjz7a7yTm7Z1PV+lRfRadIaU8pP7Dq/MJXYfSuUDwzqfo7/9Id27a1KUTXVhzcpp6P/Lwzv2ts379emt8fLxnypQpmWvXrg0fPnx44yuvvLKzsrLSlJGR4QbIyMhwV1VVmQB27dplGT16dEPz9r1793bt3LnTYrFYZGpqqrt5eUZGhqu4uHivqQBfeumlpJNPPrnWty9zWlqaq/mzPn36uHbt2tWyzf3339/n0UcfTT3++OPrX3jhhcKwsDCZkpLSLYMPDqeahaIoSrs8Ho9Yt25d+E033VS+bt26teHh4fo999yTsq/195FWXHYkRfmMGTPi8/Pzw++///6S/ewLgGeeeWbX1q1bC/Lz89dVV1cb91emQDhsahaKovQs7dUAAiUzM9PVq1cv1/jx4xsBLrzwwurHHnssJSEhwbNjxw5zRkaGe8eOHeb4+HgPQFpammvnzp0tHdRFRUWWtLQ0d2Zmptu/JrFjxw5LSkpKS01j3rx5UU899VTqTz/9tCEsLEz69uVesGBBVPM6u3btspxwwgn14K3NAISFhclrrrmm8umnn+4V6O/C32FRsxBCTBJCbBBCbBZC3B3ksrwmhCgTQhT4LYsXQnwrhNjke47z++xvvnJvEEKc6rd8hBBije+z54Tv54cQwiqEmOVb/osQItNvmyt9x9gkhLiyC84lXQjxgxBinRDiNyHErT38fGxCiKVCiHzf+dzfk8/Hb79GIcRKIcRnPf188vPzh61Zs2ZIQUHBkIKCgsEAbrfbuG7dupzVq1fnrlu3Lsftdrd0CBcWFqasXr06d/Xq1bnV1dUto5Tq6+vD16xZM2T16tW527ZtS2/+Ra/runA6nX2Tk5PNn3766WCHw2H55ptvogcOHOiYMGGC4+mnnx6yevXq3Oeffz590qRJNQAXXHBBzUcffRRvt9vF+vXrLdu3b7edeOKJjRkZGe6IiAj9+++/j9B1nXfffTfhnHPOqQFYtGhR2M0335w5Y8YMZ2Vl5cA1a9YMraurizjzzDMbFixYkPDTTz/lLl68OGfBggXR5557bi3AsmXLeq9evTo3Pz8/d86cOcmDBw+2t3Uuft9/1/7bdHYS71B/AEZgC5AFWIB8YEgQyzMOOAoo8Fv2BHC37/XdwOO+10N85bUC/XznYfR9thQYAwjgS+A03/Ibgf/4Xl8EzPK9jge2+p7jfK/jOnkuqcBRvtdRwEZfmXvq+Qgg0vfaDPwCjO6p5+N3Xn8G3gM+C/X/b6tWrdoupVy+r8eqVaucLpdrlf+y7du3l+zcubNQSrl8586dhdu3by+WUi5vaGgoWLNmTZOmaSuamppW5+fnO3RdXy6lXF5QUNBYW1u7Ttf15evWrautrKzcKKVcXlRUtGPLli1lixYt+m3IkCGO7Oxsz4QJE6p37dq1+scff3SOGjWqvm/fvo6jjz5aLywsXN1chjvvvLMwLS3NkZGR4Zg1a9bG5uULFixYm52dbU9LS3NcfvnlZZqmLZdSLh8zZkxdfHy8PmDAANfAgQObTjrppBq3271y+/btJQ8//HBVenq6Iy0tzfPwww/XNJ/LyJEjtezsbHv//v0dp512mlZdXf1rW+eyZMmS0kD82wT9Yh7oh+8/+Nd+7/8G/C3IZcpkz2CxAUj1vU4FNrRVVuBr3/mkAuv9ll8MvOS/ju+1Ce+dt8J/Hd9nLwEXd/F5fQycciicDxAO/AqM6snnA6QB3wPj2R0sQvZ8DiZY5OfnOxwOR76UcrnD4cjPz893SF/gaA4iUsrl69atq62trV3nW8fevLysrGzrli1byv3XkVIu1zRt+a+//urWdX2PdaSUy7ds2VJeVla2dX9l3d/D7Xb/umrVKmdz8OrKc/nxxx/rA/Fvczj0WfQB/Ns+C/FeAEJJLyllMYCUslgI0Tx1Vh9gid96hb5lbt/r1subt9np25dHCFELJND299CHLuKr4h6J99d4jz0fIYQRWAFkAy9KKX8RQvTY8wGeBe7EW/Nr1pPPhw0bNuSAd8hoSkpKhcfjMVmtVjeA1Wp1ezweE4Db7bZERES0jFIym80ul8tlEUJIs9nc0ndgsVhcbrfb3LyN1Wp1ARgMBoxGo+bxeEwul8tsNptdrfa118imjnI4HFaTyeTZsmVLpsPhCA8LC2vMzMzc2RXnomlaczNcl/7bHA7BQrSxrKfM+LSvsu/vnA5mm04RQkQCHwK3SSnrhGjrUAddtm49HymlBhwhhIgF/iuEyN3P6iF9PkKIM4EyKeUKIcSJHdlkH2UIifMBGDRo0Hqr1ep2uVymjRs3DggLC3Mc4C4Opgxdfr2QUgq73R6enp7+e3R0dOO2bdvSd+3adaCjm9osl9/fX5f+2xwOHdyFQLrf+zSgKEhl2ZdSIUQqgO+5zLd8X2Uv9L1uvXyPbYQQJiAGqNrPvjpFCGHGGyjelVJ+1NPPp5mUsgb4EZjUg8/nWOBsIcR24ANgvBDinR58PjT/6rZYLJ6YmJiahoaGCJPJ5HE6nWYAp9NpNplMHtj967t5W7fbbbFYLG6LxeJurkkAuFwuS/Ovc7PZ7HI6nRYAXdfRNM1oMpk03zZ77asT5+Eym82u6OjoRoD4+Phqu90e3hXnYjAYNN/bLv23ORyCxTIgRwjRTwhhwdvR80mQy9TaJ0DziIQr8bb9Ny+/yDeqoR+QAyz1NSHUCyFG+0alXNFqm+Z9TQbmS2/D5NfARCFEnG/0y0TfsoPmO/arwDop5TOHwPkk+WoUCCHCgJOB9T31fKSUf5NSpkkpM/H+v58vpbysp56PpmkGj8djaH5dX18fHRYWZo+Ojq4pLy9PACgvL0+IiYmpAYiLi6upqamJ13Vd2O12i9PptEVFRTVarVa30WjU6+rqIqSUVFZWJsTGxtYAxMTE1FRUVCQAVFZWxkVGRtYLIYiNja2tr6+PdrvdRrfbbayvr4+OjY2tPdhzsVgsHrPZ7GpqarIC1NXVRVutVkdXnIvNZmvyHaZL/20O+WYoX1vdn/B+GUbgNSnlb8EqjxDifeBEIFEIUQjcBzwGzBZCTAN+B6YASCl/E0LMBtYCHuAmXzMJwB+BN4AwvKNTvvQtfxV4WwixGe+viIt8+6oSQjyIN3gCPCClrOrk6RwLXA6sEUKs8i37ew8+n1TgTV+/hQGYLaX8TAixuIeez770yH8fl8tl2rJlS7Zv/yIuLq4yPj6+LioqqnHz5s39V69enWg2m13Z2dlbACIiIhyxsbFVBQUFQwHS09N3NDfR9O3bd8f27dv7SSlFVFRUXVxcXC1AcnJyxZYtW/qtXr0612g0allZWVsAzGazlpKSUrRu3brBACkpKUVms1lro5gd1rdv39+3bt2aJaUUFovFmZWVtR2gs+dis9mac0Z16b+N8PWGK4qiBF1+fv72vLy8Q2kukG6Xn5+fmJeXl9nV+z0cmqEURVEOSCimKNd1nZtvvrlPZmZmblZW1tCHHnqoeRRbt1DBQlEUpZVQTFH+/PPPJxQWFpq3bNlSsHXr1t+uvvrqQDVTtkkFC0VRFD+hmqJ85syZyQ8++GCx0eit0PTp06dbp4Y95Du4FUXpmb5/a1161a6GLk1RHt8nsmnCFYN7ZIrynTt3Wt9+++24zz//PC4+Pt7z4osv/j5s2DBnZ7+TjlI1C0VRFD+hmqLc5XIJm80mCwoK1k2bNq38qquuyjywM+scVbNQDlu+eyoukVLO8L3vDTwnpZwcgGOdCwyXUj6wn3WeAr6QUs7v6uP3RO3VAAIlVFOU9+rVy3XJJZdUA1x++eU1f/rTnzID/FXsQdUslMNZLN7MnABIKYsCESh87gRmtLPO83izwCpB1LdvX09KSoorPz/fCtCcovzUU0+teemllxIAXnrppYQuSFGe8fHHH2/273s499xzaxcsWBBdXl5uLC8vN/qnKD/ttNNqvvzyyyiAL774IiojI6PbmqBA1SyUw9tjQH/fDYXfAi/izcyaK4S4CjgX742cucDTeFPcXw44gdN9Nzf1922XBDQB10kp1/sfRAgxAHBKKSuEEFHAamCAlNIthIj2vc+RUu4QQiQIIVKklCWBPnll355//vnfL7300iyXyyX69u3rfP/997drmsZ5553XPyMjI7F3796uefPmbQEYOXKk49xzz60aMGDAUKPRyDPPPLPDZPJeWmfMmLFj2rRp/RwOhzjppJPqpkyZUgtwxx13pDc1NRmnTJnSH7z9HPPnz9/cq1cv7Y477igaMWLEYIA777yzqFevXhrAAw88UDJ58uR+M2bM6BUeHq6/8sor27vzO1E35SmHLV+m3M+klLmt3/uCxT/xZtK1AZuBu6SU/xFC/AvYIaV8VgjxPfAHKeUmIcQo4FEp5fhWx7kayJVS/sX3/nXgYynlPCHE9cBAv89eAb6SUn4Y8C8gBKmb8jovUDflqZqFouzbD1LKerx5kWqBT33L1wDDfdl2xwJz/DourW3sJxUo93s/E2+z1DzgauA6v8/KgN5ddQKK0lVUsFCUffNvE9b93ut4/3YMQI2U8oh29mPHm/ETACnlIiFEphDiBLwz0RX4rWvzra8oIUV1cCuHs3r2nBTogEgp64BtQogp4M3CK4TIa2PVdXgnU/L3FvA+8Hqr5QOAAhQlxKhgoRy2pJSVwCIhRIEQ4smD3M2lwDQhRD7wG3BOG+ssBI4Uew6yfxfv/MfvNy8Q3rlBsoHlB1kWRQkY1QylHNaklJe0WpTrW/4G3pTczetl+r1u+UxKuQ3vBEn7O0aTEOI7YALwnW/xccBc3yRLzc70LevWNA6K0hEqWChK93gE39zvQojngdOA01utY8I7RFdRQo5qhlKUbiClLJVSfuJ7fbOUMltKubHVOnNa1TSUIAnFFOUjRowYOGjQoCGDBg0akpycPPzkk0/u301fB6CChaIoyl5CMUX5ihUrNqxfv37t+vXr1x555JGN5557bk13ficqWCiKovgJ1RTlzaqrqw2LFy+Oas4T1V1Un4WiKCHp638/m16xc0eXpihPTM9oOvWPt/XIFOXN3n333bixY8fWxcfH6wf7PRwMVbNQFEXxE6opypvNnj07/qKLLurWWfJA1SwURQlR7dUAAiVUU5QDlJSUGFevXh0xderUzYH9FvamahaKoih+QjVFOcBbb70VP378+Jrw8PBuzwCrahaKoiithGKKcoC5c+fG33nnncXd/oWgUpQrihJCVIryzgtUinLVDKUoiqK0SwULRVEUpV0qWCiKoijtUsFCURRFaZcKFoqiKEq7VLBQFEVR2qWChaIoSiuhmKL8448/jhoyZMjgQYMGDRkxYsTAgoICazd9HYAKFoqiKHsJxRTlt956a8Y777yzbf369WunTJlSdd9996V253eigoWiKIqfUE5RXlNTYwSora01+me07Q4q3YeiKCGpau7GdHdJY5emKDenRDTFTx7QI1OU/+c//9l+/vnn51itVj0yMlJbtmzZus5+HwdC1SwURVH8hGqK8meeeabXRx99tKm0tHT1JZdcUvHHP/4x/QBPrVNUzUJRlJDUXg0gUEIxRXlRUZFp3bp1Yc1luuKKK6onTZqUE/hvYzdVs1AURfETiinKk5KSPA0NDcbVq1dbAT777LPo7OxsR3d+L6pmoSiK0koopiifPn36jsmTJ/cXQhATE6O98cYb27rzO1EpyhVFCRkqRXnnqRTliqIoStCoYKEoiqK0SwULRVEUpV0qWCiKoijtUsFCURRFaZcKFoqiKEq7VLBQFEVpJRRTlH/yySdRQ4YMGZyTkzP0/PPPz3S7uzWPoAoWiqIorYVainJN07j++uv7ffDBB1s3bdr0W9++fV0vvPBCYnd+JypYKIqi+AnFFOWlpaUmi8WiDx8+3AkwadKkunnz5sV25/ei0n0oihKS5s2bl15WVtalKcqTk5Obzj333B6XojwlJcXj8XjEwoULw8eNG9c0a9asuOYA011UzUJRFMVPKKYoNxgMvPXWW1tvv/329GHDhg2OiorSjEbjXusGkqpZKIoSktqrAQRKKKYoBzj55JMbV6xYsQHgo48+it68ebMt0N+FP1WzUBRF8ROKKcoBdu3aZQKw2+3iySefTPnDH/5Q3p3fi6pZKIqitBKKKcofeOCBlG+//TZG13VxzTXXlJ199tn13fmdqBTliqKEDJWivPNUinJFURQlaFSwUBRFUdqlgoWiKIrSLhUsFEUJJbqu66L91ZS2+L47PRD7VsFCUZRQUlBeXh6jAsaB03VdlJeXxwAFgdi/GjqrKErI8Hg815aUlMwsKSnJRf2YPVA6UODxeK4NxM7V0FlFURSlXSpyK4qiKO1SwUJRFEVplwoWiqIoSrtUsFAURVHapYKFoiiK0q7/B674n/wEiAtvAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "swiftdiff['px'].plot.line(x=\"time (y)\")" ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 8, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "
<xarray.Dataset>\n",
+       "Dimensions:   (id: 18, time (y): 5)\n",
+       "Coordinates:\n",
+       "  * id        (id) int64 2 3 4 5 6 7 ... 6000205 6000206 6000207 6000208 6000209\n",
+       "  * time (y)  (time (y)) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05\n",
+       "Data variables:\n",
+       "    Mass      (time (y), id) float64 0.0 0.0 0.0 0.0 0.0 ... nan nan nan nan nan\n",
+       "    Radius    (time (y), id) float64 0.0 0.0 0.0 0.0 0.0 ... nan nan nan nan nan\n",
+       "    px        (time (y), id) float64 0.0 0.0 0.0 0.0 ... -0.8965 0.0579 2.692\n",
+       "    py        (time (y), id) float64 0.0 0.0 0.0 0.0 ... -2.79 0.004403 -2.209\n",
+       "    pz        (time (y), id) float64 0.0 0.0 0.0 ... -0.1535 0.002265 0.06436\n",
+       "    vx        (time (y), id) float64 0.0 0.0 0.0 ... 0.02465 -6.413e-05 0.01524\n",
+       "    vy        (time (y), id) float64 0.0 0.0 0.0 0.0 ... -0.00847 0.00055 0.0203\n",
+       "    vz        (time (y), id) float64 0.0 0.0 0.0 ... -1.248e-06 0.001067
" + ], + "text/plain": [ + "\n", + "Dimensions: (id: 18, time (y): 5)\n", + "Coordinates:\n", + " * id (id) int64 2 3 4 5 6 7 ... 6000205 6000206 6000207 6000208 6000209\n", + " * time (y) (time (y)) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05\n", + "Data variables:\n", + " Mass (time (y), id) float64 0.0 0.0 0.0 0.0 0.0 ... nan nan nan nan nan\n", + " Radius (time (y), id) float64 0.0 0.0 0.0 0.0 0.0 ... nan nan nan nan nan\n", + " px (time (y), id) float64 0.0 0.0 0.0 0.0 ... -0.8965 0.0579 2.692\n", + " py (time (y), id) float64 0.0 0.0 0.0 0.0 ... -2.79 0.004403 -2.209\n", + " pz (time (y), id) float64 0.0 0.0 0.0 ... -0.1535 0.002265 0.06436\n", + " vx (time (y), id) float64 0.0 0.0 0.0 ... 0.02465 -6.413e-05 0.01524\n", + " vy (time (y), id) float64 0.0 0.0 0.0 0.0 ... -0.00847 0.00055 0.0203\n", + " vz (time (y), id) float64 0.0 0.0 0.0 ... -1.248e-06 0.001067" + ] + }, + "execution_count": 8, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ - "ax = swifterdat['px'].sel(id=6000200).plot.line(x=\"time\")\n", - "swiftestdat['px'].sel(id=6000200).plot.line(x=\"time\")" + "swiftdiff" ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 9, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "
<xarray.DataArray 'vx' (time (y): 5)>\n",
+       "array([ 0.00000000e+00, -4.35675281e-11,  2.08934259e-12, -1.73709378e-10,\n",
+       "       -2.60294819e-10])\n",
+       "Coordinates:\n",
+       "    id        int64 6000200\n",
+       "  * time (y)  (time (y)) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05
" + ], + "text/plain": [ + "\n", + "array([ 0.00000000e+00, -4.35675281e-11, 2.08934259e-12, -1.73709378e-10,\n", + " -2.60294819e-10])\n", + "Coordinates:\n", + " id int64 6000200\n", + " * time (y) (time (y)) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05" + ] + }, + "execution_count": 9, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ - "swiftdiff.sel(id=6000200)['py']" + "swiftdiff['vx'].sel(id=6000200)" ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 10, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "
<xarray.DataArray 'px' (time (y): 5)>\n",
+       "array([ 0.00000000e+00,  7.07309283e-09, -5.01222758e-09,  1.41219044e-08,\n",
+       "        4.51840986e-08])\n",
+       "Coordinates:\n",
+       "    id        int64 6000200\n",
+       "  * time (y)  (time (y)) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05
" + ], + "text/plain": [ + "\n", + "array([ 0.00000000e+00, 7.07309283e-09, -5.01222758e-09, 1.41219044e-08,\n", + " 4.51840986e-08])\n", + "Coordinates:\n", + " id int64 6000200\n", + " * time (y) (time (y)) float64 0.0 1.47e+05 2.94e+05 4.41e+05 5.88e+05" + ] + }, + "execution_count": 10, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ - "swiftestdat['px'].sel(id=6000200)" + "swiftdiff['px'].sel(id=6000200)" ] }, { @@ -110,9 +2161,7 @@ "execution_count": null, "metadata": {}, "outputs": [], - "source": [ - "swiftestdat['px'].sel(id=6000200)" - ] + "source": [] } ], "metadata": { diff --git a/examples/rmvs_swifter_comparison/mars_ejecta/tp.in b/examples/rmvs_swifter_comparison/mars_ejecta/tp.in new file mode 100644 index 000000000..f2fb549d8 --- /dev/null +++ b/examples/rmvs_swifter_comparison/mars_ejecta/tp.in @@ -0,0 +1,601 @@ + 10 + 6000200 + 0.321794727714005 -1.39371227734394 -3.555372224179648E-002 + 1.263074812292130E-002 3.554726326213714E-003 5.297001107625803E-004 + 6000201 + 0.327226184487893 -1.39016803345029 -3.472777857996821E-002 + 1.531652992677778E-002 5.307306923696688E-003 9.381180638637231E-004 + 6000202 + 0.327531202012097 -1.39478453402311 -3.768059236489047E-002 + 1.546735696093592E-002 3.024509922426058E-003 -5.220083559665117E-004 + 6000203 + 0.322202594471783 -1.39306474804903 -3.456946701305916E-002 + 1.283243270995610E-002 3.874920791552106E-003 1.016400989319078E-003 + 6000204 + 0.323275102327213 -1.39231827262319 -4.063286630727199E-002 + 1.336277331700188E-002 4.244042778378470E-003 -1.981867830972686E-003 + 6000205 + 0.325545380653237 -1.39090814848439 -3.352094137611809E-002 + 1.448539519756952E-002 4.941330064430509E-003 1.534882714168952E-003 + 6000206 + 0.322100032022230 -1.38941601200018 -3.760532734717235E-002 + 1.278171696789352E-002 5.679171344404378E-003 -4.847908241938090E-004 + 6000207 + 0.325337949472342 -1.39379280734375 -3.374182250017822E-002 + 1.438282328833895E-002 3.514905332112549E-003 1.425659990227191E-003 + 6000208 + 0.326082121124172 -1.39530189188417 -3.547938702480286E-002 + 1.475080609012429E-002 2.768683467912135E-003 5.664578686714714E-004 + 6000209 + 0.327275685332658 -1.39441408865592 -3.527833381550945E-002 + 1.534100742400883E-002 3.207690138697890E-003 6.658759573365761E-004 + 6000210 + 0.327855145958387 -1.39318452150569 -3.930697139978061E-002 + 1.562754285445135E-002 3.815694439606817E-003 -1.326230758891795E-003 + 6000211 + 0.328078710851062 -1.39251497560023 -3.504731288890989E-002 + 1.573809266562045E-002 4.146775820637334E-003 7.801126762871974E-004 + 6000212 + 0.321277816542297 -1.39195915492794 -3.578170284173104E-002 + 1.237514255024041E-002 4.421621614318133E-003 4.169667921656008E-004 + 6000213 + 0.324665498950826 -1.39490683098100 -3.970778222771389E-002 + 1.405030561400833E-002 2.964035733299919E-003 -1.524426284193035E-003 + 6000214 + 0.324929781829885 -1.39405065211182 -4.042853272331495E-002 + 1.418098991790633E-002 3.387404586890700E-003 -1.880827641470143E-003 + 6000215 + 0.327555988964697 -1.39432357646482 -3.871524764041770E-002 + 1.547961377324917E-002 3.252447191181946E-003 -1.033631374975754E-003 + 6000216 + 0.321438220244275 -1.39139099683337 -3.554066697282077E-002 + 1.245446000814854E-002 4.702568096196656E-003 5.361557644332466E-004 + 6000217 + 0.327269239795964 -1.39111145972877 -3.426641502936644E-002 + 1.533782019342066E-002 4.840795408002274E-003 1.166256090318398E-003 + 6000218 + 0.323989604933645 -1.39344910147082 -4.061548100580529E-002 + 1.371608517771453E-002 3.684863230669170E-003 -1.973271034886913E-003 + 6000219 + 0.321638711062574 -1.39087780679511 -3.541758920393877E-002 + 1.255360000208438E-002 4.956333618793669E-003 5.970160540173710E-004 + 6000220 + 0.322548442178556 -1.38969900689356 -3.516409557966647E-002 + 1.300344971633673E-002 5.539234202882858E-003 7.223652175672305E-004 + 6000221 + 0.326425905823290 -1.38912149826721 -3.516059934460978E-002 + 1.492080296716386E-002 5.824804396075300E-003 7.240940584451341E-004 + 6000222 + 0.325686803560988 -1.39394570042399 -4.039235894274407E-002 + 1.455532690983334E-002 3.439301774920795E-003 -1.862940196963350E-003 + 6000223 + 0.327419177059503 -1.39257520407836 -3.429421078579946E-002 + 1.541196213929194E-002 4.116993654016758E-003 1.152511465264746E-003 + 6000224 + 0.324490052620445 -1.39553030085443 -3.872774571953116E-002 + 1.396354978011187E-002 2.655738325987235E-003 -1.039811505806850E-003 + 6000225 + 0.321731345032329 -1.39031854790183 -3.853949894445081E-002 + 1.259940624537157E-002 5.232879566177696E-003 -9.467260253928260E-004 + 6000226 + 0.325601002353804 -1.39413056037783 -3.398693779254160E-002 + 1.451289937508628E-002 3.347891031734538E-003 1.304453798319087E-003 + 6000227 + 0.325825665573909 -1.38837078808078 -3.783072474317082E-002 + 1.462399229429218E-002 6.196020414654821E-003 -5.962467834347196E-004 + 6000228 + 0.328007765778950 -1.39171771880483 -3.491579285982123E-002 + 1.570301128843478E-002 4.541008506865058E-003 8.451475491903844E-004 + 6000229 + 0.326038941297358 -1.38837178473418 -3.696108906496336E-002 + 1.472945425064948E-002 6.195527583046642E-003 -1.662237200565085E-004 + 6000230 + 0.324865697509285 -1.38818972404806 -3.707510223923227E-002 + 1.414930108941243E-002 6.285554126260012E-003 -2.226016903873118E-004 + 6000231 + 0.326890450931510 -1.39246034816528 -3.388258567264032E-002 + 1.515051423076515E-002 4.173788347277624E-003 1.356054508127195E-003 + 6000232 + 0.325593633813168 -1.38828971482136 -3.764909195799739E-002 + 1.450925573155059E-002 6.236110043272472E-003 -5.064318314410409E-004 + 6000233 + 0.326626976582762 -1.39216397200456 -3.371070451979964E-002 + 1.502022973415419E-002 4.320342344240131E-003 1.441047410021331E-003 + 6000234 + 0.327898875134393 -1.39098636673654 -3.929119504411583E-002 + 1.564916633966053E-002 4.902652198234570E-003 -1.318429564711382E-003 + 6000235 + 0.325772175388830 -1.39276944665725 -4.081849758657961E-002 + 1.459754212231259E-002 4.020943329837438E-003 -2.073659984154833E-003 + 6000236 + 0.325970862049605 -1.39573308065969 -3.729140206480187E-002 + 1.469578998479103E-002 2.555466459007209E-003 -3.295590242806717E-004 + 6000237 + 0.323334466127931 -1.38903625723277 -3.899099401712434E-002 + 1.339212791235488E-002 5.866954932989258E-003 -1.169984223183693E-003 + 6000238 + 0.321058043138336 -1.39159433514890 -3.730357920900053E-002 + 1.226646757888280E-002 4.602020053452453E-003 -3.355804571435673E-004 + 6000239 + 0.321631711087942 -1.39236740484015 -3.509861335813708E-002 + 1.255013860944605E-002 4.219747562605665E-003 7.547452891357343E-004 + 6000240 + 0.321306936119327 -1.39062854909901 -3.740848791004924E-002 + 1.238954178664758E-002 5.079588173238877E-003 -3.874563887900095E-004 + 6000241 + 0.322359878582366 -1.39481226441940 -3.804012263745659E-002 + 1.291020757217519E-002 3.010797617075125E-003 -6.997912057971096E-004 + 6000242 + 0.323817379375385 -1.38838475633633 -3.772374286905647E-002 + 1.363092197200557E-002 6.189113301486290E-003 -5.433456957891204E-004 + 6000243 + 0.327623726544369 -1.39464240098381 -3.788413992758530E-002 + 1.551310908886966E-002 3.094792785121600E-003 -6.226598686018013E-004 + 6000244 + 0.322826324940355 -1.38886200520861 -3.642775302359531E-002 + 1.314085897803487E-002 5.953120198633154E-003 9.750372819152632E-005 + 6000245 + 0.322612433934003 -1.39373086883468 -3.975745217634516E-002 + 1.303509277261627E-002 3.545533085870808E-003 -1.548987400995745E-003 + 6000246 + 0.321924589943174 -1.39047906318775 -3.905523030161998E-002 + 1.269496323622005E-002 5.153506931519558E-003 -1.201748195765625E-003 + 6000247 + 0.325483402999438 -1.39432969787986 -3.410254622383869E-002 + 1.445474808727309E-002 3.249420234359466E-003 1.247286994996116E-003 + 6000248 + 0.324075789932299 -1.39390925319708 -3.387081182165206E-002 + 1.375870249214517E-002 3.457324434934150E-003 1.361876517960287E-003 + 6000249 + 0.325555351421543 -1.39061734515993 -3.362701221631959E-002 + 1.449032560743936E-002 5.085128369354139E-003 1.482432120453008E-003 + 6000250 + 0.326533132333845 -1.39533418173792 -3.601686458585766E-002 + 1.497382502421540E-002 2.752716572605950E-003 3.006824950365570E-004 + 6000251 + 0.328077463583129 -1.39411605620772 -3.770484336908437E-002 + 1.573747590852224E-002 3.355063147389858E-003 -5.340001490527292E-004 + 6000252 + 0.321651628711733 -1.39378675143995 -3.597380335900138E-002 + 1.255998760462024E-002 3.517899894512744E-003 3.219756884387962E-004 + 6000253 + 0.324493522176533 -1.38860008548191 -3.882967355197441E-002 + 1.396526542860112E-002 6.082635955702965E-003 -1.090213438304715E-003 + 6000254 + 0.324521711140188 -1.39586389310030 -3.745606427432058E-002 + 1.397920448930016E-002 2.490781479023323E-003 -4.109822564850102E-004 + 6000255 + 0.323857759005521 -1.39310641830337 -4.070010820644245E-002 + 1.365088915215276E-002 3.854315415219456E-003 -2.015118039299293E-003 + 6000256 + 0.324036458712472 -1.39214287853430 -4.089940117322854E-002 + 1.373925373669473E-002 4.330772779566417E-003 -2.113665711887594E-003 + 6000257 + 0.322090781783657 -1.39468966511415 -3.693872754859708E-002 + 1.277714285021657E-002 3.071421312874804E-003 -1.551662531073200E-004 + 6000258 + 0.328092592928615 -1.39374096057231 -3.586344633457473E-002 + 1.574495716493310E-002 3.540542858309311E-003 3.765457421963371E-004 + 6000259 + 0.324548049062538 -1.39233312653552 -3.331554270885758E-002 + 1.399222823514643E-002 4.236697719935751E-003 1.636449572938168E-003 + 6000260 + 0.324887975354043 -1.39268402922437 -4.094045138727274E-002 + 1.416031718188453E-002 4.063181093387097E-003 -2.133964486694073E-003 + 6000261 + 0.321944774836254 -1.39352845277138 -3.514035858345473E-002 + 1.270494439243692E-002 3.645625087381753E-003 7.341028406686787E-004 + 6000262 + 0.327218902932329 -1.38903404006224 -3.771600840082730E-002 + 1.531292929618219E-002 5.868051293781019E-003 -5.395211060156568E-004 + 6000263 + 0.325421465012036 -1.39138633273825 -4.089966505983452E-002 + 1.442412059147267E-002 4.704874428056443E-003 -2.113796200239817E-003 + 6000264 + 0.324143614425125 -1.39411979216184 -3.398376362544507E-002 + 1.379224078514325E-002 3.353215768678313E-003 1.306023380953205E-003 + 6000265 + 0.324251383572956 -1.39481065425279 -3.453953133960193E-002 + 1.384553116897805E-002 3.011593822657939E-003 1.031203774355344E-003 + 6000266 + 0.323082300628443 -1.39404098553295 -3.438431714855358E-002 + 1.326743548852136E-002 3.392184579203928E-003 1.107955089402535E-003 + 6000267 + 0.325299655987546 -1.39040136408079 -4.060268324575995E-002 + 1.436388767880525E-002 5.191928087451827E-003 -1.966942715894272E-003 + 6000268 + 0.328389849750416 -1.39205995084950 -3.873156679977102E-002 + 1.589194663687427E-002 4.371779396419775E-003 -1.041700978227703E-003 + 6000269 + 0.325793463819930 -1.38928905744628 -3.967831695053032E-002 + 1.460806896313280E-002 5.741948651668003E-003 -1.509856103742432E-003 + 6000270 + 0.326601355314686 -1.39233804164542 -4.057608377821496E-002 + 1.500756036413905E-002 4.234267264663880E-003 -1.953789639491635E-003 + 6000271 + 0.321827451960849 -1.39213071131237 -3.479524070321862E-002 + 1.264692981972560E-002 4.336789306004781E-003 9.047589577125814E-004 + 6000272 + 0.321663381504649 -1.39142897193463 -3.917019157646589E-002 + 1.256579920152172E-002 4.683789923007624E-003 -1.258594988989431E-003 + 6000273 + 0.324648930551282 -1.39538104296497 -3.904801998382525E-002 + 1.404211276485815E-002 2.729544330582880E-003 -1.198182791282209E-003 + 6000274 + 0.321223963022878 -1.39253230185508 -3.603509083528145E-002 + 1.234851271435078E-002 4.138208222307412E-003 2.916698615769088E-004 + 6000275 + 0.321349228011334 -1.39221979031449 -3.561402108543606E-002 + 1.241045455438754E-002 4.292740946056133E-003 4.998831493500657E-004 + 6000276 + 0.324379742702741 -1.38833977818834 -3.620042042538130E-002 + 1.390900301999063E-002 6.211354386424082E-003 2.099166187150847E-004 + 6000277 + 0.321457479810322 -1.39117525342277 -3.868176372657838E-002 + 1.246398360268161E-002 4.809250290424065E-003 -1.017074033132608E-003 + 6000278 + 0.325405274076550 -1.39361725382059 -3.366648162894059E-002 + 1.441611439318668E-002 3.601714171383542E-003 1.462915030537829E-003 + 6000279 + 0.321794234123648 -1.39021314640619 -3.573859188553418E-002 + 1.263050404917571E-002 5.284999178073548E-003 4.382845760531595E-004 + 6000280 + 0.326444232064343 -1.39553980455841 -3.753707767299938E-002 + 1.492986504512957E-002 2.651038873102111E-003 -4.510422847790933E-004 + 6000281 + 0.322277249464589 -1.39203923743843 -3.998141216047151E-002 + 1.286934859267273E-002 4.382021897621130E-003 -1.659732579536102E-003 + 6000282 + 0.328312462659902 -1.39188336185765 -3.539912210963739E-002 + 1.585367976884838E-002 4.459100260935347E-003 6.061477820066648E-004 + 6000283 + 0.325347229643059 -1.39100999963996 -4.082441159349795E-002 + 1.438741220705580E-002 4.890966047611239E-003 -2.076584380468826E-003 + 6000284 + 0.325755425615411 -1.38888217640230 -3.916570950307285E-002 + 1.458925958623789E-002 5.943145816580859E-003 -1.256378664407702E-003 + 6000285 + 0.325383104362951 -1.39054381044671 -4.065574588764872E-002 + 1.440515177010779E-002 5.121490288988412E-003 -1.993181473557695E-003 + 6000286 + 0.325529431096013 -1.39548980174177 -3.555264887114687E-002 + 1.447750835756422E-002 2.675764588624088E-003 5.302308780096619E-004 + 6000287 + 0.321386839414644 -1.39364212905941 -3.740484666220886E-002 + 1.242905288386607E-002 3.589413702731343E-003 -3.856558410548179E-004 + 6000288 + 0.325609144161828 -1.39058365051225 -3.365164289696592E-002 + 1.451692538887061E-002 5.101789916224030E-003 1.470252582503900E-003 + 6000289 + 0.325746594379357 -1.39018221747378 -3.388130434928219E-002 + 1.458489265963108E-002 5.300293116209953E-003 1.356688105171851E-003 + 6000290 + 0.328670099790430 -1.39148012214160 -3.754968804313675E-002 + 1.603052648558429E-002 4.658496838509829E-003 -4.572779420004985E-004 + 6000291 + 0.327050766991314 -1.39494919091028 -3.584795596588932E-002 + 1.522978835080585E-002 2.943089322051995E-003 3.842055196892361E-004 + 6000292 + 0.324537148578145 -1.39587795189979 -3.718291882528427E-002 + 1.398683809326467E-002 2.483829593105807E-003 -2.759155317796403E-004 + 6000293 + 0.325510509540737 -1.38957080345959 -3.425440936567027E-002 + 1.446815190477895E-002 5.602629064425047E-003 1.172192728395570E-003 + 6000294 + 0.325364297568245 -1.39459990856101 -3.430333353770083E-002 + 1.439585206486989E-002 3.115804712620871E-003 1.148000388020130E-003 + 6000295 + 0.323644249841663 -1.39514928477911 -3.905298580505075E-002 + 1.354531176267529E-002 2.844145614251118E-003 -1.200638322614567E-003 + 6000296 + 0.322252155329533 -1.39173521301220 -3.994150826731383E-002 + 1.285693988279573E-002 4.532357858282674E-003 -1.640000644880718E-003 + 6000297 + 0.326793296706268 -1.39091455480311 -4.029063745626598E-002 + 1.510247278236859E-002 4.938162226599772E-003 -1.812640299750235E-003 + 6000298 + 0.328394129335298 -1.39145794911307 -3.861054959678920E-002 + 1.589406283363007E-002 4.669461100775097E-003 -9.818596105701817E-004 + 6000299 + 0.325717727002889 -1.39559209734865 -3.590486599481333E-002 + 1.457061813298568E-002 2.625180796653218E-003 3.560642812509859E-004 + 6000300 + 0.321281416206461 -1.39340867720813 -3.702359920936095E-002 + 1.237692253541057E-002 3.704852481010806E-003 -1.971341397410737E-004 + 6000301 + 0.326826923545244 -1.39251344124293 -4.043303848440993E-002 + 1.511910079875579E-002 4.147534539539541E-003 -1.883055679299622E-003 + 6000302 + 0.325221626988922 -1.39387263555342 -3.377034941755341E-002 + 1.432530339591360E-002 3.475431363727050E-003 1.411553815991447E-003 + 6000303 + 0.322406645286088 -1.39162516754770 -4.006538258520391E-002 + 1.293333307369588E-002 4.586773849880030E-003 -1.701254817159826E-003 + 6000304 + 0.327187048763822 -1.39463651865951 -3.547316576113219E-002 + 1.529717784133024E-002 3.097701514809212E-003 5.695341992875920E-004 + 6000305 + 0.327480662501029 -1.39487649003504 -3.693682090078323E-002 + 1.544236585728670E-002 2.979038920099771E-003 -1.542234415895329E-004 + 6000306 + 0.328551533470727 -1.39151193343580 -3.818397457052056E-002 + 1.597189704651028E-002 4.642766584423178E-003 -7.709240381764704E-004 + 6000307 + 0.322630404860865 -1.38945264563622 -3.538899387690294E-002 + 1.304397915252749E-002 5.661056507596596E-003 6.111560559036752E-004 + 6000308 + 0.328374974048021 -1.39213603553404 -3.876134862395392E-002 + 1.588459080353605E-002 4.334156550504728E-003 -1.056427686881741E-003 + 6000309 + 0.324324666115425 -1.38831319086554 -3.634772132746308E-002 + 1.388176839359332E-002 6.224501457414479E-003 1.370783178738140E-004 + 6000310 + 0.327780384205777 -1.39040540380995 -3.520423097607128E-002 + 1.559057418045742E-002 5.189930496106032E-003 7.025188076919263E-004 + 6000311 + 0.326421010606110 -1.39541002836368 -3.608433683151303E-002 + 1.491838234857564E-002 2.715211443533721E-003 2.673183834932795E-004 + 6000312 + 0.324795021294935 -1.39575370968382 -3.816709106533256E-002 + 1.411435265874821E-002 2.545265686000613E-003 -7.625753735535234E-004 + 6000313 + 0.325813555169143 -1.39541693890006 -3.874684061627268E-002 + 1.461800386317468E-002 2.711794276901274E-003 -1.049253673599004E-003 + 6000314 + 0.322857896534282 -1.38930726976145 -3.533330948140715E-002 + 1.315647070358410E-002 5.732942908509505E-003 6.386912352132688E-004 + 6000315 + 0.326219515200231 -1.39471810186507 -3.471690125726086E-002 + 1.481874559968870E-002 3.057359724732503E-003 9.434967526056666E-004 + 6000316 + 0.321329945274699 -1.39230542391461 -3.567221135978340E-002 + 1.240091950231233E-002 4.250396290733629E-003 4.711088468913174E-004 + 6000317 + 0.324943199595914 -1.39100652024179 -4.085233942880207E-002 + 1.418762482145983E-002 4.892686562876645E-003 -2.090394316735189E-003 + 6000318 + 0.326617044426315 -1.39248593326675 -4.055213917054404E-002 + 1.501531841732557E-002 4.161136861160018E-003 -1.941949355335786E-003 + 6000319 + 0.322850100717639 -1.38930020663682 -3.535283888796843E-002 + 1.315261577785102E-002 5.736435527964740E-003 6.290342082008193E-004 + 6000320 + 0.327167129230231 -1.39467550377465 -3.550717986735132E-002 + 1.528732790178609E-002 3.078423903439586E-003 5.527146844942592E-004 + 6000321 + 0.326175366886518 -1.39179718675643 -3.352228413895547E-002 + 1.479691485656030E-002 4.501712681218141E-003 1.534218736134025E-003 + 6000322 + 0.323968917553573 -1.39007740521288 -3.395618173788408E-002 + 1.370585554848606E-002 5.352121359507512E-003 1.319662250746563E-003 + 6000323 + 0.323157023880416 -1.38871984144043 -3.622802058959314E-002 + 1.330438512447146E-002 6.023418256348253E-003 1.962687113654486E-004 + 6000324 + 0.322467171269080 -1.39052754228481 -3.454951097363686E-002 + 1.296326235244168E-002 5.129534674689875E-003 1.026268980502423E-003 + 6000325 + 0.328300367708115 -1.39034678218645 -3.766648081609209E-002 + 1.584769897901974E-002 5.218918094875141E-003 -5.150303862311877E-004 + 6000326 + 0.327360300143949 -1.39032093899113 -3.953640998478093E-002 + 1.538284830205775E-002 5.231697204908905E-003 -1.439685031354920E-003 + 6000327 + 0.326573187554760 -1.39309481142322 -4.043854310486605E-002 + 1.499363178839726E-002 3.860054860233954E-003 -1.885777639553315E-003 + 6000328 + 0.325274808255286 -1.39252110692220 -3.334038631599733E-002 + 1.435160081177285E-002 4.143743964976811E-003 1.624164745722223E-003 + 6000329 + 0.322448036910301 -1.39211775354247 -4.012807498610159E-002 + 1.295380067120641E-002 4.343196747700902E-003 -1.732255360214964E-003 + 6000330 + 0.326810673793556 -1.38893153269665 -3.834409546944328E-002 + 1.511106551665390E-002 5.918739797571082E-003 -8.501016538046097E-004 + 6000331 + 0.326775470290585 -1.38869840010547 -3.741564841183009E-002 + 1.509365786127752E-002 6.034020706053717E-003 -3.909971599293355E-004 + 6000332 + 0.326271250871605 -1.39390875709412 -4.021126897748652E-002 + 1.484432818840653E-002 3.457569751128591E-003 -1.773393662065399E-003 + 6000333 + 0.321666944784208 -1.39307414833960 -3.528651321326542E-002 + 1.256756119499773E-002 3.870272475198349E-003 6.618313559388317E-004 + 6000334 + 0.325766832048710 -1.39507466523029 -3.493921370497112E-002 + 1.459489991300033E-002 2.881043970394314E-003 8.335662585066534E-004 + 6000335 + 0.321552673105825 -1.39151724509758 -3.527455850230447E-002 + 1.251105539788419E-002 4.640140039617423E-003 6.677427985786109E-004 + 6000336 + 0.326440517111881 -1.39496814513121 -3.518250039097470E-002 + 1.492802805145723E-002 2.933716716541827E-003 7.132642876744128E-004 + 6000337 + 0.328149349552533 -1.39379829625323 -3.817632697865712E-002 + 1.577302254667368E-002 3.512191140723280E-003 -7.671424075890959E-004 + 6000338 + 0.324135308775169 -1.38971487157048 -3.416375884586458E-002 + 1.378813375374246E-002 5.531389335037621E-003 1.217018182548996E-003 + 6000339 + 0.324286880142571 -1.38901868791716 -3.482913889326955E-002 + 1.386308374183996E-002 5.875642721572888E-003 8.879967618943014E-004 + 6000340 + 0.328695490885568 -1.39186332148561 -3.764840312220126E-002 + 1.604308203819967E-002 4.469009953454875E-003 -5.060912114703549E-004 + 6000341 + 0.326127505166542 -1.39330528092113 -3.372351078671609E-002 + 1.477324788433473E-002 3.755980544394189E-003 1.434714884496157E-003 + 6000342 + 0.323411133267291 -1.38892035601416 -3.885662321348482E-002 + 1.343003877428764E-002 5.924266515672206E-003 -1.103539680879768E-003 + 6000343 + 0.323443535356482 -1.38878773000395 -3.566309895402712E-002 + 1.344606116849591E-002 5.989848281256044E-003 4.756148081073266E-004 + 6000344 + 0.326161896812670 -1.39219123381745 -4.077210626406107E-002 + 1.479025408749908E-002 4.306861747038268E-003 -2.050720103554825E-003 + 6000345 + 0.327496529611136 -1.39347379619798 -3.470387250659136E-002 + 1.545021192830949E-002 3.672652022584027E-003 9.499392933331112E-004 + 6000346 + 0.326907660359911 -1.39442320185806 -3.939126722457996E-002 + 1.515902406000240E-002 3.203183783680173E-003 -1.367913902440976E-003 + 6000347 + 0.327627930876646 -1.39427354730877 -3.562233577430311E-002 + 1.551518807423184E-002 3.277185931186497E-003 4.957716483304529E-004 + 6000348 + 0.325119576127018 -1.39376660753008 -3.370627018966749E-002 + 1.427484062701359E-002 3.527860785085645E-003 1.443240126207252E-003 + 6000349 + 0.321582888441476 -1.39198483014723 -3.515228871706261E-002 + 1.252599647208718E-002 4.408925566155315E-003 7.282035511970817E-004 + 6000350 + 0.328649590429548 -1.39230144150209 -3.789609952169283E-002 + 1.602038488443426E-002 4.252365539779506E-003 -6.285737258914244E-004 + 6000351 + 0.325091488310059 -1.39588396719579 -3.697135848001232E-002 + 1.426095158198591E-002 2.480855110709554E-003 -1.713018066956808E-004 + 6000352 + 0.323629747250521 -1.38958719080047 -3.983511560394539E-002 + 1.353814042779786E-002 5.594525746332192E-003 -1.587390913967920E-003 + 6000353 + 0.321678639713746 -1.38999980605140 -3.649592552125010E-002 + 1.257334417924268E-002 5.390493093753963E-003 6.379335151967012E-005 + 6000354 + 0.321979789899611 -1.39044531600697 -3.911106486620910E-002 + 1.272225886697657E-002 5.170194455300030E-003 -1.229357631657783E-003 + 6000355 + 0.325389900698438 -1.39279507601704 -4.088553189135913E-002 + 1.440851246594727E-002 4.008269958578387E-003 -2.106807539867057E-003 + 6000356 + 0.325843612080240 -1.38834809342622 -3.662007773836583E-002 + 1.463286659858156E-002 6.207242613925633E-003 2.401761837715437E-006 + 6000357 + 0.321369009956186 -1.39207775902686 -3.555942305231640E-002 + 1.242023645816423E-002 4.362973493929667E-003 5.268811371797942E-004 + 6000358 + 0.325504520827021 -1.39403046328794 -4.038166123959698E-002 + 1.446519056696509E-002 3.397387686834489E-003 -1.857650327660951E-003 + 6000359 + 0.321990553923499 -1.38951059135313 -3.685720838148137E-002 + 1.272758153098724E-002 5.632403135476837E-003 -1.148561291719788E-004 + 6000360 + 0.324610041911226 -1.38980909225812 -3.401440587559276E-002 + 1.402288285910944E-002 5.484798481249997E-003 1.290871203314320E-003 + 6000361 + 0.328509070048447 -1.39256561942476 -3.832567067658864E-002 + 1.595089945937337E-002 4.121733135395790E-003 -8.409908433077421E-004 + 6000362 + 0.324576175323637 -1.39320821636194 -4.080080161764887E-002 + 1.400613629028089E-002 3.803977654144154E-003 -2.064909567216185E-003 + 6000363 + 0.322653205268223 -1.39479071554500 -3.866475003855334E-002 + 1.305525364512712E-002 3.021453243581057E-003 -1.008660994860108E-003 + 6000364 + 0.325278470241920 -1.39300006526267 -4.085198724804461E-002 + 1.435341161456094E-002 3.906905553258331E-003 -2.090220168121022E-003 + 6000365 + 0.327882444238280 -1.39058964691285 -3.907161492253900E-002 + 1.564104148409440E-002 5.098824777354270E-003 -1.209850168875113E-003 + 6000366 + 0.327925534924085 -1.39112239948916 -3.496856334495498E-002 + 1.566234924454749E-002 4.835385844669870E-003 8.190532590849669E-004 + 6000367 + 0.327291588883165 -1.38906313116244 -3.674455074075982E-002 + 1.534887151431583E-002 5.853666138781197E-003 -5.914845181920610E-005 + 6000368 + 0.326403910822977 -1.39105116007004 -3.374528293181369E-002 + 1.490992673743835E-002 4.870612772465192E-003 1.423948853656114E-003 + 6000369 + 0.328034117962957 -1.39191633393057 -3.935146673723238E-002 + 1.571604208647744E-002 4.442796017492191E-003 -1.348233100558018E-003 + 6000370 + 0.323887294249176 -1.39308506929055 -3.356970859695257E-002 + 1.366549393007587E-002 3.864872212877050E-003 1.510767984034018E-003 + 6000371 + 0.324195814984085 -1.39153924083899 -4.089966051523319E-002 + 1.381805325447577E-002 4.629263443433466E-003 -2.113793952996017E-003 + 6000372 + 0.323696388851772 -1.38910872113090 -3.494303980046445E-002 + 1.357109379693476E-002 5.831122516909043E-003 8.316743061108993E-004 + 6000373 + 0.328595123870093 -1.39262055299375 -3.630079724885421E-002 + 1.599345190855066E-002 4.094569229621872E-003 1.602816391441369E-004 + 6000374 + 0.326120765455020 -1.39220930951942 -3.349891101613327E-002 + 1.476991518827874E-002 4.297923557253797E-003 1.545776428773127E-003 + 6000375 + 0.325484317791356 -1.39317547361155 -4.077391072541323E-002 + 1.445520043948558E-002 3.820168500697788E-003 -2.051612385251454E-003 + 6000376 + 0.325938226310969 -1.39515654081679 -3.513803924942247E-002 + 1.467965205409683E-002 2.840557601902009E-003 7.352497199315036E-004 + 6000377 + 0.321076313569735 -1.39158800158611 -3.754184390774611E-002 + 1.227550205973084E-002 4.605151914460645E-003 -4.533991233011673E-004 + 6000378 + 0.322061849593275 -1.38999211007948 -3.550101842230587E-002 + 1.276283627396840E-002 5.394298647622158E-003 5.557614356104747E-004 + 6000379 + 0.328448014596534 -1.39176878090750 -3.857068065062026E-002 + 1.592070836541795E-002 4.515758988744939E-003 -9.621449567273610E-004 + 6000380 + 0.327573537859497 -1.39067587171533 -3.953481342703309E-002 + 1.548829146402245E-002 5.056187780465731E-003 -1.438895555174500E-003 + 6000381 + 0.325341975557433 -1.39322163371681 -4.078006057193613E-002 + 1.438481413288207E-002 3.797342953904623E-003 -2.054653401055379E-003 + 6000382 + 0.326945407039128 -1.39080439717213 -3.413440781683218E-002 + 1.517768928158405E-002 4.992633682997575E-003 1.231531868836375E-003 + 6000383 + 0.326864905896984 -1.38950094451344 -3.503286762493725E-002 + 1.513788255720786E-002 5.637173367032906E-003 7.872556636559166E-004 + 6000384 + 0.321872038730314 -1.39444412969432 -3.725713409567794E-002 + 1.266897737328358E-002 3.192835252123807E-003 -3.126139777195664E-004 + 6000385 + 0.325656830678744 -1.38929012479731 -3.972454209455978E-002 + 1.454050572555626E-002 5.741420861043140E-003 -1.532713811330532E-003 + 6000386 + 0.321381796697484 -1.39297273615133 -3.582742182750125E-002 + 1.242655932853529E-002 3.920419428636525E-003 3.943593729806181E-004 + 6000387 + 0.325031722246969 -1.38871352884892 -3.520767103604314E-002 + 1.423139807333841E-002 6.026539747343834E-003 7.008177446325606E-004 + 6000388 + 0.321818492658012 -1.39053648944325 -3.893028585473215E-002 + 1.264249956582932E-002 5.125110426032039E-003 -1.139964859192371E-003 + 6000389 + 0.328229706719982 -1.39063530392342 -3.842585640223736E-002 + 1.581275807751381E-002 5.076248004065306E-003 -8.905313275931080E-004 + 6000390 + 0.327031282654416 -1.39046595458131 -3.436136073991577E-002 + 1.522015361013151E-002 5.159988959845626E-003 1.119306722521987E-003 + 6000391 + 0.326576412217448 -1.39384168926042 -3.418873736070407E-002 + 1.499522634041712E-002 3.490733886441843E-003 1.204666645302896E-003 + 6000392 + 0.322057300552805 -1.39399573472489 -3.539654307306134E-002 + 1.276058683507447E-002 3.414560490852700E-003 6.074230806596391E-004 + 6000393 + 0.327117831128816 -1.39014808504543 -3.464203976347065E-002 + 1.526295065839501E-002 5.317171139694695E-003 9.805147472618755E-004 + 6000394 + 0.323831937101924 -1.38978966102672 -4.008662335125708E-002 + 1.363812057058993E-002 5.494406961976909E-003 -1.711758088260099E-003 + 6000395 + 0.323291350089045 -1.39111881770597 -4.052780949033966E-002 + 1.337080761514683E-002 4.837156987941165E-003 -1.929918658028077E-003 + 6000396 + 0.322050969058998 -1.39140096523577 -3.967628168503533E-002 + 1.275745599714884E-002 4.697638856237326E-003 -1.508849692523487E-003 + 6000397 + 0.322065436123730 -1.39461505964254 -3.660470991700977E-002 + 1.276460976469764E-002 3.108312708032919E-003 1.000094133619592E-005 + 6000398 + 0.325123972979440 -1.39448560316462 -4.010971214210824E-002 + 1.427701481097967E-002 3.172327182835760E-003 -1.723175182590891E-003 + 6000399 + 0.327666613656848 -1.39046306711613 -3.928421544112214E-002 + 1.553431618507078E-002 5.161416772263869E-003 -1.314978245571975E-003 diff --git a/examples/whm_gr_test/cb.swiftest.in b/examples/whm_gr_test/cb.swiftest.in new file mode 100644 index 000000000..46a8d0257 --- /dev/null +++ b/examples/whm_gr_test/cb.swiftest.in @@ -0,0 +1,4 @@ +1.0 +0.004650467260962157 +4.7535806948127355e-12 +-2.2473967953572827e-18 diff --git a/examples/whm_gr_test/config.swiftest.in b/examples/whm_gr_test/config.swiftest.in new file mode 100644 index 000000000..b0d8ac31c --- /dev/null +++ b/examples/whm_gr_test/config.swiftest.in @@ -0,0 +1,29 @@ +! Swiftest input file generated using init_cond.py +T0 0 +TSTOP 1000.0 +DT 0.0006844626967830253 +CB_IN cb.swiftest.in +PL_IN pl.swiftest.in +TP_IN tp.swiftest.in +IN_TYPE ASCII +ISTEP_OUT 1461 +ISTEP_DUMP 1461 +BIN_OUT bin.swiftest.dat +OUT_TYPE REAL8 +OUT_FORM EL +OUT_STAT REPLACE +CHK_CLOSE yes +CHK_RMIN 0.004650467260962157 +CHK_RMAX 1000.0 +CHK_EJECT 1000.0 +CHK_QMIN 0.004650467260962157 +CHK_QMIN_COORD HELIO +CHK_QMIN_RANGE 0.004650467260962157 1000.0 +ENC_OUT enc.swiftest.dat +EXTRA_FORCE no +BIG_DISCARD no +ROTATION no +GR yes +MU2KG 1.988409870698051e+30 +DU2M 149597870700.0 +TU2S 31557600.0 diff --git a/examples/whm_gr_test/param.swifter.in b/examples/whm_gr_test/param.swifter.in new file mode 100644 index 000000000..0582bd1f7 --- /dev/null +++ b/examples/whm_gr_test/param.swifter.in @@ -0,0 +1,27 @@ +! Swifter input file generated using init_cond.py +T0 0 +TSTOP 1000.0 +DT 0.0006844626967830253 +PL_IN pl.swifter.in +TP_IN tp.swifter.in +IN_TYPE ASCII +ISTEP_OUT 1461 +ISTEP_DUMP 1461 +BIN_OUT bin.swifter.dat +OUT_TYPE REAL8 +OUT_FORM EL +OUT_STAT NEW +J2 4.7535806948127355e-12 +J4 -2.2473967953572827e-18 +CHK_CLOSE yes +CHK_RMIN 0.004650467260962157 +CHK_RMAX 1000.0 +CHK_EJECT 1000.0 +CHK_QMIN 0.004650467260962157 +CHK_QMIN_COORD HELIO +CHK_QMIN_RANGE 0.004650467260962157 1000.0 +ENC_OUT enc.swifter.dat +EXTRA_FORCE no +BIG_DISCARD no +RHILL_PRESENT yes +C 63241.07708426628 diff --git a/examples/whm_gr_test/pl.swifter.in b/examples/whm_gr_test/pl.swifter.in new file mode 100644 index 000000000..60d090453 --- /dev/null +++ b/examples/whm_gr_test/pl.swifter.in @@ -0,0 +1,40 @@ +10 ! Planet input file generated using init_cond.py using JPL Horizons data for the major planets (and Pluto) for epoch 2021-01-28 +1 39.476926408897626 +0.0 0.0 0.0 +0.0 0.0 0.0 +2 6.553709809565314e-06 0.001475122968086379 +1.6306381826061646e-05 +0.1030256860922895 0.2897796047098886 0.01422904600374035 +-11.74004209950937 3.8343124110162736 1.3902496665973592 +3 9.663313399581537e-05 0.006759127649782299 +4.0453784346544176e-05 +0.06110218027254217 -0.7245466901305982 -0.01346904300924688 +7.311995449678243 0.5941125721336201 -0.4137913843379075 +4 0.00012002693582795245 0.010044756567546644 +4.25875607065041e-05 +-0.6061796342297583 0.7761214554702035 -3.4750047790977e-05 +-5.054824314301841 -3.891667468503358 0.00019720338148272726 +5 1.2739802010675942e-05 0.0072464490746299084 +2.2657408050928896e-05 +0.2751944175855944 1.51937688993241 0.02508924593104206 +-4.835983593209577 1.344855094041679 0.14681413000004515 +6 0.037692251088985676 0.3552852357486061 +0.0004673261703049093 +3.200135438345358 -3.953498213518368 -0.05517737289975112 +2.111393749129838 1.8660266890185446 -0.05498941067210089 +7 0.011285899820091273 0.4376306456694341 +0.00038925687730393614 +5.607382165725712 -8.258649105608766 -0.07958445228024298 +1.5748468603228847 1.1414574661825514 -0.08250331331320372 +8 0.001723658947826773 0.4690969274244374 +0.00016953449859497232 +15.28225422201768 12.53905314208462 -0.1514143582550325 +-0.9198472198098231 1.0454390993472462 0.01574538863031621 +9 0.0020336100526728304 0.7807192056765467 +0.00016458790412449367 +29.47483071169769 -5.147686530859088 -0.5733441819169969 +0.19191677740340274 1.1385110364087574 -0.027844325148353527 +10 2.924216771029454e-07 0.0538346817277698 +7.943294877391593e-06 +14.14000920780611 -31.14141812522779 -0.7565722591093476 +1.073396108697069 0.23003123192799815 -0.33424529561177047 diff --git a/examples/whm_gr_test/pl.swiftest.in b/examples/whm_gr_test/pl.swiftest.in new file mode 100644 index 000000000..032262e70 --- /dev/null +++ b/examples/whm_gr_test/pl.swiftest.in @@ -0,0 +1,37 @@ +9 +2 1.6601367952719304e-07 +1.6306381826061646e-05 +0.1030256860922895 0.2897796047098886 0.01422904600374035 +-11.74004209950937 3.8343124110162736 1.3902496665973592 +3 2.4478383396645447e-06 +4.0453784346544176e-05 +0.06110218027254217 -0.7245466901305982 -0.01346904300924688 +7.311995449678243 0.5941125721336201 -0.4137913843379075 +4 3.0404326462685257e-06 +4.25875607065041e-05 +-0.6061796342297583 0.7761214554702035 -3.4750047790977e-05 +-5.054824314301841 -3.891667468503358 0.00019720338148272726 +5 3.2271514450538743e-07 +2.2657408050928896e-05 +0.2751944175855944 1.51937688993241 0.02508924593104206 +-4.835983593209577 1.344855094041679 0.14681413000004515 +6 0.0009547919384243222 +0.0004673261703049093 +3.200135438345358 -3.953498213518368 -0.05517737289975112 +2.111393749129838 1.8660266890185446 -0.05498941067210089 +7 0.0002858859806661029 +0.00038925687730393614 +5.607382165725712 -8.258649105608766 -0.07958445228024298 +1.5748468603228847 1.1414574661825514 -0.08250331331320372 +8 4.3662440433515637e-05 +0.00016953449859497232 +15.28225422201768 12.53905314208462 -0.1514143582550325 +-0.9198472198098231 1.0454390993472462 0.01574538863031621 +9 5.151389020535497e-05 +0.00016458790412449367 +29.47483071169769 -5.147686530859088 -0.5733441819169969 +0.19191677740340274 1.1385110364087574 -0.027844325148353527 +10 7.407407407407407e-09 +7.943294877391593e-06 +14.14000920780611 -31.14141812522779 -0.7565722591093476 +1.073396108697069 0.23003123192799815 -0.33424529561177047 diff --git a/examples/whm_gr_test/tp.swifter.in b/examples/whm_gr_test/tp.swifter.in new file mode 100644 index 000000000..573541ac9 --- /dev/null +++ b/examples/whm_gr_test/tp.swifter.in @@ -0,0 +1 @@ +0 diff --git a/examples/whm_gr_test/tp.swiftest.in b/examples/whm_gr_test/tp.swiftest.in new file mode 100644 index 000000000..573541ac9 --- /dev/null +++ b/examples/whm_gr_test/tp.swiftest.in @@ -0,0 +1 @@ +0 diff --git a/examples/whm_swifter_comparison/cb.swiftest.in b/examples/whm_swifter_comparison/cb.swiftest.in new file mode 100644 index 000000000..2386b53c8 Binary files /dev/null and b/examples/whm_swifter_comparison/cb.swiftest.in differ diff --git a/examples/whm_swifter_comparison/config.swiftest.in b/examples/whm_swifter_comparison/config.swiftest.in new file mode 100644 index 000000000..c32a270f5 --- /dev/null +++ b/examples/whm_swifter_comparison/config.swiftest.in @@ -0,0 +1,29 @@ +! Swiftest input file generated using init_cond.py +T0 0 +TSTOP 1.0 +DT 0.0006844626967830253251 +CB_IN cb.swiftest.in +PL_IN pl.swiftest.in +TP_IN tp.swiftest.in +IN_TYPE REAL8 +ISTEP_OUT 1 +ISTEP_DUMP 1 +BIN_OUT bin.swiftest.dat +OUT_TYPE REAL8 +OUT_FORM XV +OUT_STAT REPLACE +CHK_CLOSE yes +CHK_RMIN 0.0046504672609621575315 +CHK_RMAX 1000.0 +CHK_EJECT 1000.0 +CHK_QMIN 0.0046504672609621575315 +CHK_QMIN_COORD HELIO +CHK_QMIN_RANGE 0.0046504672609621575315 1000.0 +ENC_OUT enc.swiftest.dat +EXTRA_FORCE no +BIG_DISCARD no +ROTATION no +GR no +MU2KG 1.988409870698050917e+30 +DU2M 149597870700.0 +TU2S 31557600.0 diff --git a/examples/whm_swifter_comparison/param.swifter.in b/examples/whm_swifter_comparison/param.swifter.in new file mode 100644 index 000000000..b092c2553 --- /dev/null +++ b/examples/whm_swifter_comparison/param.swifter.in @@ -0,0 +1,26 @@ +! Swifter input file generated using init_cond.py +T0 0 +TSTOP 1.0 +DT 0.0006844626967830253251 +PL_IN pl.swifter.in +TP_IN tp.swifter.in +IN_TYPE ASCII +ISTEP_OUT 1 +ISTEP_DUMP 1 +BIN_OUT bin.swifter.dat +OUT_TYPE REAL8 +OUT_FORM XV +OUT_STAT NEW +J2 4.7535806948127356533e-12 +J4 -2.2473967953572827815e-18 +CHK_CLOSE yes +CHK_RMIN 0.0046504672609621575315 +CHK_RMAX 1000.0 +CHK_EJECT 1000.0 +CHK_QMIN 0.0046504672609621575315 +CHK_QMIN_COORD HELIO +CHK_QMIN_RANGE 0.0046504672609621575315 1000.0 +ENC_OUT enc.swifter.dat +EXTRA_FORCE no +BIG_DISCARD no +RHILL_PRESENT yes diff --git a/examples/whm_swifter_comparison/pl.swifter.in b/examples/whm_swifter_comparison/pl.swifter.in new file mode 100644 index 000000000..d0d4e7ff9 --- /dev/null +++ b/examples/whm_swifter_comparison/pl.swifter.in @@ -0,0 +1,40 @@ +10 ! Planet input file generated using init_cond.py using JPL Horizons data for the major planets (and Pluto) for epoch 2021-01-28 +1 39.47692640889762629 +0.0 0.0 0.0 +0.0 0.0 0.0 +2 6.553709809565313959502e-06 0.0014751229680863789154 +1.6306381826061645943e-05 +0.1030256860922895 0.2897796047098886 0.01422904600374035 +-11.74004209950937 3.8343124110162736 1.3902496665973592 +3 9.6633133995815387361564e-05 0.006759127649782299051 +4.0453784346544178454e-05 +0.06110218027254217 -0.7245466901305982 -0.01346904300924688 +7.311995449678243 0.5941125721336201 -0.4137913843379075 +4 0.00012002693582795246295385 0.0100447565675466429165 +4.25875607065040958e-05 +-0.6061796342297583 0.7761214554702035 -3.4750047790977e-05 +-5.054824314301841 -3.891667468503358 0.00019720338148272726 +5 1.2739802010675942316241e-05 0.0072464490746299085006 +2.265740805092889601e-05 +0.2751944175855944 1.51937688993241 0.02508924593104206 +-4.835983593209577 1.344855094041679 0.14681413000004515 +6 0.037692251088985682938581 0.3552852357486060849 +0.00046732617030490929307 +3.200135438345358 -3.953498213518368 -0.05517737289975112 +2.111393749129838 1.8660266890185446 -0.05498941067210089 +7 0.011285899820091272946487 0.43763064566943408597 +0.00038925687730393611812 +5.607382165725712 -8.258649105608766 -0.07958445228024298 +1.5748468603228847 1.1414574661825514 -0.08250331331320372 +8 0.0017236589478267728883093 0.4690969274244374022 +0.00016953449859497231466 +15.28225422201768 12.53905314208462 -0.1514143582550325 +-0.9198472198098231 1.0454390993472462 0.01574538863031621 +9 0.0020336100526728304385693 0.7807192056765467829 +0.000164587904124493665 +29.47483071169769 -5.147686530859088 -0.5733441819169969 +0.19191677740340274 1.1385110364087574 -0.027844325148353527 +10 2.9242167710294538257026e-07 0.05383468172776979939 +7.943294877391593783e-06 +14.14000920780611 -31.14141812522779 -0.7565722591093476 +1.073396108697069 0.23003123192799815 -0.33424529561177047 diff --git a/examples/whm_swifter_comparison/pl.swiftest.in b/examples/whm_swifter_comparison/pl.swiftest.in new file mode 100644 index 000000000..7bdc4a619 Binary files /dev/null and b/examples/whm_swifter_comparison/pl.swiftest.in differ diff --git a/examples/whm_swifter_comparison/tp.swifter.in b/examples/whm_swifter_comparison/tp.swifter.in new file mode 100644 index 000000000..4043c1929 --- /dev/null +++ b/examples/whm_swifter_comparison/tp.swifter.in @@ -0,0 +1,13 @@ +4 +100001 +2.894380502059049 0.2060633316227693 -0.5267473116107563 +-0.36781773526648104 3.511709550412678 0.17870205470035164 +100002 +2.402157114026988 -2.063650689101573 1.221462187067896 +1.9929454155842872 1.7670860403833373 -1.3920313116840413 +100003 +-1.76242073921186 -2.766072818539065 0.6998178789690445 +2.514640194836968 -1.5368452873925367 0.24681696899054972 +100004 +-2.13688673996785 1.023684546777855 0.2293900047351895 +-1.3918193851434775 -3.7994616097035703 0.2829665907639622 diff --git a/examples/whm_swifter_comparison/tp.swiftest.in b/examples/whm_swifter_comparison/tp.swiftest.in new file mode 100644 index 000000000..01954c3d7 Binary files /dev/null and b/examples/whm_swifter_comparison/tp.swiftest.in differ