From d9fd79bca6aa254bffa938e7d1fae1ea76975074 Mon Sep 17 00:00:00 2001 From: yan Date: Fri, 24 Jan 2025 12:41:05 +0800 Subject: [PATCH] update --- .../__pycache__/csst_mci_sim.cpython-311.pyc | Bin 0 -> 259443 bytes csst_mci_sim/csst_mci_sim.py | 186 +++++++++--------- .../__pycache__/__init__.cpython-311.pyc | Bin 0 -> 187 bytes .../MCIinstrumentModel.cpython-311.pyc | Bin 0 -> 1761 bytes .../__pycache__/__init__.cpython-311.pyc | Bin 0 -> 188 bytes .../__pycache__/cosmicrays.cpython-311.pyc | Bin 0 -> 12511 bytes .../__pycache__/logger.cpython-311.pyc | Bin 0 -> 1491 bytes .../support/__pycache__/sed.cpython-311.pyc | Bin 0 -> 9037 bytes .../support/__pycache__/shao.cpython-311.pyc | Bin 0 -> 7475 bytes csst_mci_sim/support/shao.py | 12 +- 10 files changed, 99 insertions(+), 99 deletions(-) create mode 100644 csst_mci_sim/__pycache__/csst_mci_sim.cpython-311.pyc create mode 100644 csst_mci_sim/mci_so/__pycache__/__init__.cpython-311.pyc create mode 100644 csst_mci_sim/support/__pycache__/MCIinstrumentModel.cpython-311.pyc create mode 100644 csst_mci_sim/support/__pycache__/__init__.cpython-311.pyc create mode 100644 csst_mci_sim/support/__pycache__/cosmicrays.cpython-311.pyc create mode 100644 csst_mci_sim/support/__pycache__/logger.cpython-311.pyc create mode 100644 csst_mci_sim/support/__pycache__/sed.cpython-311.pyc create mode 100644 csst_mci_sim/support/__pycache__/shao.cpython-311.pyc diff --git a/csst_mci_sim/__pycache__/csst_mci_sim.cpython-311.pyc b/csst_mci_sim/__pycache__/csst_mci_sim.cpython-311.pyc new file mode 100644 index 0000000000000000000000000000000000000000..236ae7cd03523b824e01341e941511c75f22c5d4 GIT binary patch literal 259443 zcmeFa3w%@8l`kyWk}O%W<@fs`V;f_GZTtcP7%(;ugE7xQa2}5Ek!6AW;7AVGl^gf< zjmWsIX`L_1lzdm|kkm}buWoPJ(l(??r_-jl$&3zjGI%tXTiOhhrt@+p&7Dp=)7<~s zM>;x3wgaKP^UeHza|~Pix%PYSwbxpE?JuRIBx~?}Cg=0crZ*EcZ;~VW=_8hR-3c1a z6B2a~lUuJtk(?yORf#dy<)5 z>o)gT5ZB;N8BFa-WqwI+>tI?>8nYYS>4WJ#>4OJN zU`0=bR%2p+h?(lH99-E`#X_yHuj*OF>}jx9_f#``I_xz)HO!s?`|6(6%%169Gg#YG ztJMsxF4w>(O93s5i$TlcV$ia<7_=-d1}%$|p&i!@C7w+@t~sOSvOlTmsWWL3G`loh z&L=fo?z8Bz=kQgp<|eKZ=4Ngs%vP=n<`!-h%&lBC%xxwO*Tk(xK5bkx zSBKy2EWdi>w>GYf?OX#=ba2|kn(CIfsR-Sdv|;_|kjv|-)`oTc(y#uA{e;^=9^TX3 zAbyja4sWk#bcno;xCR}U62eB#?sa%!51Yv39B!}u5@HT_%U)yLb%SlK)S_F?!Y zdK?5WV_$oxYsll}M+Y54-hIQ!pC#c2*EXGSaW4LCEo#?CM)adbr~mKnf$@aK8}-NU zNjY9n?odUr#;17}F`koSJ%Csf+KB{T!ePxNZFhBIIAPcm){WS`17V|Q#O?Ar&U)c% z?02~xLw0H;vh_+e#|~pk9yU6T!W)%sbRjpj@-RY zm$$ypN}C>T89oXLk#Q zEn;Ddpv#-=^pDM(vqHM8cM(g|M@ICEHp{!K;PZsWp`o=2>yox7kxSwdG3s=j4r5Et z894)71`bVAZ3-v1ckGMwhVEP}25NT+R(dU1`dkA?0W z=Doaq$i@$kdR;>fiap}s`-k~Kn|Hurv-f#N?e0k8KD)b*8jPB>E}3Ov+c;w9?SnRt z!;9oj&qkYTXk^rDW4;cIMNcG}cWlI=h{d`3yea?~=Jse}vQt4$IBdL~b73rVN+?p) zAW*vGB9}S7Rx~Er zC&!4dRqCP`@wJkBh|4;RZ5)e8o*f8{v~a{059}Va`1fIvjLvuj-r=ODn9I)l7??Yn z`ZBB~W`7LS9W4MdsO8rHM6p%tc$%tw0gSiF_BPE|mj>h6l>y@|xY)-mj@SRKty&+} zpB}Yyp0NI`$2IN<>&MtnA~$>@tmnC5kHgu-zBjizS<;inHfFKy>Th$7mFWXhqqG8w+MrdXGDeIR{yVy!KI)i+OvZ-hj=#J!oy6$Kmb=CJ@VXVByp2Fqvr=FYslD|Fy!sOQAIt!nnzm7G^sgzmkMGU(A# z4n11RqDM=4EEqJV6gijX3!4OEm1wLA8LI#;D`#foyd`aFZy;H;lm|O5?hUQoE>^S) zmiBiNwOKm|4otuCra5z}HPCR)Tof`F&F9&I>Flp9m~|o5-!auWUsxV&V}GZE9T&F1 zt82bsMX-VW*@Mj&>fp8az0@Uh+$L7E36?g>v7P1E{yWRjzt_KaAxUE?MzQ%y^pOhT zq`65`^W7{MINNEl@WUp33_nU7uaNeB41>cCuWfXMAvYMCK!NMzpd;OUCQ@C}@(1zb z!7$KBqlVuJ>v*nmz$ixaA5&ynPcdkW=v(S9KmJ%r+05yq;zhq~SXgm?jvbaT3QjwY zU1%ba+2Ds=PIEFE8wb+Z-l6hOT*hCwjK6*vf5S5VNz3>fm+?0(~{h^Jr26u3%?WS8FA2M-0pV*Jrg}INu(wp< zn-;alpSFBDI&CdI5>)<*9)-WW!oV)iSlHzm0K2?Sz%I`_*yVBT(=XR2Dt|@ad&{Hs z_@H$zE|y#l^{9@EB?DH+a_XPO>Ek*11TtYzwQEK+na`k3hx0{={SdO2H`}(65%-vM zz8|u2b|PlkPS`}Iu&^89_5i$aH zc8T&pNT025_z{pE00~*xZ62Kc-3}X}<{))=ZJeVYqz=w@VvL~yE;)ecV+x(=yGGWWWG2osYuHg{jccfP$V(=}*U94L z*W7;=i2Rx{GR9#%SS?qNKp%3P9(D0(6g%%^6h+TQRX8$7ezqtf2)%Wxwh@s*=vh_8 zRaK*Ly#uxr!=po-$F@q&N}UNt0+m#TLJdKg$39k*4+uvc<y);4ks zR*!~d>p=~*s0xv0{Vid5`MU?Oe}qN~sDI0L7IlDaW2+4`;}M56RMu8U=HcR2r$G@H zyJ7j*)F>1igxyuFHJ6Ke2xm)v$Toc9q@&Nfcs$Wy-iUSPPDheFRO;$%-NRnY58HM} zpMBKhuu0U!`tH%e16WpQPidS>qk}`zTDAsDki$0QaB#7i z8XJZZi4)Xg9xE0A^^dx30|11KDT_M02Xt5j4zaxCg#~0-H(OV1pg_yTje!L^VE16bbbtsAI<r z*0W`w-%a*Z^=JZR#xWF5wtMVDV{EikC$Up6qoIUN@}3yhvm%D|tbD9!VLiAA7}132 zL`I!NTgjp=Z|o2w*<S^jOPus4UW$Tw3_sSK->8gzkWU~U$WfH z$eCFk7!Wh6X0t9A&fOi3c zSlKe?ec>Z7uM^krd9~xUy`lRa6uTY@J=p_d1xu4?X$o1IK177o z&nM1yT<#K9-6O8LM<}>Quxu198$%XGzmYPzgD5r6c&Nle)q1gNy-=`Tuxt=58$y;1 zv69?DlvS*qvc<~fKm@%PKQf0bHVFlr1j}a8vN>egtjGf4<8x&fCxn7#!O|jHT0)kV zI}Wa$triNJ1WU7MX%1PMmFY%>f;EDrR$_JT^< zi>1L$yi2%}dV9IyzTRH`2$C}T6W<5hC5;4-5P!ZZYJ3zH&jT>bM%r$zzG6Z1L8g>3<=ICIS@ff$Pa^)m3vmUIv%E5~M#=8=X(b|; z`cSY;Tu|_cAf}Exj(jE(mPv7Y=yHjD%R^68L64b{XTJf=-&qDX5o<5WP>&B{8b1!Drx%Z15y; zh6(UMCrrL1uY3l^AeD2XSU#}T$8kC*JC{rle5QyG?3|gwO!mSa{Xj8DR*Qd3n7p(# zM!wj&pE$2^G0h2p?=!0h%Y=nXok&r^km30(z7#53T#8hm#+S; z0s&iTz+ndx7f$LM9wZ@<7%#|qV`Z}ct?S=xeVb?)Z_^_8Hnry4KSRH}=0SRKoh<>Q z><^x#@0(kzjp2k5E^Oohu(!|i2t!z;sgN+@;)!?=Hh?!ZvuyzLz2pZl- zVa7f?#A8s@u>Oqw5yx8KI&k2p-#v=y!hRaV+K#ZcUve~slPO}a1TvYTHT4cU>_cHQ z6Qt>N+XqkB!wEaqh7J7?br|NuX4VuQ;Ml-$(vfQf!-@Sw_((YAHxfYp^N0lTDvjso zumFcKXwtGmS=$8bw#nT$tm#u932d0nxuScyFJ#>#Soet5J(IiVEm_kCr}m1LlF1!P zI})d4Og|u4R!r`=X-S{nes0g?4tSc9C)?(YWLVN4-#)YB+`d3t;PjIng2fhW5G)nJ z)3aqydj-qtSqHF^$sMMqZ%heBI>Uo949X)gg0n zNLMUPc*R5_#;7#m=dgAzyAvdUF1^4CGvyhGXw+c z3B;Tn;T)#}bT=%5V0QDP4o4pdO&-kMc*)XQmSpGqAc54#e~c39Jg0GnC-z`i zI~q1{4kzz$grz;4#H@O~!?fu25gL$;^%mrx1iwVM02qUY(eS-Z)yeJ)c(XYi}=e+}1B01Nn2PMD~JK%6R z2LK2QF)&km)B+g>=JP?zqFI8*;_+War8{V4+yMiqNs4uzh9~#W7gs!U@Zx^Y1z}mgKCw$u3$~a|78y?UQ-I9kUHj?-ugc%#P04FP#N4o)I&tLWs5>EmkP5}h1boNy4qDMlz~g;|zHzELI}F%d(Si@hH*h9LJ3 z5xt-2+KfO<1zjRb34DZp#OS)9PelnyEKs5)4WI?96D3NBj)ZH%iJ+o+sC9t`L}(rD z?H$OXE81fFNT`zk2BJJnBc51Lzy-22W~<*j-4Q4|zw4Y&NG_agziCdL+#_wCh$DZ| zjpSW7yp)2U2^vRY3`ccB9^i{{5`yX&kU>mzCa9Oue#j+4SX2)QKm!vRHE~J!g-jq8 zzT|N7VVo$&x@gink^AJm-{0Pb`qK$0&gjD~Cx5hdvk^B38%CT;?t%$Qq9vq{i7?71 zDivx=_boMz7XhFG!PRi`aUXtr!-_kAdpRDf6vGd{j_)c0aE8157)oM~pI7tX8Hmc;g z{~xh^oFu*wD{IaYeR=nz=l_S4arMbvqCQKGY1KTJ*M|?|>m)I}SWC*e#Cjd!(c)-8 zZh{@_v)o3iP-7&tAeulGYKnx)f-9=fAZ!;Ok*SA8}xD539YSs|zyZ?mo^98S*|bEJksk zEu3^sojqsoVic16bD9w2QQp`UVGUnl&-yEz`-3$lD<1ZD0X8MKh;} zE!G*z0bgztB(4?G=s4J}Zn#{3b=`|wUSIiIy|AfA+|+ZUX6@y=t2r;03N@W#P3PsAt|xu0VaPg8jD5a%J-5rwN01zUp@NbL$B|AZNJd^ zfY|!LjWvyzyOG^iVa;A~&ED6KzV^_Kwd-E60(Q5yjVek%w-&|f1@tu!MWC08##Mlh zNL#P&e6bq^+$wCoU)+5E4OC{!)s-*S3#<2ttM?$=?i=+jmmdZ6eOlIx^m9uc3y@cL zMIeI$S_R0cTq(fRy`J-0sj%_5xbgUn>iWwiR}C*(h3eg6^=@=d^`1zYBGFiROC5~1 ztb5~luXubw*y9rSxNeggy^xYQwg0K37kZ!XzdU^Pkry9D29*aQ8RUt^vYT0XGo7K5 z4MG<6Dh5;j%<)j^dLer~oETKOGxvmwn}p0JIE|6e%|iBO7Ft?`=4g#jzeB9wfvBck)ORo!7Mbk$8rHg01PBlqq~jpX5fa0X-iPB1bo=Z= z(sG<1PRBSJNd6rDOCPvC3Q25`6R~K2DuMPM75Pv;U@i%RHgGzivJVHdyMdwq5c$#J zLq5(feM=8K89ihpQK}dr7;*u41&|XJF)Y2HXoPacHZ`51qYV;5_gKeFXAyLM1`|RwN zu!A;ZMpW9#@*H2G*i@H^^99;t$6-urZkCh=?X&jzqT*oNY};HLme1VvSc-Dj`^{Kg z({iRB4V)INHX*fKOfB~(&ZlS3MnVEb_tSExL16bj?^4I26nb_6f!s(O45Q*4)a>pVkylF3)P)sbtf&Yd)d;umo2S(wZBtf zASdUQ22Ru8%y@A1Y`?goMa*62x6J33EaiR6nlZgOSR`2ML~C8hSSOJtWm+A2Sf=wm z2H*c$*#BcK?9XB;iCs1S_gI|cbNs)yz{;ogr59NKFEKW&O}AMcUqd8mb^IwsmzG2R zP4fL~GT4fGi)?>G##=CO@?k77tb$QnV>pE^gK~h%xy(YC7hU9HN1Mn>SiZzc7$50t zv=nyUb}3xEjQu+=gjoeM`-5#lR;7@+O3Ym4ryH-dXqH?vmIS;P#&KyjVXl3?w2Icf zYQiWXj0OX1cy{5PB)yU+7`KbY?IGj#o6#$)7tF!Ea|W?uJspWs+Z$^y~>Y+#B znw0-;mam@6o6C41?@BfO2@M@$L&vMTU)}!N?$FU=q2oQFd+!&z9uT`8SkP4I3AJ8i z9@6IN>lZX?qd8aKu%KDSkjh5IQsLuR7$kgLhjOaLF^D-x0v2jK59~p@g^GKqOp-{a z)F2TP%3@TW_q@|DH)!3k7f+b9Q zOcIG%_>!0$RN{C$y2=QQ|8^k`Cp(`|B! z|L}Cyx%D}tVFu7ksAy7d1qZ(09 zy&Hs~9*goh5sGa%-LpV>pTmFYn*d7@BNd}7b_bY7X>2rdOx6uP@;Ct7EHdhpXAg-; z8oQ?_Ne8mAiSdq+Y@{(yZjc5|-EVJ48Qv!J_HEjQ7#CFPZrkN(zjsTP$xqg`NkMt< zZM$UPDG{ke7(?$4fVLGj47!F!p$Y-Y7QmDy=11=DKLv>M$H zNv(`?8>Z?f+d-3>+B5wK(W4v}+_TL>NxfjHpEssVS+FLjWP^+)X0Hs63E6AO-9WMi zss7RFW9KFWXbR9Raphj_&u3AJ|JePBQ#8s*eC}?ZRze=C-88&`3_@~wXm``_Mlzt= zO~aE86za}xll!-~@n-T$hY{BwOf;tv-9eJK3lvTp$G=B+1bPH>RY+I$0dwpSit?fQ zM5|m;eVWU0;%AZeGc4JD*iJ3fNJRaqPT;f1o=FWtjVljdjNC>V!aVfYX?{h)4YzBr zlux??heN4lLTZ@^Qz2xGZ08lDUzKzJxM$kFB$zw6xQkS=#|Muuoc9MX_j z&`EemAFZy-t*gKTAohNWMdUNN1f|GE{GDmGit2`}9~1|&+R$^c47@mL6$Vcq2QMxT z^g;mph`=j}D<9DgbbrUJ z4QtoF_36L=TJWF0Key#q+Ku2rq)^vEL8OMal5rf7_%zr+14`h$alw#n_VT?0($Apu zv&5Av1LtiV804VGPouHWi1;?3f_Q?ReEnp{JS6UpLwq;#lLU<5<71O!^J$rnw+D3R zJ7+rSrj=iRD=l-n{ z5g_i)2ZjY^v|0kFL-XJoy0+j`C~4vJr&uj`Aj~wZRbd82m?17KNfDMbRqIQL$+gC0 zn_{X4IAdhGY`bkt(2{D%QYrsc zzfJxu_Y`IBAn@L%4xA~{Mz-H3XI$Sfp|A9ul`RqmRT+{IVu&mmNvaStLXx5(oQ2Y= z!Yx?x(m1Qq4?Y$yQR%~SWGK%>iZ8_{$qq;|0t_M+sPyrblE5TsLco1r=!djdd2GA) z9kzkgYlCMWoE`fPUrJC2KpO@IJ=~y&52IM@Vry9I32V8qc91Q1{71DAN+$9f#j~kiO7+Ms#$5NtWFIl z^7ggF0>+gOEzVFLggrrRIe?B6_}mF))tetVx^-b^xn()-D?pB$gw{_#gA9}P5IXuh!F!p7M| zs2dFxuD${CN5f3>%!Zkb=N_GU)UU&e4oTFEKv^JTX7_muuGp#pERi5<`p~p}`p}G# z*8FI0?bniWLrJ+yM9)Dq3mqK}o;x!=>YorT#dN>r1AGDp&!6!>Dq2c$=t!|nYp2Vm z6Q_>hc9|B%Oq@A(t`~PHQ&wOPGpA2ip4*H&EoK~J1ao;vC*7@7>DD4~7lu|X-Yn>( z@*Mt4--PP&NL<)5ln=<>YU5>ZF##p6q+kETnHq?g>$ z^bu&v^wH8Npb$)!E&`Lma%RVEj7iGVR*Y_e#b2>e`jXU$!-R2}RKTn0=plWJQV~bf zLZy*zE}~D!p202yB;#h|;;>8>X7(lf%#3E1HQoeSBWRj{;$$yB23;2(NC+_%G7f0H z=wr$yAmuYOs0lk^E@=Z~k$sv*>51zxMfWZ}^T= z^Yv#+Upss4CFkQw5B}B3!{u9@|MoL|##h(9{E~C<;&aC)xBTH2XR>Cfvf_!IFF99a zeEjQYk6zm1OqdP+)h|E0_9f>pz7u%jk^lYtRcG_hx4-cR{@g9j#ITI*_pysFIt{P+9{AZmJ=y965o23$ z??*N}e|_fL@0DkLuGNVlnfLx@KCzji@mJCx5H6_0Kh3Kk5#>$vwY0lRdY^`Fmd%g1a92MXU2yKmLdP_kRB6t4_yX4o&oyly7nV z>a#EIoOJ)-sx$quAEZ6>XJuQQ|NQ>@>$$JZZbm&vS3LGN=UN%``iy;F zY<1dDnab@0o1NNk4ZZf%bf6W*Pq_Toh41vWZFQE`H~#TBb%Fm|XT~m=&;{@0UjE9; zAO6it&Zj^9x%r7&M>@ep&zEGb{h_=u6Ju9QqVD-h_Jm z=JyZo{o>QxUUL4+e$$>Wh5qcSb9mhHoijt|g7?=a|6r&%j0Ugy%(;_aM;H9M^m9qS zeB~Qgo$tT-7tbec_;{;x;r)MJsF~Y-)%lxK>o$i@qYE^e&AFxA+N(~(pZ@K`(|>xt z)%nZ!zvDPx^EaFM|0Q*St~>SV-={7xX3xEG@wu(valZZq?>&2|3s%j1=Wjp123_!q z^IwjwNPP1pXWh5<{#)~DbU{Mio*L(4n_qJNwe7v9%zIIE%fJ59v4Y0?UUC+dJbl&o zt1oYH`hWDb%!eD$1;07-M}L3yg=<%x4_$iw>!KE2@cyy)zSxuhmsg$d?fS`~cfO7; z&}n|K@4&>yR%g<~zgM5{daTv??)zNxx87*K>P&p(+=kAz=z@RS`ul$~@tXsv@ZWR0 z|7;iQf{#79s&|zaUGU6_@AjXiE*SW$%O@@xs0$v-(}^6q;2&O0o&07qx}amvOGBo# zKYz*jy1!aHVx?JNc_D?rgf4hJ`Tn#oHKQKd$ye_mvb4SA)a3cr{M#@7^s4idFF0R( zXOw2a?2m>`A4eA)_{anAZT>TKfpz33pL&A-`c>z@4!@99{uh(f;6LbFmC$n4sr!Q> z?acSl1@F8+wCeL6yRJGlX}XlJ|2?|kH}8M@>t9$~-YW4)vYp9)v0wNMb-|8L{QBo7 zuA&QmwEpTp?xrsI(T2Bw{ryIC!Bf?FH$wKGzT~X<_&@yUpQ#HBH-GZ28;8&ZjXyc^ z)inob79{uX`+r_T7kp>rig@OIbir7n|D{juM;9dQ{>XppIQ!>UoyUIort8%-bb&6d zajSWL{;KnzKJn>)`RTXO1^@Q`mya~3y!=0&)@RX0{p}Pgvxaatf7lC&jM%-vgg&UEk1UF9vTb`g9gaRm#EGGwQJR@05@V;!;pX^mAfE>a3w;O%Nd9!m_}n8?k4)~q zWzCFOpaym~SZgJRbV?KAhS`U8;Jd+`-0eU7v3<9q5f8G6QIjGbVi6BTA`*)(SAuz% zMU9#i^$3f4Bpzlo3Kik7ZgIrBfq8^w9yJN(y)5&4Bbn0$X|uoiT+(#Hx#a2QbJoeW zn?|$0f45=GeI-gKHWFlFl)cmJeNM# z_LZzT`&V+Vq+e;cV*h6IH*;Q2f4S|Itd~!}TK0*{cVCSnyQOpJc$)0j2hXbSws zuwaG_X_pUMM2;~KM~2}cNlNNmjHU-X2PJt9|D{hAi%QGB_*fGr#tSr=NnXq7)J-I- zpFvdyKKx6eeXFEyrT67^4wX}FX)quu;!l#r)0Nwp& zV!FNTkrqUvB=bZs30vFxNu%Ns(n~^jn!|=8_Mv0#hoy_n(ASXU9quFzkNj@R{=??>UDr$v@_ z11z#>n)F;Tt#aOy?cex#|I83bl(RKoJNm*yLQOj+XZjABoG>vtVNQ0^I%J*Bq4nnm zNUlex+s+-II{s9`1p{<7m@b(73B;EJebw)o-X2H@(RI)pDVbtQdGK&3r5cx@lhdYZ z1BRgKYdJ3zhmyAn$y>$bt%zjJ7OfTY7OTJcal=dsj=;g@=M9%rgz^@YBW*p3k+vSC zNLxR-`(|3^6XwsE1G|N^N)&#w6UCp}>)#y6zfc%b&{*XGiMmp^)Bktr2q;N-xiy2@ z39VKG1tr?pRF3DVFcy1)Bs3rmB>PqXW?Xi(sw4-UfG#T#4#SV`L=t&#ArLcs) zG0cdwag0i1wB30(`i0&w@w6Zbt|pgg^d(T^q5LF99B0?0wrBh)y0y!gc;g~bY2 z+{bK<9kqKqj$qx4AJZFI|I?b}J-fz+TlcHjrSvs6_EWvlX17VwCs>?ZG!EtXBZxnO z_!zh&&=X0z8Q0Koi#m`tnnu!AOX4kFDgxz3AIC>>s~wK0`9vB*IoZ=W7_Fp|!qFFt zKV*X%POma{q_VZ$a-{vxx7!hgvFgO|uv;0QY2zh5nsPQEJx7`?(p+`64oANo3I;aX zjzGPJ)L@a>9BVDKwBzh)RqauD=F5swZGiYPp0l8*EeulAhjnh0E1Wp)Vr1b*u)pA` z4hKJm>+cEPv-pqMGlE|{?3BRaS`n=+L{-Mm7~bb1ogD@#UwY}nTH{%5jK`DOaRPu`+t-v`r;6}aS*zpNoVQ4F#9Z=!I`Myf6^k&zRU51 zMU`58E>3MX+V~GAGL<~f57j^2i~C1C-h7!*QKlGgK^m!xKJ_Lrk)>|IiELknat!-2 z`xDP5oK=!VvwYcvK#yZf?Mv)qK9Y00=*z%%VYI&< z-7ZEv{aCSp{<6pV2*HF?K7m!Wn}maQ(KIMpP*4gMHgS%bmTSHw)CzO+mA_KGxi{k34AW< zgk?4y#M=k?OKh3SJ(iW=t&LM_fBKoR9a|J{B=L5N#2+X+^aqW2(Lta=@sqJb!Wq+B3_< zw6X-?taigVWYvNVMy0oSbV!f>d*O(@td7HpS7 z<_n4f$AjJry-@?ipkDfD9qVFW}q+F5J-9c(B+FlpCv+nH4U5|b|?vjY2=b%nm8lOWX=TB%q7FLaAue(oCRhomjcturNV?% z2F!G>p38t*r3@|;u1qcqW)^$THk)3u?aATRa)k)b#cQ>io;(gO)%N6bB`^!v>$HXJ zb=o5KI&Cp~owkI%PFu=er?s)yY0J0+*E#HpPm&De_7e)#Ya8}alGN1s2HRy|SzCY@@TQg&HA*NIV=nVg#r)> z@FpFKCmmiGC`pOf?P96!OAzWo2-UlL&{Fa0tk|Ps@RvU2DK^HzRii0T1tOnlmGX#7 zpEbR7+K+r(nG0IJl&?*4aiV3P*K)esgzIk;jyu_a0)4nJ@hT`x3~rcs8x*FMOM{ut zWx&KhhKV5zGndPQna>r#EaZw{7PFWoTA64j%~{qe}M~X$4zqA4G=89K)UoyVE}6oY>m((6>5&wUzmdk&p3m zQ4hDtwvDet&$vP$SMc?)h12j*gZqrb?e4SVq6gIFxu$!ynsD-n>#W0#*HEy%w_~i( zEdc)1TmO3Looo9}8~|L$Ssb+o9dEIIewp7EHj^e?sFovNr_}pt+t?!l{2I#0>>L>O zc-u$4!~Ok{@O9h5$;_9jiFDSs7ya!|*ClTar?Dq*AV&o~m5gQ;){&Y$5Kcir>z1U6 zjr7*@P`}HGOTbXmf#W@)XFLhbgj2C1&^rS>&_a09OVYA~??4QR6ZXGhgP8Y0^7u09 zsFaG9RA*!>SDl#h4BE&WkpLTcs;00KN(`*?xyAUg+$gLS3pWdeTg1XGq|6|X$fzm! zy>6t-w9kO~lU3FjGoJ68E4$Qxxqr5w`Q50g1@C6%<(0E57ZMTimcp&4ga%Dr)9k$q zNn|%lb`#m5ZDQ?g^MaY|7RjDM_EgRK&3Nb0O7=8OL(80QA)V|Qn$>G(D;F}!o+ZJ_ zCVP&iv3V|IA(!lVnl%lxEerW%FVHluo6A}#BzuvjVcnd5p_uF?n)RFJxP?-(uY#@) zl)D=CcQje1reBbOYAFqgA44ls^dlM?Td8isvK&8V(P|Xk3~>yHVy*#0Hkz!IUNDZHvNHnQ#C=0ieCUj&IzJNA7@NEwSD$yCk{2`O_Cq=LIFQ`RC#vZ}xjS*Gkokc#fIOgW1n z72hSK+(nQ|?h;boB1p1|#}C;~`HLXg?y^h;iy)QVC8WYdkXGC!q@qQT%I^|V@ghhS zcL}Lv5v0l}B=5ec4Xi7o5Rh&ydFR1N+yuA0f54p#x4a9$ZH8N3%i*@bEw9LMr@$?* zt8n8@nUnJB2)7k(d5wcR4Q_d*f;$~R}fm?2U zqHbW{sfxmvC>!=yfo#c#lLOL|taz4IIu49VFFB}=C+sW&t_nNDD9)y-93sZ#5?hV4>Dxy z=Z6Pb2n1t@UC)SbJbTeqdN%+s^g?LH)5p6`fY;)RNiKLsM@DdeppKCnaPy0$VaoIC zj&>eaRCrZ&ovm$XEXq5IWQ4Xw;qkgIJ@9)5PZfAk8gS0=9PT474o5$^@OIS^kv$7Q zUOb8PMCwfkyjZ(N=qiV(Rr%-8v5;CDaGkJ^dc$cWkVC?wn{G(BjnI+Px#`FIKmEke zPq#YRfbjqEws0zFkPtYC#voFQ2bvxZks=5a5rT$>T^Yz%5P3-+*RMiu<;FSUSU52nFLNV%Qz>++t#^@ENu;qx0Ayp<-kK9BO^QQUN7RG!eb)?u6{3A zL0(x}nFvEhoJW~DIJ^}C8zcqoXUD=8552QQTH=T?WMwX?0U(;<2`7PzNC)$9Ea!a@ zqSPw>>#VIo_XBtjKSwSDSeU3pIQb}zyd6A`mzj+4V?42NB5eFu5swk3=&2_>(JQ@; z*@ga)9?0Ae8@kltB{38{`NV&h68te4-y`Es$oM`PKOlof)yBaUT{;q1t-kn?>s6bsRc)FlCAPtSu)J2)f!npJZhCj%sE}0=%BuL~tz0nD zw1}@XS=ui>`o%}_MyK!Pqe4}eSk)Cvgmfffxlj8``W8%70J5YiNrOS7V5}C6)gdEk zld7pFC8v8oyZ^EMfsPAZvv@*sy=YmFBGuf3BGuf(h>L+D!MGA&EZG5_n71lay-Cd5 z95T~CB=>htbq2BpO9@%X-8Iz}pjWoYG7qwfl;$2WZ)3k6bTpxmL=a z+UyDytVRuOwWy)3R?-EFt}tKEtGkw0H|xFZyMlLTc8htt0XOA?;1<~~?7d#pc&(_B znz1WT2e8F809#x`VDoR@HbQRU%tx;0)?Ul4o#igOuXGBzJH_0c$k#G|+q550=A43= zQ`d9WT+3NAd*t#1S2hSa9b!%gz?$ca%P!o1y}0RGF?)iwD^QOzm#jvaOIA~v3qGjz zKcHhipd!Cbjvr9RY&yihZ;wZv+2uZiYU z3m46~g+>U4`wdL&_Dkh$=7d+w$=lcJU&$$N+o*fx9xd57)^4v!`fgPM?BA`{lD#I8 z?DZ7dhwVcCh-8T+U-Rx!_%28K2jUtdS1DB8s{9r5QT++Z2Qe?llT;OqOx2|bxm+r^ zzl;2+mhAg?%1^^y5#+SY$i!E#T2i8VV`}AfwBKZ3!!l3g<$E3JeBd?BM&^ zh%-omOu{FITiQ7e+;H%yqugwX)9%?MF|}#$kLDD?lx&1@N}}(Z1D* zXBLUswTQfenX~gL>0(N8z$vDzoYje|HeAUSSG9#K+h4E5#+`SNskuK$S0iAjJ>XVu z{=ZYAo5obXVKL8wkSsQJz7C6;XQ6y}3a32fz(!o@Q$9L{IPwI;;EzVZZGZW|3tE3K z<;Fh5q1@TGH2c6V2Z&wigVJC*Nx%sfz%mGhod7UR%SGSQ=Z*T zN0Kp|$#4<_aysJRk&GQbI5sobp_5WF_PAVj79)DPlTM?oTae90go{5EWaqsGxF2N@??01y~!6wXP+>7MyxeZk*n{4%;e846;vZ z02@wiXYdA|lQK8~PQ7h?%i35uJk_x&9=5QS0jD0k#(3?}|QGUBJSwHL=BZWoJ4kJ617sb{r1)J7Bg z2;i^-L~)r|v}X1pu?VjW0OteVrgGz{3gFDMFxJu5Mw(uYpG|ozC6v>0MfYNgkh4QD z?i7tXL&lx+OHV>s$R-&RNIeOcW*m~>jbb-b>{?&YGNzL}5tU4*FKD0A#ZMq5GdngN zRnQgBOso7Axh^3UuKg__D;Q-k59=HR$kR;wn{_%rHD9Y8GxQMy`O-mU=k0Pl55M6-yEtqDCM*(n6Mc2po;!#+beOx~TkJ#@J^L`}a z4M82|4G99HV%9VA2iyERX%W~?Mh^_!0pn#+gFe2CNpFJy)5mwvtrfMbK^`sL?NK^a zXRKw%D^+X8a?u)S@gSW~SUfB^nM?xMC>{p>N*V! zA=g6~g2_u50!)wwVB$D}3+H=q02M7ofgaJaYBo!(z9+P?ORU}}SoYy~m%X1-!1SAL zWn|6NERMTRtj3F}-H4lcfa1dRCjoM1ZlE)0r@z_sAa`N^T;~-#{k@!ig?nND>xW*u zFLbn5*!Qrw@8OXBgt)J7zNkFd#{NzRJ1*>-tHLF2`a5`~;)UARGhWLJ9Xc*__K2N5 zxUVgCK8RXn_EN22Qog-}_au>!1@M0fACVDT^b}1?hL4< z6Vq`rjkqx=b1RsnM65iD(-EGjTAT391$)tHV()Iq)u_X1$B3kj;XX;luqtJv?pej^xg4@w5Z`tBq>KE_+xQ0F=$VjKZX2cUYv5A zKsyu)0UxZC(x(`WtdF2a3IX1-XgFfW$=7~~+^gEE&d1`{(NU7ecA zE6H;#D02dpY1AzT39Slt4bX_IgSXAkX9R+l;WGghi<{|M%~a*>(}Q2AUFsYJooZQM zB&1L-JC{vs6Z~`H{BtAzd2#;vG!rHXxvy&R_KX9I5DjvQE~f7vZhDj)mHAVt=v~z zd6@qe5SHhy;7z1_GIHLj~eKI3a5dIJt0X*+bJvwN4V&DOo+nC)=qzI(gc zj&$ztmQ@kJ@~cE!t8KjC2pRyK*^Psj9mUE41{iIK%Riu z3Hv>bEdP7fttBY;tYa`YG`H*v8+LAJS=)7h#cxL8I-}6W%g~m+Z1_+=fCKIB1#2<- zTQ@W^b_7@AEI_!2&-9MC`cAnW<7?Du$Q7w;{a~LfVn1-~NINC(uFmIQME~(`lJPS# z7$*E}veD`vHjH}vaA*gg(#!WVo#Wr2u-}kD%XQ>2OqU1R9iVLg6#XQ)onbSM?e~s! z4LhNwMAFSHAF;!F7xalrT2FA19H(nsfP)Vh#8L%uS&#M30Vp^Jbw#1{nXm}H(YxY$#o@}!OM z_~E~c@bL$*FajS>20nh>NG}2g@`JktW36bc4H;|kQs8xC;Wc9+E>7W^6)sURpZsga ze8E^K8gT=w@Mcv#9slzJ1%k2SW=)IaC>D%WH*-pW)0*G2Xso$Y!+2avOfB^%-ptN> z;(^aS5UdfhYsKtZe=;pL)*8qX^Q#4OjcBf!b%e}K_;&+gyQhzc*_DE2 zrD$0hvaF1G1>Hh+qhM(gElnYdB5LrIkli3y8bwQE$kM2YznuR`JcuAz8bnJ&$kK2J z@l|kcrRFFaC@c4guFrLa%9@3&7BQ;@XBT6N|5V@}!B~l>CPR2_GBbN-ORz#*(HLsl zDXiEzzhdP@V(gXI&DS;1Z{vI;J@i?=d46*{{dQl=%iDASUZ3CI9Yn= z0Fre-9YA2Bvy76IcZ5y@U06FrlyYrs0>(NlJq}HT;2#!IVxTT+-z@ z8R(_u(c{abczwl07#MPtPpAJON|hWbm5f>bhbWbqN)>ssI&73S8t`_aHVDQ<{9J6G z8P8Dfr@(~j$Bn)CVcRcXiZ*BZO&QW4>3B3rn`koywUDunjP1x50y}JHqpRW)uW1zp zRl|r}1B=EWVJhWqu&|Xj%Ee2p?0&UCupADZ=*NV&4bX&#i3J1Z zj|be>2He8HkT@_TO?Yb{Q7{%^)>XTAhqq=7O2Y6&k}rwzxi(-Z$cKkr4vF`K zQDg%fD$044ju}w}7%}b@u4zSBr!k#D(4GN>Ns(`(4RVa6CPUVeBxfZJ_n|?%192%H zaC-2`AT&riBiLN9(NLEY;NTPcA_*`9B}OkUr>z0w7ZRI187?{1BArGA1Jwfn$r7)Z zW4$e=jDr!Fi5<609l}YNtB{q zE1BSTz?c6$GW=wmg8_Va5M}kysAA=XA)%;AXj3-V=}r3>f*SGt?=;7WU4=QGqshpb zDGg?b=@r4_VtT#bFb|LV;0h7KYlC7&GvwCdJ#M7YGprfGqqE26oUdv_X?ujUJ!0A( zXq~ra`7IKO6*ZAR3@vHcKZLDtIcW*$0W$UDb@U<#NQO92OeApniA2sYq2rP!^boq% zk|r5E0Rv4&Oktz+SY8qom?tq^K9Ul1#n#4|Vq{5{)~Mvps5ePQSdk|#mmEt8W!-Ee zz;ug+TU6mGk?kuMo~jC$G#1Glemt21x}Y?%E!u}m;SqS$itK}5G!GdH{B8%A1kY)x zL}8d?rHa`&g$^d4S26mC5|&rOT)GbFG5FBI6m36c zx~1!&PIu>ZP^V*c02wUfi(2T0l6z=TKc{{`m{oaYkjqq)3_CQ_+6R*+(x87y(Ob}p zW%4OXB0ZtaWhwcQ6CeaG5hq&}CsQTuLnmNq#3cAKe3^0cGs|aI!h;Sdi^`p|EUau_ zP8_UUJP?`7=5iLBpLwdZGD+E&r)~{s9!rnR5Sd+|ECZJl8x2yM#M~44syt$K&BxZ5 z7gxRlUqMhVMYL7uSG{&{S&>WKd0egvx{U5~aW!DIEf~u6-C+)FJw#leRjT>4cFIKh z*O=zd{O-5E2e!rAX;m;<%CJ??iQ%of2R?l&?3~>A*7P^u_(J!^>GrphKK6L~ub=;4 z+ni5`J9oa9>fFqKiS`ujE{2tcwTC2XmbO!R2|bDSz1|@Rs(T0c;ZX=cjCwC6#$m9M zI&WI~rKf2~V^Ob*Qoi-Micb9d?pCIcsXA3cV*i<9GQ!bzvenS`7l-|>CO^{7M0)>7 zPZv|m#CA$t=pJ%-stegx8P@f?N1;rDGQ-A;YZ-l~2Ec_SrMC_@u1wfTmDA{1w;l8% z$H>627ee$&-4TlcaX=vpC-i#}&D+bM7+@h!L*D??T4*{9&5lySRP4JH7E#;8uSWsd z79BPdEsByuKn-mr5vG*(AqOZ=ku1pHB;}&8Qj8=5PJVjUd`$Nf^g_YcLzi`;+x-&1 z5(Z!e!M6ACj*ecbz@9j)KQ8o#*sDC>Kubfc%&1ArnEFVdIhZ}G6RfqO72?eE*39Yc z=gvTNQnDs<1yF*N5-23iDEB9knxd@pYy8At%`Cv{IyHeK=j+KoKTvhPK3FAWR7WTU z$3+W0K8&}ClaNV9&iR~~9H1M(>fOO6F{jFJ@O#dgZ|4R@DmCNW(W#?UfSj3Hyd>Xr zj^QYat&H1ze&x(clo-OdzTozt=i)9gcg^f^5%1S0O(&kibL<~B;8uFJvgnP$zKiR{ zoYk{kV$MCvI((Rbo9VXc(&?=`?t>**n$}sYJvypZ&n*T0%DSX?|gn~FyTT$jBO6moqx}b{K5vJDAKWfxRn3~ktme6c^e1Mm_wTsH09MXs z&7KkR@KE%6f9FkW-e)Hso1o{Af@RmNemOWi_u|_CbGi{nS-7KmM zy1q2{#lfeCFAfVuE$}I;@plJ6`zpadKnc`^tgAxCRg0wzWt4LqL;FLaoGL*lpd0du ziPlllkPtygwMv9kMnl<4%s@*Mpr8_)h4hvVUbO|XDw%*7Lk#0l(j+A>U%jo6s?ZoO zRfUq`dAE2JkyE0iMj2g{QAkdSk{TnZqhfQ5%iDy!Wkg1S0>@A)lbnalGeN^@}P(-;kM8-jw8)6KVGUw!+RH9!#B374T|gzq5R4l;I`J8s_b^2|LdLyh9HS_f2$m(qI8OdV zq37=>gCrUF2g#6eDE?t`+Q~RU2C0YRIWinDJS42M1l|*+eWXwwi3QC{HC~{ifT+R> z(V?D$012!obqTD8iP_?|lUTb|$lK~qxhZ2oz;#qu&^)l680wQ-5j-XU$$-%e+y^*M zUiy4`HvQ%Ry%BTp$hVk7?4-O#Ke0pp3l~yiPcd95f&$?to&pIg9Ny5;7(9n#6=RRor}0>7d4|@9EcaBNKvG80ck|eg~svHP>xzdQ>FvCbc}Ii z8&h*q+Hd>|-A-CHb(+bm^NP#flqg=+GHsQ-N!wRVoit8fUi1FnxmfN62o5Qa?P-4x zBtD#b&Ue21opaBA&i($_v?BU*9tJ%xuY8j0PQwOU-NBVPv2IW-ACk(4=6fH?LlO@S zN%RlzV=p&yJVNJsXe9CDe}5o+y-Ngb2Y8m#NkQXVik{!jR^$(9g#PtxKsibp^*#LBVjS$7Y=F z1{+UWFRD|y6%*cE)(s}ux?i0STe@$y|GuSRaZ{AK2Hasf>*f3MCRO>^I42F8=44pt z=VOIrTkw<@N`xvu)mIyzE>&wv^+0`pE=Sfq>@JE@6qX-~(i6yy@MjqPE_0|=p^^)v zq@bk^oSdU$H2#l}qSjJ%sDz1<`Lxj8>u15KmhAMAsqdcd&kiYuP9k$_; zYikQWo){sQGpTvl3Hv;IK9_|}qG-Kt*hY@}8DX1FbL9*6!k;4cHG7eekBy!pQuMqY{dI%>adsT&>-le=+N+U z&|x3bp;5?*qC?|GI^55?ufLyz-2S97$NuoAyM;Ih!sBexx;NQ|Lr-n8KHv`TX@n+j zC$$J$=5E~x3Pi9b;jy;R`!AOGT6zo1`R4j895U%QiG_*|x?l z+;%l?;XX;@R@=Ei$v~gi5OCWr1zhOelp}4$Xi!}iqsCB<6IH#Px&;wh{rU0P(MQhh z&yGIqSyDcQT5lJDwnTh@X%l888|`x|r|7Y_6;IJeZap5xer>9fQxVeI-mVI%RD|rZ zjf9c|uAPA^b!oHjwv7ftSP5+nTp>z|@>ggVb|EFf*~JFq{Jr*l;P|9o=8*c1VJlN{%6kUfkQ#8tC{IF!OS0VzvZ>n$ z#p;}}Th$`&3B(Zg2pvM_c-TABQFz;Jw*v9dMs)_MJ2CrgVM*?o; zuX+ZEnake{g{pE)u`Z?ZSFAj#Q=sBHLXTYKqiESw{(TXgik42r*&o5F=x4cuG@IoH zXdgorYcL{~qOXk_Ylvd0%0z!4pM<2;u>Ba!9hI$t!d3v)KJF+xg{9;;79o8gy7Xls=XJg(&_PqxfHn;(s}c|CNwGcNKRIjt${H9!TD|a1P=- z!Ltx-1aRDbVtQaptqoyEY}ne^kE=oyEeT{30!1j6qSS}QQj}=>iSZcwNqkS?dm7&} z_zvTH7T;I!eHGtx_@2l20=^gVy@c;&e6Qep)%LYOiqV%9h1^1_x+PLyVb37t%5U38 zgb|Ff*f@zN?SmtxC}(I!)KXxj^r&!L7zsTO>I?0!!S_UjuX)2Ag4YDIdfXS>At=ME zcj)b*RWMOpM6XU`KzTf1zb2G!5F=vrpolRlT#G9GHR8L0)ZplzQMksLZ(zm_ccCz< zildaTjFTvg2IBCs4}Vi=^RHCWOfa|>$tRQl4YqpD$Hm4tlBa*7e1_pGRItd@yTMLG zsjyQ}>>;YwAKKJU3|gch$FjyWtS2tpc1+A8Ddxi?NJp@=@o#>D*3*$%Z*>O;`zcMz z)Nxr%^U5&aLgrf`v?zJEkrS2Mjv2^_nd#}^4pT5kkjIhxSQjI6;WMAQOIvC$UqfjK zwS3n??SZltuyb4DvZ%O5qsSUladA<&##LOFC|uX4M^t=cQTRTy`m5$Jr>LExkHD$x zzV%cTK1&$i{wRE7VSIH__{PKd3{m*5tND11Bktib74PPA^a-pAOxk{3PA$<%LRn!S zJkMZX26w<)HRe&zZ0$FM8!p8@K`=^;Fu4?#BFv~2oRguo98RzZW5T#_eLO*(Hv5!H zPcp#FKEi}BB}@u8qUwkq=sT0spU|Ff z)RL)NX~k)?4SMBG`f=g5yn`IpzAN7|i5wRwdttUl5bEWHTmBvCOUUC;&r9fE!4eW> zv~yi(+vvnF`s{S_500FE8d)bd-tknGFEXLPJ_1HoG$v9%YxNb@HBW6zrg=*;Pj4(! zm{zXAM`&^;YAnSHd_=6_sIe65?-8-iMu~N6Qw?8<8f!BRUyT}TGY!v0iFJEZ4bMl7 zwV8$&qQ=@x!;4X3h4myTW0#_EhV>=jyc~rytTzGYl_;DMt*>x33TIf4f_8depzeg9 zg|>gfnIA{|m7N=m*P@;J+29P4y;{|ZMxOW?z=@yQU@KB4%qYjMf}&0Pq%%SGnuj<$ z^z!kkwSKeM`GesTI&J;~2gw7Pv*`m-W&CUz>+wy5B;+|HlApLZI}ofD_8~VHunD}1?$N9KsmGpTXT$8E(CmC7XvQkuV7V~=eQwQLt~9C zopzpLy#@NY&|b*W&9JoCY`+=0r}7T!=f_4AckFl3Z(a`YJ7W2HRjyqLgs9hFT~`DC z%3sx(XFBvWmuA#=w!f~@kXMH?z3$m(ReA(>PE|9&d+L~~uvvB3(yg#unp4RSrc;%M z+o3e5{CS1{>%!cUHEgHt4V7#~$(wvFkYAi~{KG5!72~WhjRkH}5JKfn5Vl@o1>2-K zaA!f743=HC9%CB`h1n*rZNe7{d&8yJqYlKt48?Z>N3;Kll0}gleF=#02Yy zGVktilhn4b{i42MzW!jL~XVZMW&x0Bk2f!g?twmdJ-X*(5j53 zP9ys0m`O{3m;W_YxxEnZSFIl_QlXZmDr3rWWSbldNLsRm)f7AI$k{{Tl(&L?B*?Fl z8OL(S9;LKu{Gp-zP{`d<+j_^pZdWAOUnVCV@ef-aDGz_BdVA)s&=y`FkmH&_engq4 zyuR@~#r7`(@gl}P_f_+Y>j6IDWY{yh&)uVSwJ6-DBDiNa;XWO~J+}$>nF#JTHsKzQ z;1)LFJ{!Sp--P>>2=4ovaK9SC{pKdz=OVa2w+Z+82<~6ng!@7S_vbg^z8Jy%g-y6G zMR0#{6Yk3q++W&+`$`1&-`IrvY6SN$6E|v2O`ueB^HkTl2Si`;9B@9xwZAFccTEJO zDu0DHH#+~SKA#rqH-u*Se4qLpnkt-~w?nvMuBHl?&%ZUp?!%W)@K3Na!T$-7|3rjd zkoF}5#`&ye>eLLtsNsZ#1#%`4VDqTi2AlVz*1K~>!!}rz#5D}y%fNAD>5^>NS6gM< zvIG#Qu4x=7?(fBn~gJ=YW_;F@I?C~*)7(18g8k8ZdHpl2&Eqpn*qmvnXn z$WCYA#T>olg*>=5ZQCg)pIppd##{}|KaAc1=C*~I?j&?i&IrD1Qu@?VtEB|NUrHfO8=b0ydTlpj^>j^!*Eh8=RPb>Yj%2%x;t)mDJ ze&XQL(lx$M$vnwnE2}Lmu|;GCN0Q7wBg|N4G*pW4DL>1L8|ByFfn! zI=-vWv2eeK{kvM6ZL$lL=?H6kH~r11FXh(FKpW~yCotHlDf2znq;97WfBGfjyGRM| zqH~n~xW6|UIV%6m*htc|JZz6k{#UibZw^(?9 z|52N0zIkhE!ghzVysO~{Ab?*-3}51y*#?c#&QFZpu|n;4;~C%xL)-J|;46ENjG0FP zLg-5mh#VQaGd=3d3|@``Wz|A}eL>kg83y{}Hw?eg|A%v5JTzBoo3Q~=hxHanm^x*9 zjhls2vV93w8+BWW_+gzEnPo(y4@;HI(Hd~H!7*b!F9$IP1u=&lbRiG~ha3b4f*=@= z;ZDE-&|%pJpu@6bBH);t(}zHJC_(6f$7X!-Bf}>~Op~%@s4C=c2<6b1cniowW-I?W z%Jk$D7>%&egf9UYEi*j-Z|H&idkTvKLaZ;&e9h{MziBZ~%Z8(9RR(q}r#TDkbJ7iY z%nxJuA5g3=7RHR$^9&&7U%}Uxa*9U*<7WqYU-2b@979{pYlwf0$iF0Vf=D+BIz+YL z3nVCxI$diL2(a-y|gUA&NMDX{RH!x-`)0_-LLo;Jv66pv6u9+AII>Q;5vQjH_bOJef-Q##==jL<|65vK0Gls0j!;2%d{1BSBDW329_g% zI*So>K=WP(xf!avz7+H5=$)Hvq~}Wj==a$f{`W~P)oC6LCp0r#Q4VqY1or$*pOcNq zrZ8}1rX~NMDaL;xO_F4|-+MgunF$f7&sj^oB!iM06=rmg&_mE1e|a}@nN zi#~eYJUwlhLboCpYrsY%&=_z$i)r%vzdH zGcn=QOiua|#?4b00Vm>z#=mi6)_X)|*+@yYOB*+Ki~MfQ`r=@e7d6^C2FvvPKGJfW zrRX|$hZLh9XhS~}XkaIh<>UTe+?`!d_~K@+P2%ZdK#4(>FU8-A0Q?pBP1K^mI54|l zIAj=jp$__zrl-uaqoZ7CRL{SIr19S%^6NyN5cwM*z67c}JAE4K6fEJACA$#|wVjZR z(qo<;hf>Vjr&HLuqdKrth_id4@rT?ncm!XJH95i*HS|U#!m`ffE zFsFYoz}%t10J}`$@5&JqV7`~01%B`>_pU|4aZGgJI`iSIok zzd_{h6ImugrvCWv5$Rx$=gv4kF^6W6o(Su#lVO~hRU0c$nCw2y_&t1QeetNkcNyc& zTQmTd`?;BERFY9Vnz+%~oxY4=RB#W2g*-bpli^ z*OZjgHlG3$#Zq3aJA-T%lf7caMRu-~UE$g{VW@>k6EM zqON*A0hliOJc`mOWjD?zJx?uhj*F=^QfiI6W;J!En7R`dp3}2ryeUuNo|QdfVUL*B zE2Z^%(t63Ba=}8gS7(rPMpvPvYw$Dt!Q7c5;} zYdo|HqmGSTQe)TCuGPjPp2i~|>OZ{rqq4QZ%d3Nz#lfr6;8l-#baimlga7eK&y73M z_}#Vfxz+JGaa@qbg%>d~hc!JKblcZsdo;MBABF0@qQ_rpBd$bs=+uE@t*wU{Xy>Lf$m{1G%KwL7Gm-OKl+ zn!_@B9N-Hpz3{Bvvyv>;cFUd&z))_HE=Y#vxib5fEjmEedJb<&I_t49!CYz?7#Je;jiwHu->}Y;}uiJyKea zC#~mselb+eN3k@)zPTZP;h4+hDhig+pCWWUo2qR)fhqP&MeSZ0SCO13ulc(P0&h0~ zk?&^pudE(*n-@!4rP9_1C+GVeudeH|=DScc^NPRL<4Sydbn&Q|TPNk#y@(l%O-GW5 zIyBD;%f1=s>U$^kNvc@bAQd*eh#82@LNHNBH*%n)($(tXUAq@gxGw;u3pZ~n!F@|g z97onms$6H>0P8`6aZwRgTH@%#osG4w+2wj~?apP>$}3O%5Yxz-RnSKRjvl8TO+L8yw08Vq|n*@thQ~rXGOnqO|0#bYP%i} zID1@rm;OZp<4%hyt8m77OUs>cuz6clMrdJUT5!mS*5G@+XLSicxULz(Bn6q@yLhbGUuJe*~bUb{&uYv zHF}B~pVu@#nf4m1T$U#j%eS9a1FEcckf6#y@nAubI?dyLXFr;s+8u!l$SNscJdP(l zx>P3?w@Jlqjw4=E^Ri~SdpU8b)NusG&mbDmnzk%gEMHr$S!!l}IijJ~+uXi1>z-uS z1)^bx*HA!!l`aZ{X?>5O_Ia`Ke^TM=`6XWee}S|7ah#*;QP;<(HB^RQys8Ndcz_wY zd_c@O{0tpZZ(ukh>RKgTE6QeWo44%%{T`g}TTja(P)uiyyLbM05FZrKb8O)BnUSN* z@yh?9UT>e6wqHuy?@8N_$^}GBKXm6VQM+5x?)GSRqYP#0AB;U5n~!@2K++Cp*Q0$4 z`{v{3dmd;%*#o_(jCt`ar{K}63$Fr*udsNr&UM7y;@T?~Hc5p|^C^zD2btURex6(S zXl!9@KJ{5q$zq4=6?d0wNGxiWikjzB9X+Bpf4lxcfzOM|7Y{s3%XGxkurg25nZ&eu zDXsp!p5=IH=i#S$(#}D#VMuDAF1wgDPoVz0{^BB|+9w+vJpkxk(^jl%(TuswqP9-b z5+a3HtAB~0^=v)s(K527^b+TUlwJpv6I=kmcTUQvpO1SEL*rG80964#8!V7N$H<0T zm~nAZO~*>Pr)Hm&yI<5DkaP#;6KK@)GGXg^MQ_{RvU}4q9C41b3nR$${2J=XY96G` zC(PfWWc$mk->_7-eCD@eKhS>nlGNV+Vdd|ee{}u#8l=H1p0=xE z{fJagAhQ_^WEK>$K<72<*@m_3#?|b`cVVB?hWgMe>RTjz%e?kkdXB^5d__#J zl+r8l0QEVJR>$>6`xo}lYoDj*cyenzbsb`Qr|6Sd{ft7@L~0p)o2 zJ13r;@F(mtngin7;oU)Y^Q&8BC5Rn&sjl5lgl)=wheq6$}m6idt8TTE!v^6Wblyb+*O$UiJ6mmixb(_I{dJ ze?Y1~fW1bs_@GpL(2?wAn|tr=`o4bo+;@xLFBW$kly)4%7N1yhNGdtxNb$0rwf8LF z?^#Lu?y>ieiA`M?-e7-CG zS1M>iNMcNuem()s9q_br@~~nQBISL|CrJFuJufvqwC9e=kLms(e$^{J*9g@U5BLW@aQ_wPURLmUvaf8o)>d!=KH+3 z8pDVl7*cif{ScmC=A3ocd&}zFrsY>x&f!_+vn3JG2Lk{-*CR>#T*u_2CYM1fZ&|LD z%J;6ch~@jm?ENdlV)o(rp7rdUO~a907(=q&=yfT#c{xtZZClN3^JKOyTUWU6y|Mhp z)6oylO2@Bwu3DtyV_4%6k58`lO?vP@ecN+qc6EBzGd=6o<*(C%tqqVI~Z}quICgxPcK|_XuR!v94XG$g^aIR zftn>O3hu7Q_uaQWWjjS(%cg-IT??fH3lA@13epceFDP~1ew@FU{~|`4b$}5^7h!m1 zEfzYupOx)!_bls|uZd-yQdy@wfAf}CVH&4d=*J*(NOiH^x!-ZrrCB)n%^r8sJI9_J zTh4mt#B!IXYL{r(?cFTUW7z#7ixkL@$t!fkd28HR9PxM#mIE*PM885H1g z+89&gTHXHDy8SCaPU)BGAiNr|arXl_dfm_bfO_4kdX$h}_#y_x5$g@sBP{xzS_hEE% zTr6mn3IJx;CZ)A`(%QUf*)*a3S_RrL;y*8peLx^AE|n=F%^Y2Oo7v&er@tJvsWXocQtyCsJyZ(u|%oBSzkz ziVr#(eIA<&_dL%j3^dR;V$;!$5mhwN>jfoWFaBz=Yt-GgtPu;^q=Gg_qGD|-A!`Wh zT|}{(lUMADeQVA!hhe$Cf^sRq8w@mK}4(JL9q9Nwwx5}PDypAUc|)|?Sr~SZ1)J%k%XicSGfkCoRmk+kF=PI4xghURMEP;SE@J&ku}FD3u;&vfI3bZ3kK zw^IWT(zE`(;pOai&n=z%t+g;lIE(0CO0D*+z}L+4(d3L8_0RLRJh!=RI(EPLCVvX57Ncr zKB*Y9BFe5K3rA^<@8nx2Sxy>ElELILnAVN87(7=m?GTN-9Q_M}-mvR+L;0GaVbuVP zv1NhI8$9iDt!;`yvrB=-C7PVLu>hWlMo@eMEj6WR5Dqdc}!knvG%&WtFv5wdOAdPTV zQ(U{n{CX)L3p?`z57M6lpOtlN%~C-#yapd+5M=8od;Ov!_KO7#QbB`K)lbI@r3rh? znPnd1E-@1eG1+Jsb^2$R5!h(+J5|bS7xlX&{Vwv~?cGgl!n%X=@$1^uuOxpt*-_=H znokzBwUV~hqpd|(o?i)w;k;TY51ZaH>M!=AH}ycrKA&Mt!(oEo?kMTWjrpg!?Y${+ zqw(+veCQxi9KJv9hHo^CSs~1Vs6F0Br{DY^AzuBSggrqC3NIvZv7BZpF3bjIq7oZT zNDPyZ1T28y9YTn)3j9`OxHh0kDBPY5d^b9f?f(c`n<3!FDNWcpkU#7#ujJ@qT>Mf( znC2-$%EmZQTzD!{RT{)bzl)%W9>yg_XaK~eEozXid`nzNpQK}j>=Vj?u$<2bi;)=_ zLr9(6z&8^A%bsZ~55z{E2koWd{TVtA$UdBwUQCcbTrzu?2I#n4+ERMh+j%-$ZGaE5 zg3l6YG#THB4NF5-RNl-8-fR_bPT(FwmXNJ}TaLYhska+!jR79$K?e!hN6-rfLC5Kq zvcmMvfmPNV0=8F6v*&IzZ?ZSZn>-$_Na5?rup+LwL3hDh% zss~&yQ9^667YjvP9=j?LO1OODu$KaHB_?E!4HnU2MhXLWLRjgT2JZ|NnB#DCDn)rt zsSj`Jj8ELPOdHR1H5hx3^spCK&JG?w(qQB*wmbZ^5eE^CIDBQ{Z_Ut~Eu4{#&f%D) z(R#}=YUA(RM5ytY1wgR8t+A!Ev!f{`c#FMsp#3J^2-*wz#CHd8@EB*t@G{RhPGdC~ z2Yb#88^^$736F#GhcR@@kJmkgjJ|a{< zKB9K%BWfJ@&m-g`0@3-1>gz}J-rh$?K}`CvIQ$iKfBaP<6o(&y>!V6yuLI#+`fzQidMQ!|g33F2?yy!c%C2C1^5nDZm_&{(=y zDyS}f^j9S2b98p_K3w!p!RsIXi^^_O5l>x~FMSMWj3=x!{_JGuxRWUgb%Aomve(W2 z9i?QN{3u5hr?h{KF1-l##i62s%uMrdlHccue2;ECKrVK;`zPemf%vp?LmsO=U8Qa_}8DUp&s1J1;$Yvjd$ zhwiaNXC7R7cxhhqEKTRgcJ{gx zJZVj0T9cI4G#~riP(*8<3qv?zR9Lx)eeWGoVKXM8rBX|e)Y9)rnooRgD0U@&^Y)W_ z?_}Z>%X?W%`C@s?N-icCIbAft0L26YbUt}Kt8gu=b~UTkQ+HU*>XNd$<`bbqPKitN z&6Xz}@1!p4zkmAu^WQCg`l`o32b&~zwrPV{&oi?fyBGFK`YJKATFS(hADwpSTNuK$ z#htj6BGq+@hE7qpSJLgBPryB}K=_7#BkZD-R_8GdiD^TgQ?GhPMm;xgL0KJ7%7Q{! z(E0f1Ie81`J*B=4V%h;I?ZAAjS8s5%ILa9#^#X+Kn@ z>@hLznCHxSkJ;jx;i0U-M#_TXj-d0}^~^kH?87&t%<}noEVnze7E*Cu6#E&5JY;*O z!Krn1ia1i3*|;1pX134M;z)MxeB!!3&sq6sF9aBhopYkW1n!bb+8HcdFR77#GjklX zSeNr?jgfDkQDcXODdxxX?VDjTqpwD7ls~GmyV8E00zm;|v;cbyr$7l)k7B?Tl*&fZ zG1Pj>AA|+nTs*Te5i=5;d2k7MIWm#qL1WAzf4HPOF}wk|_R4!}*&Ejym0vRRQ}d;4 zz?aJWRMJ9YA(h0m2+rgvoaqsqu^~=-oJy+F7H_;NBvTP50d^-6X!LAPRPpFSJSu<1 zh?)94dy*?lu0F$M=eO}p7gJ%42Lh;+Mp zo~3+J3Va0MK$irhV73vQbqOg#s*o0e>#;F!tj(y95uAlleGSdSFazBJR|Kv{e%Nd$ zQ^<5FX5FFtas^9;VI@WfUTVD}J{$06DJ_`@$#VM9QfSono$c8|Hdl0mNd&G(5wJsP z^ymKwIw6b7zUo;ja$DWD1?PLpTn*0mxMHq^D^)*PdycIoa4*~iIF)LN_FT*bbA_C+ ze8^M9!aOo8Rz70oN5m??tWr4}%Tvw9C@g3%7_+>Xk>Y@6DB+3}1$@J;xWHD@fl!z# z#}w_4%3smKs8gWgDvM}8qtzy31gE-w1=ZsiDx~TziXMy0p}vQ!V0My}^{g^NzM}t& z8mo$8scIz+64^&c+$OGvzgJZBwAB&Pjeg-Vk`iiI zn|jwV*);F)f4+s^?*uMlfK96?&yy8_04`Tu&xAtG6rts2+y&d>uKxrriaVD+wg72cNf09@!f-O2Szl&vYV_8q?(W84EDX#8m`g4*LFJKi_*JS zUFv|yp}TU`zE5dujcW?W9b*>w>H|{vepNfI7$TH%=M(mE2J%zim3+9JMD6=+ z=keqJOa{kf8(WPO!Z2H{u@7IH?u)8T%gnKx$k&11{w+*DCBMJKpLX!Ki99 z#5LI(aEGTDn0GNDELnaaYTDa=&S}>;VjtE_M=u!+75V2_GDx54E#~RLEj+bdQJ&g( z;7Orwg`or6lwk+e&-(^j0pU%mJ7b+D^26;`$3#hyr!kB*JoNLiXbleGdstcPIfF`< zFg?|6BU*&mD_^#EB_lmTSE#L<+#S%3+Y|DQ%F_|zp^_LIGrgxW#zqUz{!bO`jj?tm zR9`kuf9Iy@4-5T#?CUZ2;x5Fz9vf9Nh4$QD*_x*H$dg(+s6Ov4`G@D?=eQR}3>0!C8o!-u#_XCQ^JL0!ci1a6w5Z&lT?Luv`? zs|b9)UhdGRChG`y_*0YB$8~*bvifaz0=b8=be(sA62c+jFfd7v03)t%vz%4M z?+%posZC&?fUSIwg&;N&@vnexHnIUppcsy!x^{}ODDtD zoKH>bHkV;P1=O0p!2OiJ5oL<)W&tZI_ThJwbu)Y|iDBiP7RrRvp)zbg1GK_o;Y`R! zwK@eT(}sS3!nMa zUBbHkVCYGOT;X=r+4FkE{_57sy#1W;ifV0zwJ#y*vP)f0gmY{~{gl6L?S%y8D$DjI zR6>Gsy=D6n2;J?A=!HWoBgsN1XYjW?gnpSc{cko|5j-!Pl2b)(g*r#rzV>-}osNYm z(&`UmYi2A=kqUp9yncoFp>}}`o<@GudCZ;m=k7t_{C0CU7}ton3XQgPUkxw>8YvIi z=(=;0yt}YT*u|HqiRxBHae7gRlTQ&S_B0ghtn!WxY3twXxhR}rZIN1@0^NH9d5V?Z z2zeJ|c@HlMDUQvC88FIz8WFFg+%*u^rs>^qTCRI^mDN24s ztm>$-6eTbsR!x*x7nJoSB35nGSet2B7d6&q8k(ZSy11!^^-*JOrs0mLu{P7NJxZ*N z#+0Ge9YyaPk-xiMO5$*+g^Q55`=ul*@YN$E?s+MR3e5Egi5)K`QGvrAA+h15B%a$O z9~)mvq5?}jLdT|;l6YYg9h+ZDq5{7?LdTsiB~gL=9wD*ir6eja<0B-tzLZ1-+I)n> zwwIE4b(1I4$&JvSL|_)lmte#cT9@5;ALBKlbHja%;A}7a{M~3qV!w=iuq&#&nAN$5`B|m?s&GYkRXGb06nP_D=3=;*FKJ{dTpd&N{4{c7 zqn;}dt03dbR_2a0GJ#IYl@i6wRjplcMqD0cMtns%>&2WmDSXZ=T$0BkKrOryaN8b- z@=Z9eN;kMerGmS*G^(EA*+)1|7!g3_F(-|hR6FFy1CqgWV#%W3NwB{rytefnS#C^t zZRA8io~$5eW}5)iCfKM2DOsWpry)6b9DRs_F7 zeiRH+qcWDdZeY$H_b|_q`tX)R_wY5<$%R`&r@Fj^?J#cKliiNell^Assj6!}@ZXg) zy2}ARprOZ=y+a1hQu}wmaOUv%0` zS;F5Z_BV-qj>s<&DJ1zOBK1V*RZ{-T6!34!^(*AUHyjL(e@Fp;lgQs9@@K^I6>>R< zJRw&pG^?RO31;NM+P6I=v`50gy2vKHa z@Ob~o&ZeEMyGB~ucJAJ}TMhxH#K`C!8w=&%z^HW4;G6GdF>2IA&Qm5LY{?u7A{ z*>)XZR(pZPxNWuDTUuN9sMV%~eL5O5NyRFS;ib~ISH^&p6M~n5%|>pGCDfFKk8l5>FoMWghfP9!Y184EJ;$JK~EQ zzBtU=aih6nziu0%_3?iQ7U1dt=gAbyzlVNE#^U_nNYD3)&^cWG`{epfA}d7x0g-=5 z1fdfWh=ovQ&X5fMR?~tf}M&x&id`RT?i2UzF{yC8!5&3-( zpB9x?#>5M_#)E!<4D^Ri1YL7-H;Rl*x^AAD8G$2(n5daU-nr5BOfAGMrx0HR~K* z`Iv)@a!3?GjX}r%5lQ-EBJ?gj{}GWN&~piyBUq%%uMugX(EmoR|4ih2ME*9BKPQq+ zjH(A5Ge)mC0gzQ zYReh|8dLxSULmf3MKTz=;61wh1G@b8ME*HR%O$`6LYMzHktDhtqX06}z>o#`HFD8@ zBG1e^WRPo}{Juc05RM?B)$xBy~~=x!?zM|8{@rk_qg8HZo~D1{s&l24?7 zNDzUKKTek?h@2#HipXgqXNcUO;zKqA_+h#{OXL+I=ZKsKu}&d3Hn6ktGxK5wJbbwu z50vIAw=kL-CJP^@@Vf2iVqo(vznHAMrAc|a=HnlxJ0gkahH&}|`Xq)w^zxM(=8!bVaU)j*1KP4g*8RPD}?g5FvMKoZ$LE&RD* zMd%14?1Ge5>o!QH&Xsz}bmXb+hi}r8InN%|`B)h}0%}EKL3oDGdoGT6%wwK2S+B0J5uv9t>Qz9A63M42>1L(X4 zIik<`UjYdMAa1vVJc*iq$|g{`#QgtQ5HE!_@u>jw>c7CcmqgD8{k$TX#05&0!VJ855Mv0>E375KTKT?B? z#gLm1=otpmz+6ORq57P8Eudo}2=1)M{qTbhZSRYixa1r!V1s%H5y3Fu^IV_r)IZw0 zfRO2gFym~3GetsaBri}EGbs389 zr%uzq4x}mN(@GO0K4*XgLA@18>{uBP3y+9teNtMVC#{bxg>^izKePk?Z{VHuC+Vx1 zjh@U#*&ICsNWF;6&ALntJ19Mc%QDI+`R6fzP}v{ZlUFd`?*$ono(E`(fKb^4@{5?f z^c^1)abUERm^B=W=jj|?I1Z#oH;`ppmpfK^-ajC<^?#TwwVe?SXFTlxC(rVX0K5Rd zl478I5Go(2~SHl2+0q2nBocseb-Z3#Fn(Ok_1YJ#OW+?g)esD!Uij*Cv zl%1r^2c;C8lv76jMa&;m_D2(xTk0HLxFY44+;xEU^^XM4ONJ(J=Qev=c1bN=Pfb$G zF}gk;>(G#uY&_(z-Ep0HJi9m>w6ngRQxcG5l5*Ppl4$TqF?V`fcN4tkUa9p2U7uu< z^54#UEAyKPu50fkKS>sK^^&gMqpODlm?l!D@#C0`^oq{{G!fdC@@mPrXC+?h82GSI z>Nx9p^@?X?+%v&@JG!6te|Y5&E{}Nk6SmcrSO#uIKkFU)i9Jv1IPDoeCk>B!xEs>& zlxLckrU9w}px%CZrZcev@>4t^KjDjXbxAWWUtOoGNzF7}O={-o>J3dlq4`W|1_;e( zDDDJZ4aHrgt4nd?^3`>^nvA5acY@IU1j3H}H4QMx0+ zao2GdcBu=Dgyd81NW{1}ujp&LzNvS;@=oECLb0epDryjO8>QUF7co8Y>1fD_Iy4li z)On2{f*xl$64ndKTpE}N?|obbAsI2nm99?rq_?ucJpp*LL}XJj!{Pvyt>}}5zC4HZ zQ3qfJN8dW=IQY$Bx1Mk`m*c*1)peEHa5Udlu5KWG_AVZjikg5~r7wDW|6BXX#Gk8c zRbTDVSHGuspMSRmM(y>hru`n%enqV~h@0jxyZfM2`VY&sg8Yk^Kd9`FT7hcb4NRNf zg@cm5W=&tes;?LI4U)dWeS1~E)1%+HJg{=+Y3vWqzJCIEa8mox4{xsm5(@ukFMF;4 z`NS+xW(M7#1_6{2gw~(uDis%IyycY@zlP*)s%1I{@l9={SjW)i%5NsG!>gyzZ9Y99)jUIg?Is-hA!mM`s=c|lWr=tWFPR!^(~&>|(L2Z<)CfLncD37n+m z;gt%p@_fR*OlfVB(4pXEwrl{1><1H*v0m3FT9f>M*T>qJ>V~A;NS!+1F+Hm-3 zzKl7EX4jGIJh)!GW39M(wYV9;1`@ChlGpXrJ04uqSEELV`daF#k^uo2?5PHyj(%`r zZU2eY{U^lzr=TDnkb}7CS^&^z(zB^|{mXyv z^1D};uC8#ajoqHc?x)rd+kSZOcZJ^;{>l9x-uIlpyxM=+(|;Mdt{0ZO${(8m>Y)Y5 zD*=)M(O5<%b!Zf*QAt^)FJeFxMSr~nm2UAIK>@q9V#&^B19g$LQp=DxpCEo(3sLb( z3sG!vQCjNgcJPn-UMiBK5A`IAfl_?!+Tx@;$9?zRVzG3$R0^>AJ(6yZN4Ll8?~eJ! zId{!+?z@d*SqE@mL|vz(>-6Y4{oOTy&Gd`+gP<1vSX=T9NZ?(Si?55iT1i*y(bWbd zY{mD}kidotF0DJ;-Sw{FNrq$uT#aNrAnFcEx`Q6wL6~Tz!Nf0|1C!^s)80x0_|!Fb z);rwdRjH&E==CK^O>-H{$aNR*E`P_iI4za7Ee}YgA&zxj;hL^|Rafq+18j-;-S#I< zQpKJXtyIw^>bfOew@26Q)fKJjDpmnB*5U4U-+s6E$w8@-5H>6Q$XOm;@24l{lkXKM z5BJ;}cgwr+i}wN5viyovwolaUmvs9*y8Y6%#o?#@ zVg<|-V_1LOV>s>y;Z@eed*V?+=Etc$&O7=MV){BZAy-HN@ z9I`B|0d*#LUk3B2*yxJ?w)Pv^cT%3DJU-;;rKx6-bJt?$qx%c@-Lb2ACQqK}d4=go z6YwZq)lW>zEh`rE!Iej-4+cda3`+gg<0L1+Tvpb~7tf1oTvm7Mlh@I~b-$mqGUh2d zA{HHyijFvXkzaK?9KBA1EAFd>qM?cgp|wxa+-FzfJt+9b11SEEK_D5e8ERJzweDDV zi;N!cF*JH>b~pwXPO7nY%H2JHAd=yATA)%{O8{Iac_~wFHe_%ewqOlJqnu8(j`cew z(*e$1$oOx9`bImh)Q#vz^4(83+ zHT2n>Fp1Wi^d_zID}_jv%{INmqm%1RoJZ^qSh8q90yq>fC+9l?jFzA$#ljLOgi3K~Pf*mp( z5v{5Yx#q-TfUdQ|>lm%5Vqx^6J5?-cmrB~mN^pzOAwoPfg7lEl}RVgB3tVtwtz6YMHhHUT6$JFP!n(tkF6#e$#xw zSk2*{z!kCJ>_+h9^Sg0`vNFmJ;UH{c&LW>9W`~pMFtI6JadxO19=>?WK7N3liM)B5 zo4LvVW7IHy2O2)Uoz7$=5`u|&N(OA%Mc%tXOmXZ?jxWx9*WycHN?EC9;_QTc${}rB z&MFJtV6H%Qwm}*6kI*7m{}vh*oM_01Nlt&L{Yv_m(>=MZqP9)awt2K|-T+?U3r$SA ze~o&OX8h@(ak0sSO->X`QlJ6AgT*eW&XmH{kK3=Zl>Tgv4Bik39}ZefSN0CWCI5Qc55= z2KMyBdy`-cT0ssl(rMGkPk66-8FruS3HC&Ll3mMbxHzC(Kwo-u!Jf>n6NA?&LJ}-m zP)Mqv^}EtIAYz? zpdaN0FdnjlZ}g=ByMS^JHIv23;ZoQ&`0Tktu1)bgLwJcom<@FHmNL-d&~B%m;044T zmsMQ<$KzjpOrPI7B%^7O^6YG_{7ywBiiP}&1p#H5lKe*niRp*%V9k>+0f)mTIA7AO znF*^ExKjTbLQMJcF~JW4T%m^1r{Qs)jmH@{z75}YpT>-9Gp^0JHn;oY@7c4McJay~U+Y+v9$;@3~dpnmV`PopEwEAy9|CzqaM>?-GK!ZVy0A zvvYryUwGYVdt_g*i`iB4NeqiL^4$WeLsWKuT!qr2#vhFsO(XZ^4{TutYgfRJew8Cp zxI?XjVkE2>X(=!h)HN1ygZJW?oi!omUTkQTs}2uIoWs%MQB7#nzlol~(P5aL!I6JR zDz`N~<7ni~^i;|Yr9)}qK-N>i7d1bFBl$2rgY$`P>lq)^lb$>KsPa=W2MNzlE(R@E z!U=ga^bAquSG?FLDuwT69@NnMJ0j1K^w5xr|0B9&L#00^7mb&AdPUtzIpiN7C7W;E znwq5+VJ=tI76kqOD>Vaj=p1Nngl^7`K+&PiVxzpXg2BqeV>q(d;IcijONF%TTu8{i zC5JFVDLF((I@!S}tg6hXJxp`PFQ&PMu?8S&O?3TPljGH9%qK?HTB=!ctwrvsXWYut zY6nJF)Y)__FoE$9Lj5Go!cohDH5;WOSVELDQ-7&Y@FAA6u7_eN=B@r%3O>XN_TZ{} zC5FnhKbDd&SZiFtx{BNge$e|PPn0uqPK#2TY~Bs;;H;$!1r0MGy6naUowef^BJ}RWYX~ejB5Dr z5>Yn;bH#Gba!dnx28x_yFD+S{PMEavl<7l?$r^}%N-owmux8;1o+D3o3viFqn|fZW>VRh1YrlVq zb1HvD&d?NtOQLB5r(H@8%X3W4ehzm@6ykplcS&MhPjFC!scJeK%7uKWYYLV?A<-Wm z@}cnH;-!$}4-fe$%!R}BP$qI5caT50#{EN!zAk#jE8@eYgm)~qAjlRb=ej@yNXiPh zxm21++NuNBVI?3f)X}MBQ!Fl>#Z}f`eL#+opyuZ?mNM0g!dQs@WbJf6BkdAD73~;U z;LjlM!UhN%=AAOv4bgJluC;B<7nXi|aumKWP3$RKD23UAcB5^OEyUCd$uAkxAf#** zlMSqL%sW&AD~>ac&RB0wjN&c)Sqv&>*bsxQ8C#9l%w!;7X4+^r-kN#c!Z%KtUk4_- zjoy&A%^Eqp$;nS#yF=^0gx_wo@#fJRSf>uGBlH9%8$-+3#yaaAm>$4j2J5_R?O0w@ zo*tjFn2cr<!=5SS^qU@NQ^D3%bTvk;Q z&+;vpdKe(diDhXWoibak1o6%WPIqqIn&EA*oiKXCJZ=fe=D&^n3N6F%zl2MGl_!q! zgY-)85{0@!PQWu^Y64@-x#D50D057#e7iH0wPO=DD=o6;H`9y5!#8G|<-p_STO&O! zBPT}WH-pKZ#aZ@#@W!02;M_Q*jY9UG4JMMFC?&6_1vzS-p0-RO^>-~i4O7R=qc&vP zgwLXUNHioOpx0r>A#hlt*~Lh4;(-eiOcA6l1d-9CRXjz<0t* z9bKT7;8$K?OGoN@Se2myGgqh{DJ#SB5a0TH7>8n?1qkBRs>A)O1$&rk0Z_L2;ISX}IMw8M@oa0h@jXPed>sZN?>bg8>-M?RrO#(eT zZ~}_yH|RW_;L9nK3?`bW4>-o1HmR(6Iaw;(`?URsdk`V#7@f=jMFi0K~O#a_c;VWX=8PCP5I6q@(gmQ-0NI6iThU<%eL6LFR zFgcw~NwlKYsGOoiyIxPtlv2w+j?pGxj>W#uN!g`*aYRh(l+rpqX`Rp5DvCiWFnQ{G z#e&}78TR5tg3jnFl5~w&P3enmMh<`?Hwv4*&AZ`OxaV2N-Vd}m^4c*Zbqsl4y@;oo zewm&oDA+-<9R>RQe`4u4-JO$1v&i$Jp$YF@%qKrHxz&1*4@Ou)DHw|fQx~p*J!Ol0sVctcDN01nU=BjLAjC6G(o`{{Ym&Z2 zzKkvMWp0sAr{as6GJBSaFM2w%x5%g8BACWrrx zJqKHPIiWcfHv5!U!R=4tl1nsO3QtvWsU?KsfJ+faUE+gteQX&j?x@y~mEsH`hut$? zcu5IV+OXNmq^N2tBr8%bLtC>AUbPm2euKD>x23oSzxDufCpn=wH`sa&ZKJAu*&~}n zdvZ^6Pr1e<&+g+!HIx?*pZ$OAy$N6($8{!*!8rg1_elUG1~)+57w-!MKoUFv3IIh? zhhY#sAV&lNs%IdPW57pj{|jmv3sk5zWYftorX9$lH*hS)RN`#pL{7X;vYBSu5i~I! zItiWCCfTLqC~>sekOWG9m|;yHjHQ4ys{nzE+s1AY;X}b%vZ{W#cU~On<>R7 z0MkHx3pi)sKfhe7QhQ3@UxYue$Pr||4)wV(rL-wu1bN6#`^hJs|Gjr>pD5TRIWI&{ zUEsc|5jwN4ZU9`-f$h>bEH@__K{Y?=a-IO`pCgOEV~w1f2kkzrCYszoIuSAdaVkG7 z`QkViUzaWUKqAYQUqFtqnbBpJ_Z0sz{7pP25ywbR;ZJOYV4~;j1n)X4*#;1&MSq3S zsRyT~PEVgD{`L58!)p>ulfeDmGr@Q-xZtss z17-rAwqV=|-pgOc7_h|R&H*Q_(w{c?L;TlV2)ZG=?A4HQ|6HLBGj+p#43v2vsFO~R5Us?xS`C?w}T*~d7Qorrejyd}sNBYIS`KJU& zndm6pYDYeA-zIOuYjK!Sr^a2&x{H&qhc;z%>~QU?TpMh+g#a& z_+@Z2kT5rT!G0$zf8OUWyEG$YRf<`aa7rtAEBTu1TF%R#fBN$9LY3ckr50+hUdXGO z>zW${bBwVIdzRs(>eE*W{bl~HYn6WfYK^d_ZmxH}Y`*JK#fANF)bw+MR~`jBjk0U| zgtd*ZMuk%qbGzrBzR+~0Xx%_(dQr^w(+iJWKKH^pfA?b5lI?o!ul5U>?Eyyz z{1Y7=H(~|H7U`Ej;%=vA&Tn3Lz1?g=0`78E#$$;uB?|nB&=C4WYh3yd><#V<>8AW0S9EbA{Ep`PmnuLrdF{24AKA_6umou_2 zo?pniIDJeTz50IIFMSr;QhjhYvYTLTy|82aTrU)Kq8lgl!IYaf{{$RG>J+kT(M`shg_>Wf|4jW$jaM3l?3%?jORX;#3)!7> z5v^R%sM2`9oQn;;LceOXZw_i!co9e+_VqQE-ObKO$nVs5p#=7Day`I`GB@^^$ zA-d3}Gy76dCw~TC3V#ZipiM)cVIqvRJ~6d=G3#dP`atS>aF+5`am&qOP&2&>Qx6=F z-!a#7@!;)srT%z4GA1@|zu^`e_X*B@!n%EsGyqzfxr1QX(sN}%aMWWEv{jK~&@V%& zKUK_VShC&BXbohv-pb6q>{u8^UEUdC?bhYB^`MMc*Rb5!di{wvhCySKbpSOs@D1t- z_y%PHeuEAp!_*Mv&*5o}2Np8~C|Q!oUa&6=9TXZ>`pnoudv{6rcN z4AtfR@~{Z!$G@DMTM zP;Vl90jkJX(8Y1Y>$1WVIyxQiuS=9h0I{?90c4d_+P(0+jEM!JFeVn6WH|P~Lv~Bg zsFH!y$=9*FPXa7-_e|;Nl0R9@STA+taj~&i$moSPRM><9Of_kJWxVn4uo*@3Ql*f; zK|T@(r{dtB)|LX>7F+#Ci|}pS8=3V8)r%B!Cw6L=dMUF;f~iFM81yGR)+0}=^Z^>g zV$q2$hul{UN}=`{!zIJ?Y6CkAXuMb8f9WGly0ihLGD+(nLE=<9$?Q=qjO!7#TnF^jAo(MQq{G_IbeQG^L5fwqw}z2zmBG7 z$}EWUDVZRtIgps1CE&z{A_1j1n=8`EoR0Xgz+>b94Um3HK)ywibIFq}eXqd^-)5hfw_`6{d z`{nHGB}?3^6GC0LSl2D&bux}FWn*#Pt%jiZG))(%# zsyE)O-gskwpn9WFy;rQ>D-`TSGhF!yQy_nvVBapkhC_bhavZ@2?u z@g|{QlVIO0+BXO6n;&SI-i02qq+#)Kv1GGQuvxHg5$#(7_AU2VroLUTycKf|Laa4D zB_L625PCfUDpTkL><|OPJZR1EM|r``*tj?@4({KJr%&s!w0?QI4%6db+6<7*50kxl zSqb;CthfkeCEmxf;^_(2s@1VqR$iid5;oG}+0_qhSR5fIlgJ4z*y>0a110%sQq_DS zD08LI?#WnN%*A{z8MZx$5e}0pHI%EFp%21m8dS`kfF%vhZh z<+=W0&xx5oH9Zl;8tHRC=R__Jgx4(J|51cjq8va&FP9hWAK`#HYy)}VuV)Ob8exsj zT#AGjAyG@bP8@$0xLxz9Oc5cbi{V?(9qDq=HBVC2{+d~JlNdY*3#1M z?DKfwD1Do=haVk>GaibXp{dDvbd;Zjjhii7oWrNUG!OVXH3hg+pmD$}o1UEloq>Ct z7#2zwERoxp@RYvm92mrbsskif=LDRyr^N0ZP?f%NkR{;`61eDrsNnG4nr1XL!zTwL zw%&J!E5u*|kNbnc#8D4= zQIK}`dNGxl6+R7FRF;QuQsxm%e-o3GJ3+|JEBi6-0K&#WO?YsvilRS$LI-P!Y%e#5 z>R>IK_#GG<9$V=|6DE$uIU?Js&9B1pjW5@1-8;uM`L0w`^}}qV{d^IkUXi0=2TWcKz{)^YY3fV4YHt-Fr^;_>vU4P>v`#Yn zNY`W8brhU&wdt>EfCrg3EBxypb5cK}zN$@qAkG@rh|C{TU-Al0PKJQ4qE?o$tvnyJ zk6gn{dAjTThtfV#@o%k8?`!+W{4w=qy?tbZFpKRYn=Ucs{p8vwD*k^0?IZKY)R)Ef z;Sxc9Px5&g`5h01clR6uxtbs4+|p5>GgBNKbqSqhl6L7%xFkV#5)LeZ%PJQ~N%A0f z!;v8n%%(u@?j~n~wtSFaeo!vTS3l$(1t(k}4dWIcSvV*Ka-(==#7WA96Fy*fMr6fb zhDjL0IBWO{24(g_xnUC;*yyiImn6wgqI{S%IR%Hr;fgUrPj1tz0cuI@JO&}K^@0OY z9@wXJNuuZvC^DXCp~E2V(2?w59N4ze`-MCL=_!|G7JH=Uf!6ZG!7Bok%jtOS#>`s< zwKofDLsxI07YDlocKpCKRiHrea(M^tdcuMz#&bH%dOrS)Rjb`Ww4TOdavx`HU& z`IkHQ{nq){KlJqvK|IKmp~8uut!BaABHCL5_7;#gR~fypH42pq|FpsfPI4$R{Ad^S z5G%rm&I<}wi{TG5A?Qw^?2DaMRnW#ffD{IeG1~2AsoC8?{i#19_pV)HTZ5!K(0}s?`fKe+PI*oF6_tk1=ZuAJE z0>p!vrTt92WIBX4!ebND(p+t)Cg!g}FH@#&Wj2&Ao1<2+4)@ir*==1|P z{r)u+0@E)Z@e9wiU4nhLXx|;M@4i(~zu0sA;9CvbZ#HZf8g_~eJB5N>V!mJB6>ZPV2c8J1GoN-#LNm`sqT0DKQoJ-&mSp)f=nEALT^JWvcB=9twfIC6p z0gp1nAm$>mp_F7ukvSwO>YNE++XuCWCMWttZc1cyFq*9uq@O3B*#gPl2bJ*nyuJocum)~ty@qpZcC zZIlvVuaNo$J{6&+sZ^*(+vhI*y_yR@}NdOvw+fb#otj z>*5Tm9Fe-8t!pQSuhv^PkE!D*YC^Ad=i@FG%-bfv71bEdXE?a8t;`hwwme5F?sYVC zDbBxt z^(DtdRbRyqxW43Zyu$h_Sz&$2e8}}B$8c*wOSR7T70q%69D;C+p}1@!36B1ZwNvQM zMY1i24PlNSeKzD5XVlqGe40CJ_kt5ck9%~oj&Y9TBxdbktj>dKYI{pl3)!Ofc)9Hi zBh)@|6TIc$iK&EDO*4ivF+?_l6G7mR_J%8AW&X}b;Y)IEjAopCJ!B5rPmb^$vdO;z zQu})?yWOxS4O`>jT!#2r@i@VI5I8U-rtcxm+(~dR1iR)$LHABRLE%FvK1%>R4S@dP zgLZdR0B#CEa@;%vD8BNc=`+YYR#d3!C!Mej9`d1s#`^D+Va77OjC?pF=5rr{D$>sf zH6sB^#{b#;A7X71A2_=KoX69mK|A{E@r;d5f}3a_>=H&NrJZWUBK4b4-~4eTxP$b) z7$Te0s2P;jqK{$R$Keak*pN+X5oy@dn*2sA4)l*;cnW5#p#UEaxPu8K$%7hpvU5$o zr0r{IZ96FS?|x#ssxM=4?8Gs>ks5s3b_Mf?!a6x@hkQGq2~?-igtW6`_24sC46J; zn`?!le!(#yItBs``79k=kzIyI=p=&{snj-U%J`7DZmUqdRd8$*9oqtqZPApv{hP6G zCcZKKjkGt@gd#AX=@%XS0Y`uM1fO21?UEO)RIzlMPz=Ynw~LPL0mt?SQ0i@cyK`AP zmJKtxN;_ubKTIED(P*`RBK%RM#hTDDE`U^Y@wCQ^vvDf751`euJ=q3YEjY9n(n6tK ztjq_|U4^w;37QnOwk9Fu!iUysS)ih|T2Ka1OP~9?-ly!gxkUQxt+|#8U1^Ld>quGu z^1ygYx@X>G{oDiN&3PW19sMp_R@2lyRPzDh&a4czq!Bl(I+m*p$i zxX8+jDBqw^QcAu+D#|KdF0V$3WD{J8X?uOo*Y|v7r&x}c$Y=Es@kPpsFSOG4H@x%+ z_L#JsxjlC~p#?$UP0lgdQ|X5e$DUSF&FIY84^St3?RzUH zQGM=(Y>+>wnP1`2{*dY-nsW5(0$0K0^T7H=@pbDDj74T5$)^b>VWQbIpC54P||?}ZMd*39e! zNX-jex>W}?XnMEmVEni7#eD*%E^-DSGxx$X7oPzaPZ>Ez0a!9+7tBAsP=2Z8a>?D8 z+~l1<#4qsWlwUB{GoN)~?{a?8LdTWuLVm6E8!o8^=OqV)0=Ty7`dXo2!;P#Pt*_?@ z`P)PJwdq3s))1}aFBUvC75H=gqt^-+N*5~^pT1tRR0)oqY8IOVYc|7MI<1&|!K*s0 z*ih||wj1Bop{($tODl$j*~NX*g7F}zaPTbime7v2N6wx4X;yq7&GA6z`zJ>(oc=kf%TAk=lNcZh!;#sGO2zE^s#acpXehXa6k zr54!@8HX*L$+6J=Q~W$+eXsO*xo6y~h- zxM0c6nP;twea}AF)!5iL+|}=^mX_%|9SqGF#mcApJGC{fXmSl zocATmmrwA@MH6uT5v6ec7NGIp$FF_xliZZhx`983N!1{!@Z>TEqNDg8X8L#BVL~r# zaKS0rSL(u?PC`00&Vdh5O9O5p)GNg3h&6 z8Ve@vM=1FGJ<(X;)=>T8I}5YZnW}!ZsnhUCeG@dPHz+6dfA_o6!k=`4cvh{TSRL zipu8fs0d&ko;_bFW|#XniP;T;qfvA;2Aa`TTlu3ZTy3@^nC+;N?QXoSd^h_OZf(YA zk(rY~Hj3Hxf}=rnGz6NC14!B9QnsUleJy#0d(MV6QWmNa6EvM{;2Y+kziKhhzfQ<) z5p!GM61S{8F{{j9FJv_?o)WS)+$a#TcFiT;wr7RzcB@|O3EhT7sF0W8W{}L*z5$hUIve;S*L=KH;6WU998cW)qn*RXPMH zbA*%$w@S~TEtKx_5Vz9FGs)mN5(YYq;B0q;jD``MNB z!f5b5c=1`wwekzbwekS_QBt_Iqo^e-{penB$1Z++2EA?AN9{0C$%+{~*#Wwj373ob zl^C6v!t-<>r@=uFqUAyRDA^l(*wcj~ugbd@=$QkqB*S}N!Slr!;7Kky$EPPI8lQDd zz(rT1aFG{t-K&hd1&Rbg6#f|&-4OEUF+K|S&4xr)$RB$BC@d<@JGSsH?=(M&*x+0~ z{Dj1|Nq&6(^H3SVq!0=*oZ|lpCfojZ{sTTl?5<1T@gN}oS2za~PZY^6cz!VA^16by z;jjS@UJoBkmY<04!t8jhGh9nHH8Fi=l9>DO-H?q(7Zj2m(Ztj!r)$8IMh>U!4sop# zz?ifW&{li?6kuV?C_5%I=lOHbox9+>=mT$|sTt32dT!H&Ef=>C=_*jSRV>_gGj&@a zb=z`I{@i|0C}gf#=)0LwC1g~ANK#(E*m1S=8tlSkXEl(oOMUZQ^Uo|7to^j_1>dEa z%QJqjShFzz*AvwoxLI%@P;g+m(79YxGT(ixvT<=-sO%IgJMYG%=K(o$e8B^^5fbLdC#>eYv>8f8sOS7q%~#R4&H8)W5I~u8sC*`MdnN3;P%LU)jG{ zE0nZh@TL9BC8hpNFAXeoFLd8t>-4vM=GelqWoIq8xNiTN&IUIpHhH=Uf_w|M;b0*fXo#DprCU_B4EA_=^0*5~Gj!cl)2d+T$-+IJ9u+ z%Av($LP;kEUpf@cg*~10M#y)DN@LwkiYfCS6x^{ZH$(mGx&bLs7NW zU1v;g;oJeTSqRe*+DgMT68^}t^}oT8|1@|*2l+9ahF%`zFlUn_V$^fMo(Pq~gl&=} zNVFk<82%_SDV>{@HvClnFddvvlC@*%Ju5g8DrW;!d6&tL`b*7+V{BiPlTOC`(7Tz6 z`&r}sZC@Moi)lEpr%H3bIn+(QJHeZH7X;ms77TJmLJ0!HOOyLqq6P{_>{w% z8kJM%T$XB7revSs9{Q+^w4xf7QCNNqYt$gWL5-p$7eAUwOL8Hj8rnj%ERim_@;N#zUazuP`EWn)-xwFEPPm51J zeM=Jq>ZK=b9yBg4Am7Pea?_$0_x&IYV&hLIUyh$rzNA_s-;xQC?hE|`@ET+dM``sI zygZ>(NFSJ(!oe7hx(fpG$jD3{N=Y({e)Sb<`YyIOP@5&3CNFhQBHSmJD24e7dYP@h zLk*lx$u9%gAe3L~oW{R9y-~-i^k9>bGPzNuC77f!raCG*#-Nh3bII%$3@G6SQ2oZG zgcOIBTq-El>?jCvpbC^ikZwmc5a%FfNQt3X_^WI*K^oA^BuQx(8#vgU3*n@OakzB2 zb38hn!?e>`Wpj!?tO*=b%0uKr^|GN3Wpj97Tkai%W9ry4{^ruPhr;2$IOuh+$UP5k z2{mr16qB@6lQf4(S{l7&s@sa&XVam^Q)Z-1&{ViTQVKH=;>Bq;b- z)6!{D+Ql7_cmpHrp4oILC*w4!L4FGe;_);h^bZOjj;bfpzoAbtskB4x*$k$C$Xkyy z8LAe`fbug>lX1xpLpUqMOaD)T7lvopIg&oQJg04Vt6m0!?fm>xOO^Qg>_$ZUq|ig8LWzC}v>q;CqJ0 z^Nz%%n;d_92d{19XYeZsU+>`B^Bru-LCs>aSkOkqK==uJtVLpB2mDOCot`bGuUmLlOs|=XTV{Twrq5SCw{>poLWMu~ zGxZDgw;Wl4><+=vDLOg>j?P;S^j`bOg##B4%pCx|p#^g{yu460c1usL0rg;V&0H7Q zfMsUQw>~#LHw}hkp+kK|PAPD^xMbmxm)iaDFM;RS@bSQ>uhjS}UkaTMgfn*3(ZMj) z+J(`VI{djW)f-YF&)FBx-Ac=x&zwiC zqT;!|^IdQW7lmArx(1&dQ!|y_5q57v+OR;N=MnRfICH0iEfUB7e_I zrBcM){&($6VFk(kAjac27#4#2-md4KnS18;nvw{o62&z;<`ZwNDZaACzX7gP6xXbu zPlVeMCGv%`g{MDLy-*Ir*+bc*`CJ zAq%lRuvQ-WJz8;jhJ2Pr;FCGQ0GwbSCm)55b5h|lNxRHf)1ora$;-sZ7^?cVPu6nf z_X~~m3OiQ;T!bHloND1JF=Zwmze8bi8>bXzN>|xFP!x6&uEK9vHpw1SV38~@m{m*G z6l=M->X{Tuhgw6HthG#B9f@1d#MNV*7c4}=AI0Yj*TB{?*9g{AD$;ptfgUQ#a%C~{!quLQDHZl~OR1}+$g4@D8 z=T;K8jfvY%;&w1`J28&ABz_nBY@C|EE`qX~LFvXg=zGdM(L;@$7Ot1!e~*g)y&CF$ z40Rtvya6Tay2A=i$#Ed0!*z2y0qeNy%p_n4~8$4&4n?g9S^75e2qV+{V z`ymzWOcb=UD%uYdTBG^$BZklZQ5DX|OmTj}5a;76oKL86m|0C_qpF%OKB>a`|1?)Joy7Elp(B%S2{PXJkcr zD`7@;J05VLaQ!0uLcd|;<82Oy$Xip1Z)M_TYL*j_NkNVE3j8m9V8f9r`!7CF>to7nmJ0pTG_5L)^b5``2a^_|EAwBuzoP?7 zRh>6zsZE{$hK**7NC-?{SCwIwrb8W}eDz3?5SSjXD&H(khdM&}>hU5WF#TRtzFC?M zb%gTOqeentdcUfCvosy*2<4+MUwXcf5STHb%3+qKLmi>~d|Dqclf}M4xGO4o;B#sT zVkTQf@yjZTU(rzHs3`m@imMumTouKyswh6Mp~zEFTvJheK|_(RqWGeU;@1cTE$)y z)y?$hbUthPbDBcZpYLFvQ>z#M9}#LltQV(|iC!Eu^K@KTGf&4!n)yFv@&!mu+c`ni zlWvk0?V@_ppQ$(#RTSSL6!d$_xX=H&iuSE2XqQ#AZ<8m`?u?)gn|}JEn2!rGanbJy z(1_Od1Zawt>_OvoTqlEAJBXJpmer>cOz2#vx?$h2nF3Qm03zHiTqea`>z^WMIv}l z#p(Z5(f*qWt$IBCKNao2o6!CTTkgD#gERQ_NQ4+$RP^YXikc&UK2Afg)R{`YfZTAx zT}7V|6+L>bV#`l5p;t+OP(FI3V(F6!eXsH^V^*G&6b{rWMxz2GDX=(sMYcMDys4U! zQlcoy5k*O9nvxu(Bx*d;s|nzZGHccJkPHnPG2S+*l}saEp`?k7^CO@YMuk=w0c}lGXlo*%6-9+sM9}EDOn%2g?H6m>y+l)2YmE_9 z{b`*EXQdikr)AtS6WnqQZiQvsN)y~F4Q{n%+!_nuL&+%c(HGKvYxZVFJ=3snJ4=+xce>h zd%%R>gBo0EJ*Yc_tL6todoh^l(PE6P2Wfq11=nz`2)JlL#^RFoqNVyBHmQdr8r(-M z*Y8mi++!Nt26xgj?vx4cSq<*fmT`F#T#p9VYZ-Ui1os&Y z?z5J0&zaz!*WhYx(v14=hfHv1G`O>t`Tej7?ngAZAGM5YxEBWV?=NU@KW-VjS@!F}E`zw;)zmo&I9 zSjPP&6Wq%h+!rn5{udKmul%b54T7w{XY6&UQ(!DA`-_}Hpm^jSC-{_}jX^cH+*Qr9oW)s|KK-WR?p5~r zxz{vLuiM_nh)~O4VKE`*Hwdc9-XG-os)p)Sq=NG$bnP7>amScAtxT@87;TGTwZ-sU z2VLfq#&x2xI>vJyG;Or!I%rz0j4eNxrra@AxtxY)V_KJrcHEejUPE%`=Q`*oaa%+* zKi5GcPFO@VKi5GcPFX}WKi5Gc@)i-z&vnp<(-slU&vnp<(Vpv|X|*!6aUc7Ts*go` zu7fW3!=x6d8Ifm&d*oJvhnb@$`Cp>Ld70wLx`VF+Qgd2z88-D#B{+(SjnXSt8 z`+B)%3jNh#B>L~n8t%`a&fZkZ6wDb9BVRg2P%`ahs3wf8=<}lXR1^V&)osm2@>CQR z8mUh^)a*Cu*1*{N?yIQSTK$3-q9W=~VVpBrn zNERhVIEGdv#D;{%Ar>T{C8pzume{}x5uHyu!EubCHPa@`*7C%X-xi+H*m;NjMfCr4 znUemW9)I$hoY5Oaqpha#qt#Z^w9#s-Y1(MD)iiCi+G?8i{%EUd#QUSIrV;Otwwgw~ zKiX;<@&0J5X~g@Zt)>y9)mGEA(Q2z{+Gw@aH0?^X)ilmZwAD0@R$EQaYUbK%+<&C8 zr2R)~4rgXCb1c++(^%1JzG>QMHQzLCw3=_4)|%#dSf>KLsV=r=WydmHjr6uW4ydS#7CL z2cn+>nwG5Im>N`QRj1U*N{yz~uJV+Vg?fGpXxiV^@T}14`6-}jzoDU3cvhcE{5{q4 z{=NyV>Qo?yehTPvze#B6lb&h{^b$Z*{DFo-&rbnO5%yC+$9+2rp4A-u5z=yh%y9K5 z828Xt*z}W_t+45Q)>hawg|x!n!91sS<*etYfJOKg!j;9z7M%rM(peEqW@TY2PMKpzo=G&L^#N zgA8r-dIK8Kn%;n>Sc%?%##xEpfW{$u1JLEkT0mV-1vJHXv9`HeOf7y7s9O#Sc~W|Br-%euL@xDxk}FHzEc4DxfL;g;40vQ)+1heHGBO z|Ei%?qyjx(1vKrynb4}o1NthUY5(1XR^L|vMGYg<78QLgTRPCkY3LP+MkQlFZuC{a zmY)z6J^CtO>61+8RZ<|7kG=|6`eZ`?V0;y@C8b1BQnbDb*peKiBx*duJTA{wPz#ZK z6|iXO5zr#}Dqztv5e<0AR<%P;C_YyepBJI@NWKc#(pT!MfJG~ez*{6=1uWW{2xyUf z6|iXc;;Vo~SWEb!XG8fdt?EzTQubq5?LM?K`YK>?otAO`y9sW&2Did8?zc^Ft2DUP zmT~{c1Q&f3u=Su_zZ%u=ADiH!uL2gg!7{&pVuFjl3Rv7`%eeo;1Q&f3u()lOasSi= zw?l*5X&Eoh5$%+yobW6|lJbE%W=92`>67U~#1tp_Tb} z*#sAT6|lI{3egJg+a|c^tANEND@IH88#KX1Uj;0#bzcR4VS{nGVU!CT=Z4I;*MCx{hkRf`YK>?wL8;B^UhzI;G(Yr7T0B&-`gg*=&OLmweG9n z`zE;PtANElZJFP{GQmY(1uX8QW!%3u!9`yMEUtB51%G3Li@pk2Tps~7K^Lhu{P=t|6qcPz6w~}=PmR5 zzfEw_R{@Lrf@R$QV}g5GgZrXo+`A^Y=&OL`cfm654^42GFI;G(Yr7WYe*`TZ9Y zTs$Mk;(pl}7v?{GUj;0}k_JK6-!onXtlAPj#`3sgUj-Dm_sLfQjrf}8>2=%N*uns6 z8GRK{m=N(H0_5;Eg0-n zaOS=WXq=Ck;+Xp?pm9E-#);Bmlg5hHVw0wQPSZQ&xmz<|s3hR1eHGB=sDdyER zrt7PKrudwOLeEzLO<~PqlctblGrC{t`6{3(KCj^_kKN1QuBrOHHH%HU46OuaJj&2l z0Zsd|rbiglqOSs)_LV4TRdSe)uL2s`++vf?skPXoQ_U?l>6{POVv{cUt5}QRtALLC z8mU=*Uj;=oB^qA^Yia9Dlk@2Mz6vO7PeWQGUjQ`3EHS$$pg?1(78ucqH zXMv&q04AO+J|0k-QYp)CaKah+Yb4ij|x;pm8Ky5_+9Y*MPZ~0vdxjI83Lo|?vwR!>dSMysc$X`|Is)3nj* zscG8#qo<}3?~k6EM!Y|IY8vtW=&5PM`=h6(5$}(lnnsLPPfgQCtEZ-Eqt#Q>v@6k5 z(>N>9Q`0zFJvBY6nR^yM3lADgvhSegaApQG$3ia!G*+}aZ<;n*oi|OZwfT_eV~xHX z{qxea)^y%9g*BZwO<_&vO;cFYdD9fublxbz-M z$?||25lnwE*Ll-8TAeqYYOeF9b3R<1H(fH(c}KSCHMSujUJ7W$*N6oEn@o)oEqc^m z3d*APQcxdW%(1ISjhBKl@KR7lt;&8|Bkd`xE%muT^in|6lGPhig9@$cd>UD)(X`rC zo^qy8&r1PK`@0&R6hp-`rGTdWeG^*Mc|hW&fTsN>p{37wswvP% z08Q}+8VWrx1vEw2O937C?I?IwbMQw<3tkFnyg$LXhql6|OEz0!)A_8euxScug}sA$ zPVLHB&r1P~5?(peab)ETUJB^A@Y>NOmjXJUw9XAOw9)GgXhdsz1DawbdIK6~C3*uI zhv*GJZzAjWbiEYN6yL?#1}_D4-1jgJde;ACVsC5ORI$NO&w}WsfJThgO94&$H&K*( zM^)|*2(2vpJ~S@{H2;66!-FhKJ!CHhG++Nf_>wt+UZCftfTp+`kpjIG&=mhjDCjqs zo|gie;@yZ8=%s+B_!mNVCj=g=v7i6l#gBtSo&l_|6sfn zuqCBLQBt&C3fPhyq$FxQ()w?xg-BisShVyAXpy`WuxOcx20Ub|+94+tpR0<`i%@zb zF9mGrEA>*qq7_EqEs~c47Hv%gv`AhGShRcbQotgtCH&B{q5PIs^{0PRb|cyKAG9-i zDPVD(mT~`|2`+jmU~wxf6tK81%lzh=;G&lT7T2nmf>ZKr-D>T7HF9j^FRWAjxoSA%!UJ6*;XD!!*nWT$e3Rqn2&b3khEjFnK z^isg$&RXWT!~_?;6tK7-wT!#g1Q)#&u(%($jJwVR7x%DPTxl2E%KM_!1Q)#&u(%g2 z^XoLh#WPwgu6D=Ts6UjM;G&lT7Wa9}{Fa;GqL%^|_XW$i6(+dorGUkK(K2qO2`+jm zU~v~L<5rpAqL%^|_p_F9t4(mxO96}f%a(C#OmNXl0gHRpGH$I2E_x|oaj#j%tuw(z zF9j^F)?UD9T&y?2MK1*`?w2g{TW^Akr{h@MFB{{6T&wS;fJIo+Ajtard-hU5m9Szj z1vHk$SvG<_0FC&X=IM3Y+ZYjQ8NC!xm=NY)tD?(T*F_(mK2dUJ7V# z+!hhdy%f-h6BZH8y%f-h_s9B^M&vE>XzrzeMx3^YXzrzeMvT__lcv?m(8hfXy%f;2 z(OQ4fv>zt5V6amW!Ak*+^D$E#b1wxn&L`A3QCe)$SkYQ+(zMTMdWSrBYvv1;1RS-O z0=isl7MnE1yt>A8y%f+CpVLt2c`2YNtXXW*6q0O4_bWXw1vJIyHC*X=DWEB=S!~i2 zS_#Z}l%babn)YQ)k1(c1F9kI1D^bv@QnUJA3QA|nG+qkIY3odrlj!DFupn8^J`xCoR4#Rvox6Wt3Y?6Wf(ipK z1*%*#9ooc6dU zPdX2uCo!XwoO5hy@);NJ^*Ft!Tu%3y(UUIcD9?|ccY3Fsts73bz0T=LkNf1Li*vds zy{?lk-Z|=#OBp?TmY+JyyGOk)=lB%we6(+guL5N?`B4n6KO(-Wg!_ta#Q3>S=@oC?PG9_!gN5eHISOkMF4 z*)hTRq2aDUDIz856h<8GI~>8FLTF-CAh=HNqA(Cv@UMZp{MB+I- zBt`5AkqmDeAsB=GdnJs{+S55PK1bE@sxwtw)xN4jRYUx*+29#%!6c6h8j17pFT?L( z;^g$16E2=V4L26?=Z(%m}3{|81TjT~$ld9D(uw`VS$iHnjc3j@a@_E#Q(MlQ+2ZjudFrt4hkM`Y0fJym z(7vbdU~hldp;(JKMfD#4dv+T9~vG!a)@B`bq@~tO2?+Ac+NeEe0t7% zysk4?4LD||PCVrr^9JpMT|K?sL&t|s0-rH4yMAR?Fsb*EBYph?P=CJcF^|VPa%Rju z0v&s#y{VMi5d!bV3itFnM zCJgKz>g(aZ0KEG$LM04#IcquBNgx{R8i9iNFJdZGIXEItns}= z-HomAgas$3-B4#Ionz3YCZ-0k3b=Dr*-%}T~_Z~ZplgO~Ie%#G_ zyuj`W_oS3#bleLwmx?w-*4?$6knwnM&et?Lc@ier6E5$wF6hzNHHM}q$@DgL7AE0y zZkX%5&M;BXK6IqNtAFpo7GJZL^3c?j8R?;c0igK`^QIPPosUV2gx~Khkw@ zSX1Qr5PLw?KgGi=eRv8Uh2i$_z|aU5$$tUgGu2Z(yc$L)oKUB3oTcIA;5y-4HR`!9&u(Oj4EhDXL);5!@ZEV(dldSFh zN0IThXQ$xBj-wQYJJsaju?<82aQ|=%|1r$GPR@*@Uz^*GGXZAz706@81-*(IZHX}2 z@V!tMX{;fgG}01b&>)Gbx0?(x$Sn`BR$n7KN|4TIfDyvb@M`s~X9*357oiFbE=_S} zgNqP`h8Dc7`4N1pLq*5%JK*cR3Vf@(QCjMO8MJxWJv!O#Y9>9qt)+EC%jUMu=J(>8 zyyv_U(Kx&S_&OxV6JPuceg*OCFW}2p86FaPgBCcUhHUG|G5ze}k)Y0gTt9nw0I0J+ zrpzFE<1>_99xcd@Jmt^?Dt0J>mnT+AlsLj~MW(c7LlH$LL0}kBKp9NvJv`98 z&zB_+QE7PMTnK}6xDVwXsY4L~)84(WxBI}*kwd|_^Da-Fjqifq1hP{?-=VI(y}oqU zIS}DcWC;miDEAHZ4fOL9*qfh+FZaG25Z)4rjM(a4dQ$$;=F3%t%Cn<92s-e-^#l_S zcMlxtANFM`sVDiV>9dfsm!yDfIs`&{uRG?W0G6=(vEIP}6z7Jz`}j{oiTEVLC?*DB zz8QYBcMuc1dXB(!cj~;yJqE9?=}B1J$g-@o8Udsu{e8nH$MFL)98kq}x4xG?c4~BT z(lz0f*a;?e@9XOC?>*>CgGChz2%&hCJ%@V-`v*Xj0GSfjRZh-5J`M{T9>ostw1=#c zV7Y@r#)SI>PXyPXy{B)_o+CrOgT9+a7J~?(S)9+uYLFB7J=M z95lh07Z!8CwHxXw^lS+|y~9KOSRvETyq6|tA`&$SfTXg0m2x~t$Y)*fcqnnf(}qPX zXYBzHxFGv*@1esjzH*cNq-RsPTa9xM_IC9Qw=&pmd^7aYplvRFIA}+V2pH+a4BcpWBk-)2p^b*fvuMl;MFoIeb8!assv z&%=wDAH>i+hK4Zoizo$^lgR?-4&F8H;$4$tE~nHoD7B-Y-^zc=TE;e*X@>#Ww+2_- z2nqA_(e5E9=v8o)1-+6#0_l+LQOs5p&W3LQRlcLh;3X`##>^^* zWiVOe;-U^aZyo;vKG-ke*Jb#UmaZ{B-s#IZISO?(;evOqv)0wv3^OOZ9QJng^&fon zU@+w`uaXh4Z)TH|sQgC_4$-gn2f0JJRO`80h!udDr<=?Eyzd2|R zwHN3dv<=YOqNZqVkttxx(V@_^aciqg(WXnGY30^7nW9~nLetKz?J`A&E`_G0TRUWm zPF)I3Teo(~6dQCYG_BpbA((t}gy$xw+#Xn_g(pk=>2k`0qD88)x?sZK01gP51eR+? zBvO@yNZJ&Vc4HE$&O#*Z3Q31CiBxGJk`9HW)0jl6wGc_CLbAb_M5?wB2~}eojY*_& z5+d2Cl(NZ~L>ej~l1&Q9W@8d*w1h}D2UFyxY%!z>jTeaqhKfRC-2Sa}`@>*SXpB2R zE8PKLxF|HnJ)o8D0We?`8sjd|N_PPmG761xA84if01O(1#<&x-(w#sWHX!_+aF2p0 z0YdJCd)y7u2K@FoPx7ww()NHdgs=w5TL^_m(9TPvNJ40{gn;$QAdo>NA+%dUz}jUH z$QY9lIxHbzT{8$|$VmvDmJqP!83Zx{C4>!@5U?E>1Tr8cgpHOEaBMONWPC~pn=B#V zP-PIvaFq}?TSCB*%OH@^Dm=Oo900v59{|qr2%WY!2_t%@TGfBO~YGc~qHj-MH_61LWa3^I z6H^=BnCc&PBI`$noL!)}(k+_DF;{914o9LNL_KlCD23rj=*AmkjYVvyDkS0_5moehr7>R!^3^UM|zmY z;q!6D@*!;6lxENA^CLvd36*FcJLMWX?U_D9Z*I1Y5N%8w|9!F^^4U+o{&6GdWcj~? zh+q;gElR=xSVl-Ous)CiAoxoGkm#iVi0V=RWN;|}!nG6tsaXo3G#m6J46d4M^)XImYLnR0=o^UL@LBc_J6@-H@r-y?S0?e}ESeP!uL74Z#L7czj zFi0747^DX|3{rs{hPS;OhPSvJCWVH^ES4k%3ENVTuq_1%+ftCQEd?nA!nPDk*p`BX zZ7C?R9SRe+Lt(;pC`{N6g$diCFkw3s-r&newHMzs?gsT3EW^p>D$EpgH)@x zB|#5r=30Q4d$0KG^GU@r=VpcjQg z(2GJL=tZFr^rBD*dQm6@mFZS`Q8)f8$YE#>ISkDqhoL#-Ff@l8CWVHY21}BHgl#EE*p`BXZ7E3DmVy)lVOxqNY)e7H zwiJ}u4uuKZp)g@P6eetk!i4Qmn6MoRZ>VeFe;f2vd=-9`qjChjVe!9#KkM;l2YwyL zIK0+_zlA@~VQ3$I9m7ej4CDS1DgFsV`50P*U!TCQA3$6%oovU9pjplc=%@HTB>r2> z-HLfYNgl(OV(2*p{63N-f?D4TZ^Jtn{SI#S{4?YV#-Es)n&2Npf@;kBJ;StCxDq?%06PVGwkwsG0NlSf@!c_ zFhVxZal3&38WbN)#W>ig!5Nj=#o-j4ejqt$~xgee| z#`xq@-1LMS3@&(not)22<#Y2TOg=YX3co#-Pyu24DYL(P1=f@UjaOU)(_Yw~^Z3@l z%o8$Nn&O-2q~vsf{`JY2n2*NX&5g+~noGHznsIUKLiNp5XCT$NoLlNo5OS;Ldf%_n zcN1eAnGfuN0DaBcD@kt^*4`|vU0f#=ZW0SO&80zfdiKSmFCO=M#hQ&blEs?cfxK=Z zuUmxgwC=gs<@~}=pMK%Aze&h%74uu?QX^1&vu`=KP|U6LKPl$6-;J~N+B$z2gWvx| zz9F36ITtHE>$ZjXD`{eCwf~%$+OpJoGqod-+HtF}c%kRY!NuI`YsHd{H|%0bmr%G{ zEZjZUH`jN&u=sP`{`hN-#o_Buiq5S!_KD73p=6I(vPUS~1L=7c{+??GguLdZI3cfn z>5!1u^+tt|*Zan4A@7mE;1MD3$afw+9(Zg-c=UtfqaPIVJ{TC~ggkE9k@v$mTYmm> zekuNS&ZWHTh{?*koFQbEiJ4_{33oDcl^@Gxd*%`@I_8fF_I1l;o0YJ$n}m-F_O-Vv zs;+HX?7F)3+SWi@f1rGTq#PIQrSR~X+0UPU?)-%hUHs5oJk-~x54~_GP`*jX-7MyA zo=d*%$eSOV=jX>R_6JIKLnS(UY=rtz!CrE^vT4z?6uaoTHn21psMtgj9uw@&<;vP? z2bSWNvX|n-%8tOM{eg-DN@B%M-SFL5TcyN5LfomSUCh6_{o3}W-GTCrgn%#{Ypd8J zaX@~i#H3}3j*^APL`Usn&SLpuj_7DyS__Pp?~<}%#lT}q*^Y^hYQJYO*6$G=^?}yI zm~T+Y_ikcLd6oaj)uL-ff##h7=Ppt%R^c5-)_kRyUGCo`W;Y0qMkpoFiuF+*st+s% z$kVc=Y@5aGM#0e}I+_A)J&;Ysbmd{GP^q#S7snR)#j$Js%a#3*mrlGJZ>!iQKmDD` z>TCONIdbQRFCP^gYemOe=pj|5(4VVHm&>;!tx~t4%sKNtm-h;e5~wB2ybd!%Pf-^{ zaxXhFFi#fd$y&}u#w$a$P1*#DSLNx&JbRU=(Zj)dstnbWWgO_UTb3P}*tkbUM-}i@ zSp$4k)-10d1iXr2sZ8tzcOF9CCq&13q1XjpDOjyycB&%=Gqt5fHlSfA3EFsjbz zv%n;y(|J3F@Hvkt&q?+1En;?);Aj>d&4JD%$kI`&8{%x7t4P}>>U3qFT z`=!Y%lXE>62bRmq5oq0=to(VeSWx%n_C@dYj|lZU#rmCZWC>Z_VpjLvn7HJ8@-^3W zI|E1gfWOQCv>!^T6Ef-+yKZJQ1~Pzwy!nK?aWR=qw^JSSw&yln*l;VaWMTNq$l@bQ zw#8v_Z5zCuM0i2%7V8=J|eWf?DZnVChb7S;IzF4vq7St>3x4YtNhM5ULxnEe{N0V8{mV)|?^RT=(stVl;Ru2LLq`a8s& z+QrK2%{R(!biH12OwtRR91H<^R zxiA~`+|DWZDTcey^R7K6e{JNJ@48w6-Kbj3u3PN5zVk*aw0wtHziYX^nf!X>_4YSx zH#Ug%yW!~@x^2rn``_%x${DgP4?T*3qony^GLW19ryEk~2#;KCU$pr*h&k&Qd#(@M z=z@;WBQ|s`H?)vnqpx?rk#%E_*w77++1O)S?mh742e1}~ZOg+)F>s7@hGSG``00j( zE&-E%Xx76iu_81j;%k+e5~E8bu$31MTs#mc#@V?H3M zqqkx7S)^yLy4GpgaoZ@`>5&1xGzfl`R_4s;uH#(Je>L z{O(Wpz0h~*z~uvi!zntP&|%8UL7XTrU#>(^t1KjH5qXmmX^xAI8pvB&hoVQ_a>GvK zqD$cd`fp{U{QP$vF^>F&yql@(0;%h6!3=lk(!k|`x!$?nJNV-Z{I&B-xv!Q8H9N$b z9f91PLhep6cjsL1U3*MrNMy!pyj5O#t#+~PYQwb#5Of1&(rN~!50at@(}bquYSXo*rM5uX2KiktfpuPOzSa!le3>M=NJ0RX(t5T28r0jKK-pdu zNhpOW{|<=gqgSi1RWFwNOO*U&m6G7TI2drY$SmD1Zz5rAEr|;PxyT2~+GN_>&V4cg zRKr7G5#7_u$H6q$GQFrOS5^$mm5ocKfr>8jd}MeyGUs!Nl(S|J4xD|8lyiq_*p}Xa zQ(9O`tDp3X`{o~!#`8Apemmq^h2B@yWW0d6U8ab9koBi&8D@{EsJYg-6uXqQ6f2gu zFFi^3a#(nUIL8XOT`sLmLxWzLh0pB zs$Gn|y7$`NK;xl6MZer5Z&lXl#$y%kTO7Q4@Y+GtT2x3{3#>VOk__mOgBK27JSa)) zg_V3 z*ij>F6tn9EM}z2S5SorcL1o9K3dUN(=8?1{S07?F(g0iDCJK&9RB%*gyLE+>@bjD( zMi%tSqEHWhTy)eyt(Dh9t(Dim-GJJaGD*8aW-EN%ujU2S?n9`DkVdc^wTI;mZ#Uu+ zp-j>NYue)HIWN2?CD{t~zAuG(a2598D(JzebFGwgu9qkYYBjH<-+Z9;5bBjx*n_K} z2S1K&@uc)tSGUE_bzWS!zdUwz=e3`aA5(9iI0 zle`G5?>?M&WzyB&%sY36!%2jVWH+pTV1r%vj_>)t(-mYknH4?BDnz!ci zL4Va^o`0Q?+al(+fF*=I?V{r?yYr^q`L^THzz0tLFfJy!4EBDK-8Ps=;rlN6rt%48 zoi|{G?-cBNMf=`>eebeG447$?)1xnK(;La8wEd!ef55)~0hhM*jm@OAL!$jqz<$Vb zX?FFqx#XHWNow*WsmYTzrY75OY$T<1i}vn-y<1(H5ym~$`q_A_eav1@F4&f`tzl)eTj@Z#!CFFOUn~O$y#k3f`?M*ceXz#Kt+P zf`zgTub(6hyGOL|3E20r4Qq^{EzQD!eB06V6{TUj$YXYq$L!KPrV-pd7TobVC}u+w z#$HnJUQ+PhdoDO&rD27#t#3#J=YVKG5U?L$>(?0LZVF=ph^A~y8*Wx^4peTwG5qz% zzqRM)_QAmR!R5R)VqVqNGic$qU&!0>`m;jbo;mpMwv;d-pFLns_HG;iGw6_H8Eh-x~G_tt=HhJHFU^aQMZB!nzf# zc6gv;UbVx+_0*~zp6+(leY4s7*9zpR?8ynjt`ex6NGqA|r;HCI8yGc<{N7$zvum;;3t*{kL%(xE@YdRR4A!5*wO z|M+Pz+xfXRL7~inBG$vIxytI{XtME*YsGs#Nd2>~l$FKyRdbd5b3$;{TtzlFA<|!z z=K~pYm6TO;mB^+ztL7?ZR!eHD221u$Z-c9qRdW@}PHWX%MX`g^_9zQ>NcV280;irS zk5v!cT;+~EO-Qd6>@}jjCSb35zge$C4dev>zs@o`y zqm1Z%Lh0FUt584%rn7sw5<~u zt1#YGzF)4mlUX%aF?OJWFIsZIOV5^N)w_!1;da%#%1^+mWz}3|)m#P7y+mXF@qWv@ zYObkr&_h(V?}Tetp%u zikY1P<2qv1yNbs93@BJto?|TSy03ayxkt}D#ugt`PhIt{a?egq1S@`)-gdn2ysP9s zTyqt0umb=1D!Aot46f^wKXwSk!T&$HUK6nu-ob6X_v$60u7XHS~#5^Tuor!zK z##eYTcj&{#eG%^ZdjazQngoT0i}$6dQte7AQ&O@L?*m{i z1;2gp6kMt=H716w3E(9;T)H_#$y#8SZekca9=Q+GLtxyJd?swVe4>7r%6xblLwr)@ukS68ja<8HWl*3jAV!FWCmJh zBuD-Q)FZj!U*7s~i1b5+Cgv$Q4;OD-lT4du;rOtGm~9a;$1-BBMZ`SIi1`)~3oIiR zT0~r98L`MBVlg7l+a?b&?G<5+!t4Z8Y>&Q|k*e;E5-!sxsnN08J1vc%pSj-YZKP(R zRJPYapC|~YGVhU?r{wsj3d1?nv`QYo@T5I%bsyaG7FJB){^gaD1XeR8Q|*j%B?lx77-s3$LUdV zT>2B@*oigb-6c0d%u{9IFs96AaAhW`nUlRBi>AT<>8AN2p_Jbbl$kW94$JYKG?UK6 zRgkz0Ca#jiWioM97{^sJaWy0#Ew5K9!=1@lrNpCzH=8w+^;CJdR+*f&Axv3X^W)9m znh#eulg6zlqoE#0!%RB!i$Wid^TB#T-NT{GsxT-I?`VH@_!s?GH5#+I2Je89&Za&S z!`A_7B)%KfoG{XClR6a;VLc~(_+#)Z{3C{I4#~sFeVrUl1z7%ObGVijbKk1wkKw*e zof@IU-7_)gY)`EZvpkc-wNw2h4Oh4$N-E}+WoO9Sp~PL0r}Y96zr z;4wQ29Kg`Qs>;J>zmn?CA_a9 zm$oIU(&Toz$I`Y&RhrVWI?pzrNF?%bwZ(ATq)`O#{y&2mH;XfU>p+bshDkX)nGv@! zab!m1*mXIX(NuBU2`!s%hfds1ow!{X$CSEL>D^3xmlBUNC6l&WgR)D>&Ezyg=|&Vr zf6x-HpP4MS|E03!o1yk-P+N^qEpyuY05JDxFb$vHGM{@9Q`KAc>BRNv#O>FKJ3!*t zUi$yHckQulUFZE?ij*kn$|5C77Dw+&Et>(DOkwrpsJdx=~4j)PYgotk!ktPJ{C zuoc6;?@;%WSE3asNP_|mdGx#Ad41>H?|kRG=cUUD3{%@)pL(8}Gdf1iS=iDR-&bD6 z2GDY|ee{0HRZ_>ap%6$?RtCm5x!!vxC>QWruR;&!anZ3_tf8yd#|i? zzO5|R@2B^wwUqW=mCfyu>!#A?qxZ1q4`;(?G*M-I#V|LIeF;`rzpCsi!2DW?{fNp_ zti(rDUe>Z`LPH=uS&)qH5yi4pP|MPwl5NI%JfaHw*m`_yJ*fU{9QmlUa@;s4;r;-_0TS+UFqO94{Qg$l|a)>QUc`Q#Ud*sZd9ZHZrbcL6e(ys*N z5L=e=c%D)okTX-BI5^()*^=kUJUjB+p z+~kX(J^N%|?iFgI^!@lyZj4N;OI4SWJ0{qo1m!vYfl5$vZ&8Bs z@zlb(VOc7)l}2P=?o}4LU6D09rcZdtGCjr9s)xr(Ln z`RFQr&Wr+LmDX2cl@N$kT26^oLR4+(-&iu_{qQAnC|43Py`@BEHr6ZDp_0gKW-T^5 z$5Kf|j?FI-IaZMf0?U6ueHt2`Th4~2PQ!Fgs7Y(Lgw4p9uv09)pMg=Xrz(LKNX`tkM?&l6n%3da-2}27H0Y6!Eseea;1* z_6OYqjWg|b4D|V@H-tvcUz&o980z!SVoLHN@*ROXB#RNdMjqws+MRbC@afqJOqyo>kz z?+AAQ7OFbKBrx^J`Ac}64TQYGsgdp}@(k4WI9UYDo1$kbZ}5((jOP3)-s_(MB#WZ{ z;M7s#P3kEDk|H;-9bUl{F@#qUK8x@V6O!VHxmi& zE|M3JDJtwj2WaCos87Tm3P+)b5#pbY+P$Dwb0j=Vyr5c7Vb{uZ4e}FI!~XzSsP_7N zBjLz}!1OWirK}hiwtMI1f|u+^!+}TydMStsx>F$HUog{%Q7<_^G7qId_m+h``@`whYHJLBHPzRznV#bCYh=z8wg?ctu3HD4mHw}^NBnO@ec(bC(bfkUH?^~$I@J8|Zqk-uN>TScLXcDQZnh6Jk;iu3S z2+cu-KIsk4`y>5?h2~tS%a`sxA91{~93X>-geuhcjf^8@AHc~&vboB_007>m7zn2gP z#{2@ojWPd)R065*qXvEl;V62KA^aF)aRNggM-VaWBqT`|_yUCd2)%!cIo6Q+8N%BL zPy>ei6!ZQcde5QvJi?3U{Z|C4?>thd^yE#XeuLCMBjrG`d=DuLQvZTf4O0Jrlpm?j zAoUgkRnvb&Y5^nu9;s<8?Vpeu!LS1elL+6(FgH>^Kq`RLPf=pcSlypMn7f+b(1t!X zk|h!PObG@e(WD-lbc#~inV@$zSr`gGHRX$DuIGGY7c%@2^6~(c)SN$`H29yM!>7m; z+Ro_IIbTvU?}KKg7Cc$BZnSd@z6M3`HuCs6LKjx58z52*Ra$@4Z@3UAMA6F82mV6U z;gJLDPyaCStW*A9!4c_5FPyKx!lX@1^*GB{je%P%f!8FO_QKD0ufw7qURp|DBYKuE zZ31_fGx8OU0@oyRO*{wQ@}-Aj592F%77wNL@N3|;34JZR_CnO#c=*xR%^@w_i&{}% zvuqXgjhSD>jn9fzo<;3aNvhnis9mczFZYVpw$($Tb@y$fVC|POKc9J0Mg1*1bWIu;A>8f%uT zU#-7dzgW0ttd3n=b}W0*n`XG;Y6^=(DPxsrY>2mt#tzA_6PQ-?-R;|RyWwjyZ_J=? zKpGgq>nOW8yk@AlTC&{p>b9%f7KhfVt*?#5*?0r}J{))5D7rSWG9gsE#cFq&DJZj0 zh-sus^K!%TaGYHpTWNW9Fz!uMHuIIulC?HADw)i&v9yjcwSd?Srr5=oFD+iWYpCL@ zJ%XWAG<5QYPRURe(?J@8ISp@k-4@ecEfq?wuT4X_L-czQUiI31Ew~a)l(q_`t*Lr@ zyaM0;dfj%*$T#$e^*ea3Hf5_{99=phIh>25B3Cap(5r2&1#a7PwLGpB?LD`Kc>7M# z&=i0|*G$z{4}5;|vy<^&(YfQ+ zAOA>K4}*KA!e-!Z!rX18E&6U}Y0VX4j**TlQ8Yp*>XKXdKWYPC?;85_?s z>52QqroLOxh)pAcb3`zWil$NCG@6f-tseT?&~?+bvvG3m{A!b6>x%8qa@rA(ip{%j zmx|4!f^$?bjfti)-ZYknQ_9{D?}>Y__g>epR$k}+^2E(21^WZB1CrhO`uL4W!QK-) z5IX?2+g2YNmjL#I#b|ITaU0U<8H+VmC;=2fx;p#og*UbdwVqh9R9mmmAlA0-Ocy^X zaLuWaD%2XqoT=7$)sT;bpXt2kvqM-vi4cPcG z;L_fGb5Axk1e6(<9(v{Q%ZGW(Ny$(hi@sXAJSsr8=2h4A_SLB1*m3L8H%_Jt*eBS1 zXy2_A%ji$l1H!WkkpqobdN$2elyw77YyBOICt>jjmM+oK6)Q?r*1leFU3WvrS2lu^ zOf^>@`uyR~9_Aen38n+0=>Ts!@S7T@3iP#8Fm;KhF5c9Y;rV=;DJ<)Z_uhEmR#U>c zn|JO`nd|PDdlTkf!Q3aB`(lNuO6zNpYx67fe5DIq4)?u+X+Sg$@TP&gCNpp85KJD? zm2`w#BODA5tr0$(>)V^hc zRw+AzjRFW(8W5Y~2)hQo^$&wdrf%tKW?RM5spauRnL{XZq)b(dW2x4TRiEg2@b)vJ z=aA5PD5I39_>$ZUQ1Mf3?OAnAxmz=y$6*8{<$6-BZMh6oc{)V*p4)EGJqj&@9imn= zctDD+ompNfH%-(ky2tLF2pY#&d9^S$|7zLQG7!DnCfOXK&9i!0wCx3k?tu)$%=5YR zW+4m}FnBvYySGc*t>X6PRZiTF!y5N`XU4Q(m)&jmoqDb1{y9Tel;uJtCVQxqRi0ca zOYU6BaE8SRP^PiCOLX+ydO&mx3blg*H;ADyH2i>(FHNLC(Kl*D*DaISJ@nl+v3pWz zoy={}b-VbXXoQq$?2k19`jAQqRGxck#G18kWw(?Wr#Fkn*40rF1`v;l?PI)Q{D*$% z7`2b0O#%eP2ZVv{#7CFur3ujbeAXtbPN|O#%ogV8{TLW9UA4A9WNk?0+A`G><33fN z=UV8yi};Emfg2XNVV)aS&AM?GXjmFdGn(Rhu!pZqzC0;#wJ=YESsu*RV{HOwgGmU_ z*;1wDnRJbi?ha>9aCU)nh@6Ax9CTjq6gZd2xp>Y6Dy$5+VazR+S1!RosGzJXWiZ__ zI1&bjU~q{B7jJNFn4V^G!V*G3aW`grX89an-U1GuYj+IpguxA?8`03l8`@~jPQlP6 z8oGExS0?9!cq5&(f(;3-LExMs=j1sjtzd`1c|^{`a~?V)5KHSYwL9}@*dsayzFUxR zjPQ<;d*pA12fE0$fJ`$pih(qvEACCfaM9Y4u)=uYX7G+@FyR>#JVT;qNU#oz)?qXs z(ub0S-b0VT?GQPbgYQ_shg!UcW?ExMf^#TsKvQRzT3~bzO#L__}G6+wC%dx zo3IV?w!!a?ilfK*M^B5RPw-xUV${!%`ajzHY8SW;k?Y{Oj?Ky=mOHcBvsnkSiyOHj zaojRc)fQoiA=+In}{g>ietC#@@FkJtxA5M(s>GsTkcak&<4~exO#zW5V;1PYoNE&1Os4LY2&%(3`^)DY_Sa__sWr% zj|g0y$kp*&9bM`(C-mlJx1e{4dKa%pn^9c0q<^LCpl|(oLWk+xW?TYeE zYsz%`G=pHX>vp6W1rD;Mx}9l80f(2>m8Y2;oG93aAquz@iiSCi7^^`05t=^-eVS@a zp+K)eSCwWIXg|#AU~M)BeVPHCEzKw}stK0pwx^jKJX1Bq>Kf8a4j$A@7h|d%JW*-Z znbV9K`i+=MfffT+PJyO!U0<3}ptWCj0;bbBc$pn$k-;WBs3}KI3b@L2yHVon;4opV z0$u&OamrZ16PjjScbZY)1Y4=wlV%j~z`X*896YG0*Hxw&1;*JXOrgLr)`@-vj}8Gm|2Whq5^*E=t8!un-ZTXLLjG(#UjuLG(b4}EgKPG`M&elV zrq5;FWjL36x!_CC0lx}=nGc@OK`zU$gC+%pQ=pH6+{=p#+y4yJn$um}i znl`%8iH5~pl*#Z{kT!D+gKimO{-lm>rzQ(v4=Yjx`PVnTG8=Lw2 zseM;AU_j8i70%le+U?6aL2DDWHePE>X-lq*Q{9B8#5bV6aKR%Z?t;w85j;|I+&xDA z7gQy)MT~FaNEUe`Q4*fJL{JCF3I>&|2YCWvy5X8ElLMk&*wQDc=>)}=)WIGyWW~+m zWTl*FIvggj>kp?sB1DfdoE(E2mZTVS(|ud8BZ2ILT#B49&qA5kDbsVB4KV z*oG;}5exuHE;1dMy99fXvv8h+7%}{N7+#uN_;)bGm`NE3!TE%_FdkeY6&PN&&H*$F z4i%9Q7G4yYpPLJlC_xn@jaV~$Zy3QXQi5+hA^3WgWZ{MBz*HnmumcenLNfwE_#DF5 z5xxnKEDWDH2PY!Pw=hJ5H&Ts2H`}d9(M|L^r0Ayk5K@x}GYCP1=Mb(T&}Sp+a%P{s&1ZG!7(XA`{S_uzta^Fb-z9E2a3An35~>D>21a=vQKFyxK{Ol~+57 zspr*BVybwxlb8x#b`sOftDVHu@oFbAEqp#FF_paRB&O^N{Ys3Jmz~6zd9~ZjDqnFf zo>93R=6OoV6elrFyxMIRMY^PbHK$A3S#A1Bwv2U3S{JW&lGeeiousYh)lSk{d9{0(ukM2;(Ce zdK@SX1x~Oowjs?Zu%M~LlnOk_8rcdM)ydGM0g-O2M6Vip_Oo~goA9(|WRhiZ__Fo& lo%Q%%Z%lXew8p@e!boxxM9^CQ-Gi=SC-bJWX{5B^{{cY%R=NNH literal 0 HcmV?d00001 diff --git a/csst_mci_sim/csst_mci_sim.py b/csst_mci_sim/csst_mci_sim.py index cef45af..c62edf0 100644 --- a/csst_mci_sim/csst_mci_sim.py +++ b/csst_mci_sim/csst_mci_sim.py @@ -821,16 +821,16 @@ def psfMaker_IDW(px, py, PSFMat, cen_col, cen_row, dn=5, IDWindex=3, OnlyNeighbo dist = np.sqrt((ref_col - cen_col[idx]) ** 2 + (ref_row - cen_row[idx])**2) - if IDWindex == 1: - psfWeight[ipsf] = dist - if IDWindex == 2: - psfWeight[ipsf] = dist**2 - if IDWindex == 3: - psfWeight[ipsf] = dist**3 - if IDWindex == 4: - psfWeight[ipsf] = dist**4 - if IDWindex == 5: - psfWeight[ipsf] = dist**5 + #if IDWindex == 1: + psfWeight[ipsf] = dist**IDWindex + # if IDWindex == 2: + # psfWeight[ipsf] = dist**2 + # if IDWindex == 3: + # psfWeight[ipsf] = dist**3 + # if IDWindex == 4: + # psfWeight[ipsf] = dist**4 + # if IDWindex == 5: + # psfWeight[ipsf] = dist**5 psfWeight[ipsf] = max(psfWeight[ipsf], minimum_psf_weight) psfWeight[ipsf] = 1./psfWeight[ipsf] @@ -1051,7 +1051,7 @@ class MCIsimulator(): (self.information['ysize'], self.information['xsize']), dtype=float) return -############################################################################### +##########################################IDWindex##################################### ############################################################################## def _loadGhostModel(self): """ @@ -1822,11 +1822,11 @@ class MCIsimulator(): # ###### use SED_code to generate star SED ####### if get_file_extension(starcat) == '.fits': - umag = self.star['gmag'][j] - gmag = self.star['gmag'][j] - rmag = self.star['rmag'][j] - imag = self.star['imag'][j] - zmag = self.star['zmag'][j] + umag = self.star['gmag'][j]; gmag = self.star['gmag'][j]; rmag = self.star['rmag'][j]; imag = self.star['imag'][j]; zmag = self.star['zmag'][j]; + + + + # SED of j-th star else: @@ -1926,7 +1926,7 @@ class MCIsimulator(): psf[ch] = ndimage.rotate( temp, theta.deg, order=1, reshape=True) else: - psf[ch] = temp + psf[ch] = tempIDWindex conv = psf[ch] conv = conv/conv.sum() @@ -2205,7 +2205,7 @@ class MCIsimulator(): # 'earthshine': earthshine_ifs}) self.zodiacal_wave = wave_mci # in A - self.zodiacal_flux = zodi_mci + self.zodiacal_flux = zodi_mciIDWindex self.earthshine_wave = wave_mci # A self.earthshine_flux = earthshine_mci @@ -3312,42 +3312,42 @@ class MCIsimulator(): return ############################################################################## - def img_fits_save(self, data, filename): - #### - - ofd_g = fits.PrimaryHDU() - # World coordinate system and related parameters ##### - hdu_g = fits.ImageHDU(data) - - hdu_g.header['WCSAXES'] = (np.int16(2), 'number of WCS axes') - hdu_g.header['CRPIX1'] = ( - round(float(self.information['CRPIX1']), 1), 'x-coordinate of reference pixel') - hdu_g.header['CRPIX2'] = ( - round(float(self.information['CRPIX2']), 1), 'y-coordinate of reference pixel') - - hdu_g.header['CRVAL1'] = ( - float(self.information['CRVAL1']), 'first axis value at reference pixel') - hdu_g.header['CRVAL2'] = ( - float(self.information['CRVAL2']), 'second axis value at reference pixel') - hdu_g.header['CTYPE1'] = ( - 'RA---TAN', 'the coordinate type for the first axis') - hdu_g.header['CTYPE2'] = ( - 'DEC--TAN', 'the coordinate type for the second axis') - hdu_g.header['CD1_1'] = ( - float(self.information['CD1_1']), 'partial of first axis coordinate w.r.t. x') - hdu_g.header['CD1_2'] = ( - float(self.information['CD1_2']), 'partial of first axis coordinate w.r.t. y') - hdu_g.header['CD2_1'] = ( - float(self.information['CD2_1']), 'partial of second axis coordinate w.r.t. x') - hdu_g.header['CD2_2'] = ( - float(self.information['CD2_2']), 'partial of second axis coordinate w.r.t. y') - - hdulist_g = fits.HDUList([ofd_g, hdu_g]) - - file_g = self.result_path+'/ori_Sky/'+filename + '.fits' - - hdulist_g.writeto(file_g, overwrite=True) - return + # def img_fits_save(self, data, filename): + # #### + + # ofd_g = fits.PrimaryHDU() + # # World coordinate system and related parameters ##### + # hdu_g = fits.ImageHDU(data) + + # hdu_g.header['WCSAXES'] = (np.int16(2), 'number of WCS axes') + # hdu_g.header['CRPIX1'] = ( + # round(float(self.information['CRPIX1']), 1), 'x-coordinate of reference pixel') + # hdu_g.header['CRPIX2'] = ( + # round(float(self.information['CRPIX2']), 1), 'y-coordinate of reference pixel') + + # hdu_g.header['CRVAL1'] = ( + # float(self.information['CRVAL1']), 'first axis value at reference pixel') + # hdu_g.header['CRVAL2'] = ( + # float(self.information['CRVAL2']), 'second axis value at reference pixel') + # hdu_g.header['CTYPE1'] = ( + # 'RA---TAN', 'the coordinate type for the first axis') + # hdu_g.header['CTYPE2'] = ( + # 'DEC--TAN', 'the coordinate type for the second axis') + # hdu_g.header['CD1_1'] = ( + # float(self.information['CD1_1']), 'partial of first axis coordinate w.r.t. x') + # hdu_g.header['CD1_2'] = ( + # float(self.information['CD1_2']), 'partial of first axis coordinate w.r.t. y') + # hdu_g.header['CD2_1'] = ( + # float(self.information['CD2_1']), 'partial of second axis coordinate w.r.t. x') + # hdu_g.header['CD2_2'] = ( + # float(self.information['CD2_2']), 'partial of second axis coordinate w.r.t. y') + + # hdulist_g = fits.HDUList([ofd_g, hdu_g]) + + # file_g = self.result_path+'/ori_Sky/'+filename + '.fits' + + # hdulist_g.writeto(file_g, overwrite=True) + # return ########### ######## ############################################################################## @@ -3370,46 +3370,46 @@ class MCIsimulator(): else: if direction == 'horizon': - + self.information['bleding_direction']='Not_horizon' # loop over each column, as bleeding is modelled column-wise - for i, column in enumerate(data): # select one solumnn - if column.max() <= self.information['fullwellcapacity']: - continue - sum = 0. - for j, value in enumerate(column): - # first round - from bottom to top (need to half the bleeding) - overload = value - self.information['fullwellcapacity'] - if overload > 0.: - overload /= 2. - # self.image[j, i] -= overload - data[i, j] -= overload - sum += overload - - elif sum > 0.: - if -overload > sum: - overload = -sum - # self.image[j, i] -= overload - data[i, j] -= overload - sum += overload - ###################################################### - for i, column in enumerate(data): - if column.max() <= self.information['fullwellcapacity']: - continue - sum = 0. - for j, value in enumerate(column[::-1]): - # second round - from top to bottom (bleeding was half'd already, so now full) - overload = value - self.information['fullwellcapacity'] - if overload > 0.: - # self.image[-j-1, i] -= overload - data[i, -j-1] -= overload - sum += overload - - elif sum > 0.: - if -overload > sum: - overload = -sum - # self.image[-j-1, i] -= overload - data[i, -j-1,] -= overload - sum += overload + # for i, column in enumerate(data): # select one solumnn + # if column.max() <= self.information['fullwellcapacity']: + # continue + # sum = 0. + # for j, value in enumerate(column): + # # first round - from bottom to top (need to half the bleeding) + # overload = value - self.information['fullwellcapacity'] + # if overload > 0.: + # overload /= 2. + # # self.image[j, i] -= overload + # data[i, j] -= overload + # sum += overload + + # elif sum > 0.: + # if -overload > sum: + # overload = -sum + # # self.image[j, i] -= overload + # data[i, j] -= overload + # sum += overload + # ###################################################### + # for i, column in enumerate(data): + # if column.max() <= self.information['fullwellcapacity']: + # continue + # sum = 0. + # for j, value in enumerate(column[::-1]): + # # second round - from top to bottom (bleeding was half'd already, so now full) + # overload = value - self.information['fullwellcapacity'] + # if overload > 0.: + # # self.image[-j-1, i] -= overload + # data[i, -j-1] -= overload + # sum += overload + + # elif sum > 0.: + # if -overload > sum: + # overload = -sum + # # self.image[-j-1, i] -= overload + # data[i, -j-1,] -= overload + # sum += overload else: diff --git a/csst_mci_sim/mci_so/__pycache__/__init__.cpython-311.pyc b/csst_mci_sim/mci_so/__pycache__/__init__.cpython-311.pyc new file mode 100644 index 0000000000000000000000000000000000000000..d0cfa51aed7d5e0c7e2f8fe2ac372f8698db5cfe GIT binary patch literal 187 zcmZ3^%ge<81fI!@(?IlN5CH>>P{wCAAY(d13PUi1CZpd?AD^CCqMuw`ToRv~oEZ<4MCE`5^Y!E7GxIV_;^XxS mDt~d<y?;IiQ#yNPJ*sWMurn03(W+fnosCFfOhD literal 0 HcmV?d00001 diff --git a/csst_mci_sim/support/__pycache__/MCIinstrumentModel.cpython-311.pyc b/csst_mci_sim/support/__pycache__/MCIinstrumentModel.cpython-311.pyc new file mode 100644 index 0000000000000000000000000000000000000000..3edfb62546b041fc24ff35024dec605978f5cc24 GIT binary patch literal 1761 zcma)6&1>976d&#St8)ASZBp7a^c{-vg7?FAVmBg=Yt|vy*dd9XmgFMFnpv%|R#Ik0 zo2*HTDTTJsLk@*pdhjW2$t9PbN}>Ni!G}Uz_0WS)zHM*=A*Ika+O@rQ1N^YkqGM^;}R%T;=oZ}7=Wr$1%j(s1QN7S!F{SY123qA zqDE6C89P;=5+d&|dXe4+9|K0kIykWq=za`b>IGx4!9^H_o)*x6fn0CE8g5qn=KvD@)@WpL?N^aKAnE`n+&#qg z2v>_fCiloz$KL}vK^R$vWmy;N#;1q%vOe}2Gg&`7!3Mq=c!*~n9iB7L3c@`||H=?z z0@C$1P!yP#;zvaT5K%{XYhK{Oa(MyD<@o}nn1{X}-qGStEIov=RR;#69(yiQaEI406()#6j zkP&y_(v9*`{?f{7{`yLODqkv1z$n(XWv}C~V4>s|OzPt;Drhr|3$Kn1F-q6MI>(0L zaCK=J!%a)6yRXBBOe2_1ZKS7P8uS=Q&M^lWrY2t=>```7S#XZ$OH;=KL)Kjj#(94H znBbxHQUURtID0(%O3+Bpl?4u~**tc9j};h>(_yLF(2eQZdP+}~3Uyf~N@sC?6zu>` z$2x_vjvt+VRiW#uin)zXy{`-PYM`-~QV4PQ^psoJjfdApOJFnR9>q>B(<@-1~hdH;`m;4wCM)N0Z?Lx2x?E zB!ew~hx?v2uBa`O);mhF*q2EbQA-S!Wa@5`kt*rZ-I`>vPjZ>0OQ)Ok#1V$fb4HR4 zVFWJ_@zcfox3|viyf?Fbt-AB} zneSeIeBAV-lkvxZ?bGyO(alV~<_#u#dn|FYE5 EKWQ%W1ONa4 literal 0 HcmV?d00001 diff --git a/csst_mci_sim/support/__pycache__/__init__.cpython-311.pyc b/csst_mci_sim/support/__pycache__/__init__.cpython-311.pyc new file mode 100644 index 0000000000000000000000000000000000000000..15e38979bd24b4bb93b9ac94441d99eae69c11fb GIT binary patch literal 188 zcmZ3^%ge<81fI!@(?IlN5CH>>P{wCAAY(d13PUi1CZpd?AD^CCqMuw`ToRv~oEZ<4MCBBh78K+cmFUOEXXa&= q#K-FuRQ}?y$<0qG%}KQ@Vg;HHazQaakodsN$jJDC0Y(%t1H}Ltxi3-x literal 0 HcmV?d00001 diff --git a/csst_mci_sim/support/__pycache__/cosmicrays.cpython-311.pyc b/csst_mci_sim/support/__pycache__/cosmicrays.cpython-311.pyc new file mode 100644 index 0000000000000000000000000000000000000000..43a7890a7f3f804c09b604251ebce60a62bcc77f GIT binary patch literal 12511 zcmeG?TWlNGl{3SKD2k*gk$PE@HS$w*ZOO0rL0Gb6*|HtL_BwXjM4N;==8U9~`SQ%r z7KPy^1yX=EgbkHtVLDI{rqHGe>;h`x0{-a!EEe7U)EERJ4B&;ZyTxw)WTXbV{n~Ty zki!{DbggYS?MFu9o%=ZV+;h)8_uPBVoS(YgP6D2P8|$15>?Vl+K!x&SEoHtr0h#** zON#!}aafSox4RnF1480EwC z8%#=c4pttg^ZFDo(jp&CL^wJeVMLJ*Gco!GM<+#&rKLE{$HI{$3ym5Efl1NaZ7wE> zs1ZQaqBqcr`6QQdZh7uo} zbfVr7E;cDm0d|(=Vw^CUqT>@7tu}H$eJvTDszOmBdlGeEB#w1X#08q0Verh7BA}+P zNsPcIMG58t6Y86iq(pFFfaf@FCJ~7XTtCY8$A!rOQG$M2+CMEOPBTI*8NI;?g9ncu zJ96ygq2qym;2c*)-9{E6qL_#nM+_A?U`-sxjsfwtJiGq6;#Nwie@YB(`ZXqap=ftLi@ZgT>{Iw3G&#NU3WGbnJklDxo)LHcqs`d2Aq z#kDXG70w{^5&=3$AM8KubbgJ*=yHHmyDr8g4#a20Zq=#H3Pem~{*wf$ zrQsRXXG)jW7HHI1)K8BBe>omRM!>RFmds5g4493u;YHZ6a0Y?a*9Vw}Is_Apz(nat zd@_iVNf?R83mE2%?RXaiwatK5;m#Bf5?Wa$kw)^*u&lPJ`WR?Sq-fZCQ3hl&7T08B ziiu5fSkmJHKgq|K2n~VFB@mXIW;CcTfW#z*k0E)|KsIhp1Ihr}w3tl9yy^$`SzJs4BP2cwIs-Ld z4%-bPWPv>Tee=&)f@F(Kg!E>agc68__Q+JBEN|s1{Qd&chk6b4K=|+@6VrNrx~yG^ zcK{1O*(y*Wog&%#!+O~Yt###Amb?bg1xLATe49rmL0g0Vr5G#HzhzGRX{8xq-LgcZ zPJc+M4KUxLWE)1^vUD5z8pQEimvBVJCQ^>7F9+(?I%q?(uUbHts8*x{R2$~7Y7?iJ z1gF|kB0meS8T_Kcg4z)hIEK{(2R7=(%0>^^RjbHFCIlpo0=h{ENKjR~#yz#czzGov zH6nwie}8pgDjwwqQcP^%(&)tjIQ^5L+~dMPKPrbNd1)Xlic%;V=0nhEE{MrQA}#;~ zRvk!J_roL{p%CZ@DHNLZ85+u!J7H`QZR8J#^>(7YCu5y+t@&EvMiD^#P*MN7_fOoQ1xlXP> z1-&8Xkew)J!l5fZP`~vh#+Y7aB`DX;!U~v2OE<_3thGd;i?zuOvI~tPHx$rnZFr>C zN>w!`*8WFQ2dK?BXQ3wTmfc1RjJ#Yqm6Q-ahDl5TG@GxC8`^@Om(BwQIQduG_UAcB5Qq@LC1XBn{Hxl z>E=ly?Sr2mek~w<-b!a!RYOUWQQ;?2rvdv*;LTd$Eb;D#=~lT_ri?hz>cJi}jDod6 zn@J`rJG;$17q;#u)+NICxMZs%f5oALSH2(SM?B;-QJl<`wvkhqZVpdyKq>{bb<$rjq&+_)WOnF1N#V+kg{G zs(F1IK;Hhml4{yZAb;kaR73u1N^|=klvER!%`5O750S2sH{l`DO(lu^8>y|LkxV+S zv0QJ~kLDFIm(m@tfXyVs=KkhVdYf#Qx5*t7)C~L?&VVMQSHqStwc+;T)oJ!s;qsa0 zR)OzwN$~inBbafj@v2wk>jBIt$&r(?x+2q)LS zq`iRI3ug=rUBQp)DkY>L(p|*4V!ZE~qo$f^oE)v6&Ma&y#fbD~;H)$%s(* ze;RQ++)}&4ASI>Kx50Q5EJ^x-7)26up919eDoFi?Q-W;AkiP)w+b(x+kiIUt30Vn7 z8JzpqtaUCd+rgdSV15l30NtH%794%Ijov^Sys8|7{5asV3v&sHPQ)QEmb_WW{43l3 z$tZZVl2JYu7kP>2#Ilc7v)^dlCnV~lm3qz`1hX1Ezyf%uz@G!1kqfuvq216MNTFvA zu;6>dzQ}51sVyT*)kYS98DZCPZE~u|PfYPaaGvroAI*6N9=OQMm{CC>a8gnLcd22@ zQ6@q6VS8TlxnX__QYyg(=@@$ed+~Gx;9Iq^sBN;ey@w}*a|i$Z?^5T6{;Ry<1tdIQ zqdQ*-Xzg^MS*2nL)i%osaZv!XmJmQ+5o+B;B+f|3j;a4xla5s5lKdl2@I& zAf{9@rjkjOoKeY?YM%)MA*@QpZ>siGDL)g2al@(=5v7uMR4NR^g;Qb=kcJrkl1K0% z^8F!g$=*QTl+r#CXM40t?9tc^U>LfJRxr~)AYOD4zLs2H{_Mi*OV&rumChsoc=%Vr zmF6*}d92txmZ<};hR6Hk9l1Su-~9J--&5S(1xj)6$WW{9#?7jqd7E?gtT6A)IrG<_ zdV4drjWw(O)(n+-?QXrEd0{6So93pn!u>lR-tkX+K(QKf2O>^mprSAmN3V z=-Ks`N!}JX89fE}IT?5fC;!~!0VkWXhwdda$!9I?pEfVF7j_h` zF7AHNv(Tfo1d1(zELF2%)!UY*@@Mk)Y-6@D*H{QD-UBE*-w6A&$p>2H(7Ze6&Ws?n zvU%Y*pwWVPW_~a?m^r%!r0mMop*>Tb0feoc`9pc8(EEU!59LCzs_qjw$&#!7mMoL+ zo#%2F^0x|23)1|RED2QSn?8MG;jO~0rKY8;U$!jHEcY&tJl^&B`%iW#drm3cgT?Ma zrDITO87j66;c0HI!;r}U4^Vb)~wVYN9>C(t@^UD5V_ zUG;b7&lITqSoTWxO76-Mx%B#%9>sqg73Z(~?ztOH)Ju5VG9IA?PJ_}Z?HK<0zrn{m z7S_KhK^r=5vE%vIaGBe> zsa->}l<;Gdd%LcT_0_&^)#$saQ&_fE(PEbn|4I+cL7x+Atz6Je^|ngA*;q-}QOKcY zH;-WIYscVsg-Q5_ZRrk?ec`Yn%&o#(M`lOcz0Yxt3=$tpLcX=DA`8oKyBESZev~GQnpP{v!EIAM;qOJ zQfMyIV|>dl^H?Bt1}(Y#_iWuY0{sGRh6Hm=x8Kf}LXC_DU@4NMW+)DhX)1?Q@i^!M zg@Wi&7}$VB-*o$RYvdA>fM_Fq)RGZi4)Pf*#rPmIv?}ZendD{EDW)771Lk-Ff>@9N zk8Sq2qy?EkbQ3Bk_!%x@F!_FtC6=O&Od3E~D58 zJ5VFiE2}cS-3&_ z23EcPTwUH#^!65xE*@Vxws`uvh4dZ$EdfbpbZvWYK3v%KfM4KI%tLnU?v=Ghe@tC!U=B`OxDbW#CM4;EY0# z7U@yNH<~%Sak76*-)Ajtxi`T|c>A;WK6&qFqrapcJD&J|?s@D{_Ma~9Kdp2R7rTd* zmf_48Snuwp6<2rB)m`v^)f;@;8`QXI2g}D_ke04L`oZ$|ABCPw0B3fLk-!>ApA)2e zjMTcla1d=B#nwFqskrxq(t5JkdNOlz7EPu}v3+6;D9_$JT<$JU|B{HV>vD zb03A2!F%HaA#apyfp`;ueTd-kgRVtzaO+2?=@>-cE+|8ooFUEyhoh|<{l_CLyFo3B zr9e_mWtYyB@+SZ8AcE**?U?%(4x+ODgXmq&)r?)s^c{$JH7uLqwANkCRi0)_cpD@~ z7>6Wa5m=OmZ#RT+6;)WVLKogb6*`!R3{Et@rgj@4boI@|^>|qhX15tsZB^UjX-6!C|vPG;9B>L~zCYziz*6A%Wb@ O)91MUw>4TZG5-q}Raqeb literal 0 HcmV?d00001 diff --git a/csst_mci_sim/support/__pycache__/logger.cpython-311.pyc b/csst_mci_sim/support/__pycache__/logger.cpython-311.pyc new file mode 100644 index 0000000000000000000000000000000000000000..a1896255c6b9a0b7d3bc2f94f8b624fb24ec08e0 GIT binary patch literal 1491 zcmZ`(&1)M+6rWj*6N9CBn%%2)`o-?IH&|S`+W$c-3IdI^+kun{l-JsFXW+WPICTg3)@% zw?i-R1+^`oK4i2bm_zFUr_I21J>R9C4_eDoL1n`zm+5yF_rX!E#vzPTPf$MyX{*x= zy*3YQCWJ~9R-Hjcb)Rto$yNGU`I_;sn)Hs$sp^c|>EOzmNq#3itTqt>Pib!=ckwfj zpHK8BT30)qu9Bvpa3UeJr)tMEn~i6qNE&qr2Yi~kXg2#vW@g4*n&G@VBk5)i=Th(V z+9k6UF@oLb(-Eos420*Yg82Eaw)b}0{KPRO%x1)uV*pKs#9}vddIo@qcJPNmf$p$w(r11iLc4_SyOdO zV*#g@3tDMj&Z+VMTA6oQ_)UyWsS6NQj?CV?_2q*tsW&X&X#$$%hxdXIUWD&%!Y02S zTV;MTekUR0asg7`YV!?CBtb%E?IyFhG%UwSY^AJRfcBGPl*V|kYE-U(E9a8NGCzKb z&_XDLs*9j+wkwUG#VU_2zq0++?TYZ=AzER;D`lnETraHHLWH%J?bX0JtB6j!9q_O+ zZdq7l-=AfxZ^$Ykhs~B{7c3 z(vU3u;tohLBE>!_zMNm|ZJ)r*sD$G)ED} Ds;_Hb literal 0 HcmV?d00001 diff --git a/csst_mci_sim/support/__pycache__/sed.cpython-311.pyc b/csst_mci_sim/support/__pycache__/sed.cpython-311.pyc new file mode 100644 index 0000000000000000000000000000000000000000..00bef1b8dc42091d553c5d64345fe59be5d1c101 GIT binary patch literal 9037 zcmcgSTWlLwb~EG*-x4396kF1+M}8#YTGk`6){fy_S*wRFCwA&s-elDz+!1F)jZBg1 z%+R)oVav^?EwvC(x=;{hp`gsxO=Uy@)FLgi{pk8LyI&fEK-d9{rY(ZNC|W4ZF0fd% zK+m~D4mq-uI3Mk`c;-INeVu#H^ZKt_T6`3w%l~;ZsUD!Hf5SpOF}2L+KY`2}lt2k| zikhL}>qt2koik1v>zpZOk)2^lo`Jk;#zpdM%Dw2B@sfN4KPsYqSQN%0X?;Pf>Y?3bpmLVC7Zp zRZAN64%EMEwb@H~r)>nBi6ofCd6kR;Ip!+qZQp2~5m><`xHX5y?5}sD$jYE3c<@->m0UFOF4R8CY?i4#?hPxW)XPYwq|C6F6x(H~h1t-HzZ|)b7c=<}Qvt`( zKmC)D`*>NT-78Qb`Nw}flMMR=fKh5~F{B^k< z=4gYT@-F}i)NTjW)?O0Ur*B@>BL@r5-0Rr4Hvg7VKK|zN#Hs_f9T(texY*i#hsU4*3Z4HkGujn-Eq3gCom(>h386RI4i=9AH*At3W zB71%hGH(Dk3iMSXo8$w{9H~(nv*p+$9axPPoC2c}S*5Xr=THKm)N*Mq;4_y0TM#Sc zVm@C#bdFDRQhF(?a>|mJP-UJ|=jBW`InVK2Qo1IlIqQCKCVzY#-s3WJ_R%e7oJ(bI zG`MvVOJ9g9&VCi=62`N0?qI()SBk&R^(w0THTeFC`tdx`80O*fBQA(aLY9a?iQAX( z2*RVj@jCJOIZ#>+lADvuAT6jEEr;M!yFrJ~B!BuQrhgiN2IQoeTMohpbRE#ra01Oh z3ILa3CxbOv+l_4(#p{Y2DhYdYoE$|2uXV~Y2jy?UPx&qYV9RdG<1a4Ub?@78@4FZ3 zUY{<{?u24`C{_u@iflE~v3{t0Xw&(2Y*W+^KBada*CWG4chw&$T$DRt2*Z7LbZRs{ zyOhA&T~B!H$xZ)nxRP(!^$pQ|h^RS)npQ&H1X>v2UTq3I;%sQcJOkw>?U-k6nwinq zJzV2bkruVIV6AYk4$r$aN(D};CylYM<0|P1I4yMySp8lc(v*Yw@*aWJ+?oe9?N^LT z`h3N>WbPbLE$?+unhT_Zw?Ps6I$>5X(eKT>w(=USxrvUZ>{PPRsDTYMvEvb0;}h&O z1sP-}&1YW=&ID4)o%`)rCVfqmRSvHwiUCvb)OiBb@FJfCgDoqFUq|>@{7;6pEQ86~ zB(>o!tXWW#4~mINDTN$K8q@sZ9`oRd+SLXrLB{za8^JUTU}ueuAp7MGpb918cbr%G zfy)!)@v)2`qP`TcBFFT@*I*Kv)DVGB%ETft3Ci4wd1YiIgONpEkddEcyrhT<>eji) zgCtm#(gq8HeJN+KU=a)r#{5K5$Xx_;iC5PP#CWULZyz?XA2Wn2GnpMTAElN*8M#de@|ij zUNF3Ry`-+>SM!B))o}Y;&h?2meH*^J;e%!My~OsB+Y4I@cb_<^hn_A>RXe*1=ZZhj zna4n8z16mU<;~8G&byIAo6LL4_VDf8R_<=^(|Y(+;e560!20nrQ=X}CgL>CcrE3T} z&X{&UOK|mEX?p$K&3IYR14k=?qwmdb&2Qhh6V>}CEB%vtV6rgru|N3w!L{D?;hO_% z1F+Z~riX|^o33{O8HfZ7JRkvqgK4zHx9E3(x8EiF4cJtk6&!-Iq1!>PJJcrMWx>iL zZL0&w{Ms1E5@b<>&Km6HY1ocfUy}#{F_V;!LRT^cI^O6GK!1`$+M*QUQzYj@7>5&;GDa)OTbogEadk4*MTpmG zmd`@-58 z`6$?PH`w!6J-3d$M{j|(^4{py_;%0jvtYCQ-RS${`oWWW@aYQvjTFx9x~NF!`Z+y( zxN!dC_WkSU_4da%qk8+X!o`pM`*)cACOgAAtW;;c&pD`igAz$JWCjeF8K>ZY_yhy4 zI|FeHR$x)-8J@G?nN5p}OF5355N5%#5>fv_6Mtiw19d=}63vW^EK9^jnurZ>w5hnY z2(<~bAq+HX!NM>NDm!aJNWx-VFi4#CR>5(V^p+26blDOZ)_{-^_Vl1_rqbTlS3$Hj zn`jFSp~5!I7TA@wCct+!S{nF(V+j?&aX$=4@HA=xA8$NKUx?piIsrajNs{uxH%~$_Hj!okyk`6F)t>r_Vo1J zhXbZeKEq|pxE8|&Dtt648qP`a94KWR}FVjOM4ol8Fr8_c0tXfs+~xniEYt z>*14PJSj_p;h7WXNe&_ij3TDyz-`Ew{A{@Ag=|VvR0U^1&NS7!Cl1PAQsePlxPFQI z>s3DjOy%bQz;p0X?fZ+&s`p-`V=V?+e06H4d*N(rd8VLT`)K8|Wp-Dh(IftkKcM1F?(rwnak$Qk;J+N$h_xkAC`L*ZDla=UE zoVXh8GHXa--+?l{;oAsqdMf*dpr)m2)v%R){a90ke}uJ(nPm$LTIw1dad6FTJY@#8 zP?14^u!#qeJ8gVuwF!FUPa;w;P3Itx~RirP?=l5!5OA<3y(tQ=9Or!DLZ0GD!ku0wn@- z0Ao%Wf_o4gs;`XX z+>0InYgV`=ZAahVe`nuL^t>KDUx}VauK53MD|(Aw!eJ1jCrlP|86Z#4wTj0LUp&5; z0hxpOKs^2jSw2;3amV9ACJ~RzFTy5e5x{NA88qj>;cdn{3S>N%&u;;Q_YYcdfD8cML$f%~YVfZ)EJ2K2))vBZ zHXaY``Zn)s4p}T%5<_$C!E>v)idx#+I!?WOGw-pDWJw0i*`!08^`mpa;#$pKkE)2I zJ@poc?y0v{)?A-t_w2w+=;pQP=F{Nc19bD*c428`TP?U7we_)@4<9`2Gucb^(97zd z+r|i3W4P`3c|TCh3hZe9hDDJNV9>QuhEh}T3D=QtfxZpuIU8;eX>HY7wASVkNo`R5 zMtv>dfLZ58Q43ZcpMoIN3ut(^1b1^O>9DWgUdo5Gkj*LzHCtC?{;(F*!sO0Iw20uT z*?$hflgF~x40g#5!JCiZ%V0ghPBGd@k`~Y+TfTST_P%S~;C!2HFJTK8Y|chcYbUT* zt!>W;QT0f@7lzKX=pOnrf`2RU4)lE2qEx+@1#o|{KX5YHUsH)m3z?vtCZ z3P>x`N$!M$lTRW*;XtfeX#yhKM7K)6gd5)++xi&e+TDow5&~4Yu{L7f$~;ye!^_CA z1}pPvFo;<~S{Z#I`G*K*5hMT@P9?i&IQdz{aPe{y;$fyf#;7@ob`vdYxXcG*;%3G= z%sQ4Fib3ZL1}SQ)a3|argN4`S#U-$-QF+q~hKpcA)C?L2Wqj5$Safd(4W}#(8V*@9 zoB{y%Xo_E)6?lUYp`Pr}L>w|dPm;AM7&u1*=D}oI0U>gc<6s1UqZo&ik$((5{uzGC zui%=470^lfTT6%6+i&*iz60=11u3jum+=HjBlckB3 zk<}51sI^2&!z(YXzO+lR>@fn0^u2J~+P?MU8&8)nZJgeWZtdS5+Unj`-(S8HeqV!N znZ04xNq3!q2e!ybFxw+1i>{A@k&?PTcQdaCj})0|unpYvp|#}7AA)5PiLO2VR%D~2 z>?{j!`}F-sEBlY?k-kc#ugDZHt+*S`I}dE4v*_y6P2zBy>+tEI$1bd4YZX;SI@6MRt|4;>VZQLi}bY=g_T~I!R;$v zdi8qYdbKrJbYki9tILJuYHPUYycdAy%QdOozB#rN7}NuUmB3(O;*&rtELm6|-!uM@P3k;_&d*O(_;IsT!cEfdz9h>2)BTycU!^7 z+o2-nUzK`%$NpEPx_0b;RjPZ({#T{GxnutWzmk4p&q%wjDBZP7*+GApZri2op!WcM paF?=!ak_)<+@mXgMiTLe*qRnP9OjP literal 0 HcmV?d00001 diff --git a/csst_mci_sim/support/__pycache__/shao.cpython-311.pyc b/csst_mci_sim/support/__pycache__/shao.cpython-311.pyc new file mode 100644 index 0000000000000000000000000000000000000000..db9b1bfcfc97aca90203f767f9ee8ea7323b2dd0 GIT binary patch literal 7475 zcmd5>T}&HS7M`)kc)-{OV{D8uj0rzB#Mp+C9|8#^BpaF#vPqiVrcGJnF~r5#?s&)s zcW7vp2PcYDQl(l)qLRC6vUJl(thAAM=_)TPtya6z44PHYNR}e4^1vHaH_|Fkd+*p| zd&WS!KMy;=&AIn{=iGD7J@?L?li%p|RS0-4y!veBT|0vO9d@Kgl_$P$1jMfp5+TtD zatR?-w~#9+sRpfz#4aLk%@zc@Q334=N-}Vy0{p4L-}h4}l&3)}Jv1E8KMe8+a2^-9 zG$wHnWdm0~Wus86no>iy9=#0ci}D0Il} zRh8R&_h+=eXn`^k3_2^u3HWU+qDXbQ9ufiTfSpTP$D;9hnxGfrAtIb0qA{9?gcGzg z>CQ1=r!zr>BB8lZl$IdVkywBxViXaM0%+gaP;`dAM#QE;A8La_5;XL|d>}&5F(N?+ zs6;0bou8WuQBF55RD~1aC>V_fLjnd35HKhxR9*t~Ns5Y5f<~N3P=jd%btDuOsxQn( z>F`{R?N$RFaKcGZ4vh+EI05}92;%h_@3q)m$h#1TddE+Uc@yEe`AC3{QQjWdpPUKP z-e4j@PtFCylYlAjBvUJC_h#fY2*CQ`Lv5}0&}KY(Tl$!cpq zZ(B}srdHn6`f!xfx_Pa88&PSFy);&HioUWbMwQnQp%5gMnw0+T*{@_99O6S8q&e+KtJAq0)U|3o<&BvZKey2vQ z`t*0))5=e_OrhF4QV^M#&B(Q8}+H%Qezu2sj-Wt(^Y}TDof)RLMsX8NW8%@6P^%J~dVs*}LSMVzV3QXRKu9cR;P|c3vzX2?T&ILrmjrRTq0*C>vohgGiIgSt{GS( zkCH3{>*YMTPDKVM}BWMIvY%D@_?hAfFNX)R=64cSnPyfm(et!>A7OLC?p1Fdp@ z{-X?RET2!2ExTu6y#l#=23qKn(;H-UNqiJOSebgVx{yirqx8zXCWOs z$~7~VGC6CBoRw^u6$7kT6@sxU<;ZFhiBv8&Fb$cKyHas=6z7tSwibDeP3|>#8#ENT znUdR6ag;$wXA*Sfl8cQrd>7^)rt$s4oocd18l_t$3cSgx#53PsI`&`DC9{`W@}uYy zyJR-B7wu*2a)0T)c>gO~O6DlFYj;a_N6lDeYf>k9M!O{ntSBH$-e668@Emk^Q4YM* zW#1U1$TE8M&x_6Uz7h<|G#AetIexK4hE(l=BxFeS9!RGQiS2=Gl_52IAot0T6?-7t zR0v~dz>S8fdIauzj|=y`#deVcZZnKt!ReqKCGj#TV$><))+x9y8P}}fx@BC&+h?&$ z#ydHyDY!j*a=lDrrf`N7 zMI4mALg9DnjA+O3r^CsrUSDs&$Jgia^(XB`w7*}1gHOD0bjG;_FU7A>&lW5qB%5W( z=@>N!zJWvVJC$(x-Oeo-nLR>P+f-eFOe} z@Y~~;N9KScGOy8rEqHH|`FnsU3(B5+5RWaGUdhf9^9~T#V)Kdsh~D=BXO2L$TZWzs zN9XB~+#iVVKYTU6;>fvSg-jw8j77%M1A>bGjQT7##{0v?y3H5-SFR2If^9{1cCaX-Wsj;q`{9 z1dX7IC4`EZ5FPv_DQMy#|7ZaVl97m@I~xm-XTwuefLahVL0}s#RK=nb)Kr+Bm`Vs* zDwGh*CwOp@jLlC)LW1`E#Mrq@Coc$ifSQ541O-!A(1mZ&@=@zNvH_K+ z@_HaLA9DKzEFPe*3Hn@Ng$m?qD~4Q`3{%56C_i?YD|K7OMO^168ma; zE-uxL;9irY#Z4qs=2U?#B4J=U3VB0w$tLx3X+*)BmVgK6DQY?lVvoT*5maQKV2I7r zusS8h$e8a-A)?^BKyrf$6y^>ErMY2zCHvxYK(2tkF8LMD2UJkS`vleXKFK)f+GJ2# zN6|UZ%|{Bj2ReDDpo5^7k6X^gVn633=&MIZE(w(dmp}^Uyn|M$s6@gKoB@TmfA68- zMOrWwv-x~15vIejDD?}#FoD0sVX%pokXJZ@n?R9`;Y@A3sqG=T>i=@)s~Jw)!)tp$ zMN2b+W^=W+LC!GD8-`hY_$6jq9?I>ew)KWD=aIn>Tf*PGpUy<-XfD?%qiI?<}Sk#meFvYf)HPJSA2Z zPl=VqQ|L=uQ~F}2cIEN|`-(mFX139h4rZJy(+|2=x>IMe4fgbj47zgafo?^Y8rw8C zXWBN*9jv)ybzt4R(LKa=4~e>uKYJ7EJ}>G%FX}#zzA(3B{2S&r*4(ygUsr9o2Uz!j z2tWEvS6<|NVRmF{KRdU4F0-(9{MR=#H$~{+dKf~-M9pKO<}vhDMJ{QIfxgPQj_|G{ zoOY1c4#GfB{GSKf-keTk{43WV99cP%I-9i#P#&_x7YrP|Lww+ zEYbdO{z>TRTi>4i;~~iETxG4%3 z-?X2r^YC?^6rMFVrQcf};T=7k*~^=~smg44)eCdU`OkXsSC^6v^CSo>AQG+S4x!|!qcwIX-@Cq z^`4X}Yi&wj;jAv+>Pl5-8yx9buAz%>=t}9c7JGV(vvly54ltt(W@}C}oNYgE+n+LI zZBB41=1-BV!Lm%>j^2%?)X=qy_r~v!=X?s@yL$g>s%leT`|0GJ$@I^`e(QPG^`}|( zo#3X!y*jks{fuGX3bLWZrh|CczN-4t{gwO4TWh^vzx~bI3Pje?nmNolJiNp6WPa@m z=N;y~!{1e6hT|YdVVhEIS0UAg<=Q(3QwOsq5I;X>>g7$nY<2G@Zh4M7HgHEeyb|N^ z9v<&u@t#fG`W$cGz?;)azQxPoJ|6e6xDOC7bhh*{R@cR1T@-wUoCE(+R1wTwgJ7vl ze7Ag!h&=HH^co+Wi;?q@&`~M~`0y(tF#?)xHHxB}h?$jsw<~lg+#n@w({UBry^ZYL Hic|axQx}iY literal 0 HcmV?d00001 diff --git a/csst_mci_sim/support/shao.py b/csst_mci_sim/support/shao.py index 575dfbd..be30f5e 100644 --- a/csst_mci_sim/support/shao.py +++ b/csst_mci_sim/support/shao.py @@ -22,12 +22,12 @@ def onOrbitObsPosition(path, input_ra_list, input_dec_list, input_pmra_list, inp if isinstance(type(input_nstars), type(1)): # type(input_nstars) != type(1): raise TypeError("Parameter 7 is not int!", input_nstars) - checkInputList(input_ra_list, input_nstars) - checkInputList(input_dec_list, input_nstars) - checkInputList(input_pmra_list, input_nstars) - checkInputList(input_pmdec_list, input_nstars) - checkInputList(input_rv_list, input_nstars) - checkInputList(input_parallax_list, input_nstars) + # checkInputList(input_ra_list, input_nstars) + # checkInputList(input_dec_list, input_nstars) + # checkInputList(input_pmra_list, input_nstars) + # checkInputList(input_pmdec_list, input_nstars) + # checkInputList(input_rv_list, input_nstars) + # checkInputList(input_parallax_list, input_nstars) # if isinstance(type(input_x), type(1.1)): # type(input_x) != type(1.1): # raise TypeError("Parameter 8 is not double!", input_x) -- GitLab