From 27a8d1c820030017a3f91d5e142d958238bd6ba3 Mon Sep 17 00:00:00 2001 From: dj44vuri <julian.sagebiel@idiv.de> Date: Sun, 26 Nov 2023 23:26:59 +0100 Subject: [PATCH] two new examples added --- Projects/Tofu/Designs/tofu_main_swe_tofu.xlsx | Bin 0 -> 53102 bytes Projects/Tofu/parameters_tofu.R | 39 ++ Projects/Tofu/simulation_SWE_after_pilot.Rmd | 332 ++++++++++++++++++ .../feedadditives/Designs/Design 10112023.ngd | 44 +++ Projects/feedadditives/parameters_feedadd.R | 25 ++ generatemd.R | 4 +- 6 files changed, 442 insertions(+), 2 deletions(-) create mode 100644 Projects/Tofu/Designs/tofu_main_swe_tofu.xlsx create mode 100644 Projects/Tofu/parameters_tofu.R create mode 100644 Projects/Tofu/simulation_SWE_after_pilot.Rmd create mode 100644 Projects/feedadditives/Designs/Design 10112023.ngd create mode 100644 Projects/feedadditives/parameters_feedadd.R diff --git a/Projects/Tofu/Designs/tofu_main_swe_tofu.xlsx b/Projects/Tofu/Designs/tofu_main_swe_tofu.xlsx new file mode 100644 index 0000000000000000000000000000000000000000..0a765f9df0df22fd2797425d8f1934dbcf03bd21 GIT binary patch literal 53102 zcmWIWW@Zs#U|`^2U|>*Wh}^v=KZKEiL7bI=L70JoA=)`VuOv0EBtE3FAhkGFuOc_+ z=#<lWhYSQ<zyIwzaa+r5!-`d12Nw0Gc8C08TxY8&)|Xq!##JA$&UhtEm%s6f@0yBt zcZ!d1vXNPzHCtD4(Z_WOwfB-GEuL;I+}yO+cmAj7X^Ni9mwx5r+Qydh_~O-qU3?9Z zS!EF$x-{-^Hm+)lQu@vqvucXnwMm!q=X5GujXCGQQ?*xm;U49JDZ!5qx*cZMczP)4 z+YYbWu@gANBFpY`N{1ax;kc&Zxl(VM-~R<qq?!EZu4U5}{$Ry9<HytsA^+HNe#Kon z>9SaH_6d*koBXsE?w#c#AE07;|NSR#pZIBRj-Mi~tZa1__Gmumys5>5t!(D^O;h^| ztZJ&=e7{dTGF8+go4exM3JE2TUQvdt)5@2zpP2Gaxa9u7n}4o7WB73R_|M(nKYR1v zUbAP~G1-IpM>TFe{j+jE`^}&8XC!~{x7;vE@Qe1*FEX)z*y_ICdnWL9x7jZi)JXU= z-}BZ>1_lOiCI$v>1_p-sqSTyXeLWDlcG5w=!v;JK?`vG2{GOHP!Le}ZZiOXx)lBRT zsN_j+?OlAN`Q7`@2@X8Hne2UqZ=R>0etp)a_f>DuHyIwC)lwp?r_Z&YJ8=46W+A^8 zH(R?hyHUXsp4rdR{`x3;tvtE6KT_1=1Y6a_9}7yi+4|1V?h`rXe>&4i^Zb>U9c2gl zzN}jJzU0@*hG#7+x1DpiwM%uyU&BkLIeF3fFDJfDT%?d5_2-mEijQn(W960ibK1h6 zEz}CQ@H{r``tCor%jbDqSh!%@oD9L>m41&-N9`|_xgE>-{CtnM)As*6WoETDDVIlb zCjT+Iq;7tph5<E1wk$ra_kn?d!Ht=LL5hKap&~~g5-8>QMcGOD`Pra^2@asWlg?%x zHV|leU(2=ed(^he(}I+q2Qls3_L7++J9p-dNf+|!|9lhW^p)jU#(sJ73{mUP8_!gI zHvRI5S-JF{rp6{2mR+}IpWDiH@7wL$?gB?{N#3n)=W<iI5n}#Y{Eg{y*-fzu6MbYf zH=KFs>9tCiZ<dk#g{AY49d|6Av~jBAb@sc9bUy!aJhqZOM0{RSy2z2anhtM{H`IUp z=6m_-k2kNsGs?zyPMB{ecS5Pi=;}t-umZV>Hg>MlVjSbzHUCR4o6{Vr%Xv6cyxy_z zlgS+06L%-^97&#ktF+DH<<1|!T_3OQ%ZRgF_R3##*<*Q$%E14xq*24EI{9*WI3oi? z7!v~nKO~%xL#Q|KbY8Q8z@G2gYxbHqu&(H2n)oY7IW(6)v#9K&&>9Kpo6A({<7ao6 z?fjZ@z9PLWeQ%6xwCBEy5eq{ZPwKEWt-2JmqLy*~?B}<w6|NZ`oq2_e+hSpBXys;^ zEtWNzwF+_Ta~OMz?<#fWTunIr#;1J6AMdrkr{}b%$=^Ju@OYbJ$=aR_ndgD0Tbec> zKJ@3f^YyJOJ_Lr>cAuK-%j5Oz&zEh2b6bldxjSzcD;(POu6JwQo5d_jGfVwrES|6~ zk$c6sp5t=)_q`JsQwrGve@cHkEYN5Adg<*S7uCx5>&*(9)qW~qtEj~d`@n`LJA}Jl zFJzeY>+C$WDJc$ZANL90F`akf@lXBlO-c86Fg>f3WB>f(Q0<-jB|@fckIxwX_il{* z{Zd#yRq@Z!HI+hN!nZEsI-T^aWZ?pFE~~ORw?cfv>}oCCvS+28l<~Oh{qw=j@~^9I z-Z1J|_T(u)Y8)Ia%cxtx#K5qeje!AFuvX;g7iT0ErKSXz6lLb67lV>&X=rrbZ3BTn z?={7K`6)#gHEr1-%d#<ALUYnJzgy=UCWQ%UcZGPK>i_+PXQzYSwW(LPFdnVroHJLJ zuYCHIFIn;2Hx`-2oQyQElV)u&3YDCGRHvfmM@5xn!$zMpAC1n<4;s{V%-TO&D!DOe z=}gw8ThF!|m}GXE8mF&v^qd&;a#m#2D<zX&w`ntlE?!gWaa-KIAuWb&!m6-akp(~P z8p1x!*S+`Uz!RGp{=D=3-`4Dn@^QAF-*$IXyV&WouP)!c{OjLt>6>?S-mkm9`Zwdg zKd+opH05_INZY^NoA87=S>8gjVWPkBr%X4yvnQ5I{x(UTpZ4Z*kL@h(G)bL`tW2R9 zb~_&GUY?`M6l>!=IqP0YWZsL9H-aJ-XHGd25EY@@(59rDx$A(|2GvU(GE;S(+}L6| z*^@RL(YPq&xLWqaHot)8Q~f7*<Oh77wPYqk(IIY5lZxXrym(&eewKV4Vy=0~aObbZ z6_<Q>ZkxL?X0_BY^}pBmnJy@4+q<Lk>^(!~We2A`w>iqQ<IUl<_Z}^1xw~yI%f&KP z|I-V@&Y9lT30^g~X6=N3(ff^eyv~~{z|K{`Y^~=Ke?hshK;_Qz8xgm-Oq~vMU0O7I zO;%oCs-m&N%j5~SRCJe1skywk9lV)Y<%VT0i#L1j|6{JLBA2zYwT)u7It6e?pNOAl znl>wTr{t=Q%@b8~3^um#vVAkY>BfEYw7%B<QzvyJtgQD&iyAkd_A}6D>s!AxMT|xB z;=&c1PrKjOyKpJ_`rL!U_gBd4R8~&&@_rSv`ekHxYND0r^s3ujr>*5fJh!oaR=d1Z zPg^zQpzqNIZmTnFVvcF8{UGxHqL%9R*!?dY#UAaCj5qP+h~=Nr<JKxIGyiZ}<?BTk zPjpx9+&EpjdZ)1d_dDHj78xG(vwr5Irf8`>QYZTv85ok885npWrAl#0WlkzMDKDKI zn0?29=lK2anmzY#-;rIjM%h!Ii@i0Mr>b@DuO!v<4NnuakEhPBE1IKLZt%i+zRrKu zzbesXKT8V3<|w*6-Z<&@?_QCH-s^0;z4pm|E5E;_V@ZJaYe6kbrz!qlryE_$c=(mY zp@~yMyHe@ZesR|vzXu_EzS}BPN6F7oy4SX@qn2T5vRBcz8@Xout-)-mvD@S{yN>aD z+&Q9iCH0T?xl@}qO?8bmyZ1|><Ku*;2*;-i?l#)leFs%GhMv1&czU;rkL=mPa|`%A zb&6UuJLjaViIuu>Qs+wQ`VNH$^URaB@h6{uI#FLadD#yBMY|Uqd;9Kn=%FnIcNg=| z(*OVKMdT`PZrwsBgS(n*pZr#`?A(*=EPb^y#7yzNaa!|6-iY6`)V?g=n;0L}*nHc% z>f+QVfy*!NpMB>-PMXg>m4BwXo9%l=r+=Haf5CDm^DnNVnfIptnL2Mutl%RH32(Ol zmko@+roJuxG|B4I>Afot@$jyd(vcNQR}KHIo?*P>xazCzmesQshy8Y+b>&;)tAxLD zU+qz2&i=-G#f_{C48hWn=3GUNeo01ZZYqQ{1QkVZBYpcXyGiW(tnIe1`SB6Ho1&3h zPO4ovbSTR_pm5F&;p!7NlF#~-om|_Lywoc+^Q!mKdsQoLDK%c<{%d+WXluwL{V!kL zCw)xYl#;kw%#p7%{m-8J@^!~|7JYR*b3svX(x)Z4meaGulseBd&uYnieQrVh##F_X z1wskOgfDGd<`NaYT`%U~Y%L{2M~8VQeD7>iT6)R%_B6Yk<jT(FDvj03dP!>%G@Yyk zesJ74c4po6EMMzmCL05qk{?u6n=G<(Z0=dQ)in8D@`tG!Gjo6V$o$ypGFRz1TeAA^ ziNQBhnPii$ZFBHD(9iayag)H|po6{JLQmUI+go~cp0MH8rm)4YpU<DPPKD#x*ExPI zhb<RgE;rYEf9B_P?tiCu#I3*n*0Q+HlHs>q#4DY|yLG2uo~g)M@Th2J`Sxd5?`Xw6 zz8`TYzr;#5e{oa3&WCqG`$DsK-TumTFZRLS)W6oZ@8)f8`Y!h4nb5xa%Ed3kK3Zxd znMmp{++1-^*|@#OVs~E4c7}ViR1YwvX5S3Gpq>@F?&uFY@mhu_2ak5CdH4n!9kP3{ z%R-K;%b2lZa>84`jZaT6UF#z@Y2SjCS<&*#CA5+lE5H6Z9)5oQqVV(cW$bIp{(PAJ z*uB3`{$Gy5p1#Qe>NP+Air3fw*8KJQWx2k5ed;l9+5f-&=l}iodV1P|OO5*b>OLLc ztggTB=b!WPa?%fOGx3-HXgKDwDEbM<+s_g?#vBUHg=SrjEt#!M?hPgk)mp0<@_h}B zA6hMWwP!bzR0C&4ZqfcluN_s?x}+M+7wurW^K?$>4ue_6%9EGQHuiMdTwr*y@VA8g ztdw&9mTl&qij!A%rkX2>ntqpXk#1K}<VyYiV}=0pO<OPRy1PsFJ($!lvdcY?J$SJ< zTgi$uF4r@83TMs|X$c6g*<3ANYGJ(bzG8c*iU>=Z+O(F?I~BbCB7MD5k7v9*F+0Ll zO_hJyq4gFqS{)@G67G8&JTe#U)LkZ8nxF8(<=BSzn=)fuxm=H{J@GxWv+vHs$9XM} zudFckRg(QZQO3C|Y0L4cft+z`C*<%J9yocb;QmsVxog+FzhGdb?opPmbZTAf%l(GF ze%Z%*oXr1S=t=*OTVXQO>ro1KTx&JE^n1Q}D%C7Sf1U6Dy-<Gt_`cn%rd@x#t(z_T zb#dW2{WwdNevM@ItWzg^XKqj9^kT2QaYo$oVvU-=@(!J&p#=*M{4^DcsdWl15qKza z#A0Vlz_c~aok4mJRT9=dnJcfZDKY<j<z{1r1siWYb=-LO!uhshv!7a1b}eu_mJplR zdpckJjAlc&SFz1zllLK++n)KZU~Adjy#MkV@kNTs2Fi}lzBe7z(RnVXB>AvmieXVt zn?<kb-D_2QubX5__imrO&2K_U;XAIvdzLOMw~44$@05uQie0FD;neDCu}_7!wz(O2 z9!#CXWIxkX;mn~Wxj$E#Jl%S<=2Aj&l-1U$Y~e0ToL0%|EN;@;Z(3bCN$1P@>bz&| zpA8RO=)HD%(Z!;F&bn_VFH~JD`0;Aj+@qWJ99Qq^Xx6ZBayAqEE%BwS`wiEMmUOq; zZ2|!c^JZ3DJ*B_;!4vISo2(4Hm1DjwnI`s|X<q#K(i7pcs=jgEy?<fKG0Rzx_U>Pn zGLt(wTloZkhyfc%&nB7N=$I89&f;FmufDVL@a@#O!0i~BqNO-(Z~nOw)%Navp&vA6 z#Yn8Il~QtCczLCVz^3J1*71kWtbBd`y5!<b=g(Z7n)%UMrM-oB)2rOHZ_ygFN+(|u zw0xIxAt$($_wc{}6`u;jee9O^vvWJXi%b&?&;E65al+%bAr_ApZ_C@Z<K5c7D|{b* zXxyZqaVO+Qj^Ta2A9gPCXVtw}Ix9AYPP<=I%we_N<groM|928wIUKWH_na^Oab?ck zh=*@3AKklivu*SGy_aKi^{1~`E1UA^>Lq#qnhN=p_Xiv^4^7dzr)Jx*IIaBKbjA-G z#g{~j|6@cg;Wz#L^LqYP28Mg?Yzz{Rx*62sFV09!Eh*LqQ{ZC0=Wq1oH+#zFa^CyC zKH$YZZ|(Cl`EHADom?@Qy@BOd$@)p#b6syvat^wA_5Z%DUkWdI^?Y>sa`ztpvsuD5 z$B+A*o5=7Z|Nr&>ulN6(zW>Mj`hV^9e{S#pzq|hX^yUAbUXTBGbpHQ$`Tzde|9N}= z@ACTZ)BnH!fBpRbzdxtff9(JN<$C>h`G3Fm|M`3W@6Y>xpZ@=T{{N5Z|6lL_zyI&; z{~y2q|4{$uas9u~^Z&op|7Usd{@=&@f1m#U$N&GA`~S|y|GjGe?`-|Q&-H&_*Z+zC z_iz8-pYQ)wKC%Dz*Z$A-{r{HN|J(on!~B1*<p1Bc|MPSG|HJ$LEdKv1{@<%Vudn}I z|L=AD|F`%5K9!gJch>&T|H+Tv|9=c}`_IoGzSn=Q|99Q~|6TjP6Q9Wc`xF1~*ZhB# zzyAN2|Nq<mzmMbp{jUEl|L@!U<8|@p<Nq9<{~zqEw_oyay)OEBdjGHQ|G&j&{jdM` z?fn18^8c^?|JPmr>-GP?|9?LJ|LOVupUdt4x7YvU{_`~c&%Q1H|GxkK@Ba^5i@(Od zzW$#7Z|a}h^Xq4Z{cYd>gTMY?fBoP3f79!K)&KeXU*i9t?f+le|G#XnzEk#U+x~yc z_y7C;|C{{3fAxQk+W+zXRR1Tx{_pSgrN7LVAGi1anqPPM|J(T=?f);nAN>FQ|9{u_ z|6Tq6-}$Tm(&T?%-Y;ie+rRh!*Bd{pbADV;p8kI0m-WZj-A}gvU-jmj?YE!TlfS>; z_%XiUcHi-3^|s&s$j$!$^~US!oS(bCz2l#}pKbSkpPy21{+YbKp74G9&3$aU?fw2r zz4>SLHv2>Uzsi5F=5MwC$87QYKf}4-s%4XZUcUJ^@?QFP`NIjBo4y>sbpGbwc9UB- zUcUKvcvf!P(PL@S(dXXTpSpT3ar2g_TUEc7o;td^h%r<z?Em+p_u8kadDrW`<NvyI z2e(vi+sjAq&fZLm$Z5+-K7I7=_0_$1b=FSXHm~kv(br8ej??Cr%Nx93{aZs?$2Yus zd&IAuw{{;r<`=HreRf04+q+50r>~w$dcGktN;lW`@1&>CQg=4$nEzJ%vv22V86979 z?bB!fdfJ&C&M@{d2Wh#%bMxH7J8N{z4yR>q-SYK|_q5&9pYm<G^ZWfu`;?u#&)+|I zsbu@6J8#=(EM0uo;AYn5-dLTp<=Ut3-u=A#EZmldsg(^n{NLw4xqa#2r7s)eb@^AH zI$XE)2G8WP2Df>-&rZnAsl1V2BfT5qte~eaLo1uK>OU0!m;6@ZzomKU%a{LZ?EQaB zroGH~`K4sO>E^ZDB4*7G6RQ?Qc>mL<O?w=s+5g!2n&+PW_bWCLv*NzxRo>iv^jM$% zZ^>^(_4{rf)7uuc?P}NA4N<!wHcfi^h5OgFcm*Gxl0$EH{O&LPVDVGMUasCo&UT*7 zyz)EtlPc%f>dCF&_vBgO*BzYyu66!ty}a)7-s_q20t@}t{d!RN=|}qUzp+~$|1<mK zxlHzc)xw-FkG#|3H@*3Bk2|$EUoO0V_JjO|2iPs<dakH_aEUv0_ltv%pM2~;F5mCo zzkKrBJI^QQ{uA&2eA4WF>hZsmXZPQ;Stn<H@V--{X~C`u0pj1&E8g#T?_}$zp3h(L z=VQf}3aguLHFlBj6F*w~RyjI7ruGZeX`FkSE>A7FCH?B3=kJu;?+-pM*_8XU|M>jL zx!c_3ZQn-x39_DiRQ=?>l#dl(D?ab2y2Db>&HC=(OE;&RX}|V9-WT}QV(y;if=z4x z9e;dLeWuMk+c=wZ>iH@adw!|RK7RGF+K*`U4Z?oqmone)d>Z)Z-M$p|{>vw;xzul8 zsQ0@U`oUAJKYf16$BQ@5FG13xVIbkhFLcW6`Fq*#XV)BmnNt03{^aeFH775s&#d|K zX0dp@(QlQb?(6No!@Qx>etB}qE%qmq-}B$wzFltl<Wk-`n}2h1&fEQZk$L=6&bzbh z^HqLNN(C`_Ir>*0_x-RvAi-%)oXVAZE5Gk&e;czaw)^tQT@2uO^uO2hK{Z`|Z^az{ zdlM`omnVPqGz07Tp1<;e*PQQiJO6+E=5Ki6w`=<CY<@ZQ{(75#6*g5iZ+GM`I(Ki{ z2aD>MvODI1b!yNc(Ac1QPF*PL{q^7d`!>IQ%-(<b=4~Ol`kF}>-Tl+EP0rYze6)CT zzn=B9`v>;wJ_IKkv94+Ab39knF32xGzb#ST#xMQ$%>NaCETdm1*i=U!nm(r%lw|r( z_x*maA-}gW#x}2Me=Ae;qZ-A9-_FIfZ2J56L7`Q3{2}-HbK54$*G@Y5_@;VfP06>5 zNBs-FUWB>pgX+0(p;Kjl-^+e~yYBeQ9P7{XH*Z(1`SN11djE3}tK^$%xBsSJFU~%G zRVe&tt?&(@=h82?@SnJ)vwK;>2aCCT^&fsb{`lhKmpQ$Yf9B}69)0|#@Mhuv3(wOx z!Q2%8s78HZPWAtLh4<f8%{g~(ngz)HGPRRVx}TTRxAXjFteL)hN)jxHxBDe{&yinp z`QC=^fcg*n?tZnHn|u3z@iy7t&vWiJTvVTV?i$xji08YHuRf0C6MwKzHpg}tRDi9y z-S)S5Ti>A<lYQjQ%jxe+`L^T7%r$mZkW~DdW!sU;_Y;4Y_SfHHe<I$u`R!Ivw3Pj| zhsDc&k4HIsd(jfVz8YKfqfgu|yKc<;Rl70w&x68Go^_DOueA9+Q%5enzx()V{hbic zmokMvDpH#Jt?^&mt{bXc&BtFpitX%<*LnEy#bh7x*D`19sy<r8@6mp!pAPqazyq(C zb1GNtF8n`IxbfV<2Zg&T!XG^n_x_=>*xx92O|ijmm2Uq@e|=8%rLNsqt_scg7wx*Z zIk)nES$E{9-XaMbKk?f${#X3j8KIxB=Rh6TZinP=7rXUm+Ukk#j!Z^Ge}SrAxKNh- zi`{p-@7^$|n4@=Z+J_&HKfd@_vN2lmS@-?qhlSsE#H>-@{A*?otWcV^<Bo^M?vjVK z8~^@-I?g=i$E3ynMn7NNe0)^>WS!d6B=Pz+b|+zZGp6(O+>%?rOK&=VV}HLpw*T_c z+d|^|Sr3B>j4~bbDfinR7VfKvev}kD?W+1ssD~PL3o9oCY<pS%e#iSVTd;p())asE zcyV$6Jezs?*KS4Bpr_3pLemRW<@Wk*>{%p!>f5d6<0l`>_uI?$%lTVRYcJ0Cu9#!1 zC;s~68su!Lr!RDBoohvWPPu)>9KCIZ|9Ad4DkRr$ck;={Hy<Are*3H^w;tw`d*0mZ zj(in=v4#JHep~qNo`eq;v3qqlK919nw~brtZcu6STcx}IoZYS;oa<o8jGJ}ak*_XJ zZ<+rtK5{jRz4`ddSGk<>e?MM)d^L*w^3`{nXVgwQs(#XL*PDoK4%N|zz}ZY!<c837 zW6rJLU;jP)^Xj`f5;lI}x4-_6$vymUXI1ne^_gXNq&L}BeN^3i{OIGWk0E&ztO=Y3 z?gxh7-oP)neD!wKb#=n+{c^QmO1^G8ay9?a2hUCVk$;}TGXb|NyIqXuioNmg*M44| zdxO3I^3~fy;n{38>&ki|`377v&G-CRfk>YT-ZA+qSI++Z{_9@h{kK(f)`1dDG01`_ z_Q}UL9|xziI<?iXO3A8Gw{Y)dzuygW7ABo-IQY2a>$aBT^H+c6aFe%<eY<bsP4&oi z*FK%F+x6pg<X<01IlDa|!8_)=+D`k*2mZes?gVVRE@$HxzCYrSd;Pk)jV=C0l{VFT zeoVUg_)$7a;dEX7<(BWSm8O8o*vZH3!*@$OeEhKj<TiPQhl|DgH_uPl`O8_QW)f1E ziF1n9x$0jzeK!LXJLZKIe|G+2R$X^*I=Bk>G3n-X9qVc7cJI9_z=bUTWTyB>-?&@4 z8vkj(O0)iU@NvoKeJ#i5&;H8i2C6aprhYoD`iEt6zn*z@$gx+*<zSi8!ZoEI_T4n) zYCe8a{pZfe{s)g2AD^B3UR~a{de0A)o9CDOI2~CFbsVU$-ERCcrTW$PvwDFNHhy*| zAq@ca{^zsHRc+46iT8)c7r~0Nbw}PhSsk4yyf;AnyQ+1?oV<HeKS1*G=GVs0X2<J* z^75TEk2fOa<*7SjJy-12e?MLR?H*7rGtU<Mb7znCLr@9uxm&g7<k819o6kA!{CQLe zo}gXXgYS8)s9o^C^Kn|aACKJf*}3oQ^X^UiaPjfk+d}5AXP<$j%%^GFmmm@v>$W4h zIh}RwU+omw_ukD}Z(A`ZZ`<qtJAYmkva9;CbNTGpiw3`69G!k<U&=?%uzquJZ3{{T zcjQ#A+`ISv+|RSy{x%<fnHK$7zTe*N<cr12XMg2O12-Mqr^}vpp8p+I6gKJ?R!<73 zKCI_oTc&vMLE)~O*Zw?DyW4R7>~A4(zMiaOXZ3rgNPqbJ6iDE4v+g_c*41h6o%-PS zvvc3G_g_9+_RQXVf9K)H7ayN|KJTxOSpR1Kpqi6c-63tW_lKt(oVMdX)7)>2l`nSR z7M5k}zkK($P<b|64X6?L?wZw%_t)wpY|hm+r5^w5JRg=7oEmit?I#A59X#+m@45z$ z-16OTf7QRM+uCyc@!i`(pjt{U-TvUCkFO>}i`%&ertO&Txx98k`MLLHy@wtY?)rUI zsC>81qvG%TTKtXf{#E0fpYqXjQ@>vM{?5Zt*X1ufFnve7hlWK|((~HguM2nG-uD0b zyZ&iAs^br-&$O9WJL#x<x_s;(mT*MFA^uUF>cVf0_a^R*Hc$9qaW^;nAIOU1i;rKv zE0=C~f3N<-!f(Pl<-d>hUq!a2UH{^G%lqsp&)?lOyWM>JWu9*nEPi(1wJ^H>SC0?W z?m{l6_BuWGbGmty<Na=EE-kO<Z+rXLeSXR(Pf&z{eN(@!1r}(a@}qu&U-kZbzun7s zcO`*bcm03ycDaMc-TTY)z3=?JDFjY(`jIu0kW0hw>YS<HJ@&lc_EO%)uY7;hp~wAT z|Lhj3Q|mr{Q(bf4lUIwODb4SJ*PZw!m+cf{r<7`gVspOk!;jDNe*d`W?q9w?>e0Jg zPLMyIz6<yM0c{NZHB#K5D!1QZg~SW<cXn^2S`R+{@?l#Gq>bLbT~!VguQwlm`mn8K zKC~_r<I3qg-r2Hg?$7GuwcB48R^2%F=lR22&iJA)o@(;8;0$2r`SIfF^UxMuDO2>L zLN!n_bDmxs`uss*RYmlphq;{n$J=w?FK@3GZ10Z)rLw0F<BvRqC7Oi?rd7mvXuLMw zBmd2JTJ!OjiPE4<!~b82Pu|uR(t-e0Q#NnBdaj;_=7HBN>y9Xgd4E~oRrP~eUw+$V zIU7H@|2BIxA1-#6x2?3X2eCGTgBTY7a}P|ZxZ|N-(f|Cg*>{jNzdy;#|E>ISRLJh+ zmk;YdfAJKXukw55l=k|i%l+Za!Qe-Q>I>ifyZI+EzTUE8j@&lG|CK+4^=m<G&iJA_ z5OdNS?`iE@Q^MfMbCG$iorgw*&c_XR|2`<Jx^e8!@yF*6?-sJF`lJf!KDbYBkJovK zC~U(YeN<Wa?bsgkbMn{LO4#`E=X=}a9R3G&+y-ggH2J;Y`s3rPkC1i>IBvRIZm~aE zc;4&{$UD}cPy;m|`2VZ*wSSct>)(9*>BIX+;MUw<#=Y&xXK@Zu=Dmu)xA6XcP*_zJ z4NYa2_jf*inac^%@)Ai4sEJ;wwD65x+xNcTdDkC*`LbFS(tNbdy>4&+SB(#~q{*Lq zU|PjIPxT7H<CoVdOW64N=X=}a9{yJeZU@Z!=PCnE`p2KXj6d=amY4$)yloz*oDS9x z(OLYWb<^2zpdRP*m*L(vd-Weq_VEX0uDF_$Z@i~1w~s#b9Nu}<|F~yTz=x&_%X2@l z_fI}<FIPYB-i!|yFJIoT@a1JMI0GPOWC`K((wwRO-4*t4uXP@L{N)R%8eDG5=_YS$ z`*-4WaMPh*&%gf3Qds)G##z((9Ml&4v9G)8*0%Z&7Pf1pKYZRI`}?@MWbLGvFU<>o zyf`}j45(&<*8vZ_ZRAz19Io})KOfu?0)@ceBu;r-J5a_y`DU_?o#(fUN8`^SGIN5r z&HN>o@7?|W@8acsR>dDIY;)89J^uJSQ?~V@yT5;S60FDN|9TTjb8wzMsF(Ki;T>=a zc$wP?3Ly7>|7=+wx${sD)iuGIHt|k-4@|50%d{|Q!GFUna9LY*bFHwwe!lmfKWcpP zwtN5F6tX+{D1B0GN#)PYs(+ALH+%e4LVk$u4Ja!<{-Chx=V8^jd(%H$e0=)b_sgf} zSH1*$%2@MXQ5bS1rTwP8uKVlhvcK=AgS#uIrFiY_em!1%eEN1(n|ZdNI4%5W9Es@N z=X73gE!kxI^>^y+9|s?Q*%Wu6|NLLjc)+KMZ^n_BO^DZCRzI%_>F?f|_`zcDpNm5J z`CfZK8DG6WJ-+DU#hcT0>`o#Dv)7)yC7185JI$NA`!{HS;{53<MtOPLeV`@;C|=O3 z*>y)=`#HTmyVu6^#$Dd#<1bTd+vZO%<+B52GX4KbeDOseE%xux{`g7%(Fb@RwLmqF zUnom{f;;EyFdn(()3@@~+5FqHNB85WZ7uqi;NDPe$<2*6pDK}>&+FAWxAK4ack%Q( zP|Ha_o9z#19Hgc)`q9(-M?fvX&!=w->Hh~+pRhK-wj-}yoqn1ZzKNI>+kN@;{EC+! zk3asHy8HLZ$2T7zef;NBx2jDXtZ}yV@=Ng>Li^eG1Q{wFKXJ_A{s$)igVj45>Kpox zhgct7@Q&vr_lKI8g5^;k{&91E$~mC@VcDL^d_7(Fnd?2idoTGcZQWO|cq(#U!+)0W zr+<vNP8d7>;rYmY<?Mvz4u80Qnm-V&2)?Jue(vad#-P0K48L5xjyHIuS2Fz*j#;|- z<730niRbojV}2yQV)nwkMjP!9ycMkXrj%c5ocXc3VaaXlU;lRNvi)+Ctg!OxD7kmp zf~&ja-c^o!|DBjB#X8Ol#5aB4Kj}b?8Q2z2_J7CTGfp~Z_A~zh^M{fUfv0;;Ggi8R zT~yO(Bl#@2X?4-Pw+wNO=1=DqJ*fC2r||Tf;M8+^K@ttil_z>UujPNB{DJe;xr3or z`&F2qMK*o6WjK-#^|H-``^+k`d-_b+e!Y{buoCJl`DeoY0PKO;iURu?&KsN;ikYv! zSD9%KQ~oLImj^_D+B#f%Yv=U;aS_j#zh)mwRD@pcF)#A*T=MxH#|OzzIs)?*_Z_ea z@3`L({=m0VUG|5{es+&tcbzWWd(hvwkClC1pNrC5`%g?iIBK}omCv>R##CYVPef^M zy*~%sRdXiaXI7b8>U{6uR{8BVvZ5}Q`){i3(0VRA<-xBhkId&d)U^E&ikTm@ryuOB z{Zp+cPA%fnetJz?T=`RI`KD<H7;IRl&wj@e)^*<4X`W_9zPj@~)=w!N`Uk7MB_Bbe z&(l5Z)XSDjZ<Y6MmSMg$p&A_7FaKQXPj*t8d)L(1`45xLVH@v)$^&02p;6%;cj{%& zrJtn=2iPCHT%mC4=T?P-p_colSbp%-H13kDX#BZGU_ayf1G7K*S#ApQ6A#HN7yGho zbJ!6E8`(}5tN%_{lD(AX+CMSgu_K|*j{C{Ufb9(NjnKI0TQlXd=+fJ3Su;=gUyFuB zUi#iv|AM`*|IC{BKx=Mq{DM5TxQ6uyVD67Q^|I~K&ChE04lb4dZo?Y>bj{L_jCOLy zqEG*52|eKc!1Gz>K=_Hx2nUsioSV-n6?3cbd~iW1-<&CvOJzU_X`bXKoeApBHJo~L z8~Y*gwp_+ib2+<KI&VvP!OLT|{u7?7@2I<Y^T)>-D>a^P-_H=Q;NSLh%7axh5}NZs z@w{9{QggYwmiT&0!Ix_a4H^7FsqKx*jxCSXvt4`+{@%M)V86p%o>K1*b@jSWAa<zB zN@^}Q2>K=N_y4@$#B)|m%pXKz?DY?NTK-(BF(I96-KECKk6x!eP~G9h?$`BPv?BPQ z?%kPlrhL|)q~~^o#X4F&iD6!!%ifZC^$q<8mOn^}|MgSnfbj>Rm^a?`^PyR8;^ZQ` zg(bf|_V?MaEs_5Yj)&9s2cI%uKJI_;D9h){X%9X%zI^OmV0iEkEXTR~o%-2x>Efy> zo41OwJ>dMo|KWUK#FFMv&Gd+=j~=6E>WNc}_AMy6R#?G4z3k3sPzL?;Q9hu&#YXtZ zd9bG|xsI%#P|pm>pih>}S-(PNkCiIeHLTZX^Mx&~uw$Cv)cyI*YrO-X3lA#It?%Z7 zWre;uQ$F)9-Mn6_l3BiC`Qj!G_U@pXKO7&0KUnV2dLF-usg(Okq+%Vy@b@d`y!~?@ zoFE+M$v=<><^3NrG4lftcz)c&5Wnyp8#p6BJ>h9-A7J@fVNbKS^}9;0b@Qh^SORiK zpSHl$1D+pkoaZr>!ZI-^F+48_id%HYPU+a^o1lnd`J=wWEB$NVO5>$)~oO*{X# zU8*~_Jhp?nrR2$yIq@2I?H2VmmTIs)c<G@6ip`&(On+YE!IFx3jng0e>S2P%YwMgT zr}ZbjOPb!Zr`g;5U8R-Qms7ESekU({$5E-)F~3oNg8JbauY$sPweuh`durmuBFlwE zyVoo5Yl>&P@Eo`t=(OZ>`$UI%f}ea2q@URAP}BBP<bn2+o6v$tQggbomw5Whu99!8 zTFl^rYfa%t&UKd>Cl-BtKeN?_`^iaz`9d)^2J_P);aQThBt}nTuC3b_-WrB^O{eqO z#0`S<t`$fqo@c2{JF(ee9?wVB9d&U{6CYKA%6E{Hb}lNqW?#WRJ@3wEk(jB?qMG6n z6CZu{IgtK<`^wn^+-K5Tf56<|U$k>!(eI8ww~9ADU(a{o@<fj%pVQYe$v65Re0qIe zw|%~71zTLx>GigCcG4fX_GrFeqOAP1$o*v*v!lKK+q_wfkn(!Y<YX@;*)^U!e=tWh zzc!b((A<~DxQ_{z$1F9cTikx4bm^w4+WW`x9YxdBAJ`wTKX}yh<Hz8Fzsoj5YaFmf zXjZn=oUX2w{(TE1Cb;zEX7vA;F&27yz~f_0+b-@S=_{fY&U1Z)IBJTr^3zXFPaizp zX1L((j;ESCyq?DwREvncOz`YlZ|!`S^;61<=%#M=XOF?f)<505Q+uYo24&$Tviqaj zlNjbT{Ac*jGX4B%Z41rt*ALEs^GvCC!NhaV5q`KX?8WZB>sZ5<$Mcpy(2}*#w5zjG zQ@M0-lF0mYNG0fP{q8fgBJZ0s<+bF}#j2X&uVZ))1Q$%@o3<oo?P5c@U{R2(`r{Ym zwe4d4v<6;~`kngPbSaYkp=``ANYQBa9hMgoiu<2SeiB)++vyKH=j@CN+<CY2!unmS zvnAOcXzuW0_YSK0Q{fZ3WOGfM4fm1s8UBp>oS^CS(!{Andl#1c?)btMo*fHu&YCIO zGE*KsUjE=$D0^7feE+r|OxL#_f)_twr%K;H&QDqzv-a@^yML`3>6J}3_6NVtJ|{29 z{ZqYTzG7YDF3C?K9r23)5pJ>z*eN0Sa?Mi%2LA)~O{+@^{JPE?gUXk4o*yfXj!b`0 zb-8gWB%OksRJ*8T_tOQJZn`Q!vO+?2|8r2a&v(xA<417Saqwu>H>8U1YwM-A)0*oN zqWZz5u99qp<~&fTs|8{_f9$V#>F4y%&_LS1boa#RMZaAy9qstSX1?1i8I*g1Po1se zmhX^n?k+L^5LEMrsT5RH?-7xQ)|ghGg7Va-<jGHDbv1t-%LC32%=$+^3;qy^nbvgo zfcHmmbsWhBO-c7W-Q!NZ?YZ<;`BU`WvsL1t+QZHu|H0A<J2_#Ymv6FWNInvOAo(n^ zX?4kir4``fx8%vvJ9-*-?Uc-J*eS^Fd=4s>4ldjb^@xns>q@q`rqv}0%b?0rmfcyW zAv<@<CVmh1ca@;#(FVPzY-ym76=%PETtC2^S^q%z14vF-qHO$hTluH&n?1HybN`x= z<lD$t0csK28zjb_Iy*~N7+dpbrsj5QFLvpgylW*A4eBR8Gyh?(VX<SJpJ>PZfd{3) z_StD0u=8&5h4s7If4`6ix9!-|y87k&S9X=~|JcJ~eDwO+fb9<RM1MdM#mf^uJ8c84 zDzXcBz>fFoDEW8xpnKv(zB3~6<=}MMp?@&dbK?B-2n~A!to|RG_u$rq=y`oyi{{j4 z@_zsqrw_P4$iz(j{Qg~|4cN2)!`}JDo%$=TS>DcTSDvj4Dt;G+zvGjxsktlo$w$DR z8QMZ9I{y=*<NXw6?Wc84Pcxd&+bR6r@%*KR;-$u^MeV7YcIR2JXB+#)CBHkoz8`P+ z@^~JoDa#&qN_N`PA2z`Sh6c5u^m9e`;L!tLpv`emd+u-3rO5OrvY-~u<v_<N%<511 zP6yTefix>Vew2LVUGR4gi~7;epyG6jvij3~+duxzn_w@i{3|qFwvn-dOK-o_oy5ql z=gJ>6cg#rG_m0Wt@E=g)paz=#U!L&Xc|X7^+hR{6{{i6-^B?>w<Xj!}?pntTg}SyM ztgq5e1S{0F{s9Nldr;X}=l=A~N%i*+F7;ObD-IE8P59c`KfmeqaU}^gJ3HqZ&UNzL z2{Rzc%_(KZ9r+a#Vy}GU)AQZ_Tyw|0=m*p9tX*&5G>`MsnF*kR?vZyvDBmBU^Pmn( zaNnOvzlEnR-kkdWn!Q8=s9e;s`}py(q2fe7KbMc7?n=ORmg~zsK6u-J+oV>Xyx%>X zzM@3GO002nsx7FTc)giJ$L?bd*Shmc8-5)5%=202;MYhN9$}CLf1-D>ORDw9tK2;- zrSSg2C!I&&4o=>#c2hCbX5;+nuxh_3WrmI4iaYOwPxyXsSe$C>@2)hFuWwpyL*k>) zJ_nyJXOssu1NR`cdU+-vmt4B}-$CmgHcKkzDV*m&Qm;^Flcsg*pm!Cxka(ggu;1w) ztkR7Cv|~Y0S@b#U9X5+AKxtGa=kdpn%18VUJb&=X_4#AtE4mYq>e%>CI~ElcD4b&E zoA%WP)RBm&_|(v<&cDIiW_}|$`x?v_g?0!u-OElsZaZ}}P*X{Ej?QaHcJF@9`2I!Z zfkTtqYh^!zyMPD03o>~S$w9yJPM}rwzLVBF{un%{<XRVBvZ4F<XHYwto!j3gz9D@^ zI-)W-HL?Fw#e$*&xl@PpvLqo5y5|dR67GW1M0PyG^MyZtfZ7z_SRX8gq=}*@i)|kL ziKtA@>&b{`P(NJ->ZV*z7x@rqGr!5(_+JBapbfZ_HmeX(kOs=>PtlVRo5<$}@^z<+ zRYb+3pB$CT9%#wx*n9-1=&5`)BG5=y|EY04UgfUcqSbrn6>M4a{H4Z4@8>_CYuxdE z{&+d4b^ICJsbiVW?}KPoPyA%Fuw>i&dT>jJ<q!L(`iA}k>JN$wLeG6dEw|h&=LPP( zaaZ}$!IPW7mA8bZf3`)>!KbU=ajKvG1!{H)#5bJ|ONP7OC}^Jel)i5tY#wGro(DCH zny0FRi|9|G-&>}t&r5%xx})am&!33u{<-MV+wM!uD-PWRSy6Spaq1&bQ&TUxcYeas zig}IVJJyLlgLZ_hHU-TypOPmp`1~hl`su29O{+^jOp=qg2UQvm9N_JV*X(K$kiLz& zY^3ISYp?Cm9FLl8ALc}!2emYr^^YF!Ki=YEWg&R<d8C4!b4_FAI#}n=>Qm6XYW4pU zT&Ec>@}F@{M500cf&IZyD}&&gd73*UX4ru$@K@6wBqjWUmK316<EP8Z7pAuU7O~ep zv;1NI@mxdaQsc}|Ev}mWAbmWg-UU<n&VVCr&*Z){Q=W4#T`U|puV5D_#uP6dTqJVd z{XuX+txuz?=DhS7>748OHK7&R{Y!gK_*DK?dKtGt<0)%9$kV5NpFg%g-c$9o2Dxnc z7z9h$PfvJP&JD19t+=Pz+Zogx?<l&K{L=waxI2RiU(|d#{gc(=qH9u>?9=N&E#LQ+ z6`uXil|Q`P!4$|Bcj<9un+@wLc=%Xu3W^gC$&;Vxs;O@geeE-75aEDV<v}$$d3)CB z%g#Lq)!7M4_q=(w1UcH`!l&ewiyb*oa}X5Gi(5R>b*F*rr-M%=j~>?tcb;S;H0|K^ znOEiAfL(VZFSLsw*%Em=vE_&64xRlqjZ+_ePFWF+)T0hNrM+g#Y!CM13oGnE?I_TQ z5ja&tY7M5E)?KiYUH$1-ho|2JYrQu19fox8Bc^@Y<aeCo=kf=a`hJ6|f(K|FqCN98 zVr8Z%T>?A*!LLfr+toJm{hXa88+1U<H$JoZK<U~Eu+GPmC41sP9Vq!31>4p<2emRm zRau3P$CBpnUw(khUlHB(Iu_P7R=xXE+(L6d`<{~j#Y-mZvHfb3sjw0%y~O%{evm?4 z%MZ;8{n!uxG$GCL@0T25_1OV0(1_05ntTD~T6aEh|M>k<25_t5o2`}c5&Z+kXEw9E zuV#h~Bv^e2s^gzf^4sG-cO}mkSxazFO6NZVsNWK5mEELs!1x0bQgfC)?bLh5NqO=w zWzR^wur7L_x#M1GT{~!4W6w;XJ?6m*=b0*v&Maqne^vx;g+VaLifhpo?C+(E9&mp8 z!*cz+cfG<X!{ESk?>{txJdaY(GV4G6eqjHf-OpX8+}3TcGraN9-^Hqa%NtOaN_~;T zA5f(1oz4WQ5|Nrn2fRU^Ke%a1vcr<k>8IsXU99%&us(7>Vb3^U@jTN<$w%C0ieWuA zP;LYb-?{(lRWiI#{si1Y-JEuu9W<16ul&iuudC1LdpA~!?N|q^QamA^fA><(>fXWM zpcv8WDEgOi^TDb)--8ys<AHX*zS~bY3~l{uvd5iz&lq%0UROHicVYdOoPxj0K8JCA z097c@`=_z)0}aVB)wDtyRsX}@xrCj1DXv-G^X1aqKQ}iWxIA&i^mCv63wFM&`NLA0 zequAr^x5Z*Pk(Sp@&~-UTr=e|=hDsU^PzozJ@dOKt+?kuD7yU9M)-r|4^YFL=`+fx z!nqZ9@^nrZpE@b8%MDH<1wnrQbZsZ3bAg78(8s0xPQ7fpbdYbp#Duy#;Jp23lgI;5 zgUTju9*8Mf!Kv3j9nq^*f4b%S$Md^A#NWRcd<pVIaKYZ!2fRLlMmbrpueaI|R0HX0 zUq$Z0>{^IA=(0*=g|*W>&PVAFxFH>RD`r?nzUax4Iqz27x%+qfA3nXiJ3uX~;)Er0 z4tRX5R6Ah|8yXC?Gzfw;lq`*c{P;uOJ>2*E{=ucv-@y%s4}p+z7s+R!kp)oq5mwq; zehT#S4>|XX9X3Gw=bapA*yYptgQ1om0&7r4QhsVI7uVYUoBdZ|=JgCvtJ!JE=iSxd z47$UNt*)H^3DTh8pXgo2lA6y!m9tWQRVQQSaa(`)&bVEhP4(<RW#03d4d*AMH~!%H ztce)tvoZ=gr#>amT<OTZJM2qjzuU;Fx>ywiox8>;saVIQ$FFGzPT1ZBl^!2I!v;qz zHJ=-MZI>4O^!dllb&Q~vNsGsG^B<s09Z=4<u6*wO2h5*!4yO8cK?excWi2(I@9nMs z8vXzO^9~P4<N4nFTMQpTLr$R5%CGCWWCdt=WiqVSIk{-xqLOPbAOmGAm30dLm_I%S zW$9z*4~U<boM?wWn6htS(KQE%a*>#+%*vq7U`b4)2D@7qdZ*mO@_yjXT9-SgH<#+M z!N=^F)t~0b&W8_gzFUgydUmhv{T|`v+q@mXo%rYeUX0+LV}fT_e|*DpPyy{7RD;xw zuUk~|yQ41Q)?{a}2KMZrnm-aT@2%>CxQ>9^k_EBa@PWHMQ$BMq{hjP0ufMhwJRWM9 zty|OdI$U;*=f|Ir;i(@&G4eqNygq_P<n~Y9Jz;9mPuEkiYnngs>AeLdQjO<7*<bY? z_o;qbgJt%?-S5=XmQz2S)4}6qT<cCXP6Rc2nC#f6SD*Wg+6MF3c|>ep-Q)V66`yN1 zwjF@fR|Y{f^8_C)fAGjA9+VS|uIM6lWtE>6Ilj#BJatrxj~QITg9d+C>{zEC_qLzU z3LcnxymUUiruEnv6KJ(t=rOqK6<o0QWrAzhb8y^#lKcT0N9X<^3Y}*O?(3Oyns@4L zUY5+``qw%k1KsSiK~*7WY_4G%tlQpQ{Ss+J=d|q7+wMp9?YUtp1#!@t#UB~%nC35d zClj;()t~AKiSxk82Pw5xE-2Z(ae@4dZ%-yLf~Q*kEt#@-g9p3csiRi<#!tYlHcma? zxd?A~?z9Qqb+_Zf@+r%IGf0B_*w6jF8h?D8vC1R;^*PAU6!()z_?Z2iDW|!Yc7N4U zn!7auI?6e<h*A3}XmGjeTjjI|TCx$E;=9jgwn9qZKhe9iBQ>wbtHk!4V1ITsuN-7a zp#ELau%$nKfFmP5E-|$0{N#msoO=Hd$+c4R`iG!(N^@<0JgNP|9eD1ONQ|BS`ho!P z{036z@aL+}@Nqf)r&VrG%T7Dine2$t2amI{guguC{jpN*0o3!TP1bedQ}VWpEivB_ zm7fS2SO*R9gj;AFeLl14bZ{4F)C_I*;cFLY@MZJCtFvSwQIZe}c8}%`o25TM?RwB~ z>O7<@G!rxuwtKt6e(9LEMN=VhvO(u57id5QTH!tmZdea*0;MdAc@Js=UAmbH8>Wgp zC3;3=emd7W{uvu|z#|@@HWPGoy(DE>jJ(F(ce*dwm)-qo%mD6vmi?#zkFU2@w#wSg z1CLbQd>sWFGY4szuQB(l#=V0}6~BWfQ4-qEU#;SU^e9)&sjufp?NLrIs$E$0&ER@@ z!Q11y{s%)XHw3LK2<kd-%!uBh{I~W_@0uym5t`@Yc1o}D+8ODdz%WnYynn&ZMV~(^ zKbenpz(D<JnbXsR&{Jov1erguRMtUTpvGs?8Rs>MKhVB%7d)H|8j{~`td$O%pK)gm z<P%%^W0&%inF@7m`yA%+eAYa8l=(As08i81?bO?*Q#XB$-#@t2Q$4Tg_41{?{~6@_ zXLc1$JAVK?Ooj;eq9=>*tXm=TQ^S1s?NkO(2kF7$I|tl%Hue`W%Y%A5;JKeii+?~G zA@7$c%Reo1e)&z%a9;kmqlY1BLFYe9_{)Uq^Jnw4XMqPSD;kkP{mJ4x?=)ifTz6mc zJK14YA*7e*s%ZymsP6CrO>)eGba)`W)IF2?YNp&~UfSFoRHk=!$8*r+TVwAd(ELLK zsFt<RIC>qW-``tQy{M?{@)_$LCZG}HU%u?KtMfMOJPxWc4~Var-6%c-X#`vP>9*Ot z`%@MlJe!wQ$$mZU&L@^X&sXT2I@l$0Uj2kU`}JzRuq8i`dNBD<-=!?RQ$MlDQ0eZE zyT5puKY)5{>O1a~l!H3v;HD!=N`A8Vj=x51Pe;i&m0O@urt-V8)0W&>TWIJsPjkl~ zi+7l_N)wCh7ngk7d+hy#OMTVgfsQpEI}dk>%ujE!0cBr*&=hD7eDu`PDCnJdNME1l z&f9lQATa<cf<djWKVIeA`A#s|FzX+D8g7N?!cQ*R530~CG|NH5<e+{wvpT4T6^U62 zV#vh2xB4%a0-wmW)ch{(wY|UNbItbKvCy(^auK`w(dRS46Y8MpjSq;yN6(%219#QB z%$ZW_@C!7k{`bHaeJ4Av?vf2kPx;)>{Q`}XcjzDdx|?qftWP^<%5UbSx81)O9$g15 zd_g0AYeeMjS+DQrTO%@G-T5w44eLJ0gr?P|pm!TP-Q$9{_rFkkvimCw8@Sr%{{R{= zHoJQFp!Y|Rms5k0%6IRb`va_Q-7DVy>=2|5jF|lBbB};M>vZ{=8c_X_aOm4S>zwz{ zhSAUMCHzLga{E1YNeHrkPTYPOG%CGf`a8z%AREv)HvjIm8ErpUE80M_h%JcP>dchm zoJ((4*GxS(OOW|P$rOQ?2Yf0cE+@9w1Q+alT?3tpe!31;R)c!~g>Emusm1G@X)hIr z2JY-n8}g2WMr}St&jppfpb54=Vb4Fq=8Nx4InKKD_UxLMujSym?t$wG$VkbIPYvRV z)$>4=qRp}&&^f0|6K8+gv8bfJr-T3722c>NhCkH-PwgZ`o?_;k`xV+B@s57}96DNJ zX%sAH9+LO|<I=Ua*_S9oT9X#qM?Zg@F~K8!yXaYg{h%I+P)xp+K`^BCY-JQIS3adL zuI0~$<r^ULTFq~5`6H%*=fS`eb3d3q#~_u!Hf8Ax>-+u$_q;Hl3yH_`fss=-m%1ab za(Mk*wE{9Avs|W9vtM0HoIA|L>a`^(PCS;tW(sy&eF(BSf55xo?W+Uc;Fgv3%M;#} zdjhPg8qXho_Zw8-PV`{s4wn1CBrk8zD$f3ykvq6%9#1JZEb;DNT0Lp{r$4Tj-YUOJ zuv>fGkip-5QJmFz`vdL|f(t6WDx>mK7(w%wpQb#51run*M}5jUIl)hCb`1XRtbyk~ zffM2%l~a+c(-`L|`Zw-^wGD&&)=YUWIrX-6P2s_L*P<jEmM?DcU=KTWIf_sD2`DdB zik+COxQ|6#{LFD!T(6n(oOS8ux&N(qe0e+%X_}y4Tli^0Y}a$vA6%v2`L%nR?tZ5} zcb#fIu5>A~8&Y}~ym=$M22^Q+8vo#-C$STo;S+?Bn&+jx+>f8LduzBJG>5X3S^MZ` z!B0K`<)EZ(rx&rLIXN9R>w0P8)K7P6TlZD|zHF$(`}1-X6SS;}J#{%sT=~gN(1^tA zpqjfp$YZ=~raWg}TD?R6|M!pjE|#Z3!{$|&4{ky&%hn*uvQ0tr{6qBOCh1(dnF?*~ z9q_7*@=gW~+}Oo8P5smooL*bQx$gY*1lY_@%2FH9e6!N(^!e8oK}N$PATw0pNw(*p ziJ;fvvXC*opW93HHwDG*_t;hIvVFq-HQ*NdQ(FW5?ARKGIvces-rqMZKXeW}<o+zU z>2&OKcuRE6l<S;JH%F(#7hQk`u2zZ2+cVBz_>SqL<SVfgo1xQS)=y7(Rn85xjGFps zQ~nu9r;0rcG-U%C^U*pAnm=;BD-;v9^ap$xQ2S}A<I8Vi@p1;=88c7V`nz*Z<n#0X zSYrd4?OY`y56x|B(qVnuk|#^|tkamgRSe=F$9eM~tdfb))ZY+<wbKxGD*HSCulOhb zM3i@lyh+t<gtWp6UaSF6dE6@ju_T|xH@)^{4+Af{xCfe8dEIyEt?WLD8T#3Q(EJ8U zBtJk^h4#_UJbOSbOGrn?`sE3)O51>4wJz^>weOz72%e<A7X5?e`g*GmL2(<79tU|F zW!g;t=~tJRYg#|>>4D~jRw$f0$X#P2-~8G){U=Wi<GjRqtUtKm4Vkc0uUk*u>@|M> z;8q`GlB(3i?%>zQLUTmsr?-JxodWR>c7;Kzt)J72rpG^6@<;B&siPehnuZs?KL(AE zOn#?v*He7Y{+9~p!DA$l=2EpfE3}I;nQzXN`^+kL?>+oJe>y}HXmLT<slbROpToR9 z?33qTcd1gF_0_Zowf20_2KN$X?Wf-vJjIoj?#eUQ+CzrY>fh;r>Z3i|=QXb`x#z{X z?owoT24X;I>Zj_4CC!h6k{0BjXar4hJ)gc5(&x6*^Ih`!bWp?e2U>Ht2Z7ss>ar1< z``Q0Zc&Bvft@5wNOk42$BJ-bO5wVwhrVH(<_i2Pys=tqm%z$;xAVXX?ugXDY|Cc{l zH79;5Q%wtGsK@5;AEcfc-<m1+nU~1Sv6e>mz>$6dd&j$$JL;}Ny!>m*qqQQ?wP|6X zu|2(Ot2;`ry#i;tx%VG3gT`&JHRltS{*hC7`Auwp{6;y*U?X_EO?8LYbN_j~>&oZG zPXmoT3B|m(dcTwbGGzrBY^?6TFFnI>f$Ywwpc&hrshXgo_`|(w(8xsmG^Wow2dm#h z*LD5Vz1w+a%4PjY?^5=5lox>V&_UOo!WEkPP*;*v?t?As(sYkI^|IwuWcMTR{6A>o z@N@nL&UL3ivj|uhcBO~#g)Kb-o0fbYzEmby=;fbPQyzTseGX~Kaq8VKv4hUezdYe( zc~rVTbl2X)*LR--m6+hwGwflfPDaTF7kr#o3vCb_t^Ni}OUl|$OI%;BfzAJOtvmIz zIej6j5B@F#HPc`-nI+CoABfr+f~QtJR4)B&PH$daQWwf~WiF)k&ADzqa$z~=oyOc! zP^U*0-0AtVNaTbuc;#52CT9Cs{plBnr{5<2WB(*E0kp<RMd;-LPs^XYb^`XS@5Psb z250vlWPc_00KU?3;^d;LMI{CAa(FzfArt42<wc;?LQLx55@UMVo%4vX?m2N9u{nPQ zU%qLZ&R79n`c?Gly;B<~m4en3CfdQvfPcEsMPPaI7vp)qK}yL^Q2p%wFwP1zfv|jv zaRhWhlFd>?+h;0hC0OF7j*@RO*C1{1n^m%)(JIK|$odGbGn-ka%b%!)4UM1hwmiOb z^WV5;iHS>j)3Z5w4g`Ow%NILR-zV8wl(6*2A5c3SG>5Af3vXjSJ>g{;4_S&KY6zM_ zXzHH&=<)FfT61GfjhsRACZO&lQXBSM_>{c-k4L?4vM&MkKd<lR+XI=%sXXwdRQ$ou z9SIZ9#a-uuCk*|kHO?;|xt>&-`_&56;yoYeuw-)$iydg76x#Ipy9(4A0ws(k%KA@h z++RLWwKD`QcK{8`rJcGtOF{XF|Aczz5`cnO?K7}=NLe;V9<)4y@7c|Hi;NgRtN7}p zyg;L|J|7`-^d0&J=S86trk?I;r+#)`x_R%vY>ZUVdEqBN8UKNXWE##NpW0ajT658` zKA;@3ma6E<(mC-f=4=&!tm$OEUJV+^-<4btdhQs;GP(PwVo!MO{G;^ro8aHZKS||# zK+9$(p8I^>aB1AG?yXA89waS=&46@IeFSe2ubJX4p{72OFKlVd8f61;I}F4I4@)qj zchA8srKdfo-pcO#^PqU+w8IQGjQ=^;na_O+nn6DoB@3JT08L22#z`YJr*nIiHy%IT zZEKwZ86=6&6pxsSt7*G5W}U{}dk<HaT<e0i^Q3EP<_UiCIrw!opV^NeyFf#v@uf(u zn3z0}hS%kDZRbEOc=O>SWJ$4}^h`z2WXX}~4>pM*w+&CrF5P^5=e}mhQod6M+2rNz zo&K=Y@PUSq&Vv_)K%)%QI(*uA>aB9sW$>7aNK7<y0jM+cI{4J#yCNVjgO?!gT?ZSk zo7^{N%4y-H-K^T6dJVKrw$DWgG^lms!DG;rKD#`|l>QpfpcwaL*Cm^u?}4=c!Ly|Y zzpned9x{0YUO-j@P8s(v#h&!ud5mrTzrP=kci1GvtyKrb+@iQ$+7+61=Q(hWR;f=p zxBt-dAL5JVJOwS$1UI*1);!l<<AuI@u~PH8vDWtOKOWh}uiXwZTxlX7pZCtD>7X?U zOHDt3+S{PP=PhiAsYCCbaRFA@mV27L?Lo^UL3P3Fz6t3~KO{k`x)uE!e~5m9G|q$j z)=YUVI`y-@GaqEwR*T*5)Kx1`9y|E-dEtD}`20uKb=lAnz&*jB<!QWAo2M&67Pf#I z|9+ssQpJhq<gK}}EWYwP^|cSOJna^!C%yOO0nd+UeIsA!G$1%pe{DQ<Q&kPT<ODo$ zAfXu^G3n7~p9!$`QPb+8geB1R1u08n=C7FZRsG(<rIL`Dj??zwc?+BQ&980Qy@P7* zGF3vGQ<ckQerjC5aZj0BI>}|e-dcHRU~U61muu>-u55j4J1>1kJ*OUQsQu^m67x@i zdHX$fB`=x&F0Hl$vL;GSKOa13#wyNz2BSSF|Fq2M<u}2fXCI`AuLfoPN{^k9&X1XE zp#A4ROg8f3I+qT9lY=*cQ<mI0rxBa;NATquNAPI7Rm7B{&27g(-RG+72N#LNn>&O2 z0bi=6{ItyR<u|tfDwl4GKZ7=D)t|C~7tKwd%4-MiCQW;=Y12WZRfyZQAp;<ryG<Z{ z^q=NGSf;D{%ui>sbDqcg$oPo(jOU1=S{pPf_Q>?4(%if6L51npb>MczeNc9Xs6Vru zAzl&IQaj<b^KPKkt^1q7i_N&!T{`%0-t4a?ZlE=`?(blWF#3>c&82tbH16Je{5rWS zmemM6+6nD;#_*h(ZR#ub05sFawhS@8?X&akwL_ub{@s5gsq+58rHP<L?=n+eE&m&( zJP>{|zxnm^rS+0VPc#F{TcJY+rzXxUs$N)B-M>C}uBj=g!tpPt_4>H}L!(V_!Ow|T zpGPWyN22OdzCx;~|I^<2rJedK4sHnF&08Vau>3*M^}C6IUC%+?LvQzYk3qHCgScJ0 zO|!$`GrmT_@6@N{z5jT$c3<@7WJt;S=EKhDU~r%BPO05_rq5?U>qKJVLlag;!SB+i zoU7t~{{U1hHLWgs02=|tvJy-GX`SQKHLV|@{d2xGQ>@>-N16XyYWxA>1;|o0*QW`A zr_L<}wU4e>^NB6Jv-Yx~(;v{HURO=X%B@q4(8-2V6DJq#Us$y5BXl_oeAxDLVd^58 zH}6naUVH4kAGp)P$#u!@-Rs|g8#|HCk3o(210FkXA2l!r%_X5HRynt>J*(d3)lP&g z$3K`_-v3<wf%2324WQB6pFg}iil)KG%|YY0^H<EV&Hf6TE124R*4D1B<qy{B-Nkp} zSImj+{>Zm(FKGR)7JGKkyM)C*esKQaDHZQPu0H;)z0;gC<+y&3p4*ehpAvKXAelLD zSNm3Dw~qM=|6vo_r!Ox>Bs#asIRU%&?w)?^L8AFzP<Ecmmjh|y6dsUESPWhY{-^=5 zoF`|>aml5(Z&&Yq{U&V|ID8I1o7etT8k7-EA7YCyZ?utog*?gbZg=W&)1|k6_v$>V zzF{244qB3%&f5#A5$cwIhBc$sRl+B=OH!8DtkaN{o$_e|-)ROL)^Jdn3a(}O_`JYN z=MH#-nmwMdAtcYrI{{W!8c!4X`9LXJOE%|yWy91@8>Sr-iK#OJS0Re~4wSBiO$u9W z3X<EX_EhcCTV<`eTUj6zhGC#QFu&<_GFOD=yyuUbT&=RNu0hOpEVYpXHNNbnY}SCP zse_?b(*HTvmCNRQsJt=RCvC-U2jpQY@VJ&u;I6&Lr%(Hn=sOJ@QjIe|ZBRSG0b1I9 zfW`RedC;zob@DTy>wxB#q%5<E(}>jx7ks&9H+ZFRLp^B3Ca>!`v@MI663Tx$@5wTo z`iUj}*?N1Lz1`nca_OC_{MWcE8N9dvycP~+JCu9csn0E^-rDY~xoZ1XY693j2ifG$ z9|vuyF)d8%h;M+*zVvz{R^`vsJa6o^z5m4dHU7I7rZIpf8xofN_+hxiEBz-pH5tr@ zuXcD38YiuEe7VN_6Q3UM&gYstwAkZL35zYomR9e9x^<tsPTjorAKa}2?RZ%ON_N5@ ze(o>|<lFYWA`Yo}nzGd9u-C-3|NVELfUF-C{viKgezMyl86B_6DEBmIWBK4yc6rzc z`l*SNKkZpuR2F>tH(D#?y`_c5(a$`UTvuXYqnt%4OKjv<+-aVt_cm=hWKd*-26&w9 z;MbQzVoNLL_1m|4u+OeWOgK;cRI{k0x_<rq>DM}-v!rg2g{e2D`x`(;IXCM6g)SWN zbay-Tx%E=y{|B-$y9-Z)yc-594xtl51+ma@|GB+L|5ISxevh5g?;3yR1I?&t?$G$o zBoB5FXs^PDnm?*lPi^3p)6)~~m2(4j?%g|m_PbQ^Ovu88X5n;O(5j!Ov){3*pR1bJ z&n@=yfcr;iAqtumDs_AL?OWVD{rGRJCm{pj4_r^^e+UFEgZpCxp7m&){@~N+9q@Vi zpBn4Ay|#b<o@Dc$89edfq9iM!sh=Iw16pjd{~L5(k?HzYC-}_##7|wy=|?>(55Ai> z-4r?vaNbZkkgsQo_?vehp@Z6|ox5OV($f<jm3ISn*1F6OTJzWd+|4O~PKJZ`7J&M6 z;N3^yZ5@{;PX4qDv=sTq^xxn<KX_7CA2fA<b)!k7=6ZH7_w1lJ@Q{e|L=X1tptucc zpgG}G->&(f7Cg#SyEbU}>w)bF{pz|?3^s^eJ2f#gH9*Ze)Un;yt>7L@!1qnhCqddR zoBtI1crIbqK6-tw;(4K%{jWeHl9rm#^)>$|z5B%t8rO;2CGk=|sOkG(iw`1u!u2Zh zmomsJo@c3H+sAPJc>JPZKExm!J7`!ZDDVA;-?k?p!(X6<RKKgi+pP9ZEsB}t@%;OD z2K|H8+N|(pl^ejz!oJ;@T>fF7EIX*-IVXN16S!wO``mgDv@vk-2-|<7_1XI`K9G1h zcZ$S|--Qo2e}H<kpKn(;_#gZ|zq)C4NnJHRX!TkByGCex<Wl40PoPCu)-_YlMM0M) z<or<w4^F98{GY%CZOnt2&~1BQ?@T(k-jnm*_QwZ;3qsFXYaZ`$v3$RT<^5a{$b^K@ zp7Njs=sp|Jy7BJ}n(6V`fB5c|Lv-Z71FZsBsqvhD4tS@KKs;#Mr3WK;TR~C66418G zHTOTj)$4mULMCG|n<C8eProy0eqVCzz`Lc8>4@jOko8e!Z1+6CGYhvQDwJ8FdqyXM z{rmg(zvmw%E|lH*1X?WlXBI?**YoR3Ak!~PKf7^(HVOTo_Rhua)Jt*C?foZa@)|0< z-SPaT#})7LO?*e$AAz=OeYfvt7k>$w@bLcd?|C6)oU{_OoM*YPm$<UhT-%u`;9|J0 znEyyWXoc%(=xP<@Ih!X-=9JBUpyqeVPHBCEC7X-YY0#qB&+6x#z-#TQuKyHHUzi7) zvAPFp?Sez8BxUIwJy7TNyS<dmf0rej!7Is5+fPUb4eNh-Y<FHLX6fhdpeAqoch8~I zPga|P{KBX3`FV9h#(6j=o?Gw5sR!EE>3A316>FUS0KDE8npLM3?OIs!LFtLaguXkU zOQs0C{Nt&409*sr3mQQeGlE;};LhjVk|#^&ywkXAx9C!2aU#RKrq}Gt;-@l!){`In zy_EOhS5DzGBJs$ZG|x=A%nV-earba-1Gq7BZoL<1(4*<}c{cU^TF}YB&9IG>>asgE zm#ce;+jV=So9}^USm`D^P-pn~IVLQt3@6Sms#;jG@hf<0Ji6#KXkzyB>7b_9<+3|8 z)n7k?Y><LauYlHw?F6k6i|z4<FW>e$l<k3L#s0=W8W+5uKR*88(~L*h>QVEjd!o7j zo?lttCHVIBxoe`}70L(fk56&2ikM#XZT<U38}EXj>pp)}-SMRqT9Uk9qHO-O#_i>r zPUx;ArhmdQ`(B;<cO&wIu_Iy=>Hcd!p_vubdDyq0=z|tCI#h&S{`p<qr`cJOu<XY! zXhT;h#;RyOVgbjT^%}Bky)7Y~03jvWHQpa9k0m=n7-H~MVl!cj%4Qbr1MRM=ONjcM z1d71e=kcvP;0S!m`1wtID^Dfr!X{(TfY>*=__}r0w?(xfvsVXvz*`6oeq|Ou1DR#M zz7;ly^YVm`rG0=^wZ{qT9VSoaf%ag6cjf-w!=iumGnUnoOXtjojB|lE<zBBox1Uq; z6=dk(2V|w|_2=FjU@Hwb1)b9m>02{p^WBe-9vt5q@aW}vk(m824|shnHU0pf^3H#$ zm$EeGkk?D||HnH@*ye#Ypuo0x*2F}vpYvR@f^nbIUD#fU$wih6i+*?LJlP#_9F(|M zmu%2G`dm`+yhu!*6?oJJw0jhJEtBNZZo!G?)-g+hmp*|efRYQQp3BpnrC0}=hf&?} z=fPv7nXuEWOK-a`G2RiO4GDVo*Fky-pw`$c(9{`tCxgvW<aMW~MVD?qZF}z^s8<Ev zjWMl=5ww5+)Qdr1l%f8#$oc6bo2ZHB_8Fu@I)phNcOC%E8gi|RFWqqT`1P{~;FGr2 zpp7|q0(b2_IQ?#<_--D^*xB5tpz&SJ9d%bDkAvoHQT7281?llmIk(QJ61-*|v~f3~ z|NPPUQ<;z#g?^n4Tb}?Lw>d35b@TF_`x^NVfMyArT{Y+RAN)G+bNW&SjEx?Xi#lC3 zA1!&K{Aty$WCu`Tr|M$$|Cihe(1!FKf9^c~tOMC|ZT}rQL|yM)uB!g@SND=`!Hq%d zz!Sp{Kobh~2cNF{3>yYJT*C!lNjDiZ^2E7x^SMLTJHQ(wSIJnwra8ca{^)xI=LPOG zwuy^ddl@vPtppx!0@csckDp#$>INAtKwXzx<@j<<BfQIe&g&yc`SpHp&;kJT=55US z6?04l?tvG~vR<#|o3;dPInq+Y5Aaz$w^Lu+PQCrR*DJlZ8!`jAK?Sn*<lFa{M9;40 zpmA@$Z?K>#zcg`T(cZ<N)-b<jJ$z1KgUZua{&QPxf(u^qfLB$_@CWT-KM$)sjDptj zPdTS&SP9>nolt#Vv<Ey{+YA~5`>~7l(;8Tp&C)1nUE-1$SIgbe&t+qF6`mLVFn{qT z4ffsDHK4K7b*B!#lZyv!u!E1Kx}Eykdg`XK$GzkG1)l2oAE*bf)>{7HR|%*4(hB5B z|3A??rGIK{k6&`9x6f<m&BdN<50JOnsIRbwRmugC)h{7iW$$^q+nss~I-x*ouBjNP z2U!zW_0N*4yCh-Bk6qCI8pyw}X6?%p9y{j-T18DN+Ppg!q#m>})kYqa-$4t77|tJO z1?z<`+ky^)#qWEaSGN5$19+uy#H6Cn@85$)cm5t|QPtd+#})^lP6am;K!afAKQ^`H z%fps;C-_3TC0`z!XZM1Z{|m+ZDunftK$}f>FDx<tZutd1KB~U7!p><P=MT{AZvXxk zjr5<0`CY$LZ~IPN^wczetplkonbn`Nf%eym#DodoJpkIwiqd*kf4c2A*Y)Wh?&?eL ztWkvyNP;_h!3DKmjG)=+&-w?yYO=0@Ei+m(<+kwB&AUsXv3Bt5XQ95iUn@bw3)$HW z$Itx&7hK!eV1x2SDNFCft+->ih+p%Y#5HK&P3`FM{)3@bH|{=$O!&c1`Kb3USJ!{~ zuY1XD_b)cvcYg(Me!t@V{IUH3_Y=mf)1Q0G?`JrFeEy>8@ACG}LX5KSnev-+>E>)7 z$TU8vaNH2|?%FHRUh^8ZeGYed_8h1|48ulhemB-ixBg#V09q41;W_&c&@KcfJNa<6 zQwJxBoag=^`YP`LWS8KaDZg2#E_P0T|E<;xTHkr@jAl;(wYvTsd-Yl8;8XAxfltz_ z5h=-YXFRAdy}tY0KBY=<={q;p<dfwx&^G>}cIR~X{!5n6@coy5r~bBqmI>5d1+P5m z69UZ>lxJIlx^C;rK|6CGJDm|#tJlu<L?vzMcO`zO?m{L;pKIt|YMct%8x3yNf%fNs z1}Ko0BK&Q<^j7)O)w^q#Lgv?A$Fhf=x*H`JT=4Q0c*7vt1lv3RDSCYYf-m1_=YsOY zz0wUjPx*dV*MQO=sB4Hiz@Y#1pA%^Pg>=ztS7`bJ?Lq$Jb8xF7yZ+J7HGec+tSVaj zp=C_|OSvaYZR#eL@K=a`IvfSwggm9X0kr)LRAPT%t+4wiw6|y;Vr(vDsg2)?JDgR` z?Y56g!2v7#=6z)Ze+p=DyAm??6tEq6QR|bXHs>^A_ozlL*?igr<OMHwedlkWl6bm& zjob&Ll=G0CL$Gqp$|zV)9A#Xz;Mtn>ki8KdpuO)ltdGJ^KzChRehQKc4*|_UvQ^(O zj$j9$)NpXs6!r>DyE-mC-ks++gVuI<e?;2R;~sZP*>Bp^x&Qfh=RX4vu7GR9zqa%G z4}eOt>;2w0?mzyFvT5F{@=l;7X#8w#2P1e%NI?Lk(5Sn7&{1<9XiHxm@&@%!6$?s! zcWh(h|88Rr?x3Hx2W?$mGlhN5lt-|!@U6Uvrm6PRk1j8hT$XGWI|>=k^MbVLj<SCK z3>u@}#rlH_GUGJ2BxNah6f8pUC1mAKC}@kxAJEn+=mvhIS<ZjDcPFoz@?3vXp8Uns z^ZvV^rh+03e(1=fCfl{_n;1YV&I*xE5Q#hWx#iN)XYv!z&HH2p9+?MC>w>0%Yxkc6 zoj%0#SqDBGQ>iXn2^t1tfB*jGPEZ^3LDl7heDd=4pb7a0e{JnRLk35tKln7|5wvE| zbPqfAx#`qf-F-5*KWupoP8ECJzscO&4H@*b*MGQ1<UDvL@2=w(M1J$GY@4n$n|I!} zul3CGpiPMfK#SU-BVuMh&_~4bSIoJ);WOk20<LwZ4$j<ccOE>pxTD25y%u~3kG<YO z=(-=!VyAncJs9?SZ?!-R?u9x_5|-LLTmfEr<_liJh_Zk@?9}H*P;v1m>%$h<5ji&= zJpK$SF7&Q#KED~VEC=}*g0NGcJ3(5$JoexH7#tBwvTr_A9$Ga;Jwnrew^SOqxBWnM z$DdQLpg{#%z4W>D67q^68JTIHSk#Xm2kjZIy4*PV(QCw-t^HH0Cr<wK%k|XWH;*Pj z2HPh*f4<_5Hn=lY`0_D)dyyP`Vf(3x6F*rkEZVR?{rRO{$XeYu?<*VoAAQDt?gwbA zWOe=eOE*nHtJ*=ULp>@RL8pX(_w;}gE9g83M67vK9<{F0Q<B}g_jxKW(xRV)z*F4! z?mdR)hU@3OVf&Y?jDkRWCFLjH=ZU)(3GUk+{JO5$Ra3tpDDGPaXotlL)RU<|LtG0> zN?yG=J?{`Sy~56=sr}g$stMXq+V~o2)%?;u?=<GdzWE5N)6*gK8)%mOz~_(C1)l04 zjU9Sc+6L_0docXv8>L=I?ticZvfbfQ-+pM74(b!YM-J6NLt<+TKfz3A&knMYZ}gu~ z&$sT>!9Yzrl+~0yQ?7GQ-OP5#ddIfM)0RH~ok@hY)~9w4i}F#>p8x+r@4j*OOu5e- zw=;Q3uwKReK4#F-Miu{E**_+Kh_l3Yw#EUEpSBLCT8}DSx+$u_{6S>T{+E4BN^^VV z73;X}nX?`IJ(~}7#0g>-DZAgP_l@>Tc6XLkH`fcj1g}Im02*e0UjE%-o@T}W9_U6y zmiO}m6<~{m4tRpLBEIFVyWc1OA2vAhehFwb*umfP_?~SBl`$XoEpLWwB-sOBqxRtY zzf_(tZ!N&<jQ30~5)pg(r+Y&4_kXHSz)4AH&v$*;D0B&EEu6y3Z*2SXd6O2rKMq-Q zCce~R9(ae#+&<Wl%CE+kb75<|{6K5G_{5eTak$^l4BAq0P9L<8?We2e!LNo#^}RqR zUm+bh16ri!v1{+)_x1Aqjr$ZK2Y%hZ4;rKV+zsCJIPJl&8IQ!5vRvo)$%n1;f{Y$k zcdtL&c-J4abAxl@x&2&{Kfs0Zyz|f^6J%kvCb(TQnQsqxjh@{i-JkNzw<bl;Yx<ro z4T>52GO-^bd+uj4l!`yF+)<YhH5s<zuHL&`$@uA)?j_bYpZ%AOu|iw3^7T|D>K==k zn#(u%{M!7ZPq6ZUO=NHa=ok~<=Res$^?=63_Uo{MmL{PdcQ#XVxv^LJ-d6Jup!tte zjnf}}?m75%+2)Us(T;gcAGM*!Ux8LJFXz@u_n%k`?|}Y!H}65!`JZ=fv7crzwP@F( zqU!eb=5x0?uz^QQqn3PL?_{T@0%|jXI+UO@l0bbS$PfUi?Y(PZiJ{UZ@HsCckl{6O zgX!eHrpmTmSdJi*UD_S~yQAcr2xyJey7SZPLb*PG)~kSKM?r_?ZSZ*B?*%;t<oyz5 z?WbQHUgkXi56-EapZ+j_Hv%7c{`i!O<$euTwDE|lMMVWKAZx<7)?NA;s(J7$co7@; z+&H8)UiUoR{Z74Xy>#=f6=I#Iy@u+Jd!W7#eE-89&{!RKL-zl$cRp^XEPJM;uaD_D zVSX2WNR`KP@WBV5el&F2dGp=vTM%;x;N!v6GkDLzccPsRI#<`s|72#9clHBVkY@5C zre!7<*)1-qZoi-Zw;wVlsWkE2{)3?AH|RXO`eyzQJbPHc`xQXvE9?Pv(m%5<-K>8M zG0yO3;b!p0p&G=XS$rwt2#uxSBf<`-gO3Q41s&fowP@Qk@LGmVTMi!A4=@L<8lGl` zylS`xwErgZR>8K_(4`3nS4}Zq<Mpv}8Emc81C$2Poc9`YV^cnXm)vUZxL5M;uk8VM z$Ql6WuBYfH9)ZpZ+bss$-4}Xp9q&)%$pfb;%<4y<L)J(wm$B4*uCDd`ev1s~7?$9I z&~x)A7oG8#!mRxC&HWbeoSr~@(`#cc3D|0WNzLbX!z%8%TK>K?tGwXtF<tNi`eXG? zuWdm)J&ixGzDj$*`V6+LdjHhj6DESTN*?X7iH%$H5H^q5T_SD}WCIydMcI6lvLwb& zV=ioo3+wd$)7)!3Kx?Xk1J6B!FRA)e2<vxiy4#(4+H)%M_XF9OS)dg+q38DR29KUM zoHv-y`Uvfek=~-63yMB=rA$2ciwm-VULJZj1ZdR_`0Nbi0}+2}OxN~G&jrooAGHPf zlrIc?@><htaoIH<=xdj5YD|~b+Ft(s(xuonxsYKS&|W0SUcrXg1Xt)@LB1LADYTm! z)2+R>Z~y+N;Vm<`J>sGy`$Ge~Eo9T#2TXQd;I;6l&6wc_ZhQ&^jfDMThb$RVIPYKZ zcUkjZTXE2yXDsu)cBfLy`RAnv^6|~F1r7a!r%E5(y8N>^sOj}{RSEE*0fYrUQ4Dm@ z-*j=W?dCr=RsTtn?SeG{!29xEKUcMYw*TXrAnPfR2Eq&vuKN~i1`5?ByAMqgIS=ka z%O9@=F9EvaUEZpSblijI&Uv79=kIo{-l_*#p9AVi9)F-EJ7)^kbNg%;mu!r1U9wq! zGiXk~>iWS+Q&xYq0u2Oob3cPLs}7W|HJA@8O{W*_T~JcC{TyT@_`xa}nJGnEuN|#8 z6EX#~X#zEwf>!^7P714asWaJeOF0g_EP*}llxW%aA6S;P-UA(H^|kHN-8<7EyF5TE zwYOT?fkq7&=Qrv@T2+6lCn!SK3q3vIwewz}Wd(l`kB2j4jctJ+s87ZD>fG_AohAQl zPaN<rsNKh+4Q~2?%WA1Gum7`lCA+-0dd<mm032_i4Y6^ju2#vy29(zEYl4<;f{qg{ zd9pNSzQ)~qkFv89Y+=y`Ij+##`Q3AHlcXT?1Y!z0?$p=TO9!hbBs)kKodzAew(a|m zAE2E=?5l%n&=yGjpZ3lz?$le*Ow|j|8uQnspzijH>F17v&i`3<%e#E@`}leA^^WdG z;Co^AOu5au^mcSj<edU=hZwvh{QGy%8ten!JHe~MIQ8zA+U<jMTtU05tzVw--Z>Yv zAFb}w&YQszkhpjaZd4w8%K4e?br9qfop}aGtABUK23Xzxy&rr4N|TK}a$jlMQk1^Z zlcjgot(ddb4Zex=0HpVM;Q529^9N^c-j~;CgS`FcHt*Ei(KUZ}Y>@+vx`9ggX+_NX zpq;j_&1H4GKY|X#01W}a*5y7~dgqY0WoJosw?6OBvsFyYA6Wi?M@puHMoK<$|KtUo zizWGL+JjHrJK+UErRH{auk`J&_w<7gvjwgCTl&LB`H26)qh7FC@$RYcJxrFHg7W-B z`ugVVjm!pZ02We`m9c#N6Vw7@*2Z!`320ahv|i=(k37&|HRy<&uj`t5?aVg><$d$G z0clE#!`6g@#>KQj<6_-{J72&1cJ=@yQ%8G07LM8X3Uo%-lFy~$CnhgMKD;ne^SiKD zx@<dam<zlq(^0d$0JOOUdiw3@@VU^z7(dIRpm*t@b7kUr-+}r)pdA(L>F0j^04*oM zzJ2>Q>(blSH5-E<>&QX9Bw4%i>{xIEZ+A@^xIeif8r&9uEcXP}=f6ck>F?*k=Xv~z zvY>;v=6)5k6R?*H7Jd0(i3oUR1-v5$u>o!8@r#=_eXW<jc<(_TWZ#Bx%zLX3L7;^Y zU)b#G7{SXwo`dQ~=%lrK+^N5Pmm;+v!PjsxYahKXt9YK}qwxpUPbo9}IoJ7XB9BqZ zuefvg@Ahw?1#aNh#k+!_Q+K04$If9K>5;PZ&ioa3?mg`0F54b<9J&<}G-y@BrMLge z-L<BW5sjZq&%u}Ir!2h_uQ9h)@;tno1=={;0NOlxz2CbasHO(_6vKaO?=;VWt`z$| zuRWFpwEod`7rO;yOIFckThJk+pyPEx2bHiv^Qb3yrP$%=?Ni$C27o#NN)!2@L%EfE zY8s|Kss!(}tbkS2puYX%9?*b*&C)l_m-?&cLDnYQ)!8i5x^(lX0jNxaEKf{@E>8p( zgN2SSzwyOiGehWE<5k%J9(B}`g)+Hgu5-bU$B5MIznfoi@>ZLKp!l0a(55M-f5jpK zpb0wAFyE6%1v~VKnd7`mH?OvZjWLH>f|4x9Pj$p`BTpGGLyI)f%+hh$soedpOE#}| zfYn@|8bUSg>Ocq5*?`VML7m0ae)`ej<%11V+CeAOg4%Bd5nVWk_%;Q}`G@rJc~{;D z1$AI!>P#v?o;{wtN=a5`8n|Wzt-@K1n6r-5>=*aizWv9q%FO$(jltt@Hy%6&7X{6S zL1#DG$k%g$7EynM%;<my$@YK-$=+q^wwD*YJZ|fM@agAIoMA9W!Omns9yFQy$!1|u zb-VuDgWv-JHFw-8nO6^*9NwVyG;=MaD+)RWVg7Waoivs6K&AM5EAWwn6FrtRZw{Jw zEeSlQIZq!vD2<pe^#dLCwY&bqq{BbH3Es{AumyZB^!43*ajkD2yPvRU71xIEW=~p* z*v+o}w9@tEH^IsSHBq2~5OhG;zj^9!K79PCj_u?G@Off<eBPC}T%|yT+Eh?mz;MUC zl7fh?=iH!`R4$elI*9f$=rl2{?d3l<F?`O{?P9EeC#4!Ly|o7~N2!5Meq@<`?Eq5V zvylf?!>8A?LB_qBR+nthdCD2r1v*MF98_a`J_%}6hl38|1f9nW9#K>WMd7FDyO*oD zLGD`GENkcg+Op?(h^tjJ=m5wkAU!CDV}iF(F1lY@=6)7Dh?)?4>T;D7@@|oXTa`h> z;9Pq&|4)188h7foxaapxa)+QTqD6POVY^qnK9-6<K<h_<j}*%egfFTG4dsB(`irz< z1|K4nV24~wUVj(;M@ajrt8`7=H6iHsiPb@I*O(*~|BJ+c_KSj!ltw;|W6hN7tV^30 z2gO}0BygIhXJu^Q&b<exr%#y<8E;UU2s%LFhu|ai71r|8)j+2ezn=CW$@VvNEdzKJ z*us+C%?kgXo{o}*1n6c+3nKOHn)&CMAjgwhzeXB(eBF2Hruir6IZNwKHBSE2Wd7sF z5AdmIb427pN2Md$m!~F9{#3P~X!phi=cixW3EoX7yT<Dyc-8yfmkFVeNzr$V=CT%$ zo)c)%&U$rFyZTe`-P3s$phoPc;JKiM!5=ht{Ftu)`ic(ZP_#KyVI#+&j`i!-OK*Ai z)g|1TYY$nJcfb=eu?yad4?30gy;T8Zfr9?acl$vjQxBGQPrO&jcklb-1Hm8uEo1)2 zc~6`DBWU1wlb<EZ8sU4Iknyi6ySs`$u+0ar>jJmXLNx{Ko&IT7-0#QS1_c@Y*>F%D za%>ssRD3VcL}y*7*N1=Z%+RBjEB-HLJOZtS!CMto=GN}~Xk7ssmOA+R`RND0Dqnv8 z&M;r$KWORP^Y7msP>x;#4@_Cszj_1TguH4_d5}UKlb!Qj^iCakmD!Zrz1)9yZ21j3 z?Lb5E(oa{-1IHiy@_jy61iX~B$<^xp5(n5SM6iw}x97I*TiM;XPaZPtwLJ(_p4HX+ z7u0%u4C6vQA^~#p+S_wZ`L@TwLu=pz#K8W(=Vhtqxa4zp5a?h*&U?EP;C+xr@Rrfv z$xHH&H}Zqmvb}%Lx$e}@#X+ExyZ`)?>j2#!u)64;Cv<5Vc*O7N0nbiG(0L{9oD<Jk zgLY2}UzrQqk?P49*YrAkZaH{|>Cf#Y>P10*yZtY`KkzHK=ft{Q$qS~RgB(UUBS`-) z=b9;#YdJrIdxQtPKm1#61z$ZDskvO(>$(42Ug#DW&=?V9{Tp=c>A|O;FCmQ(EH~EL zZvI1k;=4Ehl;s<iKloHr$)pbLHAjOsp1(>1oo)x;)&felFPknMRRo>1qXe1+1Fb8A zjv-6^_z7A|fqFKCU-%ThZ{LqJfaX&`*C2fV4(SeH+r|byNdK+wz6wqC>b6thlJ9`G zW#sw978~I!b76fD<g#FT(XPcMb)jC6s_yG&p*2$`g9f867qTvriGc18Kxv-Mkpnfj zPWzYshB*neyhizhX2t$3iaWg0K_?G?_m)n8u0aHkthu~=lhgyPPQN||*|F8i?F8mY zGxtFAL7Cqzf9Y?rKiZMA${Nut|5eEucIu?uMI+~VPj{qD<Xbak^46>{=<eYCORXnQ zE&AtrY3^F#kKlHr>(0%kK6Vq*ai3oA5Aue0aQU|1b)aadiJQ0fD(Kb$VWF3QR&Iuz z27$6VcPeDOOn?9PIcr5BdvpG4Jppw)yMxXpEQQ9M5V$;tpX_l?d`e!q%opa%>%tIg zk-&$UWEHS!AN?#7vvhOVSuXg5f%elHS8T(O4?azJR0*1OLhM6aGADn<oji*bC1q2O zL0Z21_dlO8$%B10*1@AEOXj=>En0m2d@1ONC(!xk{&Q_X^Wam9wyjQt90UJZ^B`y+ z8La6ZcIs#2rJJ?JpabDVVy1rnUER=sFx2vNI>x08MM38tE{Sp3x%aTR@S7Y+!xC}? zH0$+i3~NNrgG%FQX4rBuE6{SWDd*(2{o-5qcgIuCAN~cQ=j5%pL5)x>qhymoqhxQD zH?7nI4<)jSzqA1l&mMfbY<C2A3gf(f?srQ+Bd>My+Bqk{s`|#O-{2MPpau4z!?Hl< zQYf6~0WGan+{ZF~esv9Wi|o@Ao;&A&7ZtsOtkxD%LY&PG+Md-2zM$mv<E5}JnpIJd zp7@k=QXKEW=h`%#Ha^7)T95<Uu4uQ;XG*iFCTeH7XUb_&&|zUOdGyNepNHg4&z%v^ zlUjbf+`%%DZw|N`1#LlwR#|=^E6k_l<$pNpc@uH~9;kyox#*+%Qipk(J1oxIK{t0j zviS`^cR?L=vi>)=`1e+bbtgxU_Y2rVwuXQf5grVM)rw#%9xkzQ*>%@(LfxMvP`rT- z2TuSU;`Pb@VCvtckONzg519`;^|bR+WIJS&IcVrb9(*R^j0ExBGg8o2M}ziUa!&=H z_p$BqG^F!BnqNC}<-AAQPqcsP?n#r2{<@xuP4qo=RZ0-lE}8QAqdfF1#14%4exsmu z>LGdgAC|7%9hDC1=7X+dDK^{T#hx9s?%Kzhisy@BHi4F`uluHR7jgP(QP8^dDSG!N z>RgI^4m|+gLNh#K@*~g@!97K+>PMd=9bC8vv}^P=>(a%hn&H_SLCa*p2Ree+ESPGB zZvc<BfJzByI(m7+YiC>_c!T`AenxQj;p@8Qx3=@Z3y`x*&p^r*p_n=o*cviW-|uVN zrHhNEtiJmWx^z-(Y0TQIhHT+469T)=gN|wZiF^S>q~`Uz`8!T{)`Y{Zy<qxR93uAg zfalJp$&Wt657_>w4;uplC2u>>%;ogCk>Rf`Km$}BOTf-L06x3_b~X5@nZK)G8}XMY zt3Una^0JKGam|#~u}h$(-7mQg$T6?`k3WwD)#!*#HlQxbI(^U%@YKx1SJz#V057~u zx9x(=<IfLzQzP~Od~``L7v%IC@P^7D(1uEB*{Ma){^AwTfz6;Qnm26eod+vmLlD#B zd0{yfbQ=8JfL(W^Cu~hH&yE8{$yB~+;B{D<J8G`~{J~lQI_hQf@>}Z=8<LjYIj14J zhue}3eDDwG1{=L6pvxv^Xzl}Dp@iC;IpG5uFMADY-8+M}%$m>rYJwQOEH2pl^1+oz z*oMs_=!lsKWVsV-_|pehBI6;KrohgEEGdO;#{~^_zO6g#TK4aLM(C-tpp`-*G4^_i zp<Vs*pss*?P(o<e`q>NfIM?yRheS$Jmfm@%aW~KM#pk78m0;6l;88Z{shJ;HpY24< zvU~5e4cN8!==Ax#dEZ_@*DqY#{R1+rcys5Ui081`pr0GhL7Ot5Zq{~nukF!%n-;uC z+uI5q*<K3Sz!WcIg4xvl6qKhwh0pBA9q7^%%ZO=3o5PNS_c>!0;ghEp{dYZe_u*gL za!|PxT<~+pp2Pp-Altj`)~yLYcO30D4t4M$`;n*aT5&@bdnH7IPtE_KxufRtPn*+) zCqYGYdBimMxSthh)ND%MG^8OE9q*loH-XOGnx4OFcWH_k+OV9Z=67x{_tQb|t|fu` z997p3ZUXPDx-z%lTli%{WY_$3l+z0%HNS5Nic^x^!~Ii2;H8fL!LRE;@gS)PT4D54 zb`N;XwNQ*y5n@In=$-kLb5Y!o8Z;lYJ_odGUceq4pr}VrPc7QNxabSh{N3;|Bfq(O zqhBY17F7LR)f}8&i+&mr=$Ns&tS?`^y9PcfsOmasn;K|p;re>38xI~=f^tGL=wbrc zz?9cc`#`JPKlg)ogCmB5EqB~2*`RgwGwQj}|D2w_NnHK@!L9!2c}=T}t}Xn5m@7LN zY6)860zSF?PxMZ9OU?fHC3g;6aX_oX?jrtzAhatKK6xV#aLR#8@6cWLY@oAUS*NFq zKo6#h<ovdw?<90<^uo=s>BN#JOKi?*+`aeUw0<dnb;${E4VvJ4>adjpct5T`_{5UN z=?^YlhArC%b#-s#wyxf@>Yd~J_mvT#a|=O>`7AUc%YR<-)IhFvL9BUj_dE5t3)HK$ ziMLxL44S9noCrGKzy@@d0;rA#9s6JEjj^S2mvrs*o1jrC&{mCvX3&x0nmg*i@c~}F zv=?!dvU}L6$9<sd9UudupfgnJt~O5o#GsCK8U7L*P=ouy=~HdBpabYY+f$l(?Vu~v z!e4^y0!=4BfbN3_1=BnYSzRy4IUKU8E>;m!K5e*mR3c_+Gw5K##$Di}osf>r>KFG) zUpuk9Aom1#!jccPHf`#s23JixcvI<mo6H*UveQ4&JJlmK&&MyhW4Fj&>J9k5u|A=_ zCH^<=KYn~X*=cWy|LccnJ3IDF?%OlvIp@;N-G=v|7mgkLCU^e$gHMGo!3QHXpFY0y zIdn^3<#L%w&GUD+@0jbl^KNIr%Qc{-u%HFWeTS4}=UDyTa_uw^Xtg6~3-xoP&ADEc zaRIyT77B>J$po#;WBOM-MGSdZ6ZKsCO+oX*r}Ulks*HLKZ3KYM*aGd$TiXaeGq>vc z!AMP5u?;#2>R!OEyE`Ye{Q;W}+S9%YI*E55I*F%q>87|a;?(n$r8c1E_M_9cZZ5aA z0qsP3&8}7fJs5J$+^?0>9#3+$s%Shf9QWWce3#Itpn3eDO7!cx&C*$*ZiyHB>Y#ZK zmi|EQziFWljzx2C>Up2J7qkFfsI#Qt9pnhCx`abJH_v<iKx-~?;)?@suH=6x8zT*x z+4{Sx88Y@Wp&s07-5i8;%rkg<CFpdqo#Htkpey={&R729)RUiXR|i_8drkW}`b`y6 zKi%r>o;|1c|0c0RkW=0s*t@3&@`Zssf$b)-l%;#lt(a3P1KCZ%_@5EHH8&KrT>*NS zGAO@7O0cIVJS$@ZEUQ19wBAtyIuNt!eB?jQb>C%VrhMAub{v*6KuhA`$3pB`w?bwQ z$IsN-KE?{<B2SBbcD0TEf$$Z(A-y%rPl0jfA$j@lzjmL~105X=+VSWe6!&29575H= z7x0DovhZTe>Qi7``INlvqF>f+-UnH(`QWRq{_EnN39yqC&Y+x(;0`)bZ13Lb)83_) z?*_G9YCRfFHT7Sw2d@>({!s%uSm%RK%(V}QGOYho)uNK>>h=Db{@=eq8;G!jD*8V) zfexy8pt|FS24V@R{L`;aFTeeZt4pc9S!)k!Nx!!;2#Nz;_*KWXu3T0EdZ^I7Qt0U! z@t17t8Yh0Le!&0o*IG8iOZ9PV?{61A(5zT5^MU_^;18D4=Z{ZWJ#9!o5xn3Xtg||W zSsAqdk6j(Qpz1(9XwedAAd2greS#CXn|{t(^WblGA80q-`Xy)x&f(JAb4~Y3%j`iz zv+wu6>bs;g*Ix4}=orSI&<ISJ2OB<}3>_ca>V}ws$Y1I(&v1tpXv-k_U@KVl5}7?* zmf!(6u636x^`C0!fXE#z6VFv^f^T(13{_rgoDSNRxclex(-JSFL1+5cU4`iQ`PLRR z6!-X&i`8t<ah*t=OUOxNZ^fMU3+X_uDWCftd^rJZ%KAbi@^uX0qwpu?ZCCmwUjsQI z?Sa<Z@8BV|d)}NNH-cx}5$k#%fDaf0og@bys4NIQr|;Pa8q$3LI*4J;_MioMNb|5@ zPyaVs@4bJ6qFa@)&W>-7p_ik5RXjR>N>NOs#`E)YTEV*k?3PrhgEwT>YnLmjKmF1@ z<+gju*<v0Kd(h1L*GqHPEd9u62fAu<|1(MOwNvp;-?Qf_Lg(sVp762+wZ0=?)PLw? z1ov-3t?ENg8LEOOxF8!W(K;*3jlI}^SNGK`U5W(lim$5&T|)SIHFyOb<d)FW>_|s@ zcv;2<SnZa248E$)JNWre$ciNMx%Qf@73jBb_=SgncK5!0YzR9q40P5R=yE>rSwx`2 zEg`2Sg6=s5^|W4s7xvpE%zIM?-Mv1o9=sd)@A@?dLCdFY^b3&Bkpk_ioLeti1m2!- z@ayMG`}e8u&|(K2uL3)6A9ZKdW#*~PcWuGW(cIzny!^XNtE*K(kY8OnKV%srN^jBA zvUTz&5zr)_&xXmM(-k1MNrBEk2ki|5-T9LEPjyEf%7p<{3rqBOy<YPiv`2(<;yL-L z;H8qQraY*+{&R5<<QNmsoGkbh)BKm`o-CbHKe6O@#~gnVczr#uzVFc9lJpNj=l-n_ z0e7x<A@}j)PW^1Tba3^AesHXTF1BNkZ(d#UFZ2Yc0}I-I^Dp!Synpp%>74T#vNG1v z%%GE4nocL5lH4<8^H<OiAM5n_k&xLep_sZf#N4I+(;CN@di%^jfG#BT4lmy}?-=`| z;}2HN0Uf(120k`)`n?12dBr_bK8sG>d<{~efVM7zPBY(?T<~|t7HkU=BQ>8(Yi+mY z{<Qj_odPfPdbCwjHpluJv4y|X!M3z&DR^}3(d*|+Z7o5yHF!D+e2o5@>xQr^RA0VA zET5ZRbdR^Ux?Uu{=Y?<jPtciLutUQRew74m(&0why0T`<XJODmskyJ#fOdPhD@{C? z{{g()_O&nQdIHe0`@Nv%3UmV5%W{8!)$5Wy&Dqv@m7qiO{@Oxs;)Hc7?cj@yw96IM zpMC;oWa&>I{i|NufR1_H#hvp3a`A+?@R{eJ1++h>JV<JSU1A36;1@Z)EK7D=x>+A| zq$jAdTMd~7I_kMI!aa>~UbA;_IivbfXe59}LQYFAy)9i+>9sQgI^rUm^I_+KRp5p3 zeqGR^qIKD=$fHr<De5|f{WoC86E(|%TCd9=R9!#FCx0IGjurQ?Q%}1tUF4fz>mXHh z9<=ag+R~T<-XA|+F62}KO%HyQeB@nF=>t927`)+Ny0lk%cs1x6Gte46c`I)CnUL(^ zPeC)zkkL%g(6Z{9DeW5R;WHr<cTE47{|Lvtw=4+KyH*KWJ_j~@?{m-%qo7NIL8B%% z0Xz2|p1$r~(pqSr*z6zZUOU$5bvAmwOE!b2%3dCajkAD;e?TXh9KL?-KoqRo8hYyJ zIw|nP?rtm40c+qD53$e#Tr1UOBQ>Xsdp*B@zZiNR{%p{ZPmQ3~H|ph_kV_BiwaZnt zpZ@BevO6mMA^4U8;U_;WwZ5PF`NMF>onkv|m$#`ut#W$$jV=Dx$JMs8Ac-4vj01Q~ zYWI53$<Lr!iU+&GV4G>Iz@ub&`5#g(FWP~&-h&R@104ni8bED2Z9W&gr{SY2s3(q? zRs=07+pzz6mj7bVT9QvaMXivXG$-uAySe_!&0p?t7k1_z=p^^AU8mmeYr0nw{$2hd zWCS~6%A=nWG4}eopgr;6@uF__XP?2n%srF)dZuJcXqHdRGZ%RB^HY{A=*U;0&XNrp zpuNC#TzcB+>#t1&rF7&iGEbI(j^4lbuzPCJR?yLSpk=5VAnE+;I?x0+#v<3(%%GJt zANALR76XBjX(L!G$46o4us~?n`sk+D%*fYjU03&FcV8N_))_ofW?2xlZi5DBUq<>2 z&>bbHN9lU*+!tUOF{NnhxkC&#tk?U!Bc?op9?1sTSB+ROb^p}e6P`Q&g3b?&wXELU z?EvZ$c$9CxcJ%p-ApN<Zv-%Z57ND*YU3$k4)B)Ui6EuwvI%OdNbY|5D%^fyN?yO1O z#Bd%ohZTPv){X<s(3H7>>Y!`!?;sU~i;`^4d+=21^!wF$8}=Rtos2vEel_A|;jmM0 zn?MyrGAN&dcJ?>_wKb2J0^TVCIf#PkGvZ|ICrj@n1)lO<0^b%0YUM)?^1QR=s}bbL zG4Q@1SWgJFhRZyJ@7j_(;9mH_P|JwPpsU~19vFXMt!Vox4!;l=v`%-swbpiN?n_d4 z;I;aLC3oH|nh09&eU#-h=(;4d9oDu1RyXbyZ)Q7+d{BuD<Zc$^Q4r;)-&{ahXU_B4 zcl{tMdCpmhfKIFe9dP}T`-4o((&o3eh@PV7&b{EPy3QTYgU$GWZaS8TSz3MYxvCCm zBRtZ+zWkT(QkLEU*Us~97Hn6a3>unUF&%nBr1hQy&-45<W`g%XS8(c`+YDb42U$=i zUJaj;Klr=)-L;nz&8v(4`JOm{aWjDW(>gbhKk7_&+=eBK1D-n%e*@oC1}bcdK%?b; zUHjJ|kCngM?fb!^vqZeYYiBfQ^rYdx%u-OVwg0>^>-6Q|Q=AWmTK!g9jM)D<XUcEU zrHkL%mS<anIyh7L)`-kcht6@Vfy{A;Yac~g5^^VBV{UE!S#V!*`GZxUS;sn_b@8QU z@X-fwF%GU7K?A7w11(=?K#P06xclIJR!u)7A$f--{3&wrZ4WBG#h?GogB*G90Xkg< zd=YEi%H7RVK{IMir}dG}h0uQb&*kMBXZUU)_<kL5=h!&C7G)ydQuDjG*7rr{4p@mp zhA1EWwFNKf6N#~AuBbmRJ4Zx5z9~DM0e(9)cr6)d2j0}pkRy1M1NG-Vpa0-f<;%y* zA0#!|u2BZ<7DJrj4mu8?a!!C%RsH$HdQ6}}%A)JGYp#9-ZD2k4PwqU{7I**M@QOLE zR#l>(4ugi6!9_rQKj;Q|&L5aPtiCni^Y^W53=QP36@fQdxo3hGk$_IdTlcvV+8ldt z1s~Pi6eJfuCGY)*uj@AJZvi<9bUO@xc4^Oq^d`u07~ne^5DT+61<9FD$=fcI!uafP zl{jS8r^2O3(0wM}(a&G&9VjlC$`|+GF?`8M9B7Q}=2bc9k$A5`vj%o`Ek8JafR;qC z|E@+o&9PFmpIwW6b&y=y?$f9_8nUks`ND9oN}E7StCDYt-~WPgo5zx7XXx=y-{-MI zug*RF9_g}b{ih$nTWR+^Z?AR#1S+%JML?^_K;4W>pb=N_2%oJl?5x3h@9!%5Pk(kV zsXi2_>0iz9n!$#3`s{a^Yy0@G$b#2B9szBAK-oeA8uQAp_$TyqjWJ}i80cn@gPxl6 z`X|&g*)fBzaRMER04rynp6~?SMtQejLjURIrAnYp_~4rp<v}Y`cKnzQo*Aq_+RS9B zd7c}zz~)zFra0`>tHnW}%bq=!fc89qs?DEM9z0rtSUsi=T0K^#>@-h*{<W9jWN`qz zbyaMs&D%gV$Tpz6OqF5qqbSv%Ryw~dd+hY*^zvA6E0+^=f^CHzTR8YoI>@FQP}d1| ztU*!GymD|dGCaCZ8MHYYywTc!Ll9^c?H_L}2jqdbTHX})fXuc)HkX6$u9~_T>)!J7 z>Y(*$kDA^>)>H~9fsUX%v<bY)5<D3PT6g7LfYS1=+!tW=x@S*wwtZful~za52Cb)@ zeqH^@hXQKZK@SC}*Z!`g{&eea7XO=B9~OQ8z6U&7dGIL9=Rd_hE2p1(J|Dc`R%ZhE z&UfU^H&Z`ZEiBov|9L8?T?smWCU2KG=*CCTF;<v+<kuT(mA9RHkm?UQvu|;W2Rr!Q zj@Q2F-|vHG6hR%bYpWuWdg+z10lW4d78ie$3R;;5nbQStQ#jqvRtLJ0`LoW!r|!rZ z6tvOu{fDJ1t8f0*0}UD38|eNA^<5*Td}?yl)X$y?UYa5l16rL38H@uBi{1^eve0^( z3d%~LliKEpoIn2HS0(7w2-FMy-h(bFE_HuecGxkjY44k)si0kD^QV7MOOXd(SPGly zM7dDo$&x*C8h78hWMqGU95Jgr(IXu?9^kKejVZF>`2#K4oDYbT+Q3JUedCLl`ZT?E zK4^R?9K7C4`2**xa{=YxBdI_|6?8X8q~`j&5fz<NKP}oE6n9M+GH1!24xRjd3Ox}R zakk7K?_JOJK?k{4yi9TtTXWpu{zPW~1ND8E=9aiW^s_o)%=&&W=ol{0C4Uv~mpZ^E zg`opRaW4N`B$+{%Q?kC74{oUMJESx>etJ<%lgIP3BH-mt-}%(wM_4B;106D-5PJ%A z1tjF~q}M?T=RpU-{$*DKl_1D>4=@{p&g|d*{m0*~ui)D$!AsT371^)vKDXYRbDjK* zf3YWyt9F;v#d;O|J@6$Be8|@OB@lh=)}V2zlZ~L(?oU}e(CJC}%Rm=df<^{V&&`;j zxu4y`DqCq!^LJUx4<a%7R{6_t-lp0(^HX)hl+W!;SwPocYVOc_9=;T+7<RAruL+O# zLN;i^Zj6jleh<2h0Ca^F=qUQC^FMj*j*E5_)rE2$N$=1<7-|XIY-tHvZujT(p1QAL zUOgvPfiAz?Y{Pn8-)g@W`0h(66V$nBdM%E6{eXuRXc`REy+(`)OM*A2!u#!iqIU{g zYA%moa_6L6?I-AkiyF^=s)PGkT%bX4P~-d2>5t;CrhbO5v0g4CskvNOYy0=_moCj+ zBM6#7QkuxuGeuYeya-x%?)l9JY>bage{ks%;-t`?DVHUeDl4o0>-0zmU)K!1)N$E^ zq$S|<!BO_SEt#XIard3}iSWB8rJxNDP(W$PN@&i*dW|Az3E90?8OX{O7bVagh&4Rx zF8wT5gbZwemWV^vKi>nLzWB23(!oGY=nz&NXiNz_LZ&s>9(pPm=xjk~7P8b_ZtcbX z-l-353t`i);)1{HK0_{kK<?`<zpL@1ROn?{_CkAiLpIP+@t__VXbUQM1^#gn5EJz> z+qhFNJ1-s7u1RtLZFB`)RUZ31{$a6?$CA(OOCA1z*Q<W0`>68-vT+YIm%Uuvi#_ht z$$7Fy4E~5iMcA*O_x=!sH0GcG6f!tgln%Qem`m^b9O$}Jjpv|^K@(8=te|l*&;<*S zsXNe;sPv^Q)2q)tp9x=9gmiYQ<))x>=^=ULN>3ht%88!=T7UlZ@u!@%ec-FpKx>_z zU|DegnRV&z9m}^tHr|5P$fno7&;j2IzG{wri5OaE(lS2KGHO!Mw)GxzpbZ?dI-Zuz zlR@W!Kb>v`-M}$D{uN^Ngs0{GfL(hJz8AMW-T*pl30|aH={YXheC{;(dab2~AE3z@ zv^4j5LC`#<yLmb%PTkyG3OklEFk%U_`cs~6=ys^-kNrU>s=y92n_N`4prpE4-~Pv( zwcuG)CC~z=^X?D8<2aywax>CFC&D5v0Di9SwOyL~S8V#+R#4plTDJ_EP*j@x6twvP z?Oy7lpmX|D&b^ym$$q`=&L>bjEI;S;hvTF02TAm6P4l1XJy{Y{H?c^(qWL_0c+Bd& z{XtjIB@5!ue=>!E_Z}j(mLfH$3xn632j4wfr30!Bmx4~b*4umVXq6yrCFPmrEYp$a z*1dPm3ACzaI|*GvzG{xWHtP?G81O23&}n*NNL^OY@%vio{Tz_3v?4M4UfTQzbzQtZ zVr|Dv*9Ns?F7?O1QGzb$NC59%>+=S0-2hz=w$vCAWs#cGHw4Kk-M#m)`)kcc&}IS9 zg#h!XgGNebHg$79gB<n4si!yNJoNY!t4%?AkhU1}Wzcbypjm0i^>vy%?3Vo43>{XM zhh6gtx*cbyO@P(w0C*3y5!6Ai;R0>DTY_?C6!@fg(1~NWqabs9tl>`&c<gKdofH8- z%Dr#@1n_mJpn-bO*wE9~Q*U|qd8PY<HYS1Az507Gg3elcbq+KFeemn-bNlzRL$0|2 z-9>4sc|Bg`?o}y4$kxxM?xJu5wDkv%^{;5du0IOy>zVSJbL!@OhoJL_t7O)|&!Ptn z7oZ*pY<KEw&#AZn8o!kM24|WFe;00sG=_VoXy?2KpR)|Uo!|$oW#q9l4wTKr`QK!n z=LPkUK23P^6MPf;<f5&8=VfB-^%Da@<4E5A&*7aOzf+~#6@RV!?Qgh9`{253f{+Hu zbN;!!>*7lvEUCbD*1G4;dmxLhFBbu=>w+!EdkI=mxb1rd_$b?iK;+AAuUmVi^L9Zt zqJfT+0xeo(3xAs63tDM5kLTmlAKaZNr|QeExO4A8_vE6jen&ukpV@qVh#QS9G>=|C zd%(K@W98Ty=a2B?RwDm_&RT}t@&exQ1v-EOI-vv_W8DW{nRa)=?}l%Wr-8S+?XCu$ zl5sHA{E2AUcksE|o1`y6hYaptf{cm1%qw~J<}CP*U*$mkyS(S-eEpbQ@N>@|@aZm0 zl}4b0W@aOA1==p`#U6I*9QZgx@Q?!Z?g7mB4m<U>=Tu}m<N#$*yOz&uXCq{Zpzm|Y zep;D&#X!&$HSB1r+wYcq@t;`VC20HjSsHw3Qtc^QHh6E`j07xQvhCVl?%6?k@Ev}T z!$FR#c9d*T!*aFEoji@XriiO<AQ#y>{<GXs!Zr_eW>#7GX_@=e2VZT$=ca-OZ=tt! z{rhSQzGxYB$(ZL(+W;#obwvAPgW6NJ=SzP40G&8^a29x36a4f>KdYjkJbzHX?dkV< zvt!plr=&sG8-k9N_t_u~xwHm!Rm2b2j^n9C)eB3i+x2f(gZEg1V>#^9+5P+2A1!}y ztG{__(dYH=8-EDJ{4Rv95k?vsdzZ5J_T>!FURRIjKc8#d@qYez`GZT>enM`EIzE5J zbd>FV{u*<0|K0~*gNitWSSDs_^JL_@-_LSW(7VK?cU&!Dw_||1Vf$X9E?Ro5uL8O< z1AHj!sfkmI_Ae-^Zr9hlaW_hY8FHO4<SfDI^1OdQD;z60_3WjPPm2643O=IkAT-rI zI1vfj`wiX*U0e{#cMaa*26d3?oS(i~HVv9LLoILIdHneJgG&-WU<XAM{b_Zvge?yS zZ4xaH={x7S^ET)zVUd`n&BZ~WF+jyS(E1ch@Z$Z(*HMY^mH6sU|2e!YJM37Oycecr zvmfZ%_sV4te)*<91|2MsF!7v~DQqncsP_M^4jMPmnp<o50WwB^;~vOWzZx^=_Q5Zt ze!BxcG6@<K^K$n)Ww+q%j;EkrNaJMC<(kJqckmPhy#pPXBLumalMB4}%F*xCV{y;z z{V&3r+rcYM-515}V%PEd_%XR);yXrm*abNX=i$d*Jy~K?roQyQ>>q!_MYoU3gZ5N` zHtmLjr&oLq9<BNY%F+SbS;Sx8gx`7*so8JrmHz!K^4!n9SKt#7^}@RPk*@~;Z$o(8 zbn4>2dF}8`RiHB*!JDes!$A{upfzA18ayls>IugL?z-zVLH>mNwbzg>PoTv<Hfk!D zA|b^#>XjP$PYa!2mYsLp|7y)n$idwa(>`r@b^^S760=o!X3BBdrJME5>e#2Z-T4f< zLtOtUBY4+b_jKqfV&DR44|tK*p5VScQ;suFT^y>(|5_I^3L6IM)q?hfv0tAJT8o1+ zN3zsr9yHI$f44!-GoS^FNbO}%o{`fBO}Re}zXsm235wu@GdF_{Wq2&)Vr5Z$?K2PP ze)X4YAPM5<_7eY1LGyM`y=340=cUBUsnt=Caby<pr=T)^LOq|}sY>J%Cc&+R=bTe- z_p;u+|9MN_VFnx4>7UJN{(#p3=<S#Kv5Olt>b}|iz#8~TaYaG%(x>#@YuUpOS;){? zvf<!!PS9R`(5we$7iy`EK4>pQ_t!sLK*y|sx*!qYZNs2ea5Lzbs0XS$YOW%BXltfC z7hSse`}ZUEf{=Be&AfL0;F~}X9%Y40FnAXf9z?XWy(;eoSnW=HE*tX;eDVfg-2Kli zxQ~dfbbtDdFMfW|o9;$-@X>XVuOA7=EM?Yy%9#$jWA@2O&>hrusbAs!H1JT^TiJaP znvjh+p!w8IBIl2TE-APh0a~yNo}B!<3N*<DxkF3)=}*M{grLo~8?<0Mu|eYr&gdhO z8-h|Y4^QTsTN~w_z%Z}jKTG(_2Nxper-K8@Io<X=c*O*CO1$3tyOQ?Pub}fHx2&#? z1YNz|@SkP+@+snTtX_ZiIk<E+B$$v-x&`f{+;w-?1oku2!56<Ae7dgs;3~+%I&po_ zR(;UWsx0XKBv`}ByK-)zRYmQ2=+?!2D~r-=pAC0tv9Evs5!)S$O{Wgdo)8b&hFEp} z;35&ol4Cx-xsB*~-tW}wu2XN_T5Oi=S+g0mxB_&&(%f3G9@gvmRu(1KF!lyq=Uuwk zRTDCn4~l_BQ=r3m;CmKk_%rTvf^DpuIQ7#m(1APa<!6*_Jq0OCBc^_8g6xhC>pH)g zW%{~v$EQC?LOE*qI`h)a(q)hVPf)|z6I3w%-Ls_<vSSC*UW2XvEPAqJPnr7DPu^3S zozvq%El5zo`1P?6=uTX<j`@xHprOc9AH{``r%5XJ23YQneGb2GVcU1e=?bBipt9ch z$aHWmI1${rtJnPgi`DJadq&XtL`91vUhFP>@Y3T7_!dggFu`LXtCvmR&1FA84hDp8 za|9h$rgHb*{qOrX9&-R4R|Z-GF@FkE<*^4^bMN0{29?;Ly?ExZvq9BmB|t+>puTFo zi3NOQ1hfZnDtKAhac}T;%Fmhye}kBiMSNBT;M4oBNmjChSFr8SdOjC)#60IG(DmA< zjZv-y0&jJkvf11A-of8_G8Le87OI*Dzs>_4L63PLb^`cp{BQsE#dC-;*P27_#oVvP z`h%kelp#K62Q`3Bf3#bI<;Z*UW{DTKK>?+5>8B}Z|L7l)82O+-vqV4((m{*;??F%C zF9F%|L;U4A{fT_rzC%tbJNVo6ygMUkza!|H!nx@yb~~Wnp8{D`_tRAK0H}=wS{Mgf zAXgXn%T!YUv|hQA>&W^E(38oFAp0mEESVy40o-+7`Z>(2VCuPO&4XKwKZJ3e0pGMT zw-ItzCuk?-rXWA?{jy?a;F}`VAAk-JIsQQD?#X!<KsN!}$cx|cF5i?Lhqwf3Q;?tf zly@xrGvH&Fph8Fxw2tsW(e+BuvaC-HFJW`NptJGQr<~irUmbKQFzC3L`}c%n?DT8@ zi0s)8IyvQ|BxsBcCCOUG23YQ1{S;QZo_qg}bDjM32TLljY`~r~<+AA1+txLiuf9Q+ ztbLt*ZoMaHTW-SAIe*^CLBi9!AQNer#LMz-fYtwB^B|kMRu|QIa)ItYjd|l8UyixQ zDOks|va6`N8Fa-7sMq#DYi@ll{}V~j&6S>9kab@|F;<A-_M#v^{waCpGGEqh2KTDM zU;ddj<pFqH4Kjj<);_ni1#K#_?=1ONY6%{WVAg)h5_XDnjYvFrlj*bn@v+aR-8SQK zxwuv&)<HPf^)^FpozAtYK&9yq7M)|RwJ6{z`CS^^+wkJt&&_>m#x40Nnd>rzUa38< zTFtC|^toh3;J!(HYo?s#^)k2SzIA11m<z8B(|@j-?INIM)HiqBW-ArA?|hf>mD&ko zruhx&4_tp{zhsKkd^l^#?DHq8q)*RZv-ESNiQXmF>$`3+<~M9VcGT<R>u(3xAAf3; zJmB^5vHAm;9UB~9J_8NWKR9=3MWg#;(15%hNCRm3U5y=s|1mIQUZeO^?UyIKI2o4~ zNB>`2n0caBu6H%)y2?F__n6-ItS<VuX!8N#4}!1!4j6w}RI^kjQgflx(-`ff?3iDL z^W7h>nsd55V7nps`qSgz1Iin0I3Ky6uxC7fp!?I*qLigA-&fpuW6?2H=uT~I)`N>X zyx2DfHD#BY2;6u2$M@-rfW1Ra)IIRQ9rHoH^|IO&)b#ega<KF><6pJ@1-6Y+6){4c zCI8s$4oIK4FLm8|?*2yp4~)lV_Pn0Uw?_o3=~K`l_9^FLf22TNm1O(-qjg8V!amk} zlGl%WZwUIci}O>D0Ftx%SKPVt@U?T-?g9^98>#Evd}~DP_eHX)T-sS`-SK|mJBHH9 zg)fi){UiK=XOBj_r+b*!kCsb2x##EfdRxEyoDnk>G@sX77<w+gSgeA-rlC;nf%21? z;286=EDAcrt7R_D{j28Xzk5sB<fgkQ&E5Y@=L1_scFg;df3r3pV1N9p@}-rj@)3Vf zu>RRvB5J95)zs5`{>6ByGwMrZzk?$AXH9y{x?QJAm7-SMcAO{oiN|1GO3Xgf51)NN znrc*KEj6#QEuG!V`cYr|dgO{m_s5^61pVY*ssmmO`Dd5x5m1DHe0~6yF6K<pT2NBl zt{<GPUB1l=<nz_%;x$<-HuHEaDNgS#5xCDRpD^d-jx*|?;I>3+hPGZhxMs?!yOlfG zK&it(dudaoX5Boo3idtiyO=+)JeE|fV*(|m|5JBY)mg}#`}%Ls{=a{=?6+U?_vOEj z-#^cu-@g0v^zGN*>&YH|`SIVuxcUA5`s?-X{@(uhXU(5fcmDRj-w(gMfAMMjgLf4+ zw+~)jdTjpVm&d*5-rgVJ&CDXgz`(%4z`&r$pd@c0aE+OPAxngTL4tvSp(01WJijQr zI3qQ+q*xzJ8R=Ez=Jd{R?7wUv(E9#QmqxkS#^Ys^LZgqUITbiKaqMWnmaBYgiPs^$ zKh<;AX>9v0`7=QCV@cebi_a?y!}9Cq`k!1Zw><B~b<N$UFS~rdop}9gd+D#Qr*~gZ zo#j(}`~9uU7haeCiuQjw`|Rqgmks{S=!-hvDfnV}Y1N$Pr)#JA-d=3C{qo&{i+elY zw12#F?pc{l%H`B6Z{N)-u}if%A6k~T+pacjhkUY0+uwb+U!~pnXY;zQPE7La%DHF1 zEKe?<dFgD{^_e#|%}crRyY$zzyXLJvo6o-b^&)?6-K;gQQzB1(PyDg*PR)(Amur4z ze|xn5{kqLROw_FR=GV`AcK-0?>btUin^)el|GI4M?Xy`j^Mp3bq<NKRPJgtaI#F#~ zEZhBAfs;4;d?`HK75Xh>-ZoSF$t9;%ru=j$ORGG+hb>z@A~BLNwRP&u>Q!Yu&m%$| zLsQQm58j$QC-JJ#s%d4*mcKe;5teE>CAKYp=aZLKj(&`ZeXTpEYS$!MpKx`M*c323 zG2mv~O^E`@2`w^Ld=_jnJCo?Zs-ik=%8Es2KfSM?)M$I4YqrXE?Wl;Im)5V4eUj?1 z!{yLc$u`l_4O}@oVTmWX!nbV7<bI<y&)L(sN9&Yq19Nm&aYSHV@`^oOR*qb6wEQa8 z%?R`p6n)C%o#eK&W95n?lV$|!C2Ms_J~eA=jSsZ;Xe>KqxLC&cRIZ!m1rLj~Damd- z4`(cxDA0Ub{r01{A5&zPI9MxsWaYM9{B$zWBxU7{K&2fun<SSyC@Y@P;N%do>|xuW zen!JHMS=HbPZz_tB^#DF$!<}+J=>L`y6MLCG>I(}+HxMaE@%~Co9DHqU20RKc%!Y0 zM8ztB>AW`&bRUotlBh^Dj;%akeQ<(0C(oV(2e{=fI<m6GB|KntEeU9B6ql$_RI*yp zP-4gN{OiY~(JQ+QR2vmVS@kjngM)ZhY7}~&o-rxO!?o2r%5_RelE>2G4<F)B&doS6 zRb^R7l9y_KfoiAb)EPmZg6mWg`D-VXZBF_3@L*?Kj_0E=Ez!9BT^m61ta0^%iQYoN zTHYN~bd+;Ele|({b3MZZw=Mbbd+n)(x+`rCO%q|+Z+7^^G?DHNs;!;jGPfo@65O`L z;A7VUUCpiys;R>M=8vX2@2Ou?x=1lTh)eVODfN54S|3YtF3qTyC|k5}-a0Ltxv%Ok zYf8-Vd-%oqze?nhFNd2BzZ9vnFPY3TukY=&Gxlwk7Xzxb6xd1?AC?%sZoazI_C^2O ziWk2$Enn`f45~|yI{V{K?dvLA{^_r$Rv!DFIIsQl{oRZ6%By1D@vQotIOG34o%t8@ z=5w0YGoY3QSDwvTqr%9*5WvL1AjrVLkdmJqP?TR#te>1;lnN>Y=1w}Db=W|__4{A0 z#d~#6FI>uTyK!5}RfmO|x_3D<Cv`;et8DfAx9^j(A5&ba<*xE`cO*mf_*Z%UJ5ZTf z(&{m7F^h@0)2rB4JG0ea1Xb+}bR&{YLSAsl>Av1QDf(1zvw6DC4uym2f=g_lZgcuA zpmX%bXIqVDn_pZAyd$+dqGRdp#M3|cZ!_wB`opNr;j`gD0lUvz7Z>HA51g6xv6q=y zj$deClC@EjJeZ!(u>4Ps^ZFQzsXXr7o3Bl_*vTK~cWhhk=CX^?$#qTpudRRoch$-( zcUM|0*d;Qjh+Rh_eEN?~9qrjO_+M8%TFdgWr03U}+dEENU{~s7d?E0CNsZpTcr)z~ z@k+-_bCx`*bLN_*pAlG7uDt%&c?VOzoii(#cy_M0cl^tb8l(%mEJNxU85ne#85jhR zgEX<A01~2(`L_*tTHpT`ZF|>yCqbgkAuMo-Y}Jho44$`*)h}-hh!k?U=JD;dx^>6n zHSLU_Uv}11{PGFBD{B@0WyT4%_|Ds|np|Jhu3USs{c+Flvd^I+M;>Gs986a2^LXS} zA2eG;No|{ltKWx^&S{qGmZ)4#V&8HicfPf+hirk}vb{XM$-H(fpYJV9->Gw2Ux>eF z-?3k6evOZ&?0EIF-uZO`n_bmS+2nJVU)i~(^A`N`{Cjwoh=xgwjU11a>T}<z&D%2W zO@6vqZ|WTRtJfM1m-1cXwQy3^XZ_k@a3+V@Y+{)0`YIEdAjfaz^KJ!+>~APIu+oM9 zWK(+Krn!8J%O+a9yguRH3fGjNBCWoncXP9UN9wyfi(PlKY)q-`-Ld)2ggP^6k3ag4 zCOV#x$WEAXIacZ#PjO4hABOs4KVCm7eWCQjCf0K0`)m6z7T;-QPp~=vYuAcxx+-fv zfA+~Y6ulDM=9l+jee{RwHRtE`_a++0`Q2WhKfT+i{;|`XY_B-6*#C?H-i%Bl%!t}B za`&G65Jm<DaaINfVFn0bWME)e(g<QlJLl(>q~?{xhg24%7RQ3(58Z%I^F439WME+M zW@2FA2I*m7V1O9V%E-VFUzD0ttgi<m(KT*ad|K}V0|SE_GXsMZipJ?o;941?7gRSV z<>zOE?8a(Nb@Ju%a7G4(FeU~DeiU=|GlR`RHUZrc2g@?*7BDd|EN5e200jz2Kg9LU ztYFQ>8Hq)yDZwR0nR)5OAUB|!Ahk#8WIrPVLozc1120HB-2b8+5EDu&b5gNtwZHLR zaU&}OL$EY@*q`SFYc0u0%}vz@lZM#r-t_m+>-k$581A{VftxZQ>mjb+$PYFE<aKb1 z2HLQ}Y6^OT2w_UwQLLtbn@8y3fL@d$Oqp~NOTmh+9lh8<Xn%SNNjtO<3GilR1F7I< P;AS|+#K17=42TB+CA9ou literal 0 HcmV?d00001 diff --git a/Projects/Tofu/parameters_tofu.R b/Projects/Tofu/parameters_tofu.R new file mode 100644 index 0000000..e5a0bbf --- /dev/null +++ b/Projects/Tofu/parameters_tofu.R @@ -0,0 +1,39 @@ + + +designpath<- "Projects/Tofu/Designs/" + +resps =360 # number of respondents +nosim=2 # number of simulations to run (about 500 is minimum) + +#betacoefficients should not include "-" +basc = -2.5 +bcost = -0.05 +bmanufEU = -0.4 +bmanufnonEU = -1 +bcultEU = -0.2 +bcultnonEU = -0.4 +bCONV = -0.4 + + + +manipulations = list(alt1.professional= expr(alt1.initiator==1), + alt2.professional= expr(alt2.initiator==1), + alt1.expert = expr(alt1.initiator==2), + alt2.expert = expr(alt2.initiator==2), + alt1.domestic = expr(alt1.funding==1), + alt2.domestic = expr(alt2.funding==1), + alt1.foreign = expr(alt1.funding==2), + alt2.foreign = expr(alt2.funding==2)) + + + + + + +#place your utility functions here +u<- list(u1= + list( + v1 =V.1 ~ bprof*alt1.professional+ bexp * alt1.expert + bdomestic * alt1.domestic + bforeign * alt1.foreign + bdamage*alt1.damage + bprice * alt1.compensation, + v2 =V.2 ~ bprof*alt2.professional + bexp * alt2.expert + bdomestic * alt2.domestic + bforeign * alt2.foreign + bdamage*alt2.damage + bprice * alt2.compensation, + v3 =V.3 ~ basc) +) diff --git a/Projects/Tofu/simulation_SWE_after_pilot.Rmd b/Projects/Tofu/simulation_SWE_after_pilot.Rmd new file mode 100644 index 0000000..bc2f390 --- /dev/null +++ b/Projects/Tofu/simulation_SWE_after_pilot.Rmd @@ -0,0 +1,332 @@ +--- +title: "Pre-Registration and Simulation after pretest" +author: "Julian" +date: "11 Feb 2021" +output: html_document +--- + +```{r setup, include=FALSE} +t<- knitr::opts_chunk$set(echo = TRUE) +source("functions/help_functions.R") +``` + +# Introduction + +This document has two purposes. First, it simulates the data for the preregistration at OFS with ID XXX, and second, it describes and exemplifies the analysis of the data as described in the preregistration. The document contains the r code for the simulation and the analysis, and provides details on each step. The user is free to change some parameters and values to inspect the limits of the design. Here, we simulate the data with a rather small sample size and realistic parameter values. While this html file is static, all required files to replicate and adapt this simulation are stored on github under the link + + +# Preparation and base dataset + +In the following R chunk, we read in all packages needed (including noting down their versions). It also reads in functions written by the authors to automate the simulation and read in the data. The functions are stored under "functions/help_functions.R". Further, it generates objects which will be used to store runs from each run in the simulation. + +```{r prep} + +rm(list = ls()) + + +library("dplyr") +library("evd") +library("boot") +library("apollo") +library("sjlabelled") +library("readxl") +library("httr") +library("ggplot2") +library("tidyr") +library("magrittr") +library("psych") +library("stringr") +dce_data <- read_excel("simulated_data/simulation_SWE_after_pilot/tofu_main_swe_sim02_tofu.xlsx") +dce_dictionary <- read_excel("simulated_data/simulation_SWE_after_pilot/tofu_main_swe_sim02_tofu.xlsx", sheet = "dictionary") + + +source("functions/help_functions.R") + + +models=list() +results <- data.frame() + +sessionInfo() +``` + + + +Next we define the utility weights and the utility function. In this case, we simulate data for a conditional logit model, with no observed or unobserved heterogeneity. All attributes but cost are dummy coded and parameter values for the simulation are derived from a previous study. We also define the sample size (i.e. how many replications for the simulation) and create further objects containing the number of blocks, choice sets etc. + +```{r} +results <- data.frame() + +# Setting the beta for utility function +basc = -2.5 +bcost = -0.05 +bmanufEU = -0.4 +bmanufnonEU = -1 +bcultEU = -0.2 +bcultnonEU = -0.4 +bCONV = -0.4 + + + +# Simulate reveaved preferences +respondents <-length(unique(dce_data$RID)) +no_sim <- 100 + +nsets <- length(unique(dce_data$SCENARIO)) +nblocks <- 1 +#max(dce_data$Block) + +setpp <- nsets/nblocks + +``` + +We will now create the base dataset. It contains all factors as well as the deterministic utility values of each alternative. This dataset can be used as the basis for the simulation. It does not contain choices, as choices are subject to a random variable. For each run of the simulation, new random variables will be generated. + +``` {r} + + database <- dce_data%>% + #slice(rep(row_number(), 2)) %>% + mutate(RID = rep(1:respondents, each=nsets)) %>% + #group_by(RID) %>% + mutate(across(matches("._x3$") , ~ .x - 1 ) ) %>% + mutate(a1_manuf_EU = a1_x1==2, a1_manuf_nonEU = a1_x1==3 , a2_manuf_EU = a2_x1==2, a2_manuf_nonEU = a2_x1==3 , + a1_cult_EU = a1_x2==2, a1_cult_nonEU = a1_x2==3 , a2_cult_EU = a2_x2==2, a2_cult_nonEU = a2_x2==3 , + across(matches("._x4$") , ~ recode( .x ,`1` = 15,`2` = 20 , `3` = 25 , `4` = 30 , `5` = 35 , `6` = 40 , `7` = 50 , `8` = 60 , .default =99999)), + V.1 = bcost*a1_x4 + bmanufEU*a1_manuf_EU + bmanufnonEU*a1_manuf_nonEU + bcultEU*a1_cult_EU + bcultnonEU*a1_cult_nonEU + bCONV*a1_x3 , + V.2 = bcost*a2_x4 + bmanufEU*a2_manuf_EU + bmanufnonEU*a2_manuf_nonEU + bcultEU*a2_cult_EU + bcultnonEU*a2_cult_nonEU + bCONV*a2_x3 , + V.3 = basc , + ) %>% + as.data.frame() + + + + rm(list = grep(x= ls() , pattern = "^b" , value = TRUE)) + +``` + + +# Inspect the results + + +```{r desc} + + +table(database$DESIGN_ROW) + + ## For all respondents + table(database$a1_x1) + table(database$a1_x2) + table(database$a1_x3) + table(database$a1_x4) + table(database$a1_x5) + + #Only for the first respondents + table(database$a1_x1[unique(database$DESIGN_ROW)], useNA = "always") + table(database$a1_x2[unique(database$DESIGN_ROW)], useNA = "always") + table(database$a1_x3[unique(database$DESIGN_ROW)], useNA = "always") + table(database$a1_x4[unique(database$DESIGN_ROW)], useNA = "always") + table(database$a1_x5[unique(database$DESIGN_ROW)], useNA = "always") + + + ## Over all observations + cor(database[, 6:13]) + + ## Only design + cor(database[(unique(database$DESIGN_ROW)), 6:13]) + +``` + +# Simulation +The next steps are a loop, which replicates the simulation the desired number of times. Replicating the simulation allows us to investigate if our estimates are unbiased. +Within the loop we are generating the choices and estimate a conditional logit model with the package apollo. + + +## Descriptives + +Here, we plot some relevant desciptive statistics. The first table shows the frequency of choices for each alternative. Then we create a new dataframe in long format, so that we can plot conditional choice probabilites, i.e. the probability to choose an alternative, given a certain attribute. + +```{r} + +database<- simulate_choices() + +table(database$pref1) + + barplot(prop.table(table(database$pref1)), ylim = c(0,1)) + + +long= database %>% rename_with( ~ paste0(.,"_",str_extract(.,"._" )), starts_with("a") )%>% + rename_with( ~ sub("^..." , "", .), starts_with("a")) %>% + pivot_longer(matches("^x"), names_to = c(".value", "set"), names_sep = "_") %>% mutate(set=as.numeric(set), choice=pref1==set, across(starts_with("x"), ~if_else(set==3,0,.) ) ) + + + +``` + +```{r} +for (att in c("x1","x2","x3","x4")) { + + + +plot(aggregate(as.numeric(choice)==1 ~ get(att) , data=long , mean), type="o", main=att, ylim=0:1) + +} + + +``` + + + +```{r} + +for (i in 1:no_sim) { + # database <- database%>% + # group_by(RID) %>% + # mutate( + # e.1 = rgumbel(setpp,loc=0, scale=1) , + # e.2 = rgumbel(setpp,loc=0, scale=1) , + # e.3 = rgumbel(setpp,loc=0, scale=1) , + # U.1 = V.1 + e.1 , + # U.2 = V.2 + e.2 , + # U.3 = V.3 + e.3 + # ) %>% + # as.data.frame() + # + # # Preferences derived from utility + # database$pref1 <- max.col(database[,c("U.1" , "U.2" , "U.3" )]) + + database<-simulate_choices() + +apollo_initialise() + + + modelOutput_settings = list(printPVal=T) + + ### Set core controls + apollo_control = list( + modelName ="Simulated Data Tofu", + modelDescr ="Simple MNL model", + indivID ="RID" + ) + + + apollo_beta=c(basc = 1.2, +bcost = 0.2, +bmanufEU = -0.4, +bmanufnonEU = -1, +bcultEU = -0.2, +bcultnonEU = -0.4, +bCONV = -0.4) + + + ### keine Parameter fix halten + apollo_fixed = c() + + ### validieren + apollo_inputs = apollo_validateInputs() + + + apollo_probabilities=function(apollo_beta, apollo_inputs, functionality="estimate"){ + + ### Function initialisation: do not change the following three commands + ### Attach inputs and detach after function exit + apollo_attach(apollo_beta, apollo_inputs) + on.exit(apollo_detach(apollo_beta, apollo_inputs)) + + ### Create list of probabilities P + P = list() + + ### List of utilities (later integrated in mnl_settings below) + V = list() + V[['alt1']] = bcost*a1_x4 + bmanufEU*a1_manuf_EU + bmanufnonEU*a1_manuf_nonEU + bcultEU*a1_cult_EU + bcultnonEU*a1_cult_nonEU + bCONV*a1_x3 + V[['alt2']] = bcost*a2_x4 + bmanufEU*a2_manuf_EU + bmanufnonEU*a2_manuf_nonEU + bcultEU*a2_cult_EU + bcultnonEU*a2_cult_nonEU + bCONV*a2_x3 + V[['alt3']] = basc + + + + ### Define settings for MNL model component + mnl_settings = list( + alternatives = c(alt1=1, alt2=2, alt3=3) , + avail = 1, # all alternatives are available in every choice + choiceVar = pref1, + V = V # tell function to use list vector defined above + + ) + + ### Compute probabilities using MNL model + P[['model']] = apollo_mnl(mnl_settings, functionality) + + ### Take product across observation for same individual + P = apollo_panelProd(P, apollo_inputs, functionality) + + ### Average across inter-individual draws - nur bei Mixed Logit! + ### P = apollo_avgInterDraws(P, apollo_inputs, functionality) + + ### Prepare and return outputs of function + P = apollo_prepareProb(P, apollo_inputs, functionality) + return(P) + } + + + model = apollo_estimate(apollo_beta, apollo_fixed, + apollo_probabilities, apollo_inputs, + estimate_settings=list(hessianRoutine="maxLik")) + + apollo_modelOutput(model) + + models[[i]]=model + +paras <- length(model$estimate)*2 + + results[i,1:paras]<-apollo_modelOutput(model)[,1:2] + + if (i==1) { + names(results) <- c(rownames(apollo_modelOutput(model)) ,paste0("se_",rownames(apollo_modelOutput(model)))) + } + +} # end of loop +``` + + + +# Simulation results and Post estimation + +The results are stored in the object `results`. Looking at the summary table allows us to inspect the accuracy and biasedness. + + +```{r} + +#inspect any abritrary model +modelOutput_settings = list(printPVal=T) + +apollo_modelOutput(models[[3]] , modelOutput_settings) + + +describe(results, fast = TRUE) + + +# basc = 1.2 +# bcost = 0.05 +# bmanufEU = -0.4 +# bmanufnonEU = -1 +# bcultEU = -0.2 +# bcultnonEU = -0.4 +# bCONV = -0.4 +``` + + + +Next, we can estimate WTP values. These are not relevant to inspect the model, but help to interpret the parameters better. We do that only for the last model + + + +```{r} + +wtp("bcost", names(model$estimate), model=model) + + + +``` + + + + + diff --git a/Projects/feedadditives/Designs/Design 10112023.ngd b/Projects/feedadditives/Designs/Design 10112023.ngd new file mode 100644 index 0000000..d60e4d6 --- /dev/null +++ b/Projects/feedadditives/Designs/Design 10112023.ngd @@ -0,0 +1,44 @@ +Design Choice situation alt1.cow alt1.adv alt1.vet alt1.far alt1.met alt1.bon alt2.cow alt2.adv alt2.vet alt2.far alt2.met alt2.bon Block +1 1 0 1 1 1 3 6 1 0 1 0 1 0 2 +1 2 1 0 0 1 0 7 0 1 0 0 1 2 2 +1 3 1 0 0 1 1 1 1 1 1 0 3 5 2 +1 4 0 0 1 1 3 0 0 1 0 0 2 3 1 +1 5 0 0 1 0 2 4 0 1 0 1 1 7 2 +1 6 1 1 0 1 2 0 0 0 1 0 1 5 2 +1 7 1 0 1 0 2 3 0 1 0 1 0 4 1 +1 8 0 0 0 1 0 5 0 1 1 0 2 7 1 +1 9 1 1 0 0 0 2 0 0 0 1 2 0 1 +1 10 1 1 1 0 1 4 0 0 0 1 0 2 2 +1 11 0 1 1 0 0 1 1 0 0 1 3 4 2 +1 12 1 1 0 0 1 6 1 0 1 1 3 1 1 +1 13 0 0 1 1 1 2 1 1 0 0 3 1 1 +1 14 1 1 0 1 2 5 1 0 1 0 0 6 1 +1 15 0 0 0 0 3 7 1 1 1 1 0 3 2 +1 16 0 1 1 0 3 3 1 0 1 1 2 6 1 +|||||||||| +design + ;alts = alt1*, alt2*, alt3 + ;eff = (mnl, d) + ;alg = swap + ;rows = 16 + ;block = 2 + ;model: + U(alt1) = b1[0.3] * COW[0,1] + + b2[0.3] * ADV[0,1] + + b3[0.3] * VET[0,1] + + b4[0.3] * FAR[0,1] + + b5.dummy[0.3|0.3|0.3] * MET[1,2,3,0] + + b6.dummy[0.6|0.3|0.3|0.3|0.3|0.3|0.3] * BON[1,2,3,4,5,6,7,0] + + i1[0] * COW.dummy[0] * VET.dummy[1] + / + U(alt2) = b1 * COW + + b2 * ADV + + b3 * VET + + b4 * FAR + + b5 * MET + + b6 * BON + + i1 * COW.dummy[0] * VET.dummy[1] +/ + U(alt3) = asc3[0] +; +$ \ No newline at end of file diff --git a/Projects/feedadditives/parameters_feedadd.R b/Projects/feedadditives/parameters_feedadd.R new file mode 100644 index 0000000..33b0ae5 --- /dev/null +++ b/Projects/feedadditives/parameters_feedadd.R @@ -0,0 +1,25 @@ + + +designpath<- "Projects/feedadditives/Designs/" + +resps =360 # number of respondents +nosim=500 # number of simulations to run (about 500 is minimum) + +#betacoefficients should not include "-" +basc = 0.2 +bcow = 0.3 +badv = 0.3 +bvet = 0.3 +bfar = 0.3 +bmet = 0.3 +bbon = 0.3 + + + +#place your utility functions here +u<- list(u1= + list( + v1 =V.1 ~ bcow*alt1.cow + badv * alt1.adv + bvet * alt1.vet + bfar * alt1.far + bmet*alt1.met + bbon * alt1.bon, + v2 =V.2 ~ bcow*alt2.cow + badv * alt2.adv + bvet * alt2.vet + bfar * alt2.far + bmet*alt2.met + bbon * alt2.bon, + v3 =V.3 ~ basc) +) diff --git a/generatemd.R b/generatemd.R index 5808e8b..a03acb8 100644 --- a/generatemd.R +++ b/generatemd.R @@ -2,13 +2,13 @@ rm(list=ls()) #file <- "Projects/ValuGaps/parameters_valugaps.R" - file <- "Projects/SE_AGRI/parameters_SE Design-Agri.R" + file <- "Projects/feedadditives/parameters_feedadd.R" rmarkdown::render("simulation_output.rmd", output_file = paste0( stringr::str_remove_all( - file,"parameters_|.R$"),".html"), + file,"parameters_|.R$"),"_noheur.html"), params = list(file=file) ) -- GitLab