From 47ff8692d65467cc256847872dc9c596d2e5258c Mon Sep 17 00:00:00 2001
From: march3 <gang.wei@qq.com>
Date: Fri, 24 Nov 2023 16:19:11 +0800
Subject: [PATCH] =?UTF-8?q?Python=E8=B6=85=E4=BA=BA-=E5=AE=87=E5=AE=99?=
 =?UTF-8?q?=E6=A8=A1=E6=8B=9F=E5=99=A8?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 ...1\345\222\214\346\227\266\351\227\264.jpg" | Bin 0 -> 37470 bytes
 .../halley_comet_sim_02.py                    |   0
 .../halley_comet_sim_03.py                    |   0
 sim_lab/oumuamua_sim.py                       | 488 ++++++++++++++++++
 4 files changed, 488 insertions(+)
 create mode 100644 "docs/oumuamua/\350\277\220\350\241\214\350\275\250\350\277\271\345\222\214\346\227\266\351\227\264.jpg"
 rename {sim_scenes/solar_system => sim_lab}/halley_comet_sim_02.py (100%)
 rename {sim_scenes/solar_system => sim_lab}/halley_comet_sim_03.py (100%)
 create mode 100644 sim_lab/oumuamua_sim.py

diff --git "a/docs/oumuamua/\350\277\220\350\241\214\350\275\250\350\277\271\345\222\214\346\227\266\351\227\264.jpg" "b/docs/oumuamua/\350\277\220\350\241\214\350\275\250\350\277\271\345\222\214\346\227\266\351\227\264.jpg"
new file mode 100644
index 0000000000000000000000000000000000000000..66b41d21280d52b93fc622dd1b89f9fdb4a113bb
GIT binary patch
literal 37470
zcmcG$1z1#J*Dkz=5eDg&Mg*ipknWTg5R@({DQQO9Ap~J4Nu?V^TDrTtrMq*Gj^T{&
z`{lXLIsf;6|M`dO8HSn7=2`oB*1hg^ueG`VeZLG4D9I_x0U!_nfKXq+{Tv_zU_3y_
zK!1RNfsO%zU|?bqU|~Ich(&^jk3;Z?gq-{l2^kqB4HF$D6$3RH+2d!A8J;q;vawRo
zad2_4a51s4viv;>2m*m%VPX+uVG*-Xl2Nk!zkb}e0EC!m&S)cG&=UZS5CkR!-FE<V
zsClA;{_z0*<pDwiqvi|2e29gOI-rsOKm&onXb-^X=no#Cj`l%)4m==4C!&8Lg+Z)t
z0(s&@!tEE8fyp3U)=H`|dd$dU>g@jzi|i3O1trr{W)@a9UOs*SK_Ow8m$Gv53W`c^
z-fC)T>*(s4nOnTKw6eBwafN<#bNBEJ_#7DYB{&2g9TOWD|LuE1VrEu$PHx_h{GY$d
zD=Mq1YijG-+B-VCx_f&2#>OWmr>1|;%&x4it#52@ZSU+NPX3&ponKsDUH`QU1c3i(
z7V7z*hW$@=5u)rud+-4K0P@!^5SklmfC(R<)4#wVl2V75I1xYL_QNERj>;%&eaOJ0
zaZGCJJc>oe$h*RX_-oofEc>q+=KtTa?B5Lgce`c*954vAcwj<60=UwEv!tW_|1|Kj
z+T-g!B+Aauri2wy9vW9(bt>33Kau->l_Z*cYOGZ6k-0PK?v!G;E<iA&Z<LhAwJO6N
z$Qt2`v+UNFmOeP1J$2CVf^Q){pQ$M^hFp)Nj?Jy^&*zKR&R%|R*C}VWUsJ#;Fd|Co
zI9hLgXGBYmw}kF|w!!J4vg6N1)pFXZ4oJh!AG)zoD5Hz(nkt)#ICDGlgsqk}SeV(F
z6z-R)RGqYv>{W)cF8W4){!`7kCc0W^lt!lLl-KLYmz$#cY*y;&vLhZr+6PAwO{dJg
z>&~eh-uUWjfsYwRdxJRURNq;KJcb#LuhrjixZBDHDQtYxEKSv6`c+e39W7LazremF
z!Mf(WR=uY^*RsLEyOURMWn^3Lr5v09vU=6>D;dobU6h(OE)d&h^wRmrZec>$GPz8w
zX^XIDi}6_8AUx~wpYicf|D#U-v5tpl&E<-WF!QUN-$OzwDIqjv!#QYWkdO)m-nh#6
zq0@!AAq&Y(1`1p%J9}<-G1+%hop(+3e1pYlFud#wA|YEZ>!^Mnh56Vwz0Vc{VsyE?
z5K5zcVbSz#eeLK9%l=|k6jkx@&hEmtWgz|E_G1j3-&Hp(%)%QnUQ88S>a}sa=^*3R
z#%ZI-!J1QjNaH9r^x|?Vr82ki9>BZ@;&882xb>|$aZBskzE7T)vF@cVs93aT?qTy~
zd{JvAGN75-FmW`0Y~SZJOpQ28$BT0<;>e^d`HdN7UF30&&sP;sT%kL>I2SnY^Jj>i
zM`xt)0GdXi|LlU0$(IAm14E<RIXt5M@ZF26NT-N19OL7d!F(wg#pL3(xfyBORgk>q
z@_~Ie-j~uJAM9?WnCHgYI5ZpfTir_wzYJ_SDFS;vL-<aKz$~(Te8|6GgL~d5i2aX8
z!y2bxSngtyq(iaxaW=v}G;dZeXpDlxJ4ISWEX0gHPV@EEcjU$3h*M9YZj7ce{%rQQ
zN5=0`SK`?<j7^-vO6gorI-G`FG;Fxs!;1^)C!S1(_grxhna8-dZu|ae?N$Ot$4`C4
z7ylS|Xs9SrK_UdRq9_-#@gJ+~H+wjWc$Qw%NUfxBL+UE5))Urin**QYeCNV;eb#{Y
zQ1GgJE?0P((+j$a%cJErH*=<JsI*Q&>5(v|D?_VOw!R+GG|^9YP(5~(*nJ~gvAMb5
z+0fy2d%3ohC@Gab*X~xdb3RQUMr<HY6=;`d%LRc{>+Nw=*x{lV80#(;sKLY%(`&jE
ztlGBNyTJTvKd7<Rm*gf^b!k<^5o9c5+8o%p!AsCcYVhxiBC0P4we!j>gOqo>PWk+s
zYCMCQG)k6fMzJjgt5&DKD~W02s1c0svh&%G+P|2zzcYWwx-*g3U2xrJgAHS_#dfWZ
zHkqFf!kgJvJ+xWK`@(tXnma)t`SnAzEYYhUV{ASrKUMW<mjVIPX~W}Rt0&W4RXj7s
zYNloOiK|w|)&x(!#fZ<nJuBqgw^JLKxl40v9we|k(HA$N@;9S%kC5Dh@9^HRH{)H(
z8#FanrB3bM^pD$3y=Cdtxr#1@B-_cpxd#l{EJ%6Rm|WYo+*XIfQn_Q=X8AgE^Umd(
z1)j=~C{Uw+!^$Y@f>2OT!x!dN+*~+1Cw{4jh8eC$h4K(jM0+Ogm$>eQYq!%9S5;xr
z^*6;O%-Liz1}22;0~n?6N1oFk-!&CID=mEVAb5Op{<mFTvK@6mmzNf~H!U45%mkh4
zToJdo|D)-&=A+vf_Ou&p*?S;UFMs|Ly5c|OI>Eb>UXwg`6~c-Y1~8UA_D<*8tU5~(
zHe_B|l|E_Z9?pN3|BR!SQH^ogUs9B&t?s4cL}gy?_8+5}dBcRFgCdiwS20Z^gJGXv
zgBaJ+?9W6~sT(8ZRiqMe1BtXk2kuto1}hCB`V^2x6(Cr9+6j|;0dhg&)Vf%!5v%)U
z;&2;fXLhxg?RkEoxo$^u-S0&csw5h2i5CX;mANX}i+8DRaFx;Qd`6X88VVV9@FBh@
zFkGD0&LSV!U`XG!^h+`=q%JseRK=t{&uGwsX<q9azAZ)jKv{I?nHWU<YFaaW=RymK
zk1$=3Vyc;YADFA*r<kBH<mog#o(QZCG8s&s-h{Jj)RZ+y#<bj!gQLzMH4G2vWZrIr
z3yu*p%aNqqcdHy`%y%)Y%=ds4X`l9+OfF+I_wQ#Me9=DQ*7(SNgGzvOa{Hb9GAB`w
zv+}StnV7ichkIbHxyZS4VyBSd5?VwYvuQ&X73mQDet1Wcx{~BI3phW{-<!LdeACcK
zCQYSO{`LYH*RQSO{_$3DYy#&hF=w3agDilo92fzFr!`mAH!fw`oX(zNhw#kdzF0dg
zehV$m>)e0z@;I9z**f>tJ<uqc{3#aHuc0*@>?-6MLhqXJ{PZS{TW|3trs6k2QCdyE
zmDm&fWm`Jx*J$dpB*6du!2!6|;g#pi?4S*+kLT4(d{8d0g1&TwOG;zK1k`LL*4_gT
z;Y>B|9}JAi-&LDlpr_9eYSe^M6PI}{JIfnAch|b+5NNrmR-tw_kaaJ@S|wwC*b`+y
zX4TvBSo4hu36gr0sp*Bj-d*MFdb=Y<Jo5|WpZ;f#PEjq?EIq5fl&LJ^<8Spd@(lYT
z8Os%E27kwCE-1{yd#l-#H~fZWR3}6YEnx+Q!~AKE@zoz+&ni;8u@&R!bM>_XTm2MN
zUWsE9`@A<CUh)q3T~*U)+of^#qfx652=0MB@_?`t1;~hLjM{kN8iF>DuT4ko<2OVv
zk&4%PSJlSvil2Rwj3!>m;AGksoh^<tJ8lkQjj|h>ZQEKc*(V!H-Q^!g<d%|-oZWZy
z5SO#0H4$6;#8)SZAEp@&bfxJZN3ts2em<a1Ui&IMA<j@;;g$Ai=<Dx3jXNCT&=wl?
zlzU*ZOE8u&EZMt(mdBttZQZA7?jYet;Cqw3@P5zmP9gzMMgxLZz)kpN`P?Gcuex{-
ztc8PUSr71vF@+qBfXlF8vGbv9;;yy2q58si_NH?n`D9_Z+4Hc$$JQ?2c8bfo8nK`1
z=D#RZoXORRxBSGX53M7cuFlhpG10O;KYjcRd(Z0+<s4(#UFN1p@Pvoema|29Eje^9
z+EGv_9ZRd}xsIrCqIL{){VvE>TU`)y=mdga?eGB<fy^x2MY7d}lU*;V)|b~b)hhqt
z>nKF<WC)%}fw7~#QpJr@EI!6dKo*lcIGD!sbaBXX)@YSH8Wh%!DBMBvuWQ8wMTd%3
z#pE@|3kw~Xw5lsWW>|iTZJ*JDccp)pjJ(D;6&HvQxoM?MS>Ji1UsUn2_es5U#_d25
z^zlYN5$2Xy+3Pei8D72PT9udmKZRdJ*~bzP$Kiek)t=L!YaWY5BzQqR`36Jsy6yE8
zKYx!Wt})4c75SpxPV`+&7rR8K7DCS*evR#Jw;X7GH~;ZZ!Wwo2*Lw84+dXjh8ojkp
z2On?lH<fKw4Seh1sf_UcvycA#2R-D*M!8k@zzQFHPyK6S#f?4Zb4CunC01$#W`f(n
zN)S&y`FXAix52GM5LLdJzlasg(WDV~tH*-bwmmGUo!0q91+Ax+YG}^F>|UK(lvl1i
z-+-}ogk_=-8DY@P@;ZAM5y_uquh-yPuOE6#yIP87?8>+~Vf;B^xstk|A#%T{;dI($
z&1QVtlUZwmx}9%sxH$>5LK#u<q<;kOP6CA}E-2_&SL>*-TSQg7FK{UhlN5;f^?zig
zzNT>#Q4W&akeY9Z#^Y&aj&4LSEqt!qOv@RDoI@8nvOSz|>vg^=S;nC|ipToKPj|VI
zRSCa!oXfBOmC8zM>D=ZTQY}lxnHjxs)pJVfez{0@(8zYGzl8b0!YeTJfY6`ugL96V
zfwK@Z<Bc}1;wx;=ls1DU<!5p(#MC8ePcLrOKDiiTB(2)oa^i_6Kg)|+WMzOxZ$azJ
zi9<v)o7q`f&=t1vr(<){@pk9Bs;tic@FzjfGoaQ`-cS$Ut_Q4vhRx5VhMXN`1y;rW
z6n*=0YEA3H5m!Gy!00n)UhBYH@`GN9^^ASCw1GJHX~}B1_?WlP_ufraH7O$8;d>w=
z|7X*aiR>x9xjQ-IMktSRq@)`7n4ARkyY**#&UO>kjnO?&Zn&d8<;;yc;IEhi@;wR&
zI30YacNs)wj6bi_jXtg>`V?cyRS&=pdiJCnx8*48&r~)-fR^ky9=0{+a_N(+x5)xF
zsdRrIpj25~VOd^|-$-Ux#4ci9<@hl*AY;d{IImdO#)zn_VWQCfurOmnG)TU_<4<69
z3&->Ep1ykkA23^N=&H-V2i6nyV88Vo?ucT6Bdp6d6@10imUhOk6{{OsLf)mrrMWX0
zc&Cp~8a`E$4u@o=N#D?yKRH;A7*p&>V}<k2*BCf1fP}V%T?{%&Ks#Ub#)?Jpomk2m
z3EG`3-uIrz;k0VJ$06tJb3f(laxc(DE=+!!-ss7G8ATSa#;wR-EKvW>B0@}z?JVs!
zCj1SuEx!4h%@bu@+3`IP`mCZ6!xl2u-_J`t*DzYg>cC4_6CL=g0S5yN-b<1CIctMN
zo2FIvSgzF?&rDS3KNa!r?0JPINU=B7uv=%x(dR9Jz*t&sDw8>#bA%<rfxI;uEe^{L
zZMUlJUa*})``tmXA4!@s3GrE>-Wx?p8r9f(aXNRu-j$<gBR-O`_ki9;K|vYSCn}r5
zdB{XjoKH`LG)$U3HXW@l)S|oH_MG`>cWN0wn}_(e)A?*tj{o*WY_x<&LGF@IsW|ME
zUXr?qgNY>_FHc*Yu6+srNAHrC2mP>Fur}T0cClVv3Bige{)$ewoPcx8p#^qiudH(g
zfLA%;J{>B7^LV@NrRAS!4cl+U1>NbaBkFRbyrCaH(SJu~kDyE7?Ul{QJWu&KS9uuf
z+?rdBTe;@_PMH{n6L7*15T$+@Dec3Qk-BKnFjnkmtKVVp5ZKLzVVB6g8J{OmWpZ$L
zFf<f<S~_7OkjOVmXH~V}kGaEWd1-Fc-7NFf3XDw>H*@)=#m)0$xJluR;oOq7Ufl4^
zz8+6VHj~Sr?psLJ!@ZuK)}PnSew7@YvJ-r|o8?|2BBN&CH~N*X-sx&n%wiG7grz2I
z0n|C9^UduAIdBACc{YveECn-iNY`L!j}j7-d7y#eZi;OuI85QQ)vn*q|K)&kDul8y
z6MNHRLH2VH1B<S6L-fy<4-4@cjaJ}*%=$TpplRYg*1YHSS?u$jI91#{{h^~QO#Cb=
zec~GAHLvsXiz%%VDU1ZV)|piCJZMirG$iVh=>L}nb|R0WEb!cFkckAS&Ek!HjkXy!
zrY?z{%95Sw`-t~FTnR83O;ln5xyYi9_$NnV*^{KP3~FB0UEzRfiz*8jf3tU5d{vQm
z9F9ss{oos*t3aqY=RFW&AF^^s=*S44Uz(4+wsQ&s?w;s;nda*&QJ18V2!*YKRyc1U
z63AF(gGQ+Mb97hEhwM0X`<sH|L8}_<@u&@1-@0+jH~uy`0V?u^wnd)2Sfj2{4XHR3
zat|0!8(-tNg8EQX0!_5QR%aG&c2WMd9by|_oKM2s{z%Mw!4QU9DsT@xy4GfwM0cHH
zFJGA|UEdm**dM`HDKntkGZ=(&u^I=t4a`F(i`s5hDO|4X=~4ja)7ac-U*Q@QkR+)3
zK&|e9Q)O*dLuFVmui6cE+8qIEZs^57)xNC+HXY(5Rs{N-x7FTO;}6v;?Q*(bI%Hdh
z&!)`KUc#|}jU8F0zn4oQ<oi9)#=OFKaBlnMR%IaI9<U(opMJACLfnodlHIvWJ-AE?
zQ7qQ;C96vdwxlfjYzAT&e>N$(FB&H?khl@(x<vR;oT-%=+u}WNN@oH06@sgy7S0pj
zx53SNW$@{eTH|=?$-Ir}wK=-vyRtZy?10M2($;1)%+;cHMo<2v#{>qpdRis<vkRS2
z0BDXi=fo}%u@HaHf2Yv%csdd5g3eD*Eh@v%#(ubN_jT~++jRNSvb+6g&@P6Liiq!x
z3O5plv%5cwQEpdk^GlgyaKc3&%vNr%P<jstcX0`T{1+|(Eu-=^8Z0rBOavsTQfF{s
z-!r1b?%Cd-L2A_>QT|yM7TNEGdKm!E>3Vm!mq|LcY|Zu6GB>?1#(#TF86V)H)=#DB
z(t|(846PHEm6!YrPfwxzPXX;cut@vaYx+y*@bPy811@K`*p!MLpJHoDjaZLJ+zP}{
z03__`G*e`U^L##BfjMHAg>Hk~3hU<wf~X20od82<(4#BQl@m%;kmJt|P1OD>L(nEi
z9#biwJ_H#nNW~LltvpBw{J-Q;{}s<9n$z!xiaP4Q{G-P3Kw(W|V89&!O^mPos*M&y
z%nkHzVpg*e38aRe35sf~9Ovo#fV;10mK&5=K1~xDsQzw~HaEPDItc?u3w7CQze1o{
zy}q<?t-s8y5jZiP@A)pSY9N=;)AG(Q!ZW-D{t+BnJ;q0sl+RJL|2(HQLW;N}>V<`~
zBZ)6Jo_nX=rE%+ZZqDj<MQ95DIdkCTk$T3vspAO4?`D40Vm<}6LAyhH#uFT;N5pzN
zDg}GfZp1d&j>LVHDy^SBL>Ld)HH58Iy3x{2jgj*usFs`L78TZ+lU1k9+H&8ezIb}(
zp?zoG3a)CsTLxN4e6Q|-jte)+&GP&`qgQF~6W4E7R3{cvLQ=bJj>QVNm35buJ>CZ?
ztR*{y%qS4CFcRg}^f2tQk1lo0YGte(sP=2rM7c_>3fgD<@@4Dem^;9blj_pnl!|u0
zlDBc_JaCYie%8RjAL=D~gj-lzK_nT4ZPjr?zBXIv6etKDhvf5uhTPYXqPH*Zfp;jE
zJQXO`s%xf3*sjEha*aJ<^t4Eh9II6V*ndua&-3->5A`i^az8k9AxCQVO}{Dj7sXx)
zyNExHyxA27Hu6!B2%6$T1$7fTx<8K`P{0sVoO^^+ton8jm?)Rrk>WZ}0NUF30B^@V
z@cSKX8@|DwpB!~OGyra$-0V5)AT+y58}b$MRQemJ-`;xjRDgY<)*%-itndxnf(a+O
z_HdfFIEp_E>ghCyxFnP|SU|QplTvi}w-QHQKKDTS+%zvL8!o;FdIDEq*Lz;~Kt{JD
zbq?Fz%2{?E^CF!duq*2Dx$L+>i7PZa3%@{NANIze2r4cj!S5Szm>X&5cc!7XtAd(l
zW?f0h3@uM1<c9)?_Pf70`?uitvfv4A6$I~Ax93O?qXDVlz9bMi@EXD)OY#~Ge$hO$
z_c~Oo+!*gFpJr8cn6m6<w3;!}fWZBNbQ#LjFs45gBgN3kYtAkFLQ6xKKwd_7GAsu|
zE+S%egKPiktM$2~j-I^3hk7%q%b58z7bDvn`O8Ea$V$;&{=9vq)VW|R9$>E%Gq(Yu
z%|nztJRa`62dD&s&dy@=JZSrdgy`cHqEd8*iJ%N=adXAC{?&-o!r9Kr^8&}&pNl2m
zJkZx3DSsZ4PJ0E_sA2OUJr;6<n9_&qBwNTbrkpm;{SnR!RuwIi=k^@amS5!{!z?nJ
zJe7}LlTE2YyY01m%Cn&yap$-|juRmTPAM6QkYdS%1OMMgye|=4j|Ie%lliYNI!S<7
zjZw*mvK;^E#fTM%ZSk80f6k9NU45KErzc$rJJC*>wzZc`P-|<1a{DhE5N^G}l_l&|
zU7Ya*e5g{*LRFK)>nU<2$IWa;^u_W4L;KW^`17(s7r#w|0l$S|Z_C3-6K&`){!Zdz
z^11f;4-;nR?=5n2T>?06-9?^8Ma#u{D$gY$w}{x&Yk|vYK4|8yR^iL3H;4so$J<U1
zcj^mpI6O`mewP|ns}bua--F~B$vfVxb;>jNdc)QGczfGJ>$X|1$GCxg1SxDF!x=V3
zeh*CR>lgAw?B0Ao5ML%Va!+&DE={!A9@x3{K$c701EpdmNPkq+&gxE9&cNg5^?DWg
zDi;<-+CPF8f#Qt;$d_qY_dq@L9>|h)Uad<b3{S;7&LoaDt&KV6F{8@?oGa920oUt*
z50`G4>2!O1RNc(ugmaO52A!_6=?p9Pz?T!u!2<*;jBqr%@B)?NMl{WMyp`ai-Rf6r
zV=Mp?Ixwx4pAosT8B9gtBaP@DH^;u67=dtEvs$x&e@@{I1MX1aquPeAX|FuYNdE$X
zuPE+{(|qF^5`6?F^AEkt9mT9WX;$b{MMGmEVHGIap;v2S$P7o>-zK&~+B<ib#@bz5
zpC-}s=`CnbBcdX<07VNTFB`|SFSQlR!N}Wy+umX$loyg46Ub2d5{Y2Hd@TY(b22Cb
z#Hh^Jam=YFi_BgpX-zutVW+41Oce1ZVLWo7PR46b%#j61P?vzDBK9>zb2+&xZ8e)Z
ze$3$DtfMUE3Y6AX`~Ie@E)qnzx%b1fDbg*nRRG4sHz2Vnp$AMHpI3{Ax-iH^4ciKP
zX!Dn&nWjn{mQvsNb4daIm<7nkqjuq6r(DTATtmniqiPqJlD{Bzdw20gjI%(${n-<X
zRKR+*qnxj%R)wBK;k=pa6Oa|P{3&>dPrK|Z^#CkM5gRvdTgO*e2RmQ-(l9;NkDmIF
z%C(A{Hh8O!Zp9nt6cyt5;D8cKrZ3iidE}%&$c&F=6$9mFF19c+Dbc%3fX)PSP9}VI
zAS;75rDp;*i*;DrLGYvtnVas@tO2PDGrEwx{mdp-&VlE$aOW%mS>_$pKS@=P5e*8%
zx-wq9=-s6rJ`%}$V4-_r>`LcmC;;)}m)EGg7*oZk3#3#&%S#$_I^Qfp{8h74Qc@<`
zDR!&$GR{LE>d`vRL%N~N%x~2`L^u8&b89pVw}Z4c=ixS#Ea!htZocY#8>IU)$nr=2
z_K#X$W@e;=_=>sm<~=aZd;8<>n@VuFSflFS_exMd7{7pRB~I|FA8g$*l099byZv)#
zEg*Z>el_t{{a|Cpr{|fHXxe$_^WTy347+ToKha<KS|`brHkxxyK3Z+TLPF;C82r;%
z(xOBtHa=JH-T0-8f<RYGSONH}`~wa9$Hyw7%8$1<Bj=m;4q@w*X{xs@mkaKTlu_QA
zHW$)E3LF~e{Ah}aa%~0R)^{Mrh;+{}Z^gTerjFU5n3WV^ouk#LWr1&urrdbgY<PJQ
zJVu<{xDsl|DM9ASp}iPc;#Ab4Nl(_&={3JoyMEPT4*SR)3eA2uq;+o6n6&ATeB!gW
zm4;=(Q*v+xU7VdoCQ;sa#+OC9Iq@wFR6pVtE${1Gloxt$>Z~YUNwBA}%CPWdMLrgP
z0V)x}zD-7}b@AbjB5o+4I}LG_jf$yH`|p95MJ{~f?#m&ZuA>x78$0n)wXfYxu>2*!
z2BQz4KVq3^Mh?M6P$SlRAWFN<_-sc28ELceMQ>Ns=;QV$l3;+qItOt$-Wzj`<twD8
zr-TZb(BUfrZxll?EP<V}6vP#mbTB*ki?Wh$w;gJz<;Ol3|1NcvxWNJxm*_w6mN;JY
znqPzo?!^@sbfEP7OSQYX<~tfx6lU3zN{6Bn<qpY`08y5rcJf1g+Ae$iZClhCJB?4+
zFO`uE8~4EIfk6(yvS%m|_Q{a54bDYk{n%GbrDJ`XyHOS@E(CgaJjT*i9Eky$%r?6&
zOUKiS^SWbYl?9-5#1HR*&s_(wKn6b&PY~8SZ4Uh={)wV97>V)XOg#U|uPcjh)-A6g
zG<8~}|3rwtD1jPKI{csm?bl(aG0Ke^J`#>w*V<eQw)vM<0`B={46IB=eOza6!0X!2
z7GPNaa82F{j`w>nn;!97kaK_E2+~}j9!+bcDuIfp1-AlI#noB5^E<;umo5Xs?vD*-
zMeL(Hhb7eza%w2x_G{f6(9XRy-DIbGV104@JaTh)r!XgZz3l5!c9akoDL8ltnB`U!
z8glZ*=V7yIE%8DG1!OBk@0M;Sx;{#@2<!i$qt~G&XwL|ZAZY;vr3Y}bYCsVyF_EiO
z*bl}<dGc+LkE*K@#gvq(=2(hm+nKR67*e_<n_pOqUpyx6;RUUI7s|3OF^-W5>hG9c
z8NU^={%J68dtj1heIn;w$}oTA`{Jb76q5Wd8Id>}(XAp-Eud!^<T|LW5n(5w)({~=
zcEdv<x+JV)+A{m_rQ|XwIj(C(k`^~0Y!rSj=Jjg)v;^CBWPV$Va3_Ru>RM&Ok;llO
z=Nmasf=7o^#@DS4W&q7uo&iBQrxG#Vc0nnZnCLB_rxTZjXPEH0xIUqW>ib$Tu#XNn
z@&3#D`~%JZ^>syT1)%XZluGDIiuZ_B`=|yo;-c5}dA-CA-Bm_a^O<9MjD|z_r(Mg1
zaC#K$7o?rn^b&Tu%z;VPlSnUh2GVYZm;Ja(uD;34VxAaNKcJ@zlK7PIUU`${CMvQ!
zQHCdNwdIuh?OP4#t{8E3uAjMJnV4esAQNef-Qf$X?98v&e^dX8*#Da7fUQSMH_kK#
zt*7B9u2RmACivZB*i98#N+A(wscpFIcb8^XY?~MKtCCAU?b<=y%(z#y{Spg!Uz9Ij
z;4OWl^7EkL#^CiF3zFQ7Eh@5=)aFY#Di3#_hV5ND+~JISQ^Qt^MqzvB_;=)$#vNaf
zY&35~(<+rw#5WNb$|0FmpYj!M4nAz};93d32P7*{{`JtHJGuMFq$#t`A;g$7<g@uL
zLU%UHtc>8PZEO*PmeUNpqjgQNj}W<29NB5Iq;fn0)oB@*8%|^THWU-HGvt0=Yyc%i
zDOgA`Z$4Ds5ZrufUioSE9w2#hi`LCMQgaZs%lbW!Jx-K}byefjNzgT%LrI^%nP6f^
zRJl9gmN7fLYjLG)31h`qBI6c?CrKSPswlqBH$h(-(~s5g<L2!=ajHk}Fg4*>v3J|T
z+aTJB#?Pr`)8p4Q2JpH_qB4}5ziCsvm3UtR>zME1xdIt<ji<U7Cw0nLei3C%QTuJ-
zjVgG#uofm-S#l}~!<E#hzS78t!PXzG7~i1d-_db6PflD;pbDMt7V=WrvVLbM4Gtw@
zpm^F31fsJmgq-e~r`UJ8Ft5nnu3D$5I%7J4Ta`kyZbIviSokH^yq4D^5FQlm2mMb2
z^5(c3DPP10_8;SpbXNBz?XNa&t5XXac+~AlVIOn=+3IweCj9>1{kOk5JX$P*1x@&G
zkgvARH<}6V%$Yw9-XKnokghiegSv#dXm)c*Kx53|w+asQHEBAHwN^nR-?$UhW6$KY
zR$?#Sw3f(^-9Gnroc&bo0sASH`^+8Qr>i?l_Qt=op7^Gi8-J0!bx44Yl=Q4R=m2c(
zX;_>_d&c!LitU0*j?IGU!3MLrZ)ZRYzxE@iss!;eysXAWq2)%G$8O{wDHzR{%`dw7
z*6#r#J5h_}$VdHn1%<mDP>!aa5{-@4L^QSjdNEKB`|s%)RNd{#7)NuEavohn%r(9Z
zLTwI@-E8FuU+L!4@8{V%cS<Da&9CDSzm3>+2F3)12^$5VpNi&SGP!SsOa?B*9mJnb
zOZ<xXc2@$sODXTa5lCG$T2iU)LM6an`CAGa9@$YQYR|d$6=Zt_X#*H1<f5AfHvc0z
zB7?KGKBN%E_sk(X2UV;a4nT_NST3@hP;;py@4RW_n4ig}Qda!x1L1>w4{ehlE^yI$
z2emSh7mXKWF3VFSdC6nOn5ybYJw2k?)>d-l7!UtAx;vFNEyKia{UoOEU@}6CO?aHJ
zE9@S44Ec)T67`DUr2r8IPq&yHdvva*o?f4IQ~kC^;=bjm9ySi6rM)xgF{-qJvFkKe
z<D?^G676igRHh!+k(~)&WBKpqOzcHd_rR*|vK_JNj@V689xBnc{Ju&T-2)>uIFHO5
z{JC4@>GK0)b=TaV3)V1zqVf4_jnXvoodUz=Gi8B+MaB5Euap?n2H(=5(tbz)Y~Ai(
zfddjEv2KELhrM0nJDkZ9@DCWpHj5{JhCTk>onA(f!u!dioOd9PbyLWPvUIGs1=Gj8
zv8bnXysTBr+X(HnvmyLDB2;<$4MG6i^=)r~j6_162#m@DSWsCE_Afd-dLi~~skF{)
zhxLC*-HonvHxKk-Pf@|!Zz?(#irVFu_3LAz%#?SdIb^qbo6fe04+1(}kvtlSut3ck
zlWcm1ah`$p)5l1@{t;9{Y93bu!97F4lU0BAr3rHsNZ85zb<OWwYs}tZp}jO|_U~#9
z;*}GZYbc1~`6omP6vcUnF_ru%bm1TdMg757S6_|O$?su&n`-6Il+oCYk&>@wC-E4C
zQ6Kpu@iXrx$2iI^56!M)3@TBqLl<{50mq-u-L+qnN++cyci97+ETFVWp_I(^B=c<!
z=_fxY6r2W^T^TVEI(rvTv7ERM&FKWels2u$`Wc6tv@N|swGnDy@M&Qe_y3}?4R;Am
zdHr%&#u3!HFG}Zz_QN(R7#g}g5l!SlLcD|{tFmHgb*q>oOwGN}CUq*(Dbfjl)Ja+A
z=gavN@ce}M0e<=~Cq<{Roo{;RZB%+Z8r0s&e^No%z253YFmGxsr@B@pFSu{pE{j%M
z=B{JL(if-}q%_nXmTWkBp<09pZ593)ikAy@sAv^}2EGfNS-0rg0Q~#~I1QjBNuSu=
z<6&@an7wtbY@|w9cH*0_MwsQuMunSQ_k{&j;bP1+D4@S`7~Eo^t={<HnwX}Hu4n4$
zHyAW$sol#&q!o&7ulVSP!~|h0#2*jWgF+4W-AP?U5qqthj~FhSR_tWL+q2E+x&yhc
z>5h*)Hg(j5b4!xB%6I(qCxJBYKf_#ZdT(ju`!75gT1(QoiT;`!@5}9;cK|-3Yx~Ld
z0mCfu<Gjxr9DD6B?2U?cJslTa4St;j@hisET)5IV{!<_m&L*h`dVFlnlAENLPCcL-
z`l44fMiHgvJ3DBhZNimU|65HD8YAo_zeZg9hXF;HEz?BZ^2EUb1y*1nk!BJAeb&*1
z=MK&jvjTt3Y<%?gc=$`^n5Y?3Z1peW9iny;7N-rr&1d$Vb!8m&s8-HHTh;b9BzeTz
zn{kKzN@Aw)^yr(nH%i%LYU_vO2HyK)1W6BjyCBYhkOUemUAAy9dwFLntTg!^==Uh6
z{3KcL^H!VLXlb<>%j%l5{rl-<c70`C=ssm?hND8P3U#{`y<peafTl~&gP1|{)|E7j
z8TWGS>ZW>~{U^VE#dnD*?JRRQjCQTpzYLQD6Z>-UOr!Y<a-zX{j!>$>4dxuGo-kR6
zsyIjs7@tj|)ItVBx$s=IoLil;@!YB#(et_iCtmuCE3fwyRQ0|PqdVJg6Tt!(jMM?_
zFdj+LeEHMtyvOq`czJ%3Gy)yNSG0V_7t<(GqPMa4<1!)gQWffm_py)UPWp#LU#%yg
zx7$ST^usy|^f^+wLX)=Z@e)e;mxUS=O0093;;ZNCBJP0?TE0R`>Ju6o5t+WiX*Bmh
z!SU9e-@!eQG|g_P78$`j%0SI#b?zlUCDwgrhC-yjQXmSHJ|B(9xl@W2L-P7mocE^c
zV03h+jWm)gyrEnY<bKEsQ2*GN1{~NGf5L^Fuaf*^b$Y|+P%<-e52zc=_O-mbZ6BGI
zMGz{vXG{iB*gczwl;$vfpv|NGun!MmVyeOLDUHGFG`HbHiac#apQl3q2x8k)EBQFB
ztMI1ege;4d&0(Pz+5RcETe?x6wvs!0Y*R12rjI4_7x{nay8*e#Lshf}=U24M8Lx5M
zDsH{Q>B1{V)*uJvbb}$JA>iOmc?NnNY9Cxz7*Urfw8-VFJ@%(}bbFH2o&ToZYm7kD
zx#KS_-}n17Oy=Pp)RCj;Q+o9V3f2^4ZeM)GE8E?xfi;57Pz%A={~A4o|9VTMDd?oj
z&3~5+pE}s_Eqv$b*@5>PYyOPPxo*1x=^*7nVfw&0$$k~IYdmE+VW{ypdiU#cGbblA
zDDGtTu!9r!mk6ET1?+gA>GV(`8Lf)-M)Ibcc5<pQ+vu79?!_ye)#*2`<-cXgWosn4
z%|GL>k~4kv(e{;G8*cpTxYWkKj#(=--u9{L`>M;5KK@Vi|6kL9|HgOHvMl#6ug4w}
z{O$`OuMfJw@oV-XaQu>{_M)cHQxDi_VB)WbdAQ|TwWvK!S3d;*cpSV=gkkCl%%>T}
z9VEI1ow<%4?w{crwXGevaT^W%vQ{4u-6DUFUOR*bI+*?OX{d?JYYA(2w)4BLJ#vC0
zQ8c`jRL1EtHMOmNOHSS#%C_tHHjMC6iv;oZ=g`{~_C8_0N!-qykM0*;S_o=NSF$iu
z2i&f&(vc)KjN)&-7>Jt$_{B}R9R;YhiSyk57?e&}Hbm%K8@ZeKzpSZ5U1MoKQxP4)
zTIn#_=7f^WN-xYxmFI_D%kQaz6ZvmF=~|Yjdn5SEjAh{Ss7}<-{&1O9)ff*io?)od
zkL+K4Th9)Mqt$Ak)+P;KUMJn1BcYAQYe@M<2(sDq9vCUW_=6PBy$5`nQvewaDP*n0
zBHd6{{(ctsLObPGij=pO^Kf8WoC&42UlhTn`(gAbg|V4%!5lGFqn%y{-P!3rEw({b
z@YI5hPq@>rPk!A4V>$&%q~o{phN-Mk4||pMUrMOk5%8kCtwdc>8HEaw-(gnP)@cd6
zNd$`>a~a5&+Aj?0A|=-OY*A6XAODUKUfTz#8_IVuot)`oDxznOKN~F@e?p4NNQw(q
z)o2+d-JkSL5boT3PRG-T)liqDOnow$q(7O{6wMTv{nik!2)1E!y9eAP!YxMsnf!Bz
z3t=H9+%U2-fAToGg`#;g+>KW@;wyzu<hPRNPK2_8ZE^dzZH=tnURL;A`?URzfX8-`
zzZdtZhX=z$#vSL~P{F9M<u9$%JUT7TT6GW?nE35f(}k+|hASeW{UxWUsl6M9>-Jn_
zS5`U08eXvu0?x7jwt<ycz<Cw)O;o$eBe?}1mJHejk#Pe*quq;RROY=O7QHn2nz=Zu
zWd6X!v7F@rt;fRLSZ&vbV?o-uOnAgEf@ZeU^>IbUKZ>%<{HGKgMIBSsglcx|ouaR3
zg&0f=3-P2LxISsbaKJIX9K)vS$9=93BQ$1sRn<?faK`)QOCsS2*ZQyHym+J4hIPu*
z1rgox%HIt%4#k(zn>{l*Jq@+!13&6pBlt{`%>mFl%Rv_os8=AsTe_}c*7#gy+_;T_
z6YGcBm&f!Es0rL~2os-sQ>|kAc#4s!u&wWE$Zuw2g1@~YZXjhOZ4up~raTL`xlu;M
zD4W?8>2?1!P{vZcOa0YTWQ|oOg;k?57G$dH@~xGPr*$Z_lM)&ke#;gw{5_gVgui@(
zalB|fV)o?%0~@2xy!;ThlYd7tlozy_2>E>0=xsm7yHjj7-)E?RnIosF^A!Lh6uTLw
z^3)!<M;)B$ev@lBIjy!j?0rR^AgA*#IU2%D(Xx&(j1RZYuFQ?fy9_$96H%$GOgvn9
z54X}d+*QP&x71{O-Vyo?wCyy^_bZiijH&6%+sVo~%ce0-&SD(?suK^lRd;oD)1E8I
znbM;!61)n~=7kYDg0DjFL(FzHye=$+wDdQ&NX>jf1$92dS4;}()50Q{%!IeZ4=N6$
z$pib3U$8m|`aHS6Dm=dSm-U>ADZF%R+sk4~3oa_*$cZzbWPY1atf%F#FxkQeL~y<+
z1<gCFrBl;&U@wTRjPN!-zk`yD6{<9R3l0lms3>_#z>hbHhoL}&3pkVf7xjt%V|kV(
zJ0~6A<SwO&an7Axr*^zDmo;k_U0Ee8n2{R7aFb0<_X=#R-ykek|Nfg`WkWR-FQ%*2
zk2a%q=W=?>F0P^4%~4-J-#XWcl4VyX8S8Yx`IMtNoZ~|_yKSD{to2B&Ai7zcX7f>=
z)`y$v{>=s+9D9bnw`{^LyVJ(~Z>+R742leClMwBlNb2Lw1pb7T`Ye}X16FC%?}{`(
zQ?t{K%pT()k`u2!edC?IaMZCjvHIzplyrE(6I+`yl(mu=uvSDJl~V3Nafa3LDP3ED
zhb->*i87IG&Lt|^F$~3_Buy$**t6MWyg4{4{v|ir$c&VZPdl5DIHlKQI}b{IUuNpN
zT!9F?_HOg#ZGfSEsj42;Yk>l?ZMoX(1;m8EF~wd)y063gz5R8RsK^18AeLBnth)z_
zweEpyH5%u=e$WvLuGS){c9k$lDHzoE%hKu6QHf71G&_gJ3Y>Qtfp(b@C5aXP(JK8T
z+U#R==d+31S5)3+P+Dixpv9-xg@pNcQm9S{>PEG(I#@{BiO>?P7J>{>+UWM&5q~mF
zGaOdQ`Z$Uy_}2sT$N3-<9KzFdJe&VTjsZcF!^{U4AebRHZJC)Bf%3S*A4xxF!r5Q^
zt@UFGPpc!qNq4N>Y_pmPkhxy<J6xBC7hh1^h$buzzXu;bYtCuMSNr-O1<wi~w&G2k
zOYgF#9LGKyTeN$gIp@Oq^%#$nUrX+0L4wT6QcXlb`#sQbECE+<_B&vLh~2g|_l7%H
zj$c}Y9H)i832(NnhIJTpLAJYU19ZdwQSo@~glghYvwG_FLu<Y{)$pZ-D&9UK(tCQM
zj_hllCw6FHkfEu$ag`CHiJk|r5qWtbnyScmavtYc<(ZPw@2^D>Ee&xbr`s{Mx)e~3
zG)}P(fQ+{g*}sW;>H*Mg8Pv>MJR^OmU3fqqM_w!P#17onJQu$|o4KXpRME92?s&RY
z-MnG^+}rk8>GDz6pWpgC{3(w9n<Hhtgk=oY7aEw?@1_^vKZ8jnHNgPL7u`GXMtiTS
zE-jMd83zrcuw8W=52s+%_cDVU+-7NolHeY}0E)!-(Q20;JzWjc^2fDeP4Ia-E~uf|
zqGRgD>B)ZolYI1_<)nYsKT)+HRwq@RnFtax@%a%jOF@)%bvbB5&p_NOc<z>M?9xlq
z;mRr5M~-#o>#mg5EVorr)XJdz<ywTq#fJKwwM6N2*w!j)CD3AStu{>0CvLu<YBirT
zL{!gONt~U>ru~k*2PCc?kmwv_X6I3d=A(LRGGC9YXRXE4Y}61rl<2+7YNl%Gb!v#L
zxX{`5+DO?Cba*z#tysOtbvUd@{%OEZIqs=P<mU$W0Gp8UyYo_{xy7hf_W0r2gI-%<
zzlAK5JGNR4eBMwMv7k|9_zTI7%x@&7D7yZ|S8vzNbHew-=z$S}Tad@`Yvj%4_9wXY
zA}vz`&FZQmW2as3(anmo6H0|D`SK9~4mnFI7g<K7JDA)dSM>_8<tIriC^}1G84BAA
zuSQ0Ypc<f*QN06d?bP^G)`hBQ6x4SQ?A1<+O@?&R-WUqB8(&d#NbI2!rt2q4o45dz
zmAe;GX4dV_8&yotfDac*;4mkws<4T!`OLGYy{bun==~Qd?5j~uKaO(xu0qc#gl&lE
zn?gZyTn`>I<Sd{02sq}bWU4MQq#em#q~rES9_d>-(28u5xcKs?_{%D)q6@$1pgX%0
zKz_CYC*~M6`~Ac5^Doo>Vo(9{1h!p?>bii145&q(?I*7uOn8A>=)M-4`jQ2pyy5?e
z{IZ>qj$Hn{p1tmG?|LXRh>={XnkB)&-UPY(6L2{#qAdAPS&m0-MQrX$vrAw)9eCYP
zz06wJk!t<0caVuFK3*S{^pXNlaiXnuB*z<6->79qcTG(&d}9t(B!X-XcVbVf#?UjA
zwP(;3IS*fi1N{H0G?YkQLJ|5MsTgYI*5lnSdBkaD3uvPOyETtXmTiag<@LG=zvP7E
z?r6bdjpcY}-Q-8o2%ALQ6Fz+4aQ*S@^r%`#)_J^(VqIf-oz=`Q4BBNoKxU2@FlDb+
zxP6;NXz&87E*QD~B|Yh;=+*~c@M&Z?Y^}d-=^PWaVdYf}`l9iqc663cZuv2u3tF)(
zm2L<rubSnBEXS=TPzn*ke{r>H2>XjrM%1&|Ga!;#t$!#)xkmx|wzZg=O*ENtbq__J
zeJNS6?{L=EI+>kUMDI`|87-3!pFkxpn2%5Sek+pxTN3qG>qqfBC046WfQi((U$vp>
zgPy*Mz)i4<F|tG1<sKjwxzRUU<a6Lx>vg9KAsPD+@;1ycX`(X+%vZj-9p6sa>gJ97
zRT?636xVxajo0i%FqbzsJA<!dO4Upz#JgYj%VQl$es@h;<Z36k{T`cj2G&>({=;W^
zL-0~!97=`lEL@VXS-tK%<uiG=Fmn$m)h0Sz%vb;UIm2ruqH=g}2Q$P?V|W>GXR)x+
z?<tEJ_UVJ-9jo|4jyy6?(OPV)7PXPn4+=e54UGD3s*tb^Gs_#t?t6eX(dRRyjNpsO
zWtDZ&Du-sBk#OWoO@mL;W!qT^;p$&ZPiKp6JFdP};~s2MYsy8<-wx>9B#*P?I8SKI
z;0>2%Jx!7Y)5OTB@WX+RDipb+l6b>o*%2P5LU3)h3mte69JFg3sy*>UNO>>n9<Y0e
zgqpjDs*5UfzY=2#(Z7OsXQuJcj|VI;nA<Ad0|-Q}^@(LoXurheA<K<X6u;<!uoZq$
zW0(AT3#zlOJ2T=vrv+`I_e|Nj*wYVJl(D60%AqcigNn$v+YJ(89Y5j)wA@tJ*$Q2N
zj?^^vcgVO<L$&rw9cX@~^`ssR?X2<%T?UlO^vzjZ7G3xw!@ypqil&q94)>>#`pSYP
z-}d!B5Y9;7;8DH1niD_NnV@_s@7tgjFr`m_4@4^JMA=+_=H`+0a?FlaJHmS?8ZEC{
zIk`ibHGj$3p|RA8JJiboLjU-vUvZ)sPu`(O4B}CQ5=u(;^4nIUtD0ApQ7z{iF{!t&
z0<P#L#%Qpt{p^IjyjG?Ono(Wl-qZp+%5kUoHs3n)@-vR7k@DCQ*IN?E;OJTxLkC}u
ztxAOHGlU|E|I`s^A|S22*B~KV&+WIJt-CT(u80cRY=SE8`2Q+!XeIPOz8GZ~;_<>4
zL`fu+7+MI4+c3|(JD<OGq00AKLCwrb#WqA{RELZoJS_$-M;=VUjdn?Cl>M%!jymM<
z*M<q=cJ49+Z<8Z;bQ0ZeVRNu)O{8gM1FE5lp*tA1&w$$TS$e)))(BP)Qh6;wx^V6k
zDiWy41oa3+9>Emge$&DDu`v$z^hxc{#v&V{-Py*&EXN<^h68^Dr2^3-RK;q4^8T09
zwl}dvjx4VgSCbXM*q4bl-37mauMK#OgIpultV-G-YnIIHzc&iw<+bSkfsP2L&>o=(
zU9kxgaXqbae5brC9L>^S73WFf#@uVq{+$>2x3nuWI^I`Ujf`kk;PXs`+PN;p`6F#W
zo6k@lPC|#{F82ymK{8F?L3Ic++yl}m)fo-GUsUoK;nVPuvSPHoST5t)(FVQNGTsB8
z>Jp2ur`z-IN+^_A&}{c~_w6?;>jRf-0{axHLo4IHyp}h6t!-xEW4l_QQ?X#Lb8u!g
zA#ahd;8+a7#@QPG)jR19*;seJjtd`3@exZo#YZ1_uAwmhE}teG(8`3X|80db9zRX$
z->~%OoG)2SlqTZ4E?Z61{;d%WueZQToSOB@k7u_^)z$a?ksWiXRUMP911pI;qv$@C
zR(dI3=po@tiaS*A#>jzQ{mt@&`wRW&E%LRyFk;nSd&^-~Z|oyvvGQFgO^I;PEm;RN
zM8b>n?xX>H_9WoP&Xl+HwPRfOZB>KF@WkV`$f>f2gG`UATKDHJ7w&-vR!g>zrVD;~
zeS7gS{L@>yV?saO^)YqD`a^5-9|epYB(WO*QPBF|tT+b>Q%wpFQFo}=dc-SeyTV`l
z2TVk?M@{)zJ2hKIJ--Jy2j@z+-V4}w=L!wy08)l#u#UWyv9KtoqJ=|=vm*%j*c#PX
zw3MOKE;AXFWoPQzKiOCVw1XvE*C83242*O(vzPvb(53x)S&w(4b{Ye*2EMO9AtHpq
z3@>i2z|njY+sW>bj*irFU6$*^*tHXk(x>Ak@5+KsY?Xr#m>Fe7I@*LVcS(a=Z`o~8
zP80<}Mf4jB6r#?*h6pH;!|g~wbKWA8oODLY9(f+=H+;P(8zNr_gj$%R#;e%zhsK__
z5YG)-&{H>XY%<@Sc-)3tc?+rEHS=a8a0BO--R|1MK2`QpdJhcG$V2loO0P*y@_mFA
zO4hnyj!S~H+n*dpB0bBOdCd_+LB%;bMx9^f88%6Wv$#twOflC9XY18~-D-g*&QZG~
z?|Ogx81rAbI|Ko(xD!aBysSsfbP~h2%?xoUC|YMlsUL<Y>LJaiT7Agzp;B=pzH3{`
z@la=Ofy*FT@ya;k)=5{ON}EjM%FVGs?KfAK{Est|L^eL+kGm@AUFEiku4*YIE|kgC
zjyjsdR4-z>g!hsyN)E}_7PWpNspt^u*08g)N5~*H2MSov97-5EMX{|H`DAAQLD~G@
zY=nPy)NA-Ewxr`EV|fu|ff0O8k=mu09U#BUL^}+)MqCT)zj*vVyuM_m7$5$?qg~P3
zIrOmWQA@FAy$>78i(K*yb0%YpaJDhE4b066VA3{ao2dT6fB&Zfcm*B45&@#j%+h5G
zvXR*v<tTJV#`U3kN53M`j_X-(el5y<%2;z)Z>!Y<Sgl#D0BUW^tOfncUsQ&V^b)D;
zr1?$?p^EF6KhvE^P({*y@z(1SRe5k&fu2uX=UH2)Ew<8~0@znfaT`*bVe@oM+fu#2
zooG*G>Q*YZ-^*TKA)4h?_qt+sRNQN{-vEJbH7!?d{*gHC{H(2ZmBucAnt4Te#LgKy
zj@{7Z=%kEU@)1CLImBKtEozeO9(h1%|7)Zh>g37{Uvu(u|JX$&_|y?B`tDGO&c%&_
zBA0X;x`1jPNNSH@cw>e;#y9<q08T>j()-;*`~j?WK#u!`rg=5@W-Jo-k#kVYauJ)D
z@X}iqrNuo}D|_a|b9crW%)n7h9Vrayu9vOz5ay9_2k%GEGxCjgsux+xg=5(n<F%q%
z=_7Qs2)SY&qlf5XcF7xT(b)b%shFSL+_Wv*iB+=I%%`-4y(xb-#{uuQt5(ctb)<R1
z`1h07qisWLE<<t*XZ94K>GN|7HiE%}7Csws)e2LBMZM;XnV-?i!naCHlrGqK$J{*2
zsXA!U?U<$F&-etNJACM@op?6&RIBpOv@7hz$_+-hw7gMeXX2r^@><WC5792BYnffX
zjfACe+)}Ec<;kOmWiL+5SCV3)1zU>AqS7gq{(Tqx-9!J!QDptC+<JOsVb0@|E3-2N
zgDweew=Ej}Dq9f;;MpBfS+4^%ezY4eFjeiQ7oYJ!L8j1eRhFopa3mipyA{xre`aXs
zz_%cy^F^y1`MZAHwVwRFwWjnSqxw~`Q`=|E7DON;)(I(vAkE&k1aJvf9j(~&Dw3Od
z4pn?Dk`T+?nfmc(Dz`|?wfZMtTdRw5U0ZHFl|q^P%G7T@&nU}xM-?pUJm{b4wx>OL
zO0O<Mpt>oh<JY-I%Ub<S2eI&T_oIufTNay&b3J6venh%Z6A_Od@D&1yfuHs-JZ(E4
z^3qQdar1{heDVOh1{~G1<qx7=ZSuH`Trs|st*nHj9ho^cS=3j*=JC(I4!G?EJRf!9
z8p}^K;0#*&*p~{HkrC+c;xqJCJ^A&lYW*OXeh3pBzdIKpI{j}D`ImIgYYP~PCkB3^
z-79wb(p;{2LC?&51qU{^QH**QNdxQAy9bc1)v#-pX*S!Rm#+Bd#pzzS-1{**H(3Vx
zsM}tQgPKpot&Go3<4}#M#{=NCM)I_kWP$aE$XL*j!l#ediOm*;2D6>(FLV#QS_K!=
zeD{62a4m)UEvFOhc=umrW)}`Mw+&`Dhg5`iBwKG%^;kr+2Ql$~c@PIn--ZddBS?s4
z3KtTxT>mh+X|x})no2?bM{Qpj*5;aai?pSs6p9v?;_fbiQoIy*D^j$$Yk(GaFYfLd
zoZ=3}DemrWNzdCoXV164nb~`0uIrp1x$=gjB?*tLd);fTr~bNRZoFzL4-#^4UbLuj
z?W-$%tu!T94>jBAtI4h_zT1ZxTT)u4^3BLUt_i}G{Rt-&^AgyD>K6%i+e?BBT~edP
zj2wCBIv<hbQk7^M=jJQ&(b|f=2KW{lpcDD8ZvY{b9}#JF<OY;#7jPib-c#qRJo#Th
z=i={3=U<N3Uu04SYO>00EAfp>_?xAr!Ydwcpd8<#&u=9;wQ?{F2i27G{)DS@x*P&j
z+t|iEBkuct{?lC1qO)#pF|+PTQy4#XyrEI!=Zb``H$(F2SM1lx(ClxAX*J{Qv4>ex
zQS-acERl6A6gg`c*dhy{7umBUji5P@zc~TaKq`Z)&bBl)9z5CQb=?F}1?dXH6}`e7
zw56932%`wM8``FTabojdVdOiz?g321Nb;wdIls~@qEdmvm(-O8NSEcs7yC3d`_C;1
z<ll0Q1!g0~_`n3g__ur=n<dAVM~3$&eO1g;lH|8)b9&FLI7YF=pyieVN(}nyBK5Dr
z`nKIJad&rhy*`Hf>-opOb}N8UorEM;@##Ig`Di)C?D(|vz0Fw9LLl_nHB9B4`Vh1#
zqle~lM8C*K=|8XXWr4$-rI?1t1NGa9SI*4r_&{XW+0b582E#}%yzr5nz8)TILWt55
z+`EsWP&GR{U!r<O9OHi*TI62X2N#~9)^5`SB1$x_>XL1h;lI*}zD-rh#(mRwMlL07
zMt@Nq`rk-CXaR>*#`8vw(|PkNrIDu0V27yQBL7D)-`3+P4_51<9?0F>%HHjQSZ4JF
z^hj|&Fx?G`TC2yio0dY;X(iLSE=JGnSuH8MYOAig(~Ob!G+az?_?~C`;q=?p3tU3d
zxS_4O56q5>?B<H$MefA;?%C!zZ$m<)gkfi417LJr`g>t*VbXDi(Fad|*WjpuOgI#%
zB$R&A!|eXetQk(EM9$8&($d15m?K_;qvUtVCVv|#3X@3I%kl&Psj@v{@@XqaS%cLw
zW+nprh8;mTB5y)?txepNPIR2pnpno|;#*7w5f7#<{j$7Km{eIIU=887hQ4C7tDVyL
z3D?g#1v;R%1e95u<)%ZMroWa5k7w7l$M=tPz&o#jLP)|#=&A*128EPEAE_~!!qXSr
z062_%8mKLG<Y2y2fQ#&;3++a><_dH+Btv>v@=4+>168LHh80^nbAlBMR|9l~kWo>^
zEese{>ATc9zYYyNnW~q&>=SD15{y_8GgO>--!*`&mBi<tNLgcHi&@&8oWgJlUARMk
zpq`z!Eq3#KxMIVBd_Y9uV_7H@#3g=Oe5Rqb7x-~JVXz*N+4#L&10WQ7d9^3FJ9m#h
zw3>5u3^Qh=<nhga^kE!7M&Zv!gImO!AT91jw3$xzw;etu3#poxNiyKYL;8P36H4>_
z)~L??2^Vc>&i)gQ)eL+gt;&>;+wVebMr3)g25PCdJu9!cv=*k%njjv!fnS@<i6WPM
z_bxE|1<jd;<QtcL2sH$IWG^WN(v4{@YWWIXtG9y`mOq&$s}j91+~0rJc}(3lrD8cZ
zC?gK30(qog#@sg$wf!JO^Qw(XqTxtxPzVnG4j3n~f|1*@_wM||Tp*gVhuv`lCzQT#
z!ft~zfj1ee9G6q(eW4fR#d!K8-!M~ad7t~dT_^#is%wHMd^4Rv;*AToU*_Hwq|wgQ
z>$M$&2-Wxr({`o(;Zzp&F|~CwZy)`-FaAs?Q$F3C0M7$OoZQ|EyB`Jl(pU+MJUPuZ
z1eQmp7D-}DAs&|2Dt5|Jgi#ePo1xGdaYDAVE->F{S;CvCLz4=8>cHo9U)5_{2Fi9B
zj0ANgPRU>t$CLVg#`|-O(?KqzO=pWlZzD&aTL-(0N{AojpCgg3_|l8uIVdJc|5t_=
z1)YzHuz?$QXMPlI3+^{lY`kY=o&C2>Z$T?m4R7&2!1Z60<*z-%aSIC%!Xw}ad0T6P
zDv@GO(F=p|qX7foBqE&oB8Q|s=LdqPJ#zXl|3)k5CtPG1fZjhekwB@Xw@l&6%1|Uu
zMo!W^4pp&<=MC(@7j&ES#Uy9%^k5zPz%-&kj0JWZ#KIHad6K}dsA`hwBtc@Zw=+y`
zNxKi2u<Wk%tEb67;Y7n5K2Lr}P!0Df8qC^@($mVZrZC^bdAGi$n4e!+vPT0iOZh)?
zM})0Zf-1f$R#l})k@hEvGQtD~YcD=LJx(NT+evt<ZzL2A5CyfjccIC}W|k%{R^+Ce
zqW-Ff3R|_M)>%zH43}>_Y6%P)5Sh^`VQ2SFhdUOYLnW&xhfqtHxEQ5YTUAmIn8}8h
zATD-O$DJ*0ttctx?S$GJpnyR6$H}AlAHccCrJj5pxQ8?Hll@xnq6mdokDeX~zM5)-
zAErt$d<wA?RYYNks>?$3PX+ex{EAD#fQF)^xC_WNm4=fy+}pRf$_iu3P6)-88ZlPZ
z@gP=5{@dDT5B<SK1@nZaES2RS6=qfVXgH~Y>FD}^HiplycFJSG6b=|&nL|Mv1c2!D
zQ8_GDPJ7q^#zas9Mwa9wVki@w-m-;Sd<Cqt5uJBaO&Ub6^e|*{B)hYtp6q_w0G*A4
z<Xsk0aot|Nuy>Gu0ts}G3Bx!fPt8tx7BdWK@Ph!pl58^QQy$WdmrqNwbw~mn`^alN
zNpvLY$BaFfso6C9T+R(f!74l7KAgd0J&0uvXDx{1OCEhXdNowr^ks{8f@dJD%Yn{S
zsT~VGLX45}-%`zgqoRR!b>5AuB~~ky|Ix0!Sw8vl=uzd1Z(a;oUYtH)d&a}d?x{|z
zBv*V1>L$#XhyC0N@^?4xVEb?BYK>Yv=QI5w7*Y|%9HQrCaxHNGG@>o?PJhBZ_&ET2
zqsOy9;dc4d*o7ZEg3o^+&oAuarg-O123SB5YM*|>efyvd-Vq0k3V@83Fa0qsP8yY0
zaO1m-m_4IV<jwiQvi$ZOu+dF%`mH|+l!(Y~7#QULst#2Q28<396>XvISAxQL4E+NY
zcrH>fN?3O0-MPvo0d-Xhd7OiOBQIv0qu|q`xWathx8I=w^8Gc2agiuBDBwVMw{XDX
zRC-}7%Ohl(=n!WP8_}YV!CEna|22r^#F6SbRddETroR8HLHv({`j3w_$HxvZc5~8p
zhg7*%IL1+12U<(q&_*+g<{VQX*S1>1@XmsA@Dhqr2XOsD9cBA!6;-{4ALM2&`{!oc
z<L>*sS_1Tga~BgwWLqVV)TSJf+P)I`(vTC7`FR~!?0^c7%)_u0tu}l9Hg+RB{!kL3
zrgwwSm;rXNeH^^|ZfZzPWzF%$s@yY~KI^Jc9DaQ(!fk**2E_r12EV>jHPO{!D|&ZV
zjur#zhA}Gw`=0cDQZHZApq&**B3v7H^ie=jbzEy5fH`XpBd!fmscf#BR@e$ZT1W_{
z84F+Qy&uS^rrC5;b-#IipzdPja&$-)qt(E5Du#K?l9EBJPl(|B8B#!Gq0|X%uLhOh
z4=%L!pBYbNSJl_V48|g<Cib&SM4i=(19-rckeFUgjxRj*Nq}LN*<rCL!P0`+^$8b-
zD1$ycz>_k}vG*(5q4Ix-e<EAnSx_q@cN2Nz->c9n6G#^SIB4vpENwXZ7V)!@=QM<f
z?K~WXC#%+u!#FdNGZ*`%Amyl4nOjleLJF&*J%v^SyB%ee_ys=u2+_K+S{?SzJ3>vN
z;P{z2<%-1v{G_jCRxf?WrIX<eI_aJq(CjXx3r#<$;KZ%zAkfPq=&(%)w{WK$q-_k6
z4~&H0dZ4ldXCr2e^5#xH<p2lR?P;<!z$Sd@_hTHd2-u%j*5j<RMX!H3_I6u`V7KF6
z$f0;hc_Mbs#l@ADF3*+e8jHUwAtg3Hq+IKzK(O_N?^{59?E2$(@&j;M(wcYYhE1)$
zo?g<ZoH2+-n!GAb|M7&$+Uk(8Ve?(`{FyeamdRFA8>Nv|Gua`~24AKr{zyfY6KZD$
zs$p0McoI)G!n>hE>KkO<CQ^3VcJ&KIFW;K?)i$8XZVR;g_+6+-ULB+$W$jGmC!Eb;
z@22$G0bADG0xSH7*qB+ee5mg^%oE7UYBPDD-BxZRPoB>P+5%b+<{*g9Pq@*cpK#Zw
zZ9NelkKFlV@2}38JmrJKvK@JhK6v)fIZRs1dq5KCFX93&3zgI-l|Se%pzn-m-HshU
z%?ySsWj;)<12wY@fOd0ez<CxKC+zLu@*GIKYO$pu!LsaYvPTS}Hhm+g>U8TQ-S5w|
zGUy+KK?aQz7uc%~Au7J|s1kDo*D+pXv-x9>?XpWhE$hziIae?0KCdq`iOP-Yd<JjK
zOPxZTmvR3HBXn-LV;gBB7WXU6!>h*r8c}HtId+;M%Ad$#kWzogYa|8GE`hs>UF8Db
ztx>ZOXnn+Z%x^~~>{=#)H6ix_qcI4)k)#Efwshs7=tTC|lR(i-F}@7z02+dUyFP=z
zrc_2$3dS#7^{0H9;<)1XCq1bg&eQvOi}gxB;g}}4AF2%>D1O2@zyi`W2h(Ut6W0=@
zCBGEY?mazxheWh*HhR_l8KfirgKeo=od3QGmc6Ue5OYm*$Bq#lJEmw(`M<hef9`TY
zrWMY?uFVuBYF6BXm=Y5hJW{Aw(~w&v5xI`2Ky2YCd$-wFCgmLVN|wv4C)2@uf%vq?
z6&)(~2!=RaAuDAokbv3Gaipt(Tg|#BwP(&pP@hSxQg#4rpi%z_c%S%P6$@@9=Bp1>
z29~Hyvt>sUK?q@->3{fQvRcD$WOF6*T$sk+_QwP9IZ848{m39VkmhZd+Bg`#mR}BI
zcg`t<Vi$We>+42n1<I+f<VT{r_3%^+y36c<T)OHxx&t$v*-N>nH}3Z$ev5nQm4&y;
z0JSFyBN?e@p6l#zioUTL?Q%NFQFE%wEh#AuJLPiv4BI7iCsb?eB(&*@dYg$3ry|b>
zKiBF8EpC=CU%eMZ?YutaoZfH9CX414i3^MxgIUfB25jEOYb6`9f7-*=cYN&0M{|3Z
z8k3LjFlRETthkl~@zn~7ry+NxH1rN+lUvm*{A7zhS$861OV@3UIY{oHtSI&H*v!*+
zq%-wBm`wrrqfc*Opk;$ZZj2_?@N4p@e8(_GW=|B+eLJCf_ro_S?#lMuCm9s5f$;g%
z_eb=0d|m0yF=vEmUE`ju1G8&A?RWuFRuNy9b=n-fYU7e<Hru+sHMkb)>b`)y!eB!U
zM8)#NvY|^e%_)*%eEA$0kH4wiXwB%GZ!A3tHRVDMmY%a1kL81f#E4CM7G75hAf`k<
zc`rzV_@x{7b<2g>UFVI2Xw~Wa`Bc6&3z-xj3@x+xN^(|ZI`DL#_Ys_%QM`chpnH)y
z7{x=>u}sz)n(X9z)D^U9iXm#-kcG}f_}-a&$Q7Bz?&&wAV8f(+03O`Ld4MLZf)LhK
z>$PYaw)-P6q$Mg1am!|unx(cTYENH;GFz_tp=-?-WigD)GSK+2r(p+eQ-`KCl)pJe
zG_Ok_CY*4Tyx!-EPayf&lByq0bNK{9`N8{hC@s$;@60H;-6UKRCV8Z@a{G9w{%PD|
zT}C{{-ucbP1my;s^4E(DsOm@~-#0#?h)a=%ZgaH0=BoA}_T`t_tQ$8Rw<IFe?wCGb
zvL^oY=0(%APWNno$_n8_jyzJThe%U>?Wq=L^SHQdZkbaZ`~8>weRV{6e=p|J*FalM
z@HKQEmO?U6=S2w!E60~FcP#)IkI^8Vfc0(S{?RLeS?yxhpK$R?`XncUXxmzDp=|RW
zd=<f_bKHET?H~0&cPYSY$6I9M4_BO|riMm6J2GsFm%>q@aQ*m_=0(cq$VHH*YV^E3
z?|}5gJ!Q)>7Y0@=Fvy$pV_^FZp!Y9Knk@-p|E5%ZX&xH`AFH?){DBgy=Fd|Bm9G*o
zG_`vLP*f#n*`nkXzrWBEH@DzAjUs|39oE*54)C7kz$jqZv-kLeoQ0VS5?g|C69&`!
zPrEL}VXq$0c4QSAePt{XNfVv=39vu7qs#6WrVzFCyaD{?wD^B1Q2TAEs(}~A(0X%?
z*^45rDTPbtvye1(O4v^rxm!1czto4d016kqo<Co4hu4|$YbMcu)fuvklb4sT<V!11
zUq6qD7HiT9v*Sk)bNvz}nW;18TR#xZ^7N{d66sO-<XUD>eo!GLTE*cxIY;bajb+Km
zSdCI;HznL(EtAM}GOjt0L`q+YaAp3evA~Z)8NuPf_^m5$9!o-}MnB<R*fR&#ih~bH
zsuS)-typh)=d=p(uFI%}XN<?UQ;mbHIqoZ~j>c8_O|-Xp_sY=d)+0T|MZP0~UH}1B
z<v-zS=mjlHFLE*T3X-nzv_!9!zW0Hz6mHgY4UruRcHh(LRaG1<ZN(R^CYxqPs8qEW
zpFAwvxnf@TJ-{|2ZQDNkqK-ci*?S_^Q@GhReP?7mz37iuo%oUDkUPKP+uda`pIT7)
zr{k@p70S)p0arEE7;sS%=xD`r!BLS#`HW8QrKydCX;SN#Gi}pquF>=-S{i6V%Js)P
zi@0_yc8?2|fqc2G&muX+FMn-ff4nJv4hF)L?&onS*7tReJC+}dt(oODPHOa|erlH{
zYme(2S-668<ou>KrC040r(gQOoQ69qQ6S<;=8VV#=qS;5eYAP~h0+kBOVw5z<p!YA
zU$j%Bupxc|pifaHe2m^_Cn}+2$>f>+{iE1hZLb7g)}0QQoB6A%Mz(BNZJwxa4?Vi-
z-*6gsFDuMm+6@<*WWeOXV5R}%Z+e(9u)0WSxce)-`k{^CW)up4yVgi!Kx;r+q!=CT
z`DH?reX}IxJKALq2eaO{=?=d1An+Ln$=yZmZ=T9hM?K7l2UW=Kn4!9YL{?xVJdjj?
z)$l7eigVfpn`79y8(%kZtvJ&<Z;(`qY4EX259*WSH#WQ@?AmdBStWuw+Dg^Ovu3jw
z$00xAP-8k9a}v_9{IBjig%lf_5r@OMDg5<VY=lYekAu(E8Mo?V#jFVi$sYNM^?w7r
zZyzjXIrPxY>XvXdm}@<S*5GqJ!!q^w2zq$FK};ua?ZnX}1P<H#$gwAY{~cvEHx3pA
z&cj`+cM<wt;N6Z@S4D$C*TY2?V?t3ms)#Cby0cFI-8WZ}N!Dz*_HcO!g0(nyo=4N|
z&Xnl7O2cf>R^KqT#Jkk}<;7Df5Yytc&OQVVkL&b78^34r@{stxli7l0Lg0EVUe&1f
zZOmbt>F5$Nsz-CEL$|ZUD_VXUez302vrU~*HV%N53y+PPll<c-`QtT!!t}shz0$!q
z_O?$*>o4^q0H%+u28vW_I2|6~5UKO^2<?tx&mj4S9QfM`3~R=aFXWJMEX9nj8(RHi
zaoC#XcHJO<d8v4uk-s<t;6nbgA*~c3NK3aBA9D4X;WI6AA};aG4oE=hjn%-m+-hBf
zQ+^J?SIZ?&!OYO1RrI{22;W|xj>u<MrA^s+8aJMnOt}+Qm(K|(LO*sD67=&j)`a7J
z5m@IbXOEkc!kgXv8ZFntO;w}9P+~iZG)jN2WpXYER;)f5&TUJ2J)WBcZ$pkefTQx#
zn7xLT;?-{bksy?=MAJFIn#s~wB)7{~C4LhnQ&_bQ2@Z2&2+8Mr(pgCoo7ugN6Pt!d
zG^EBE);zM`Gu`;c1uU2RUW>+ffHs^8^+=#Het{sWiFYBx*R~J4V=P8WV?=vo`y=xE
zSYEAqK>S73{wMq*!bF@5i^`))9yrvRF)aXRbj447B!^Mi)3AnL?SMm|YkY2~Zc$3K
z>Zm>fP|VebRZVgSPM$O5jrj(@JHmA*zSCZ!JDZX2uhkgpDbene%b)7+$u)1sJyUZq
zzMh{Yz=h3Twu4gDvRLcNZE$J@1j}?x6Ra5_2^4Jk3AUDOM*sz<?06)2_qmUBS_NC!
zq46rjnlfs-8MGyn!~Qslia)n95>O$e6@mTkxT_K9lR2|9_O3M<ZyIZ8?DBXGqMRi-
zOpe9(AEe!ZG%R{G&9PFei|gqJ-iDe>@<XmrXBmUu48sS7n2xrUSwUgmU|PAqhtb#-
z8u~p-7g3GwaMBs_<-CmagE>rU%>srNF+RGTT0ijgqC~H}cO^_6YKE-A{7{?N8i;gR
zv|)H<UHCq@0p_>Md)<fiwMuYXuWzh>i!jZ!SO?NE8pxN>3FPR^W-ZK(AOP6_FRKa!
zK3wU~-p|l(%t03Ff#TJZ^&zPACmf>sd<6Qu?qajOhB~4<TRHz1L%Txa(0feKp(^J|
zMc0HwQMWX+1Lrepmvl|N2a}Vju)Kon$mM4>Gr<G}7^L?`w%3RzNRQakZS^=FbZ%0)
z!;j(#_@7=ZzNT7;EaM(nUWyDytN6v@#-SvzyMVgYUUi2JUB_27fm98~qrvRjA6;(2
zXdlkz!gH~Xx0ce$CJTAZ$Y|BdS4=Z(b$>)ndb|g;5X*`{V+njmOCMj#{Bp0#!mXtz
zlfw)v#&Uqlh9%P-O3(GTC=-CHga5(PW4;wvsBM;D;BJgWTHLcv`#P<tvz}8tK1QmT
zg$uB-=q{MMy<gAoo7a79H(lmByOdoveS}o-y0TgmpSa-e<Ct9tXHD_uC}MiNvLJxZ
zS=<Tn!GR7_6UQQT_KczX+gk6y;q3cIPe~)>&0jC*k>8>W0><0#B7ZW<;eRpc-=DM=
z!xHi}^c+1rlayfNKu$;@$#s88@XQ!!ytJGFjY|34faGF~7xyl#^Ry&J_L@CEcYW3p
zOF{DcAK99JyjD<4xaLy;bqjp`3Flo78j2Vw2yQBb9Uq*4M``@M^0z=ENv}0)g2kib
zqyMSwf&Zs>%qVAQL3y~K`W>a8{xRN;;AriOA38qnQocO;p3<k-D}jf>Q~)*{ztX#w
zP`W|AOt%++u>W3V9m(oUA7v#uKx#kzm}icAg}b=7d_o>`t7`cc8@j~`y!Fc8z?SKb
zX6;f)0p(_hUnD|F2_i+Hi;<%MTU7g#>Oxb*qN$5ySv9!?t&;^xvQXn|xkpVJ>Bofv
zHPqFPN}W#Bam9_|+0Hq8FPC=DIk;H~i%~rTNl!%Z9uYsKz`vMnX3n`O@5JG5B^fI_
ztkrO&48P?rkHfDje;Imo*3Lg7LWibpFDq7+jWbFQfyCRrmDD1V)4$|29H<3ExSCte
zumca82XIGINGj-U+lg1BgF4e@bKMjv1QM0<iKR*UC0%$)yK=PGXY5j;gw<pV5Qd0t
zmIP6F8f7hPHhQ<S-zmy3jP!lzyWpYD^td2O8+K=fdBuioqFRBL@*IUg><@BRNT77h
zlA|LjCu3{-q|0)=8~aXG9L*%Kv78C9Pr>}`pzWjX+*XOO5Z2{5v<7E{CkPNxAFU+r
z&kxkELYrl9@;^S}Gjj9GOE+TLip@9|VZ{FmGo(c%OWwfcJyhIS+Gv0@tW2s5n}44*
zpQ#HY==zYlMoacM!2;b^0&l%mv5-QOmL>y9O~V=xKBosHV8n!knQsWEwY~+h!yf-Y
zdO&DUH+207B;}24(D_Wp0KR}gyBxkQJ`NNVe2(XakuGdZ@A6M0>qj%N;Gne4=>d0L
zdx+B3rCE>tox-YEduKvt__y1c?{1tR*X!#SrS~MaRE_1(FQsi4+HnI@v}}`e-)pV{
z^j&aeaJM)2gtCf~Jgmm*jSNcoUfQa9^p({oqYwY+6}|ikhwr2~%ZO*0O{<_djQMs`
z`$KSeZt=^h8UtVx_x^o4DCppLi<jmPXAao2Qn_QeOCt@FRM-Ms(&OtgyIe@zN+M;@
zfMh71hMO<%lRVokVOD=s+S>0fWgKZ*mYR=M0_5dbr`#yMt%+T2M;&xj(0pl#`zB3m
zM)&vM@cVE2bq~Z;Wjxok9Mo7Au>ox+`}Y<cf#?#n>a=7P*C@7F=7BAv%9SLuh`iK%
zAJkeiVHKeEdBE@}!1`?dc0voUg(}?>RNmE}1x(3|KSGgzxQ6G<@p2pjwnbNOsM*~Y
zU#%Z04otQeg$`cKaiecl!F~97W*5|lFS3GZvHqQko2E1b?Sc{i=sIU2y{g0ATil5J
zgbchQUWJ9>f>#PgPk$|A0q3P9>03c0(0R`LFakmj+M5(NOaszT;tY^~g-FokW+3>8
z=~RUoiZ^~m2($?XhqVZOZu<$RTg2K_sSv?lu`~8*P2;2dM8yEvPWB{_R<C-|hGoI0
z)G_CkE%4$+?_6Xi2}3Ry&MtH!;+^W-ae*rKJ}F~fiV6NPT0OCfP5PZWL77gciJSau
z*sDf-&Xs!qb@@8wE@A@=MVI>ZN9l{FzAk+$ZdM^H_qk&UZ$4{IS&=0S5z9+HrINe+
zV)aJcJ>_a&o#c__T1dOG>9=#^2759hFL^km6PV+LiMNB6I){G2kh<;=>UXdodk<2Z
z_DuTMHCSOloaB5Wwb5c3jYQ(?t;9OmuUUUJS(a^Hd&+q70_&EMVidhX9sxvHq2T@A
z1s8W$n=5s_9N>GxE}?0aj@XMw%Qccuh4l3ryU@crD+$AUTDx!VEiLU&61qU<&jT|x
z45@2;dKTQzBYIk~t!lSyrY)5?JNRnVnv>rzy1{a%^D^nl#+ziw1wuy-WY}5zDK4oJ
z#hB=K-hBbQ4DiE%ejYB!((zgCew>m;c8`_#YUijgsa(%HLN|IfR*ndca*C2s^5~%<
zJ1Ie_vT9wx8YuDEA&kZY@~Zd#{f_gTh@>(Uiyi!JKQaTgY{>7%qR?vHE8ptXLf{5f
zgJv_QR8QcX`m-zJ`d1X|TMJkP+Y^tbkl{4p0oliM`GGQDJlOFG#1TkvP{!;UlT-*G
zzsw7+MZbUhgdA4+ME)*(t$gi!?XlTpp_H`Sn7r`cRj|Jz(nRt|ATEYV)dKyds&`EM
zMiie<AEP}WRG?|jpBM$VlGgd_U2=<*8g1VgDdyj9Rc?o6C9c~Bh|_g0*jVuFQnebA
z112JoKP%Lc>D8dt{w*twP$KX>@<6E$`}Vk_dJ6`x#3T<%e;9{}lF`FhtQoB=CAlyF
z&;LvugK(AHp$p1E0SO@Im>)*qB=o`(;cqL7)bp-h^H<X3=H#gV!~FbTg(vVFeD><i
z<B{5rn4Z&d&Ytl7Fcr<al||iAeiS^NbmBzbVi-Opvq`LD@U%11_gfQ^!jhQjOpi~V
zkAUESuHPnbR&=a9Vw6n13LeFZN{ehqUwogJr1)~NGbj}w(5~kO=ttovg2<FNC{yBj
zsbkb?JY|J&-2!E?VwLgILQb-Krk4nlEQme~LKFFsmhq2Rx!>3>pg;I|>X@CXeXBAX
zw0#wLfbm%2tbuNTb;62W@5xr^OhC#8hrO4a1yzH)Ym!pY?tv%Xsr8z1;wjFY-t!x`
z%R};W9Lenl1fFtjhrpO?B%5Ae7oN3Nx7M2a%BH00Gd$9UpkOJiV3pYLjL+klqntCd
z=J9unVJ-PViz=BPgXfh$XiCPA*va!ov-W_3+>JbHmwmpC4HnR0VMjlRE)~;P#|Qsd
zbgq0xa-E#3;oB!0g{-ZOKG{uzP<#0Te_gVHWX`tNULYTe{@UjBT}uv1i_%1XD~mb$
z*3x@MQElB3K2c!mdpiQ9=;f|XnAX|@%A!#^BPUqzj8-aN3O_!2<Wkv@pHC$M)r{*d
zuhznv;-IbxvQY8nolW^L`3E-|Du@SXju;ihl=J1)70_CCr>tG;lXv%tMk*wt&p91#
zP0ZBn$BA-Gm@z;vX}7O}6MboPMD*=?u0@hFp)R6U*-^Im)^nS`mu6boT-4D$7XaC)
zuA?EG+|wQ+?jEPOs%p>$37AIDlS#SuHxV6GbW7@$z8m_nagQ#Qc07xU3svddJ_pbA
zr`iSAu2;;_WF`K>W(X5m_5{VuD!f<&3c`%fUj;szPdMUP0O6OJ59NEN$S$gvUl!FW
zE)*UCqeqJD|Dw6ashmTael^~aiYgx@DzTlrZEtlcre1}5=DlG(>H;v%&=;JjH=)tO
z*wyQbXhD}Kglq3e3gqdzAH;TmHT!~P=dZx5<{Wz8Bsr+tjO@Yxcvpk7cy-B|3MsSh
zE6M%JMEI7Y*OQMBO=+d+f3vl@^(ONJvfZuxpKv7Dl2eQ9vCVqwxG;52$;-TV&cF8d
z;J<A;e;!lF`5-t0Y>wr$;5ehSyxL|q?j>zuzF`Xv!~Ke=J0A^+(zjmJXYe%&V(Nce
zIuyS@UIYbN&DV#J6ckC6>W;%q(wL9CX77B`ca`Nd8$4*n^-$gdlM&@tlkzdP)hB51
zGEyjHyIZL^ek1}@Yv4^WVKK=nXjNKg*v!J^nr{(PD|+VrM5zh-g6X(_xFqF@1Pk?+
zT91LjQ^3$7FI;+zr0HTluR5?7|6SxN&^W$VW9Xd<Zp?GZjN)vQV-&k0>9q#5M|)05
zXhD`nGo_lx(n~M8>)V(3U!pw<ugzcAd%4S39>jFS!S7Hzf$M|yip9T|NSqxa6sJO4
z6lo{t$zY?KOV^yz=qxzV#V6q{U=8+Mqc0Voc%LBgg*JHP4C`F@@n>Fn`}BH*bDs<B
z^Q{JV&%&^|5KOQ;Ppm|!t7NesE(|>!?m{!K(O2wTj4~Dl6i)D|{tkoZVb8C7rz^y|
zxjB8GDsppqepZc(;5nfn&v-;9@zKG!Vw_qfPBz_efgFVi-H4=6>m~N`-PQy9Ww+Ih
z)AD=Cj|pn<b|B~asHg>RdIY;t4!>t~>XCBVs=@gJb9KGhO2|9<tK5{w4tu4`lfl6>
zrKrCx&$(E!tHL;}FkgEp3`AF|8{=AJ3S`7zFr&?xR?O{<(lUBBMe`C~0~vZx8&&rE
zuo*WDoYoP`j+<_0#59?>Hm8VGDvp3VEL$;4jT@({t@S?BIofABO_M9z0von-MMP9L
zq~K~v`MPhyd-f6A8uxf}SH{+8D)eVnew1>~){HsBFFJ5Hbc@`AeKU(qhmy>G_yKj8
z*zjwdhABo;y%vhT>{ESJ0a=#%Q`B6#bqj4N2yt>C$w0{s<D31Hg|@pU4)?tEDq8EU
z#|FF<(XN8~oTYMag+819#l%9G&JVb0Xz=hH)QEubSTLyOh4?`xH(-aAB-ay<VKwUw
z*uc~R;8-azwe<g~aQh$a%ws4?a-Q7nE<whF@)V}3L+{z;bEgABlUrqHg<CbYUg$LW
zW&U@+Q#Lwjn*U@=yX`K67_vjB5}iHSKdW3mw>NXQwOaU^W7>a5sElkQ!BZN6fhKt9
z`1l@T^Q&=?7jjnr&;!N~$eH<yTLWmCDt?>Dzv34#f5>iV_tKC!Y(Qz}zPx_k@X~W=
zQ?5XeuGW<K?qm4!t+f>!Xbvd(6VC2bqo(YVB$h&g(#2@=rj395#4<^M0cJ?;U14QM
zKQ0XC@ZusrzS>KA5br0Pkw3P)f=nEKQI)2kr;Prm0LpD^%M7rlO@lufECOy3x{JC7
z%2mutV79$CZpDhc&5br;CXY8<m~_aCF6ggYoCTJgr@smD3oZ36LIT%B;+>G`ST$`T
z4(_>9>+2Jpx>}~>R3F^Xddkq7W(;<Fj&-6kpOl-JxM|OGlz$S2IC|F4v5po{r)xIk
z+yqa`nO1j=v9K~kC)^6ZaYe*C*7o6uYqoWKeZi+6UX0}67!RZJofr%z(7?w{Rkw=7
z;W3aiH@`{6;Xx2(VYPgR1h0W;lEw#1+-lIIx!cq!105n?RT<6L-Cko1PPpT=UvmN)
zAbW0SC)@r8{)$|zI0|%u2?_;TS0dk{@SHT?$gn2(E0X7Wi|9etMp`^=3%IUmun>}<
z63l6chsJ9$=O@XCPXgp3`uYUex7oO%7jm2GDyfWW+_n`Qu?m_`|5)+-@mfK%5qxx^
z^IjAPprAv)WNZSm^-D4}AtVhamst+%BWENw88lF5C2>d2s`H<4z8~>{xb44YF;w=|
zY9o?1sfeh;j2>)9_mi^d9@<*R^=+>+><K&xqVA6J0^`Hj{3vPvJqm-%s;g<%+vV1@
z-Euj1v_2BPojFyAI%>TPQu?vbZ{~>Ig{o-!l@F4M%^o)GEh=VVQ8#8xe^2q;#Njzz
zcoffl{|#|km#Qq!V4^mF%xJxTp+4ohvu2EvsJa=YNOL+(oX_v)vY<n>0^(bf9~`4u
zr_jpB3nt^Q(X2^bl94WiddhGOS_lrEdO<E|;y@P21HBCG9)zQ+uxMc<C%=)yN@cPb
z^(lESS2d6`^@|V8`aPxkS8-{Y<&B~^krkF>Q;_@h1nyXq_$@wo*hy8(XB>!K0D3dv
zjN(wOnM~Wy@hwm0<|$a87k~3;t?UI_G@VO@izC<UHC2SwUZS&#3ky>Aw9%>zy=TOo
zbD12{($nFcv2hKC6#zUJ`TFZ~DTpz`AvC@~#Qzw2H~AHVj71|tBXIL4+zPV&CB4QS
zamrlWS(mmL@@I+(;wW*UbCUD^Z~6O;!OPm>@jC{U5Iwpj@(21o@Emotr>qIc^^7_F
za*Lg_-<TpLMm9I6^wDpy_U}m~e{>tm2*d4bwvs|8r|a6zj}d4|4$&Lh6wP0tE?PNC
zV#r+zAICta0Ot#xKWw$^m?{-^5w4jo%3T6Qo&**Zai70f0uiR5G(p4`O5+Z54X~3<
zutG#n^`2yEc9hP0y!R5*UufxzQ9g)}hLT82w-on0o=L@pPOcB@^%gv!n#XaWdkm}r
z-;MvDzxBU;4D+4q$P&BR<@kH-9&Q#qR-XGR_`~nvE9vHi1wUl;T`8MatVO5Ziz4=Y
za0k+zuZw<^7SjjT61<(be1^jv0s7bgzBTy?_YP2@*`xuf)t*5GUi({Gj(9<0a?Tc<
zsY8+4J!1m_DEW}v%L4O>)`B;G`-@@0Dfx7B#d|a7o(=dgXP<ko@bEa203}}-ZKag|
zr#ubNZt_F<EimeYV!2Qyu=53BjIx5cFi!B0M@9d0?07D5qOS&t!4RBWd;CXC6e43r
z@B#aB_bNHr*VFo1(>U{YgkyuRs!`Pi)0gDCu`@S!&(PvE)<A~^)|--+jB<;u@7#OD
zh=(uGjP5wEV_#sqAIM*~ekKM?!^jF|rb3rzF!UIS>V%8$?zCW8ELvN&f!+S8BMQ)g
zclKrWVf@%nxZajnu!fX9l-7=^L@-IKM#XnYV}Svn5e~9Btq)k{&pr7HNzO)I<F$QL
zW;a(;EbF9IhOK!N-y#>hH(PpcDd8gV(O&El9v}C>@a*M@9<S3~`bo(ZZ%$U8qTND}
z?4T{MVnx6CO=|{ttFTSvs?uOb{1)ni1KU0mhz-;gL{fq@6D3TX&vPI(q5bssKHuOR
zDK3E93BBUo(Xwn+%=_`&^LHsaIzAEH&@V?XN$z*<4YCx$D9-SlT78TRBXPO)E_V1^
zE4EhwL~wx0Tux-+%-s;x*R&|qRmB^yk$wHMiWOM51+l9A<fW+LP%q#MKf;2t7Z*v|
zc<(vD@5{gpXgAnKNmDbkq6}P5?c#$fE5G--+TWgDAXF8$_BdQ>FLwMQbN`At0%}}9
zZulKI{Y!dc6QDKYBV#zc)^7sNKt3hs)44UNrDf!aLT7|f&mvFhN5M8@;DT7|(0k(L
z0BQHVH`gDnHtB$^5a-{`07g8Z`^Z(wba;H2T{WdQ$Qo96?b06#whd5BFS=R}&Mery
z#`o-vvrCpDi1Qg}y!t5k^cq`@Ge6%v4_Bzy!+&4SFrA_f#pJ#Z`~?<J&zbz5uh*)R
zotEDCSYp{~#!WfzY)o7a5u+1!=X4hu+)IwF5FOtEteMgof3FUH0gGRf5jce;nhZn^
zMa6#*V9Ln(<55K!`r%9AuBtqd67mMZZ|P;E#;GQ;x0s1qQY0`IhHQck_vG%$smfol
zjY_~@_OLC2UwmOgM7!Hr=a;A$^rYuknJRN6#~-ZXckw5K1UdD-1U5ZzA}@6VXF7j+
z`*kR!z&x_k{LU*h3_IDDVFcuECvhF|@j|4xOxEBsn<GvkIJ^UR%LNF3kU#$KYg)i+
zp!HVl^T^~!xZDZ?co;1pbqdhys3sSW7JI8W(hG;0bEleCvyUK}|F%7To=RLdjP|~F
zhXu<MtyJ;IM9^~y=#?arzcM!F1w>?%69C~M^_PR>_vIw3=9Apq;ZlFwv{5q>v~2$-
zViHYrB%d7X!+_)myQg_U(T(`xXdH35Ud4`R^I34L{3Z%NB2GQg&WBAD^t%_IiDCUG
zE?g{l!P24Wn4IbHX2rsRMy{u&;5xI$On|D2r5fM4Y^>9ZiYV@V7+Bx2c#^=$<{RH4
zJlQ(NRGP_+{-v?{(B8^9e|=G|vioCHaJ8YQqajQDh$~sl6an7d-E_iLZrGKYYK0ge
zaAD^9AbxvEd#!aLbjh9KtcHjhya4MA25Q@`e1*CO>KwLoD!n`cU#~CIo>!@^2D8nz
zNocY=!k{oto)#t?+A(P{N}yu=B}kDG17sM#52)WBT)r8?7GZFMpRPT8;4SFKjdzJ3
zSpeiM7Io;%`~44(-{s;Sed{j#;4ihIz8ba>6^<>BA@?2m2fC>fF1xHyAq(uOh(Ag9
zS1XV=iRPzP@@DlUB0;i=(!O>Xv02G$FLfAIei|Uj^7>57K~o|Qbvv}=hbA#d@q@qF
z`vKwJ=ofi|WzaP0&!P}~4+gH9Cder1bqpfn4`C9`E@O&lDT*^Qnk(mNH<*a66vf76
z-k6vvY1W^!UcgGB2`8;dx0PseHJDnqq)o#;NJe6dqy$B3rOo7uZerZcP}k(jC57%%
zGR9t89bag}j8vs~H&W6r?UkGFKlUhi)RPlyD{LQ;?g+JvHgL*t2z9lCf{r?(ga7>s
z@@s{m$%+n!t2MO0GOHSFr6w?vJEr{!C*VH5tmT;|tN{78?qCfCAPI?Pv06!fST+|#
zHjr>ruRMNv>CP?Kq7q~5{Iwr8=OmgXyBJL;^M+^p5gnaaz9&$~|0&l8eCq#sF82^f
zL18qV-Pa6@{sd`?11yDUpPT^`wH-ZIuq>_Vw*}ab+Iqkj5gubf4<42A_y&l>jsYu^
z6ht`XgH*u}NO95ow!DdleT?7;<M@%zhB^*+DXDv=T7#ttq%Xihtb-BLKC6Sfph>~F
zUZ`zTSj1*yrARi1-7p|9N4xpJMbN3)H~JuBC3jz(zWELapLSORY>k%7&^ylf5)~gM
zU;c19y_z9DS!(Eo*_`Lh<@t-}gR4-D9yQq}a##bt+l{pbc^KRNX=w0h^x2m9!ue3|
zWCLhA+g10SwbLiC%T@%q5!R}o2`GD>jud3Gi(&(xC+|E5JTxBm4Bg(g*h~S<Y$*Rs
zcq%A~&C^WNQ*RX@u%0ac;NU?!uH{euylZXYbK8WrDb0LgVSRHsZdCul#bVG`0(#~n
zv7OYB6=C8DL^#CxZmphj-o2;;n^ellsH?Kyu7gUA7PMt;$)a72;zg%YLr8`uIaI{a
z!~5O<d8q%oyYk0c<JUE?6MZxB0VG%**O(OYnXZI>Du~gtrsnJt7S_;o)78nBCUI_(
z>lHu|=yWM?rvzZeAi?tZW^;Dt5%-QSo?|;*ju(W^CFZ4}Gp)zpfi7_&$KQ94-x?c#
z&N2OZ9nP-<+VEA!DL0;5@IQ_MP4@Nrz}kh<()<gC3A)K8=&Bl1UW&7r)n;P_fic=`
z1`?22v4(7(!uClv+&yWJBs$AQxuSEDsQtymRcL^O#&9u~apqqD1sMAO{q0xh?tkwb
z|H)(E`~uw4z)90nhAMAIB)l**AxYRN$<hT;<DH^-P>y)jNYUMAlwM4fKnY`g=&{Ir
zVM&9LnUP6U2tQ%GA)qS}3<bW{-P^#Uk%JK_;J65&YxAF}=l)sX`|Aumz~mv0x7HM7
zc`T`_5SZTsu-jEj#ScZfi#2r<xSQJC@*;iYlf-M8*lyE@C?MnJ$A@iSB3jTd_4fED
z7D9OT>S&^*)R-BT78rcu)k%(|Qwin){a=9@7I4&!^zR;-|4?cCr#S6T|GnaOcrk3_
zff0-TW}st&QjO`SuqUT4oGhgiDM29fI>%gNc7+CyV!B@4!8YMJSMCqD9Hix`mPRx;
zqXe7Xo$tG)txmYn`?x}@NczW<_o`AHSFcCh5BM~QrSLWkMtUh?qGh>016A{XVn6zq
zqahwxR5fX7O%bM2H{!d!kxxNrn4!g<;NDU>%2;RV@!jrZL`FoyMkuHm^lP)PvBJXz
z&MXmW$2wj}AhLqsi@(6mzY1}OF@x=G=T;?N!0sZpK$^wQCXB_*U*W|!b>6E4RAzn4
z9XnFZ0h%N7zoUu&*Y`E>!*XU8f(u#nE<cvbi2@1*{Y(u8-|cTMEru`Z1Eca_`{Hgw
zGDwyYrzChnc8K2IowTk-5YJyz&Il-qbYlk(MTbHTrh_&#KBszbG{DZ&wl9c^S{8?9
z>}wzx?0s%A^RsTq-eX=J3YJSENwhRMrvDE3{-@7l+l}05eKh^1L1x{8^_gxt!0?E-
z$f4$0&QxOlm4(*Qrxa0uF^k85=v`~IeLue}!fdC&z}*Uq>Q&MItT1)sh=b0+ZCg5o
znALEb0G-|g#R=kZ#X3!4(b6j&$S<Df%bcg*>eIplTAr~s(L-fwe@Vieg1+?v{m)|l
zIcZ-UfP)pIME;^uhW{^~Ll~0_=8o$SFBx%rj7L@DA79t$RLD2)1DIVBw&U7X1Y7g5
z#v9m@!nZz|j`J1K8dJl;M_WjmRs0CMtMMZ!fP{V~ySvisO*faA5{P5njSNvHfAa1q
zE}+TjfmOpV2pIG3V_DQo3DFNS*KV=V9UmZKb%6FU4IuenH4&G1)ux$d#SJPbi~)dC
zp5>2Q_&;1DNGMz<WI}0UH$D*(2b7}bdZpEtj7M+lv7!KtL^6QP{jbjRM5cS-OBuus
zUSS*}zoFbF5uH>)W|!rE)-5>9Ox3XueLHMH27A^Z*B&6(<;D4UDc#;2Q3m~NR6;LN
z?t3_02B*RF)!Jl?HG7_$w+l$&`&gzxl@_SET+Ga~Zh9LAvV;RrD<k7b=S{eO#KBGe
zR>u*q+{1aywe{Tz7*VVVNwAgDkS(~<-puyXY`wqa=)JU<NU#`HS^=O?_-nc0znYAS
zW5MxDa~28DtTFp~t7th8ofd+FCt>|8D0{D1fKxUq-e3QOYgNl%4;Grp2ci>tYpb>U
z0o^#0+~5r<EN22@4>C>fEPOvQuyuMh|K^N%Qvz|%U1hoQhbFl9a5?EG-1)%Fi`$0v
zvtiYy$)hqLD`t&soP$8w3o$r9Wz4~MJpQOBC8;LRpHYajI`G7!`W!OVIBkIe;8lR%
zR<O&)Ch{Iv5PR(#Oka3omb-h`7RE8_I2dSI6PSMZs_NaUvMBCCXIjp>DPl|-4Pty&
zOYs`l=aYvWAg6BbzKfMRP;9X$-8`O!L4m^#>*)&jKyt+sHrr8My3-r?lPxsDT^L%^
zq=y!?+VRThoIHHF?*)hXdd&F;B%(GxxRYQO&lgU_Xf@;Or^{EW`T1IVs^{;gqZk+W
z=U%ry(ANZ_@_pDaOu8gNZSxmeTqgo4+?eG3&~pLHq}}w|It%3Hgstn8Mhhvm@e)tM
zpI_0^&j~23J&Re5m{)7pMH^$dOt8>W4aR%1o^>|fTOaun9pWU?zt(Kvjm)YWkY#$O
z7a+i&Z=rdWGS^dDjNwp9ZvS*hzwUOz^7bpvV-<M3c8c|~ql{}Oe)5~3lfBXRW0a0T
zM+YS@Qm*E!eZx`l{GZK!+s5U|KqJm-^aQ^R=n&#>9Cyr#z8fFd#&E@Ol@*XcDJc<|
zZ;(z4z$+_%YmI0zm%`^1v0{z6q&s92*;!J|dI{QdK`f;95hm6Wa|`GWWB=C37?v9l
zVAu@6(?G8MDXq&UNZp=5ISVN&lpW#|YBbg6=Y>-hBGE9_ERkQjzBeYNpFP_$<aDOZ
z#ZEC+<6S<$Qry09=+l>hMv+^d#X<7rj^!rtqZ!Ja^N;B6u&tOPbJmtg=NF%jYhL>?
zyo{w-_bqW(`|2S-=FmdtNu@*QNtu&2rFz83EiY3aZ*;5SqpQr)`Qe%|oCeV^D|Q^{
zmX7%lXyqI3zRN0N+jNH4kWRt3l+NLqTXe?{^Fcx?RbIc!uJ{6|ZmMzr)x7w3Kg&ws
z0q73K`8{r@!8+On>T5LmTVgseWFE-oAaPQN28Ux-wlT#ei}6<xzP{vW9n|OMy3Qkq
z4t^mAT#VWPlW8{2)Hj5j@O&RJ0Ez_4HQ53Vh!}_#iq%14i2udSVPeM-Jw52eCejU2
zZ3%ONW-fx7s0h7elVa}0@j&Zq9Y02kjM?vO?-Vk0WgLvZTgZ4NP%j1VVZsIj5YO>4
ze*BPv6)UI+MRyI}g^qpL)=^@oDP=wg!u8r(8`t~z6%i0b86?edW$<YZbWS~4y)zDB
z0TMY)##AY)8X}+{n~UeS-q2xuaeSmm%W*{@#zd20pzteD^BHj2ehA$0&<cjh5rs`=
zxH_E+xmmt3tw@%!tHsVp;>zd*874{;myO*fxr}Y;T~5q(!!Gu~7lael)N`T@$zo<u
z{{hhDT&&TZC;{lao1o_F9plD^TU7>VXiVQ-Y>b^f$&Q}2qtI~YiG{?y_#EiwYNB-V
zZs1+O9kPXm+U#rI=BQm{;EW*^xUr1eoU1ny*-KSO(PWQIvfp>GcqR6Lkz)>2Sx^@B
ziqns@$gC$<(59<%;x-M&`T2=+!=>^R9GEl&mlw=E8wXGsKVkR5`2Cr{WIhVo<})#w
z;;GTVUx{J;_SG47x(wi3=(UwvOO}z6toY;JQVJZ&tE#D5lAswMWm!C5RB_KZ6du?W
z{)|RseP06`vidv>e?qUe{*<zGcq{qY($(|<b{=_Fv>UCU`<l1l@|xvvOF$>uGwuEp
zfSU8<**y?_o)|chwwgK#M`F`=TeXe{#BL`R+up0*c$mPn%yhu}bNJ`O;592qvBeT)
z8x>D!?I5^H)+B6<g5C#LVb2?dyG+hr>-ZF+puw?!!|?*n#ea!2Tt`uhGg=A(b906G
zWzmWi3G26QIApjj22M28S(g|sA(}C3PRXU!hxDpt&%|u9dEQ$s8v##N-Rc!NY!?Y_
zFo7=D=dzZL8#_Zs#--|Qf$#Jkt?HCdCa!H3R+!&gY=<-jG&Ux!Y2#5kIeqJ;07@>t
z<BMbb8+Wl2uR@c|8>h8pohb%z^-(!h%eetMj^ZY!DQVr*-g_E^t;%&Y1asBOUfWgh
zOIgO2wjYP@?oJ*3#GML=^%fs(V{cwy#sEnOcL=nlhpkk*7zD)@Qf;d*eH^5=BPH;)
zli!Iu6Js5-eMk#C*0?4{pQ<JCn{)Crao2Zs#;NAqDu2D4-K*coiTu%X-q!m-I-0#l
zYOOw%x;shiI>6F_8Gm<%JARtXFYGe5m}t1FvOi+$u~7s?B<fR&m2-_W@|X^KgNhrg
zWli^NQsah6iZJ1Ko@4lElj$Aa+#V`aQR{NS_&MqGPK~;e=&ytC4vFH%pGH%XDdSN3
z8lgCs9j|gFuISL~O!*}Q*2>x)I=CORXw5>CFkeN_pB!~VEO*Z?_jM7IE<n$2Af38+
zV@#CQW3{{jQnvbnp)vg~ZN|jru?*zTnvT+v8^1roj|h4o!}jYg^IHiVmolV*xP@m@
zc2!2RHoQ5Mm-s#-v7-0#$Bjit`rH!rWC`ao8)t;<cS;eP(|$~%Xvr?)M~YI&Xd}3%
zz7~D#I$v68t05<8Az5ZljQE|7S5R>?r)r(|sOS=e;)N-Qla*B7{FVe=hTmAVU3#wN
zYo<Kx%h1@Co^EKwkK^HTM8A`OUqC2I6{bj+UKgU;ZG6X9qS+R}%(yx3=sClx89?lq
zY_jbyw#QXTO8f=Re57=pi0)30u$MSuQ)9m>cQ`4`r({El8|9_YuBs!cVZlYoGJ#9{
zOxM*M=VW{YPQzo!kMbiC0;W1!!ms;n#)!l>K;A9@34qhVe=EduNO;=sJ2MM*o|Lg|
zUcixtqiNC!i=HgIR~F2^v=M)S5mu`O7poxlU!l8yRs22w1&fpHnrLs;L*10@#Z3_=
z-TXE_JZemJw<W`7nRbMqJ_7Q$J&s_G>A-TY{rXxw;;MU<utxUM4-Gab$?2Z9Zza*v
z41SO`8CP&(HbL!4xB_-g02uv*gP(Yyfz71DTmfvg!Uy&<Q4pJT=W*fy*B*iV`9FAy
zHG=&}=DD+vnk`@lG@!_mgR#H80JXMYJ)}S2zBFP}zP~_MCY^21X;fz%Q-#*F3IB^~
zpL^et-H_bm1bj|bl-Bd47nHyHdw+gY{B`d9Z^ZfkebKMrnqS|>b8$f5`Zyanoc{0S
pfei2x<Y>LdbZMH>)dyaWFW0XDBn|)HBI$qiUH||78OhHd{|D-&jLiT5

literal 0
HcmV?d00001

diff --git a/sim_scenes/solar_system/halley_comet_sim_02.py b/sim_lab/halley_comet_sim_02.py
similarity index 100%
rename from sim_scenes/solar_system/halley_comet_sim_02.py
rename to sim_lab/halley_comet_sim_02.py
diff --git a/sim_scenes/solar_system/halley_comet_sim_03.py b/sim_lab/halley_comet_sim_03.py
similarity index 100%
rename from sim_scenes/solar_system/halley_comet_sim_03.py
rename to sim_lab/halley_comet_sim_03.py
diff --git a/sim_lab/oumuamua_sim.py b/sim_lab/oumuamua_sim.py
new file mode 100644
index 0000000..c86314a
--- /dev/null
+++ b/sim_lab/oumuamua_sim.py
@@ -0,0 +1,488 @@
+# -*- coding:utf-8 -*-
+# title           :哈雷彗星场景半真实模拟
+# description     :哈雷彗星运行轨道使用了万有引力，其他天体使用 astropy 包的真实数据
+# author          :Python超人
+# date            :2023-10-28
+# link            :https://gitcode.net/pythoncr/
+# python_version  :3.9
+# ==============================================================================
+from dataclasses import dataclass
+from dataclasses import field
+from ursina import camera, application, lerp, Vec3
+
+from common.celestial_data_service import init_bodies_reality_pos_vels, conv_to_astropy_time, \
+    set_solar_system_celestial_position
+from common.consts import SECONDS_PER_YEAR, AU
+from common.func import calculate_distance
+from bodies import Earth
+from sim_scenes.func import create_text_panel, camera_look_at, get_run_speed_factor, \
+    camera_move_update, camera_move_to_target_update, camera_move_control
+from sim_scenes.func import ursina_run, create_sphere_sky
+from simulators.ursina.entities.body_timer import TimeData
+from simulators.ursina.entities.entity_utils import get_value_direction_vectors
+from simulators.ursina.ui.control_ui import ControlUI
+from simulators.ursina.ursina_config import UrsinaConfig
+from simulators.ursina.ursina_event import UrsinaEvent
+from simulators.ursina.ursina_mesh import create_label
+
+
+@dataclass(order=True)
+class OumuamuaParams:
+    start_time: str = field(default='1983-03-20 00:00:00')
+    init_velocity: list[float] = field(default_factory=[3.34, 0, 10.718])
+    init_position: list[float] = field(default_factory=[0, 0.5 * AU, -10 * AU])
+
+
+class OumuamuaSim:
+    """
+    星际访客奥陌陌
+    重要时间点：
+
+    """
+
+    def __init__(self, _params=None):
+        # super(OumuamuaParams, self).__init__()
+        if _params is None:
+            self.params = OumuamuaParams()
+        else:
+            self.params = _params
+
+        if isinstance(_params.start_time, str):
+            self.start_time = conv_to_astropy_time(_params.start_time)
+        else:
+            self.start_time = _params.start_time
+
+        # print("北京时间：", dt.to_datetime(timezone=pytz.timezone('Asia/Shanghai')))
+
+    def build(self):
+        """
+        构建太阳系系统以及哈雷彗星
+        @return:
+        """
+        self.build_solar_system(ignore_gravity=True, start_time=self.start_time)
+        self.sun.glows = (2, 1.005, 0.01)
+        # 创建哈雷彗星创建哈雷彗星
+        self.halley_comet = create_halley_comet(self.params.init_velocity, self.params.init_position)
+
+        self.bodies.append(self.halley_comet)
+
+        # from objs import CameraTarget
+        #
+        # self.bodies.append(CameraTarget(size_scale=2e8, init_position=[0, 0, 0], color=(255, 0, 0)))
+        # self.bodies.append(CameraTarget(size_scale=2e8, init_position=[5 * AU, 0, 0], color=(255, 255, 0)))
+
+        # print(CameraTarget.targets)
+
+    def init_settings(self):
+        """
+        初始化设置
+        @return:
+        """
+        from ursina import color
+        # 创建天空
+        create_sphere_sky(scale=50000)
+        # UrsinaConfig.trail_type = "curve_line"
+        # UrsinaConfig.trail_length = 300
+        UrsinaConfig.trail_type = "line"
+        # UrsinaConfig.trail_length = 152  # 尾巴数量刚刚好
+        UrsinaConfig.trail_length = 130
+        UrsinaConfig.trail_thickness_factor = 3
+
+        # UrsinaConfig.trail_length = 180
+        UrsinaConfig.trail_factor = 3
+
+        # camera.clip_plane_near = 0.1
+        camera.clip_plane_far = 51000
+        # camera.fov = 60
+
+        # self.sun.planet.color = color.white
+
+        # application.time_scale = 0.01
+        # 摄像机移动 update
+        # camera_move_update()
+        # camera_move_to_target_update()
+
+    def create_orbit_lines(self):
+        """
+        创建太阳系天体的真实轨迹（太阳和哈雷彗星除外）
+        @return:
+        """
+        self.orbit_lines = []
+        for body in self.bodies[1:]:
+            if isinstance(body, Earth):
+                alpha = 0.5
+            else:
+                alpha = 0.2
+            orbit_line = create_orbit_line(self.sun, body, self.start_time, alpha=alpha)
+            if orbit_line is not None:
+                self.orbit_lines.append(orbit_line)
+
+    def set_bodies_position(self, time_data: TimeData):
+        """
+        设置天体的位置（包含速度和加速度的信息）
+        @param time_data:
+        @return:
+        """
+        t = self.start_time + time_data.total_days
+        set_solar_system_celestial_position(self.bodies, t, False)
+
+    def create_year_label(self, trail, year,
+                          halley_comet_pos=None,
+                          pos=None,
+                          label_color=None,
+                          scale=40,
+                          background=False):
+        """
+        在界面上创建年份的标签
+        @param trail:
+        @param year:
+        @param halley_comet_pos:
+        @return:
+        """
+        # 为了不影响 2023年12月9日的显示，附近的 Label 就不显示
+        # if year in ["1986", "2061", "2023"]:
+        #     return
+
+        if trail is None:
+            _pos = halley_comet_pos
+        else:
+            _pos = pos
+        if _pos is None:
+            _pos = (0, 0, 0)
+
+        if label_color is None:
+            label_color = (255, 255, 255, 255)
+
+        label = create_label(trail, label=year, pos=_pos,
+                             label_color=label_color,
+                             font="fonts/DroidSansFallback.ttf",
+                             scale=scale, alpha=1.0, background=background
+                             )
+        label.set_light_off()
+        # 只记录年
+        if len(str(year)) == 4:
+            self.last_year = year
+
+    def update_comet_trail_alpha(self, distance_sun):
+        """
+        根据彗哈雷星和太阳的距离，设置彗星尾巴的透明度来模仿接近太阳有慧尾，离开太阳到一定距离就渐渐消失
+        @param distance_sun: 彗哈雷星和太阳的距离
+        @return:
+        """
+        # 距离转为天文单位
+        d_au = distance_sun / AU
+        # 渐渐消失的距离范围（开始消失距离, 完全消失距离）
+        HIDE_DISTANCE = 3, 12
+        # 彗星最大的透明度
+        MAX_ALPHA = 0.8
+
+        # 大于完全消失距离
+        if d_au >= HIDE_DISTANCE[1]:
+            alpha = 0
+        elif HIDE_DISTANCE[1] > d_au > HIDE_DISTANCE[0]:  # 渐渐消失的距离范围内，通过距离值大小确定透明度（慢慢消失的效果）
+            alpha = MAX_ALPHA - (d_au - HIDE_DISTANCE[0]) / (HIDE_DISTANCE[1] - HIDE_DISTANCE[0]) * MAX_ALPHA
+        else:
+            alpha = MAX_ALPHA
+        # 修改彗星尾巴的透明度
+        self.halley_comet.planet.comet_trail.set_alpha(alpha)
+        # self.halley_comet.planet.comet_sphere.alpha = alpha
+
+    def show_milestone_lable(self, last_trail, dt):
+        """
+        远日点: 35.1 AU（2023年12月9日）
+        近日点: 0.586 AU 上次通过近日点：1986年2月9日 下次通过近日点：2061年7月28日
+        @param last_trail:
+        @param dt:
+        @return:
+        """
+        milestones = [("1986-02-09", (0, 2, 0)), ("2023-12-09", (0, 0, -3)), ("2061-07-28", (0, 3, 0))]
+        for milestone, pos in milestones:
+            prop_name = f"milestone_{milestone}"
+            if not hasattr(self, prop_name) and dt.strftime("%Y-%m-%d") >= milestone:
+                setattr(self, prop_name, milestone)
+                if pos is None:
+                    pos = (0, 2, 0)
+                self.create_year_label(last_trail, milestone, label_color=(0, 255, 0),
+                                       pos=pos, scale=60,
+                                       background=True)
+                # application.paused = True
+                # UrsinaEvent.on_pause()
+                # ControlUI.current_ui.on_off_switch.on_value_changed()
+
+    def s_f(self, value=1):
+        if value == 0:
+            return 0
+        return get_run_speed_factor() * value
+
+    def camera_move(self, dt):
+        """
+        摄像机移动控制
+        @param dt:
+        @return:
+        """
+        # 摄像机移动控制数据
+        camera_move_infos = [
+            # 条件：年份
+            # 移动的信息：
+            # 按坐标系方向移动 x:右+左-, y:升+降-, z:前+(接近太阳)后-(远离太阳)
+            # 以摄像机视角移动 f:前 b:后 l:左 r:右 u:上 d:下
+            (1982, {"x": 2, "y": 1}),
+            # (1983, {"to": {"ct_id": 1, "t": 10}}),
+            (1986, {"x": 2, "y": -2, "z": -15}),
+            (1987, {"y": -6, "z": -12}),
+            (1988, {"y": -3, "z": -12}),
+            (1989, {"z": -8, "f": -5}),
+            (1993, {"z": -8, "f": -3}),
+            (1995, {"z": -8}),
+            (2000, {"z": -8, "y": -0.2}),
+            (2013, {}),
+            (2048, {"f": 3}),
+            (2062, {"y": -3}),
+            (2063, {"y": -10, "z": 2}),
+            (2181, {}),
+            # (2082, {"exit": True})
+        ]
+
+        camera_move_control(camera_move_infos,
+                            cond_cb=lambda ps: ps["next_cond"] > dt.year >= ps["cond"],
+                            value_conv=self.s_f, smooth=10)
+
+    def check_create_year_label(self, last_trail, dt):
+        """
+        检测并创建年标签
+        @param last_trail:
+        @param dt:
+        @return:
+        """
+        if last_trail is None:
+            return
+
+        import copy
+        year = dt.strftime("%Y")
+        pos = self.halley_comet.planet.position
+        # 里程碑标签
+        if self.show_milestone_lable(last_trail, dt):
+            pass
+        elif not hasattr(self, "last_year"):
+            # 第一次运行，则创建“年标签”
+            self.create_year_label(last_trail, year, pos)
+            self.last_label_pos = copy.deepcopy(self.halley_comet.position)
+        elif self.last_year != year:
+            if hasattr(self, "last_label_pos"):
+                # 防止“年标签”显示非常紧密
+                d = calculate_distance(self.halley_comet.position, self.last_label_pos)
+                if d > 2 * AU:
+                    self.create_year_label(last_trail, year, pos)
+                    self.last_label_pos = copy.deepcopy(self.halley_comet.position)
+
+    def update_halley_comet_info(self, dt):
+        """
+        更新哈雷彗星的信息
+        @param dt:
+        @return:
+        """
+        # 哈雷彗星面向太阳，这样彗尾就远离太阳的方向
+        self.halley_comet.planet.look_at(self.sun.planet)
+
+        # 计算哈雷彗星与太阳、地球的距离
+        d_sun = calculate_distance(self.halley_comet.position, self.sun.position)
+        d_earth = calculate_distance(self.halley_comet.position, self.earth.position)
+        # 哈雷彗星所有轨迹线（注意不是彗尾）
+        trail_keys = self.halley_comet.planet.trails.keys()
+        # 哈雷彗星有轨迹线
+        if len(trail_keys) > 0:
+            # 找的最后的轨迹线
+            last_trail = list(trail_keys)[-1]
+            if hasattr(last_trail, "entity_infos"):
+                # 如果是轨迹球，则轨迹球会记录当前位置与太阳、地球的距离，并记录了当前的日期
+                last_trail.entity_infos["distance_from_sun"] = d_sun
+                last_trail.entity_infos["distance_from_earth"] = d_earth
+                last_trail.entity_infos["time"] = dt.strftime("%Y-%m-%d")
+
+            # 检测并再轨迹线上创建年标签
+            self.check_create_year_label(last_trail, dt)
+
+        # 更新彗星尾巴的透明度
+        self.update_comet_trail_alpha(d_sun)
+
+        # 计算和更新近日点数据
+        self.update_comet_peri(d_sun, dt)
+        # 计算和更新远日点数据
+        self.update_comet_aphel(d_sun, dt)
+        # 更新文字信息面板
+        self.update_text_panel(d_sun)
+
+    def update_text_panel(self, d_sun):
+        """
+        更新文字信息面板
+        @param d_sun:
+        @return:
+        """
+        panel_text = "哈雷彗星:\n\n当前日距：%s AU" % "{:.3f}".format(d_sun / AU).rjust(6, "0")
+        panel_text += "\n\n最大日距：%s AU" % "{:.3f}".format(self.comet_aphel / AU).rjust(6, "0")
+        panel_text += "\n\n最小日距：%s AU" % "{:.3f}".format(self.comet_peri / AU).rjust(6, "0")
+        velocity, _ = get_value_direction_vectors(self.halley_comet.velocity)
+        panel_text += "\n\n当前速度：%s km/s" % "{:.3f}".format(velocity).rjust(6, "0")
+
+        self.text_panel.text = panel_text
+
+    def update_comet_aphel(self, d_sun, dt):
+        """
+        计算和更新远日点数据
+        @param d_sun:
+        @param dt:
+        @return:
+        """
+        # 哈雷彗星离太阳最远的点称为 "aphelion of Halley's Comet"（远日点）
+        if not hasattr(self, "comet_aphel"):
+            self.comet_aphel = d_sun
+            self.comet_aphel_dt = dt.strftime("%Y-%m-%d")
+        elif d_sun > self.comet_aphel:
+            self.comet_aphel = d_sun
+            self.comet_aphel_dt = dt.strftime("%Y-%m-%d")
+
+    def update_comet_peri(self, d_sun, dt):
+        """
+        计算和更新近日点数据
+        @param d_sun:
+        @param dt:
+        @return:
+        """
+        # 哈雷彗星离太阳最近的点称为 "perihelion of Halley's Comet"（近日点：comet_peri），
+        if not hasattr(self, "comet_peri"):
+            self.comet_peri = d_sun
+            self.comet_peri_dt = dt.strftime("%Y-%m-%d")
+        elif d_sun < self.comet_peri:
+            self.comet_peri = d_sun
+            self.comet_peri_dt = dt.strftime("%Y-%m-%d")
+
+    def update_ui(self, time_data, dt):
+        """
+        更新UI，包含天体的位置、时钟、轨道的位置
+        @param time_data:
+        @param dt:
+        @return:
+        """
+        if dt.year > 2041:
+            UrsinaConfig.trail_length = 210
+        # 更新天体的位置
+        self.set_bodies_position(time_data)
+        # # 更新时钟
+        # self.show_clock(dt)
+        # 更新轨道的位置（TODO:短时间可以忽略）
+        # for i, orbit_line in enumerate(self.orbit_lines):
+        #     orbit_line.position = self.sun.planet.position
+
+    def on_ready(self):
+        """
+        事件绑定后，模拟器运行前会触发
+        @return:
+        """
+        # 初始化设置
+        self.init_settings()
+
+        # self.set_window_size((int(1920 * r), int(1080 * r)), False)
+        # 最大分辨率的高度-1，保证不全屏
+        self.set_window_size((1920, 1079), False)
+        # self.set_window_size((1920 , 1080 ), False)
+        # 显示网格以及坐标线
+        # self.show_grid_axises()
+
+        # 创建太阳系天体的真实轨迹（太阳和哈雷彗星除外）
+        self.create_orbit_lines()
+        # 创建信息显示面板
+        # self.text_panel = create_text_panel(font="fonts/sanjixiaozhuanti.ttf", font_scale=1.5)
+        self.text_panel = create_text_panel(font="fonts/DroidSansFallback.ttf", font_scale=1.3)
+
+    def on_timer_changed(self, time_data):
+        """
+
+        @param time_data:
+        @return:
+        """
+        dt = time_data.get_datetime(self.start_time)
+        # 摄像机看向哈雷彗星
+        camera_look_at(self.halley_comet, rotation_z=0)
+        # camera.look_at_2d(self.halley_comet.planet.position, axis='y')
+
+        time_total_hours = time_data.total_hours
+
+        # if not hasattr(self, "last_total_hours"):
+        #     self.last_total_hours = time_total_hours
+        #     self.update_halley_comet_info(dt)
+        # print("time_total_hours", time_total_hours)
+        # # 摄像机移动每10小时更新一次
+        # self.interval_run(fun=self.camera_move, fun_args=[dt],
+        #                   total_times=time_total_hours,
+        #                   interval=20)
+        self.camera_move(dt)
+
+        if 2040 > dt.year > 1987:
+            update_halley_comet_info_interval = 400
+            update_clock_interval = 400
+        else:
+            update_halley_comet_info_interval = 33
+            update_clock_interval = 60
+
+        # 更新哈雷彗星信息是每20个小时更新一次（哈雷彗星姿态、哈雷彗星距离信息）
+        self.interval_run(fun=self.update_halley_comet_info, fun_args=[dt],
+                          total_times=time_total_hours,
+                          interval=update_halley_comet_info_interval)
+
+        # 每50个小时更新一次界面（天体位置）
+        self.interval_run(fun=self.update_ui, fun_args=[time_data, dt],
+                          total_times=time_total_hours,
+                          interval=60)
+
+        self.interval_run(fun=self.show_clock, fun_args=[dt],
+                          total_times=time_total_hours,
+                          interval=update_clock_interval)
+
+        # # 每隔一段时间运行一次更新（不要太频繁更新，会导致摄像机抖动）
+        # if time_total_hours - self.last_total_hours > 50:
+        #     self.update_halley_comet_info(dt)
+        #     self.update_ui(time_data, dt)
+        #     # 记录最后更新的总耗时（小时）
+        #     self.last_total_hours = time_total_hours
+
+
+if __name__ == '__main__':
+    """
+    哈雷彗星场景模拟
+    """
+    #     远日点: 35.1 AU（2023年12月9日）
+    #     近日点: 0.586 AU 上次通过近日点：1986年2月9日 下次通过近日点：2061年7月28日
+    #               2019年5月6日     34.772
+    params = OumuamuaParams(
+        start_time='1982-09-24 00:00:00',
+        # init_velocity=[-2.80, 5.10, 8.65],  # 1/8
+        init_velocity=[-2.774, 5.126, 8.65],  # 1/8
+        init_position=[0, -5 * AU, -10 * AU]
+    )
+
+    sim = OumuamuaSim(params)
+    sim.build()
+
+    # 订阅事件后，上面2个函数功能才会起作用
+    # 运行中，每时每刻都会触发 on_timer_changed
+    UrsinaEvent.on_timer_changed_subscription(sim.on_timer_changed)
+    # 运行前会触发 on_ready
+    UrsinaEvent.on_ready_subscription(sim.on_ready)
+    # 使用 ursina 查看的运行效果
+    # 常用快捷键： P：运行和暂停  O：重新开始  I：显示天体轨迹
+    # position = 左-右+、上+下-、前+后-
+    ursina_run(sim.bodies,
+               SECONDS_PER_YEAR * get_run_speed_factor(),
+               # position=(0, 2 * AU, -11 * AU),
+               # position=(0, 0.5 * AU, -5 * AU),
+               # position=(2 * AU, -5 * AU, -20 * AU),
+               position=(2 * AU, -6 * AU, -20 * AU),
+               cosmic_bg='',
+               show_trail=True,
+               # bg_music='sounds/no_glory.mp3',
+               show_camera_info=False,
+               # video_recoder=True,
+               show_control_info=False,
+               timer_enabled=True,
+               show_grid=False
+               )
-- 
GitLab