From b236bcedd5c335fa40a82bd0121275c963398d08 Mon Sep 17 00:00:00 2001 From: Deaglan Bartlett Date: Sat, 8 Feb 2025 20:40:21 +0100 Subject: [PATCH] First attempt at SN inference test --- tests/sn_corner.png | Bin 0 -> 1436410 bytes tests/sn_inference.py | 276 +++++++++++++++++- tests/sn_trace.png | Bin 0 -> 530275 bytes tests/sn_true_predicted_c.png | Bin 0 -> 53238 bytes tests/sn_true_predicted_m.png | Bin 0 -> 45057 bytes tests/sn_true_predicted_stretch.png | Bin 0 -> 48279 bytes tests/tfr_corner.png | Bin 0 -> 1221391 bytes tests/tfr_inference.py | 18 +- ..._TFR.png => tfr_likelihood_scan_a_TFR.png} | Bin ...lpha.png => tfr_likelihood_scan_alpha.png} | Bin ..._TFR.png => tfr_likelihood_scan_b_TFR.png} | Bin ...sh.png => tfr_likelihood_scan_mthresh.png} | Bin ....png => tfr_likelihood_scan_sigma_TFR.png} | Bin ..._v.png => tfr_likelihood_scan_sigma_v.png} | Bin tests/tfr_trace.png | Bin 0 -> 396381 bytes tests/tfr_true_predicted_eta.png | Bin 0 -> 46960 bytes tests/tfr_true_predicted_m.png | Bin 0 -> 47587 bytes 17 files changed, 274 insertions(+), 20 deletions(-) create mode 100644 tests/sn_corner.png create mode 100644 tests/sn_trace.png create mode 100644 tests/sn_true_predicted_c.png create mode 100644 tests/sn_true_predicted_m.png create mode 100644 tests/sn_true_predicted_stretch.png create mode 100644 tests/tfr_corner.png rename tests/{likelihood_scan_a_TFR.png => tfr_likelihood_scan_a_TFR.png} (100%) rename tests/{likelihood_scan_alpha.png => tfr_likelihood_scan_alpha.png} (100%) rename tests/{likelihood_scan_b_TFR.png => tfr_likelihood_scan_b_TFR.png} (100%) rename tests/{likelihood_scan_mthresh.png => tfr_likelihood_scan_mthresh.png} (100%) rename tests/{likelihood_scan_sigma_TFR.png => tfr_likelihood_scan_sigma_TFR.png} (100%) rename tests/{likelihood_scan_sigma_v.png => tfr_likelihood_scan_sigma_v.png} (100%) create mode 100644 tests/tfr_trace.png create mode 100644 tests/tfr_true_predicted_eta.png create mode 100644 tests/tfr_true_predicted_m.png diff --git a/tests/sn_corner.png b/tests/sn_corner.png new file mode 100644 index 0000000000000000000000000000000000000000..daa6c22e69e4bb8f270d9c303311a78acc8962c7 GIT binary patch literal 1436410 zcmeFabyU~ex;Oem6chtH5W&DiKt)nYMHEp4=|+)GL8L=i$TAQG5tUR>kVa`q3uy#t zDUt4OxS#1-N7(ng?|bif|G49hJMJEbu{ToRnDco)^~`feO7iS#ItDrlg|b@woaki= zg|?J@ET_T0Jm>vnga6{YaYpfmjH%8It1IT(l#5qxm>HSgFw(!e^QN}Bg}$i?Cp#be zvBNv{Zrm`l;OF2l{__p&rsldFyCkV=@Jm*hol~-)Q0T9a52|qCFntOYg(5C`Qr6nP zr(V}uR(o!KaQxb3-42!);uI>byF0($y{Zs*FLF*=DZaVZN^b4SyBX2ph9A|XM5W_2 zm&sB4PH4nwmIFGP1g z67gaY|M&kwg;PV>1gZY(Q~Y}p4Nb?w*IBMOvoCycz2^VtKBMHR*eZu(t^a{I0=FQ%PANxR2 z`|~v~mhAt{&nRBWBTF4QhK7d^9y$~e&&J8=?eFi;c6`zGT~<3L66-cIV{dL@aVNKh zE9va+bwfi#tG+E6_t zy?fWi*v7_Yq`xuf<;&gFHx}OHRKQV@`Az2L=Bx4-KhoX$5|9guOtGi=u zd@U#l2@hW$9wpk}m?eDX%$dbSJHdxjBvy-$kMDNt!u}aK@i^ht{+$bduB`GutHi3( zh2#X(Mq-)6gd<=#sSbIttdUiZ%DbuMRCLqB!h8(a0+*+RNlgTZc{X?AbFF zGx^2elZWrYvj01`P5QTOTa{G|O(9?0h-xBiDa!4*Zt(KW#0u z@v{eU^J@R=O;FTC&rB{-Mpkxdl-jq9#QeOxT^t-7PkCgEGw`6&#~bWUPEXdpHD81i z9BY{QMR4c)d3b!S1(=NPiSEX0Fv(G>ki>sHO6?rvIuqC35w&i1+wwI(zczf>U_-jX zJhht#-=&53m&EP#lp^@(90hQO(a4Wecg>)+zdMud-@VC=OwO+_`jv%Chcu=5MCKn* z*+yTi=96aWxYvV$>uPO+{>(`ErR|pLs4D21?GI0{sp&A(xpwV{!T0ddxtx;~)qD5u zH4GRZ%w;ae^J=n-u@zeQ!wHFWd*^oMwarWp>({4^J_)>Csnw8fWAv_wTg8@P(f6t2 z0y7pvJSh){&!0ahnAAVF8E*R)eA1Oaqpyu|dR)k{qMB2uP;EsQ&qTW%4DmGJihvVwn@!)RFc5L&j8= zrM|_3Sv%$|gTiGiH%rB;XE4ZhY`0afHD7;cX@M_Cn0CRFo2&mopt7K1W}Vd|lf&)g z-q{3KN9U&73>zd=kJ}j6SK6!p`4eS9In1w@WH`=FiD(Zb#~a$js>+FFo!Ccp#J1^H zXjsf_;t%Xi)tfciIqlm`Qf}Q|z5TMhyqJ>Gwhj9we?=!ZOno-A?5*|LA7nptX#>-i zEn9v+>`B>4u99xg*W6iIjrmi73y)$%cL3lvveubLuT(RHy5_r_HcQx9m*J{aN9Jey zooJhF4s86HG2Wn6A~+s(y|5|Ab#koEG<9aEz@{&?y)sSA-27Ox+0zoC>H7}j^}f1+ zHn-I?Y>mGQIF(cymCqRlRxuZG%+4h(t?T*&06Rzr?9Yq;uLZ03|fu(>0%!^k>8 zyRn(i7B}mRd0t<5_F-bMrL*kS%1xU#8NKgKB4nYg>@xY{p+kpKTA!$JkB=;Bs)-|h zG##9}#i6<-fl-R_E;oB?2XkndkC-)y#H0_M*K9d8f*-$o@19{gxBJ_-$0%JDk!L-b z1PvV2@j$UTnu|~LFWo69C>VaeQOi@1qL^qPS1-A!%O>PGzWe-S!=0t{f+K5bcidPT zBIdC~b8co*05_)UINAF6$&+J7)p5$P8TP|PJQlr)RdW%Olus?!9j8)LzKlgU+!1v_ zlh$lZYfCEPu8ZqUs-p?=<2KHJdQ{7$M2>1EEhgX`%Yg$2b~1RZ_YmWmWQwrpHCR;M zj%O(?vpul!>R3+=*OFG_q?+61mxKfAOx}B(Wp;RcqAw++vSszQV}D_(FN?+IGcHTN z8o#yG@UI^z8XrEKI)DDW!O{-P==5}dJb|n6)v-*G2E~>~pDfhc*=@7?Hg(`ar4D=& zn(p~{k0bMIl+wFullQ%paE&aX+#Obo7Hx%BexzI<>8|nwJ)8ge6t6$@)vMXA)H5N{ z@-ghO@oby#Ez7yd%eW@VILgb*E4bECt+%tFpx_`I+p8-dZ~O61&mTN+fN;R?u^K^m zRkK}#Be{=%Lm=q^zW>5+3(QErd(cldSRrv!BO0AzirKD9f}}4iVJC@87ym5g+CxT36o=*y zZWo>(zk@xpZ~14t(XP+gx2RJGzihW@xjU`sHnrojvd>2 z?Ar1j{5Qf!Q(5JsnwaM&XtAufYdDN3555|X;w@i1#Lu7K&BynQ&`+c4BkJ~p zXBzB>FH&ySTPg(F4~gHbzC!D~QzTZIHq)}Zl0w>%pX?oVhX&l3-tjc)_U@8{XpScJ zX~v!U3k5?ZfIn7bzVCfM@65o*&Z=1E*z_Sf6r3HeUcDk)E$65!&FY_-(|s+|{aNiP zEe~4aKCaP1WlBe3NznWDGA3*!(z5Tzv*F?s924o+^42RS8k}KSze`Xs0PDeaq$6jz zL|7xFVtS;ryyklDj?RXR6Yq3OUIzzL!d18L2v&Nh!&RxG5o9eOdF4v;zUYsO;7&vU zf3ye(8(ZrwS|-gPhq30N0*?J`Y#!*u<~7>B*Ngeu2e2u>g-h@G(0qjS-O|(a+TUKD ziBiva{)9~hZ{l5=fAPM1xh*HF0|EjrvaD||li+BL&t+Qqq2FV$Gc4$sac#m*)snAY zDZz1ZYr|(6&PzyCsko^NUQ5njXoOzep16`rEbBu^5qW^ za*_Bf$F||3ANMPnFv%OewU&`9N?sg`Q_r;LLp!98$Slau*DHU$n&K%u$J?Jiyv8(R zbn`r_yX){8fr0GAz7U7$F@f2U@}BqonG}}&`z`QSC={*J9|s=uW-hiq&E9@X3J45viy8fEopE|+F=1V zua4JS?Ppvw3^26Br=lImf*HN9KXdv4s$cxaH7Axs3zL{OqaWz0%{)Yii}-imVYp8Eu9-)!wsxO z$pHTN`kQr+9zPa6u#v703Px$lbuFw#Dp~m#PGLc95*3}P zBbUDD^1gkyaaMn2YiTYnz5)P?y!U#`D)La694D@T@9|W%@|;hvO>1*IbB4?1XT)f( zFNa1k_UZ}|M@J#={g-sGU|&Q<-S4Y^U;aiWO!I!8_x^Xlt`Jju3`B-%|EtNB8 zKJ$Beh=mws&5n_V_~S!Q6JIH--uQQETHNzT0!x@Q6_%j z*)g}UKwCOz{7^-TqZ}81Tn%-%hmuOYW%a@p&VDm4W9R_vY^pVk)kfjI(NTKKbq)XD z0CePK*x2wjo27gij8uqi!^>_bB_+gD=%ry*jODS(qUNw5$K70vFVEcDjK#C>Qs9O8 zx!Hs4?62_;*lI)A4S2`A?fw=iZO;w2@|x|``f|Uy-l~z1>-5PsKap04toga=9`Mn#%3($kJa zMc?1ZBzFxk9#V<;H;wkU--PqQzQc|8`HF^N74z>DffJ5xTy%h?M|M*y?;Z!7Z zN$@k^(PUND+~=R4?>*zQ&^gjT%WR)5!I7~vdzU-o2e&QiUtN~1)Wc4-s<+&N2FB&m z86n?E!T%Bskv{c8D!{E=PRqFPsbQd&)I$C2^wbMkG9tP==j)?Ao%Qc|E(uMtD5zO( zYRUC%M|;pkV`6lg8Ki{1e!b+#)lbXY3ZI=GZY$a*XnXubPu%JH#*;h3)m!VaKEQ82 z;)RtQ&1F0d1bTLdgEFt<0D;%u`<&A8@}7%Dt|xN|XNiD}Ru8^1 zfL=ybik7Hu-MkFd(z@e?``6-P^>)=H(Z~@rQmx)`55d^DxQ8w-ORCN-!shQ2rKt!2 z%W+wDSUnx9_>Sn-Bbt;gs%BJ;j?-+XZ?B?i`baPQ>$RWB)sE5^k9~X|d3jw<=8Qf4 z=>Gkt*)?^hC9d>B=NdB|fP`8*LOk0Zon^X$RY~bdZ4V^0*9E_wTQo8tawlZ^(FdOq&i*?-j7Ks(IdC zRj1ie6nxW7##t-r=q9GRV}aI|CvDoi*-N2Uy_mOURZLpfAxg3EynUd>Tjyh(b4IZg|3iNx~MKCB?6+NmHXgT zdOD+XR8-XJVhfp%JzR2f?o~Se7o+WC`8ilw4e!Yo|6wWOm(|`_I;(GTW(KgZN3CDA zo%?=kfiLa!_TzflWJTsOI~D%Q^S+gTBu(x8O-q4R58rbv=UCQn+`Rd+vhw*WS9YO* zxpWf_PE9XxM9))QE`G(oMCla zW*aN;Fk5+KQuG_!wON_hx%XYuz?GdsT~$(2%6Z=V2VZ)P*d<7_{h+E%HZ4J)v!3b< z+XmPO4KgS5`*+#2o4x1DUtjFL;bpm(mGxdi!g2Hs#Z>dX_y)A356xy~7tR|T0yJxh zdRM%ABVD}p;8iA@U+bQ*tDznFuTRj6MVAFLB{a!?veoxKBcB&^t9=)qir}XX9zMLZ z>A7J{^LLp`d-2oet>-CT-rlak_JHb#H8TCNJP&Mq863QG;^SZESl%|+y`jr(FMRg# z{x(jRZnsk}n&kWy$@!Fttlf$Zq;5NE8O%m0xzAK|#eD~j2=t8UzV{qJ zY;AX0W16tne!8q=54IccBboy#X`W6qt|(Uz@ypb0S~mqcPU(L^0eBM@n2+D^rS!#V zsskJxx}Z@MLMyqAYuEk*KeG2FN_u|23d#oX=H}*>Rn6ddxVsaFnoc@VMH1vo{6yESH9{DvQ6$EuaF{zeUytGMjo^C zeDsca6tJUqV@3>I2A}X&Pk{zdNGS%132%yNpvj1fJ7Udz?=z~xNjX~}-#vTw@S3PK z+l+4R>yFg%6y3c&>c*ndAbjuQno=t0AYTeNvI?=p72~fiEAr#n&&g@w+HODic}vgR z0(Mo4BCaeY=ggSV9L{6MzJC3xSsiy}CdK+_K*SCeY3UqiDr#$}l!pG^Po8{nnj2;Y zS5MZgdvK!fD;mpbpq5-li&Y$%6T3;@3pwk#54#Rth!knotl!!=g7`78U7W zJq8J~8X9ZWGa?*|jz-R%%Z^witWIest||4!ibkV8uddGYLhA1wH>+<^QPY%$if;g5 zk%~MjG(TfuJJyr;d5cD&oPt_np7-s1xsoEFAj4i0M|YRN8dDrTH6& zbe7F~&Bk$I_V&SCrsj&|utyC??)enZbyh}i4|1G32=dXBP%dddJ6iQUQo%DaHuE{| zF^Dv(fuHv{%mH(T@!vzqS+gGBLe9e4D};15+4HXC-1+lAAInH8_1iYfb-z3t@-6X0 zzvZ&#tb{RZIm^^;?}Y+a#t~Vuy|sxgDEB%LKn^Rt{mZT|)tv|$fE(YtPGRNi zRH<2Mqr?0#FJL3gg2*+7F{!tqGI^n=UJKdh2nnJ*?#joeQ3>OgNWT+S?C;vxuYR?# zxiR5%x~BJjor}X35C44MynyN!|M%6yGqikTreQ;nObDu^8W^q8;Itml1kF0Zp4PD& zZPLX8!iUuFo+%?^1JyERq)dX!8h?ihS3(U^yOYUNMp5yc08{MIST2XOJCTX(457*JS@~a zr>THEe8R7?=4Uxdgy#kD8oO037!OrRl@t|`e>GiuQ$a2lDgupMTem@8=8GrYiAp7l zcveEE!KPp*5DkI(nc6Q!4kQQl+0qkGrW5}GREqIeCDk@M0ZPvoSB-X7syuN)G15R=k`?- zdW>J9R~|#%{P?(|vhoF9t^l9N0;L4Ab*5RvYuuR~d3%%VH#>*^lI`NTu(~f~D#fS$ z^Eru&^f9Zs>dFUw(O1OaWA$#nxJax<&=J6olxrYxZPKWEqR+B<I z@caztwr$&>Pr)0MyS8+fEtm~Ju7s&Qm(fgxRbjLRaVuZLM~@6 zVj8F;wg%#@IGy`s)9iRBK)lDSpzVm}m;2j7b+4IySh3%1OiA~VpI-sMdodpWN`lNk z6$hNy;F&8fadEG|{{?xSHrsZ2MOjfBrERYw<)agy*s|=|GohMMG;JDdI#XS5%HI&~4$3 z8m~6+F-g)de=YC2cl`<`>n19K>44Y^$p-)zrIPL&Rp#uNO3eE{ek3of0W)34`c?cO z9Nk{E?Zd-er_V`9?0Ix<=`>VmviJcimNfnHzNBvGboui76TJx=@F5e&wvp~S7}+4o zXgP>J!C;qVWDsq5jW!3GkL&#V^HpFruv*A5v~#f^1W2FOiq0Pf$O3Fkk_(ckXy>9H zW!o0F@Se&lUskr}G};fZh7v}PP5IZUT|{S_n=Em9J3CfKn~Q${i#mAf3BRAWo)&73tNy-OHw&5I+6B&iVEu*~=fqqu65=`;WFpDoD$} zB?@t0S~t5+v40*&N*<=wlw(Od&t_7G=*K*jX%+BR4vvC^4`s zZZHhGCCqZVSlB zde2L*g zFH}}m_6Q&#fsEjUS<1=ojf`)%>%`4#+`KyXlG#Sv#IZF=@9&YvEDr>&nQ!wL7Dj2i z5ux~t$M)fT*4QKFzs|Ll1YMdQ>s89Iwe%P?n0BWK&yH+B(*>q2jnxpAk5bI};I=aF z>(}s#$%@E=$YjIVmxm;$#3eo zw}qqa_x`YhbfT1QJ!o69(*5IOUte!1gageVp_-~~<-c(NfKUe&Smn4tte20^nLdbj zL6N6AvGITbXhg7{JHCa`LY@%x6q-6fA#ea2`J{?vc+}xo59n);Q1z}~KZ+%>mvK$` zTWtv>C7R8wGhSHxbajWsLMq*b+7GY~El_@szykhdiH?@1ngYkpMf*_DQh_Q?Fc?0u zii>`EwRiLr+B}#eu-^JL?f}{LsC`gWM-?PztHfE>bE(T6D=Wui`n{MoZK8UJ5Yz5Q zMD~O>b@a-IrTZlPzv4;q4P;ZhTFegz&BtB+e3zj`URE~BG8()N%Wi3Qf_^z9K$u8a z9^xo6Admo4hc~MjI*P0vn7>@A&NYsXbv>FQxL^k4HNIFVZ?liSN4VLcm=a|6Ovr&B z%I%4jn-4b;m7(uPYEN2s%qk{W59(0$K~zxEx8zIRAf7nEDAl zm0tOU=|%JcTss5=&FJAna&`2)O%w{~`TK^^!bsWZqmXC@?imIv$jSoob|MHh-smJu zT-t&sY%O?lRKJ<>EXBh~d|A5}YTA^HlIhBEEI>aHS8srbzRYP8{7Eba zXf~WqbM`SQKj{F8X|jnEg);XHx~9zQ3s=B-MtW+zCD>E{!n27^+&5LZb0Zz`=ssVC z=DGxys>FF0&8bjt$)8{lx@Do_gd^>{_&mxY;Ci^~{vE;L;o)D7cB@c36ltMG6XTii zm8Vd-2x24@6fZGHxuc?r19U&hdZ8AC@+UOVjlfFXv)rpcrDaWD`mkgl{u@+p{$Ws1 z38dX3!WziG2b<6Kt>j352m|FQlyqvdFWYpDERE|TKeY_*{40h@v?Sn!@6eDr97uxR z;8&1Eiwd?QJW3TRXD(k?cnB{(ob!)lf)-&3#!5fSJc>3(p7onIZ-{tqKR4Olt-4Z> z!NV&cz<6;+*itWpYKmTXMucptaXI4n4pyp|#@mh>JmCGg9G28(fA8nd&sW8$xeh1k zRMZ74n7Rc%c>>2?f32j#ELAqBJ|(3+(2j_h#AU;)4M8;4|N570J#eG|Pe}V`jZJ@> zH(Uz~xZ0vb(ACgHk8M$;@^by;XRcC(SoNZ$I+Fq-1R3K-gdy#`-y?c2Zt@T#&Tw2@+!uz9(eA1+32;09?0d&Y`c_=`kBo5pgg;*w{j*$m4cL-YVNuQai zKEt9#ccfoFbJcmxK#9|ww~9%scJ$p&*r>!RK^^KJ0OK9R)+kRjtQvj5Y0-9pNFjGO zvFcUD>|nMZpn}};3pciun_U!MCftCOFNUWhl8s{1m7ksK#56hV*o4$=cfAtlMFN;qW#A; zgq0C@hk!#ht9qaQta%|`ozboek=C@KZe{94AW(MP%g~~fh)4LHpSGK(=iqkI<%6d` ztmK&axfPF!=&ppjNr#A09YGaACE{bxnlXn-x{ceYTG^+s+X9@Mq%SJF6vK^bh|~;& zKxA#GuCA{4V}H8Br3le+^6rNyw!9Sd(1JpAG_@ECmm8#$6O1vMCc zQ0sbv_wr{W%Kj0LPtw5Pkb=4mf!~^t%M%$)7YP~ZIhBJTGkV9*1BgI`%wgz(VLr85 z#?QP(Xv`vtY2-h%_`>t{Q z^CqHMk)Q~aMnLhB)jU5$2u^@QmHP88D`l&I^$|R1PHk;%0%bkKToru4v%stnMY$tt zFq*MY*}W5+6o~~b40YvFZIR!eq3O|=R7;A?S7=n&ft0YyINF`c`o1>FUz zhkRtWyjwNy1{6nLUS8@JlZP16DI`eI8Yz)iT--*=MQgrKP9`MJclML!xY$c#Vq$S@ z$Qj^4dxE;}a7HlcZ-@*_4TBLux{e?e)CT5(uX{X+OZ?k{huH@aj-KEypkwC0hKX9h z(`e6$>%i8C7lN8+A2vOW3T~9pOCRuZ%u_Jx;ew*2p&9KeZ(YYz6-8i^K z-Fg?PsSS4hdk-GGV>{XvkH6a0x-4bZ#y;J4yN8qIP2)PlNRydVvAT$}Q}x z9`sko?@3Kf^)sj6juPeBDv+zLu1P!Cw6V9%K(= zn#7FID}qlsIt&9$oy;~FT@d40W^K0{HP1n4TiUN%_%sJ@rE3X#`-P`Ru0%iFwv&Uy zQ$~`{E#v|JUG%j_fU4M2uzhA)nCCU&tPxGS5M5fxVXE|v%>LW-%#x5A(SD-ezgKbR zstW{te%{~qi<*_9c*IGFZ2R^AD<%Yrn$RCx?(7e20s`8FfG|Rb11RW2*n8AwKkw{6 zrd!NNLFs4ss}=z|K?=%U2t1dQ$U!SVha+7H@EmGsK3BNwzP9M3|FSjICRYf{y~w7$;asataV6uJ=68_^*xI+MP3&A z9eEQ(|7JlH2%C*!hNzQygH8KU{W2Cjwc@{4Y#^2j9?qIW-(mgPfuLh=?8OTp>-h(? zT{xys<)f9a<9ENs(nZ(YhGGjm2-CqIq5=9T41tqo4RKdaJ8RCN)}#3tsWjSl$mK4i zgS;wDUH?Hn7$SnfpbgG=OyT!EheFT+pvq`W4uM6~x&R5;`Aj#YXoUo#4LVy51 z?IXmnP{;r+SEDqXTsZ~lEBNhC0U~c~K`R_z#I$PbG?z`q%0Wt2R%Wnjbl@S4g?P!( z2MPlP{YVfKjFl(cyipkf88zmV&Vn-{mBbp8;sQrij8%|eAlZPVi94ZmA!4-bg*Douhx{-o94ybu@ME#c zoQoDo!eCdz5{HZUyC5KC<>lLi?0Mh!yj=;6UK?=yf**G;&&m>nS@H8oIRMn&5xp8+ z*(Fz-Xe@!pgwzyZS|E}daPe{x$fl$Wd{}9CAI)z)4L#E;f*!u+<*gWljtaa^h}b;UGu>E1e@_((WQ}!lq}Xf#JLJ|0pr{@83`S5)&b*G0eCC zF3N5jn>(um8-(nbocinI6#^$OUWA|nU20T+*?<0OV|NA? zBVdliPf2~mCB`cr4j0m3u1j2ywlCLYWQ@AA#^;6Y<+AxX_A5mt6cTSHsa&fDoBeR`j}d%>&s(19 z2JQj^MFD&UY=uaK;nKNodu-=XH&m%$rZEbiU3A%RUMnY>Wi`31BrH*QZekAv2;__0h^LIi;4qqj zHvklr8j{@lKPffcG4D+xSCuq7b-z>?Z$8URRhnq1K(mH%dA0{n!#OCwRdgRBmy95q zX%Z4+SE;1??Ulq95=kcEA(Y`uFeUDYg=D3S*^FA2R3*o}kFH;geTX}Syi^L5edKR; z@7?Q$U49JUgM|A@5GnZBg-leaN#v)WVb@c1A!5i^BspA!=dNISYp%X0|A+u|Q`iI~ zTytiN+QuG@B%^I)zb4n;+)kuY3IW-s>4Rd?u_{thYtw}hO(8~{>u61+`;^>C8_SB1 zgQL;inVFi4ORXckQb-1bwM%U~tYOpYLt?W;q$hI#IoHW^!~_Vavq@l0<*f+1KD-T0 z71(js!n%Eo6d)LCVidY?6(W;McKPrfa%2Lo6#gSH4zA`x#@zrp{gIEs5MRGy3M)Z;qXyk};Xi!SQj7U&BCSy#Q zWMi9xJ2qboOBrIZ*3S5L{IO;0zBI4qA7Ps)?v^&j}aj@6yAE_238HEiHjI0-# zzagPBaL%I$oGB-pWaJDIc!e2;cfoJ%EFeeZMA|O1l2gpk%yDU|K6m{0)P=#I8Se$m ze?Tm02drh;NH++EcO7YSDzpBK@v<1T6LUz1nM1uA_RX4E=1jCc0FFTzuWwPak$Wvd z4w59*$>xRna|KN(XDjI>sK;MsQ0hQOPHwUtovwbm=!JCpRvgL^LQ&&ud#XjLMw5Mp zCW06jkunL?B!01|{xoX8Rdn)j&l;yn8XCvo#sqg@k}HDk0R5gboRxFnIji_0LI4 zwgGatz@Gn!yKIFHrG-vKv;5Ss-f@HXHy6;x)an^-!SA$%U5AGW+&oj7a0we>m?gL=Pb;Q1OnGp2i^qS3x>q zZh}U9B>5V!RmI4o?ywkv=UjfjfK4wIqE|zNQA$k{s1~v+X2je^-cFH{!opP#c0d#c z9`<*f>IA_$ia-Olv9wmeARd^Ngw>Q%5JKCp`4`WW6R^L0;a&7z&TY+F%|GVSUL&h}p+9YxW7^ z<+q{QXMD+o+#3fvJA~_|9w4#lBCaY7h>D^bPhhAmwV02@wDARIXOgA*`Sa>OF1_E4 zOYw?53IU4(jZR`&;>-9&@T~dpZV$4uw&2Iw|HmPiOAyB#BlV;AvIzzdOpSWC)wxRD zE?Vr~O=n~Do`7~tWj;huSWBg%wg(~$fKJSpOmr$H`+<`{_g^4u+{*7ixMzgEr~pk zDC9}n9iN~GESHVh2G3p;kq?Y4y#fpJGd!R87Q+(*ow@U*pu>bhrJ@#^>ez$ex@XhV zdA!TxhLurYYCEO21!rOMX{4*-9<~4+SAv0|tslYL;d&4NONY>8?B6BJGyV2F;B`T1 z0`^i2DVqKQx6|qG%3cVVnCK;*xDk?xzgk}(-$_DfL(l|H9Lm_ho z!N+=>Py)dvScx$M5C=PYn54B(CU7BSoWllDO=37h1k*AL0*lW{M88YHy(%T>ts~L` zaa= zg{F-j2+vtJ*mkpJ66A}k!+zUO#vBRQN1}JFatR0yp)MrlKw_=Hz~4t7J6_LptN2ij z4DwoM(kyOK(=PF8p6jW3hp=`N0_3!snAat~2@)aXZxGu7tPSe?OT=!8wWMm6vkZ4m z?L$KCBrR-Ia|{@BCBbZDeBdI)G6a-Sv4~)U^fvdDdc+rv;ZZwZ@FGZK9Voq6uZjHo z!SsmPYu(!3ZrY}F{?9U;%7`4^+y~Vab*TKkmFaJ-9i8?9atD0Pe{CRk>jXlQaFq#< zLlu626FeGV#s^PLCvUv>SuxYzW?+&Ce1{@~>JkleU$`^n!>O@>#-@iMInrsPRg!hc z941yLX?57~G|ptS1q~WQN1N^|$lX#GTpI-T!m1FnRFm*uaBk=`Nf4@l2D9Yfj|<=9 zsycze_qt!#D~)E=eT^G=+S7V~|8*zd?)B3)T4);E<=d%j;S|GU)WrK|bQt+2qOE-c z5QgYkKu)csWg(rSi|LjoNV0nnJAjz`eN@(9FN4Q1r*ED)?Y>_;_Ass)#&8&nP4XbN z3OfprV9i-38x{&l?xCe=;(QlkHIu0~l3~ZNEh^5fd|?)FVr>lB?7<8OrqxKw=v}$w ziFha|7F}<`h5&+xF^v?%6$ER#JlV9du}`9!&PXp~7-tM<2v`+0WkihHfmtCy^v@|N zFEw_En)X@Nyzi~`*i`Vu0wOd5lrYD!j53H?RY_?L%c@-SYZkoej+IvxU<&anIJvHb3{Uk z4~xyb{2~t-pp)e9WPx`|Mr|ojfmT$WbVBfgbW8}0C$LqS+!@W6V_I(4=6gpqvng2M zp8%I|y3kH!Ct&*=-aK8D71WrYR(S#oe~N=9SMJsrl3fhxQ{(-(p(}uR4dF@se*7gG z@7V1iqyn$l45WVIT8JZrvb~mMEsdwa4N&$ zfMS}>yNzyX?{})3Kdzs}KVtR(cTF@>V!FK5E^vABlT746U?w>pVu-H7ye?|(eE^#4 zka&h++bkKt?j`!rw`5Zdb%#+o;&BnM{!jXXcmo*x!MFl84N)<%^UZ-+2eFurLUo3f zdG_$;5X^b~4wB6pXN5N^R@k;@g?1@LDH*H3im45Zjg7W9Zl9FUYHw@1oR;QK9c$0z ze1CrBTr~FrNOiQcZ1{k)2gsoE5Pta}$u@&Rl0QN?2LRwR%<-|QB!S9=#BgNf%~=j- z?Hm@NKyJ;!WICcUjF#V}q2P1y;=VGe5s6vij3?&_rO^kX=@9VTbqJ@!fYVDJBEN%z zLli17B7|+AK<8q?6av+P-76U~5jC4x$RC&UfcWkx&90xD0tch z$SM(CjO;wz$sR~&P!XYOV<-0@xP6j!Xs{(0CqdN3VK@|j?$BbF3) zEyfMo;8kg3i?zYLU84{?-UC#*sYCvS!tIXj8&<`# z-iW1pK@TyKAvb2*Q_~B}Gj!bS)EKi~3}hJuZcp>oP(HorkcWAr8_q(8m4}H>j?{uU znTbLMkkDc10W^EW?(bw=^Uk2+5L|{^`IIA=P=_Q^HTlbsjP_8Cq$!7m&1pEVE{$)w zdEEGi^Mf5HE@LDDK5YTGMu>qEDeTKKY;`HU}x>u$e?!nw?g8sgN7XI;A))|c|fhvK0uxN1l z52n7aEee-Zl!+}P#{fLqbG8g^n~ZZ)p!*g=ce94u$_J6c2ISs#!2Z+U>j?V8byUq3 z@|^mw%9GSbnVoN2InWtaE-tZR1DzXc;YksU{vxJDI{_(2G2Qwo>2m;#A|6_qQ(dg2 zxd6daXX8kd34nB>CX;1A25YxM5h2;m47>3*F!a?fBLoqlXnu>=Itw3a2k3(Q9!3t> zB@iPK1ZLndNpG-|?lm^rJAM0$lEPQD?Kn2+8jtA;dq@LL`hndiUfcRif06&X05+Yt zaVQ*;W}6QyW)mmT7#^M5IJurY~YK8X&N+WFd#h|6Sl8uO#Y`!!VL6bi?-Je{qGkx37NYC`q|)25FQ-#Z&^<0 z4hZxDk7dr!*^$Y53P~KnOPAdJDEwPm2093gdWsZA%s`TyDy&|njT(NKz&JpcHTQOt zSy(B7rZ{5_I z$Ye-{gM-uUdo2KlDHD)Jioq(s`k7{m11NfE6(FEp1~UVF9l(tb4S;ta_6ni|#I4Am z#A!=drQik_FH~`eSsh{AZMA@`=DlU>s3RlZ(x2!#LR4--eUL46odu?ELJ|4Z;wh{% z^ec}-z(NjxlS)X0GmP*G7289%P@rI9NGb6JVM0W@rb1edbCkdcuglmHDndic&9L$_g`dPS|^-o3BB#8*@aYFPM3 z7$YBnq0VT_zUR3owG(YDJ3uGZ4M20*P?d0xFG+vh58Sp5O$w(#h>axy@4AJBM z5q#uiL_Cx-q>6~71_uL1$c!~t?7j(oUwD?!mObXI3U188w^SO z_M$VW5|$ZVc3xxl>(SXE_@l@4N_Voq?>Po3kZ4#l6NB(3Y%|&+D_qB@5g9fkGs|gK z{l(H}cbqw}v3B<^{_GX|{j)1JI4y0$V{E;X%k%;_fU`%Waq=P(Hy8yb({<3Bh=fE2 z@F2Jsk}+a@ti4~6WUPWsT6mgX@Y9+TJ=tWg09S*#vYkM+WQHH3G_1UtV7#gbuuw=m zBB^@esFF~&YuaClEv?&4uUC z@7xnaEHXKio;Oa8w#fc;hhjfhql+FFYovQ7g0cg$3(ax8(xxq*z%v5S zXh@J2^pMUOv*E)?oLvVil~GpC$J~~I+AT7|_0K6Xsj(96>}pR3Tb%$Cu9%quS&)N5 zNa6yA+B`X|v_3&Y!&mKQjT;$1LO%L$oXH`BBipo`@FIw?O94tFIXu7@AxlS&g)0lK zovxs=HFi6-3%2(~VAeQeKjWlkAkh7I4T%1Ep*i6coorsOq_qP%Yb4w2g^e{582&JbP z{`f2GhwTm?c<`oXy_UvG#Oc9}uMnC}PcktSTF`nrOPY9aAQKtpC#hTlUkGo8m(Nl? z4kk+u-xz>5@(zPwMkCdLNaVF3z7(6h+8%jyLQX-%6AeyisQ4NFXzZ@8$c;u?{@*jH zzz>$?es8m>s9@}xvYTf6(s;QFt+M++kD$`&-t>F=0DqLMaxgdSp(CUjA_A#7(3FGu zrt4VXWr$f52mG9gN#8yV4ji*px#?Y^TC#(dam;n6SSbZ^U zs}!q-7nZ}^DsBnqISe1fV+WSO!@&@Y-uLiz=$&L3h`?@$l4LrjtvEmetg7Stn+;?V z18)*}CzA&U7s|O;!pZ+sISA@;oSW1oNp1q{@oVx)6bV#8zc%3@>5jyJLvP!=`Q9LO z?ou4Qa52rqSS3YaJ~_tdpIczIa$w7Gy0^KnVtzKmBkNU6d`x>{YLE}NNAK95vfHH| z5PK{0M(5?RdMgdWogN5H=>;F_cqVAaO~y}&vjWJuG#H1B!2gEidlA%HzcQ+3pN18* z1eCut$-a}8h;iF+DJPLdQ4!=@V$hw`HotI0pK%(FsxgJyjkk{TDE2NdJ^f-WIVX#N zqTV`R@L6~0UqM*oZ{NM^{*{KtPkZQNsD~PdZEM6aOvbHnZmDD2|4w`l`dDJI7O``1 zgpB{{2fkijWUhC@X}*8{I{FhpNT&4gYwg7-BT zAAk?mhqTr-QAwc7S|JhPyjs5D$>YZ)OEw6waUH;bUjhaYR=pB!V`5U_ zWEVoqbk4R*M@jBQW#wqPW6^x{A(c9Zxph(f^9!YyTP@(VBETWM6M`K;$YZnS_`z!w zoFnFL{{|zElAaVZ5fZsa%b6|FCCv($t@wVi_j6r8t_-oQM>tc|;X}B6-te zmSo3~m}6n`tZIGbeXzstA|2td$7S0v$BguV*kF^i(N@A$DWFcR8$kDr4YM-Y(%)<{2^={0rI zLXjl#O`@Npxxlp~ygDyW<Z2bxCjxVYtGN{9U5#BI8d)e%f5* zY{`ivmoe50^ri%5h=_Egz+wjFS9|+A97;eo+|8TZaP~>Qk4ez>p||%tw?6akjfK+aDgc_b7J*kxj7g&3Z+hVRBVT*yh4f_H$s22Ud* zKhaUlkg6o-+sGh=i33MgxPdpXz{{en-p-8>=U##0x)5-OG1r;GEURc zO_l@3yuN7&^Z{qRvg&W}M;eVxBNKO?9KMAj)%3KhA7sRu5H)C?yBht*^w7fNIKE;2dBO4xpC=@=T8Dmdkze5Zm z3!Io+Pe47t&Z zKy7x_(u!B~xE|S{Aedz>kXe~zS!AP@)!Q|!Sk*qNWbvccEMP(8?esMsVd#Ky^Agky zQWW99go%;CMsk2p@Ub7NIA~|>1sxroGQ_n=cOh&Aa|yralUUXhHJcd4!2Q%dyENeJ zMz>uDp!+?OUjbU;k>R5}36|S!6or=~1I&DQbL6|DP8@SM=FzmWT@E+g-q?-NCFuFV zI6`iD8LnNaDoFzl1!93WJS_jTT7A}cIG6@>0FhcE50c0R$y}ivTHpu~PKT)*!~t*) zuC+pEB8f%pY%*(z-<}FQOvU^2?>!{6MG(bPtp_VAw<0SSAt%=h&a0SYq%;uGkYuN2 z++mBOMU`2<$3_S>KE!zzF^O!9YXDfuF(;?K(FvCOXG9h_g4$SFomkTxQy8B|1Qs>pSy&8Ep2Mo{N4-;nhPfC_k=0XjKN|9%`w zLQDp@Fz}GCK@B5khqj~UoB(HP)Z}IQZ4dAv7OjoKo-<~uQ&ij3|DuaW#!ojfaVGQ6 zh7N7?_AqE>LXtX8^Q1@Oi#?B)Bf^Wr9s@9%;JfSe?cZbiAPYo#K58>`4RG}+DV&6b zj{UH>^bctz=kO3K6rLfOyCl*TUb=kb&2lhWATZ;Z9kJxNA`Iia(8rQ{hEo?v%)kD< z6)nO3$n9N+9%f!9zCW5XZW$2{5vw#iLb8M;K@WURLJX)6SFT)%l{+2}vru(?-RdR+?p^Y?FByqI#0Up5n3q^DzSQh`3v%v6DG!BtC(Z z&VbW7=RX_)o!bXf6i1vr16RWlaJ3}2PR>ci@e!xinwStZj6}U4y>oB zG`GxZF3enTl=@g^^Usr8NK~8Pq@9Y(wqnwbjEIrI0c2ZN;$)&(M8(rgg^GK3AdZRH z0x7$9HthGpI2hIm&P?SaIu}NYila1zBe%58hh?(XQ5&0Q1Z%PnuI8E@?WK$ zBgO;pdajaTFbECg=ph}v2Cx8%cUhiqI8nv(}~JQrcDh!#)fMqC&zGrR}BCj0ZwE#Bn9mg&JjCVlTP!~jY>E* zh*+~^pcY_wH=*=IXkzPqD9ay@&iC``9IvH!7O!HKnC!}LyaHL2B!weI>E}?hJpxWRUl;Jt|3`zWh`|3c zc;|!Ctl9rysfag~Y(FBGq!1s3SoZK3Fx)Vt=>*Z2#7-B7y2RE7Qs01x1!L(hSsZ}Y z!EfDe+S>BN2zr4@Cb6w*3hVED;@6OqOCzKBW<5Cn5p@+YxyD`6k zV?0>U+sUj7^8ZNVp2q70nil|sV!Zxo9kDI(1tzDN>m>S37z;Vf2Jvqq(Bn=1gq6M< z=j25HA7q_(T+jRe{@>YT9Xl(VGNQ;Tbq>c!$SR7AN>)bELfPXevyf3%R#HllN@gfY zk|Zf3m8PV$fA@#;z5PDlug@ROaYOZfy`InKV_etuxSn4#rp%qWVSbs)-bc%mj>{DV zcy-S8iy!`(>KtLaRb-{Iva&dB(Jp0i^qb%RPDz2pARD3i|JYI8i@bRL{9y20Q88Ow zm2yBZK^cajmmy;Rg|RX1p^Q&Y?3JB(v4sp%iBwD&>SMD9Cdd>y{S_M9MaLxp zK#QCeU{MCdWbE2QZDt2NO!!n0yxyuz4cA27gKz=vm zk?9rsNSEvGtAM#EIx#o7@YYu^L7wj4zn>BDqq^p%F^LZQBurHw5z1r+nI@)3*Zl#r zJ^zEjCHcBf^-KLf7FQG9_4|W2pmch9W7*7w`OA^#j>~#}N~DHQpL({vG}M^K?ZN7qzHaKYH!`b5q68}o7*WSFYB`>^dapL+2r*MX0)_Wa4^Ui-c zln}a>NYJE=T!@)2F{SsVx7VB>EgmAl9hu{`jErw2p=G40yJ^E$Zn9koC1VBQd4C&`2Ew1C!bOL;tSZ8p#+Jg@)E$ULhfw4qE8QD zc`C{agL0i(=*4s<$sqLQo6!K{FcR*2@&CKzxbynzwcRpkkSsCsQwVDuf!-sc;;Wzhx)(|spqlzpzx}+sLPHC7NICXUKt-`7 zWM4bdfC^wSJZ)WTAnEnas}q*43l3OeZZaq$FvzNg!HyLx@ov9JasjK9Au3?KNg_v} zq82qC+&YbzP*YDayn!<`#o>@`ti|q9xu`QSt$VmQcS#YBO3kOaJsTtE1*ky40vD%9D;LbFksq z1u2CO?Rw*%fL41S^`qa1z9gA9i1Ha0-k5n)gct&@gCzu9AM?L#agJ+;EvR6pJ>_*i zVYK%aT$+RkaoGbKkWz;o^11zEp|KeOwWL~P`0AjT3e3D1JpIZQhveA@M}NHjhp1{q z@r`u>?`6R79_Xn&YRj~E%CBo3c(${k^Ez!xy2^%m(YlrOd-Cc@VdP& z=U2#^GNoC~#M57)LhWo|O4DwyH}h-jjK_v(b@1BAcaAzY%;%A~q0vjT4&9lgc{~g~ zmtf(@CIdlu0}+JaZL2ThSJBn+GFnsR6|b*Z`(#<&23Nno<%9q+)LFw!A(I>JZwwYf zsomb|1Yb_f@+#vi-W!tqDk5#N9sbifEo+THK#q*mAtEwEY$aeS*vMRvQdq1OQd97m z>#C=&I(qbIa-%;+V>AL76DKwwXca4S3oLuiYB4^!>9j zLNhAZwQSL?KW_Rse%zDs4p16Qek$1y&A(Lsl$$$b_E?}Dezm`JwaJ+L8;K+MwUjW^ za_jgW2Y(pv3s@2NRj7uTDCb8dh22Li@cEjQ#Okp9gyrKrDP~i}?RJm;&GoYAVt8ZHcLE-S*62qcG!2hW;M6NZ%Aou`~i7gpmlTC5ebK zMkZ}YCNf0*Ka)h1liR-j{Q5fLXkK1k#)kfb1~ur$XWM(T@INAkQ0)`*wwU`;I^VEK zJt8)|q}nxa7VTQg=Kjv>|0exJBjqjKbY%7yWt0^#@FMz&6S_JJ?BWWW`D>X|H_KQR zffM_~kAS|PiGf^VAagsSjb#r$7IhgB^q2qlDrO^T3COywMd{~`;l?fz@c2OAJtDmg^RcvN>^FrY6sr@Q0OzH9M|Z(AY$*oC zpc|SMom3+?KkpSu4(QG1-d?wCu>HSG1wuS<{H^)9SE>TcA1xSGDqx^^6oPQ7X=wBX zQ;|vwk=!>*6C6tK1M}ZUycsvC@2PQU^@VuFV&kJ|En!ZU;Bh6jCG#*^pMzqz=FVTA zG+U?mNzVx49#+wHBe;Fw9q)h^NVG2EI8JdqjVYeMLC&)`-}&O4zv4aDdx%jskCoy0 z{~Omy!>YgHTtUqR(+&o!IdoO}3)Kcy-m^mH2)fIxL*uufhvO?t?Wff!VH$O1$a6oW zPQ@)+6V@20=!4>D)~q(OtOLsDn{0BRoBK#1a-Pu3q}_$+H?tR2h!{X5#e`e2)p0`i zGLoMz0|ugGSBRiV#0Ey?4S$W@%6=CsqyA!>hz#Rg{PYR+lb5ESsZ0Ic9i=+wZSAum zkBh|alCsQ=Nd3*HwrZEuEU+6&f&Ys3URTt-uIe&t_J>^$J~Rux`k$jwljkyvnMd-tSmjGAb%@U&2f%0tWmoy!Wk%ralChROgBkMG}$-s~H^Kj%$YXO#uDFp0*FV|#xnvmU|JbMOudR6Q0 zyalH{N1Cds5ftyPQBYl;jE+1os!g^~Fs< z8a{<^O`-$7HYoBb=mW{V@5wgJ@!)MM{uhMy$Lt}Q-GAd+TaAO$Fkb3AQ_pO&yhG6) zU~&#)0&EbiL&Li=Ys5^~YP@71{<%VC<0-q#<{zpi|2r&(k^mXfc*aKpu9>GF64PaL zP1*3%Oi5m#=9cQTQT&$E^GDsKIcT}DkmkqKz0fA^Ywnf8VWUTneq$ceOys80XMpIR zAPsmR8x~khF@|)`)eH?YT+_=dycc+0xl%T4cKz>?*ywT)Ejei<7h@3WOxSvgwG+Zi^OyK zXQ^;CB1S*^)4$ErJ@1a; z7lsOZAHA-)BXt|17a8)Hf7@eXui|iM2li=jMhe3z&nCt(oi8U^Jb6+ zHKCF|Gbw0H|9&y3O7va}YAf1{D8%<@joyYbm9BrM3|mtb6+THjWy+3TuiA_{W03v# z0dtQ&U2?kk?rqYjQJopHjFb+V*YhnK6nXQ-_U+y+O`r7f=oP-~L|WNc?=Hc;^o@#= zvPMSb9*@n#27OkQ?U$3md#L#jdw4vglpeA6V?B0X+^Z9*rDc0HcSS}-tXZ?hXlrz> zwv4y1ffRJlX!sZ3OndZb2me&2^A_mO$oipo4ng`U2_c!)S+t4UyZZt6E$9dD5!%J8{%^BJcfH zQu1u51NE=4QaPnt?K)t>Ld9rn|N&G2%7_c3rPb({MM1k@Ll>1SnS-JHx!n?DQP%2JjN zJ)5v}x0_q9QKOC@JJ#Xays)i4KB^t3I8WZpy3OQnm%ql2O~_e?;;sKo?Fh9Y)4jBY zsO=fpDOn2z4Y2o{C|3c|myoY@HdS~~GM!6WuBE%-T*9S|`ilVE%1{0NhCkNEuIT?^ z)u-(2iJ;y6W?pD|S5N!y=hf4uA8 z2ZmZSuDU8618loeQAy6sLGx}Y^dJ4x*7K`^pVWDzZJ@>D=0u6{1YoE*7=8hquC{Znt7JLsV7p3y4<^^Yh|{0aS!}%gd7^W>Rks_Rnt3F9xfKy=5qP!WrvN=O&QJVfv4zD*r{Z3tPVIvO`IM96bvZ^q`pyr27EnT+mhIX_ z^Ib2`@NIHp_-avpQNtVOd{||%VVdTh7V-Q3_-(WPWf{Ms#p$}2l2%*3O*RNf?jLzc ztMp9hdEZ6X=Lp&^gK3y1Jbe6EjG17D_M7@Ip*Y~ePR7YV`0ud-mkL3p-3~kttdRCgZS=#r`V9 z)?LNNB09@y^hfQ!Q%@SX#CvHqoLKur{3}L`9BEaausuITuhYr+CrgL6>e8h;z{TW$ zhYMK{$HbcQg9XETS%f-74{-7i)?yl9?8U_EuC*@AON?bRiILKDza*K(hM;KGrp-fy z9!8gET}oQc-d?q8mBFNZ%7?%H`s*0YAGR$6<^*oW;5`8v4yp;f`Ydt_g9|m>*_OBU zc<7y^(B2j~rEa+^@4azd;gc)7*DAK1@9LZP78-J`Oh$eulOfV^UR2v6>VWheK z>(SGvdtm9-e*ZE>s?962E^jsU|3cGxKQ?wR&cd2pE2U`Kpsufh=UG-qf%V|pcFVI5yhLR@z$&f zw7nC5$DyLB*T<-;3gyi&$1DWuoJV8~i`VmW_?)tTz)JHr^VMnCZz(2abjOA4>@Mo2juDjJ!e2(DyHbdEB{vUTme@xC7tq@|!AqK{d;WtciDxS|$T~ixL*?prEDZ9c5spRO8n)>A znN-$!#2Tpi2K_8pg+7NHZw8jClT(&(ffTcS`*!gZ?cBL@-k12peOo{jHC%Xqdt)8P zA_XCUP$4d#g2y_3`zTjkV5ekT@NQ&jUdjXVr^p0o71*^(Ig#0nEUvY)J0}*v1Pe;C z>1Y0`|L)y8kNw@GDk*p$IeFp2g`t!415Cr3Z01f6BMv9xMcG)vn^bKI>z*%LwSV_+ zqk@#ZTXTPZkGt13*HA5X-mz;}utTAFt5&W0O+QxZ0spWZAL^#sJNf?%;Aw6wm*vrRI} zo3Y;|JM{1WoB01~^lhW-3Pm+*^Xkg%y%f)Xv}?EX;>Et?-Nn=>**_{wnw;(1w{L3w zAfIOTQB7tpK3iP2t+Gsn^9ugaS=t%KBDUjs5V$IBl)A~aXOI>90s=IIR|kEFxdY2E zNN)CWj`SuL#vIj%7YmvX_jja(J7^4F_3reV5BJ3q2BcGqi~Q{K-YDJW^E7VMs9TR7 zb@bJetM7l;_8Q#L-QYlp&w!C5Qwmgc3!PNl$A_;_pK-pCqHEW#Cr7NU<7u$t(c(p* z82ZkaH>_U0nrd_$k+K7L53Tzb;`Y2d8@;G3uCGoXqcHGPe?+OU)jYr?#w&N-sLU~M z_gN@wuTBZtAyXX8!E&(Ao;+vl$HI6K$}LCL_W1jszzt=qWXgx{DBwTJq6dkItwEZv z8uvH#Z?|Be%$n2Gu+Ei)H4)D*RA8bN;w5ZviIJ?(=3+G#-B!w_bn0I zx^G2AM(yi)+ce`w^JoSMC0#kquu>(&h&e{M+P`Wn~cV>UU5dMAH(UY~mH z!SWE%y6@|%`;U1*7J-y`mp%VF^+r=eG7r!|=IX$;jkab{ z4;?vjq29}84eAn;I4PoR;OTx%C)#B-2&`Q5VEnlRqZ>=7i%Stj-J4r$d$KpALi!(C zoV_=VWKB!7t<}c(pBq_;y&TB1#(ruh*E%!Qeijx<@M;Q#>GUDuc}4ck_rC7Fe}4<` z65ew|sdVbpS+jO+`QwPXzXzSF#A+L*b!H80$@E)F4yB0(iQg^_G3|Xg)%ba}m6ZC7 zjCy90h+Qk*X)ejRW)`(isyiyJXG^pgFtZ|9Ql4;#HwOo6U7Y4~fzbFa!l||U_y;6X z8Gn^q+{TO>~0TmPx$M~1Q+EM%yvO+TZ}(wdX?&(63k6d;dUPhBu&}8mfE7Tyuyp@ck^xlV*%{3FtI=~OewYvRY~BAzsi&zl=>P62SE-=2178|Vm8!OML; zovF*;wdcWj3?Qw(4h8RGMlUk^y`pKUQD)Hl9cYpGn>JEbrj?r>;kZ`#SdzB+?EtDW zie{>9@ty9byKoM*EzK4i#zwqV33xHacc=L!z5SsY-`s5fQg(`Oe296#gshKe-G+2@ z|NP`cr&CiW{PkBIu8D0-~hioY+lU`v_xC*O6yYN$lFS z3&thR0-=C?CB{9T{Nlphy0%}QUG`{|oU^W8g9aaSbNx=H?j>3qT4