From c387ed10edb88145c822faecefae32dc0f84f7f1 Mon Sep 17 00:00:00 2001 From: jaimerzp Date: Thu, 17 Oct 2024 09:51:14 +0100 Subject: [PATCH] trying something --- src/tracers.jl | 8 +++----- test/test_output.npz | Bin 69358 -> 25674 bytes 2 files changed, 3 insertions(+), 5 deletions(-) diff --git a/src/tracers.jl b/src/tracers.jl index d590c00..bbf8d03 100755 --- a/src/tracers.jl +++ b/src/tracers.jl @@ -28,7 +28,6 @@ Returns: - `NumberCountsTracer::NumberCountsTracer` : Number counts tracer structure. """ NumberCountsTracer(cosmo::Cosmology, z_n, nz; b=1.0, res=1000) = begin - z_w, nz_w = nz_interpolate(z_n, nz, res) nz_norm = integrate(z_w, nz_w, SimpsonEven()) @@ -62,13 +61,11 @@ end WeakLensingTracer(cosmo::Cosmology, z_n, nz; IA_params = [0.0, 0.0], m=0.0, res=350, kwargs...) = begin - cosmo_type = cosmo.settings.cosmo_type z_w, nz_w = nz_interpolate(z_n, nz, res) - chi = cosmo.chi(z_w) - - #nz_norm = sum(0.5 .* (nz_w[1:res-1] .+ nz_w[2:res]) .* dz_w) + res = length(z_w) nz_norm = integrate(z_w, nz_w, SimpsonEven()) + chi = cosmo.chi(z_w) # Calculate chis at which to precalculate the lensing kernel # OPT: perhaps we don't need to sample the lensing kernel @@ -164,6 +161,7 @@ function nz_interpolate(z, nz, res=nothing) else L = z[end] - z[1] dzz = L/res + println("dzz: ", dzz) end zz_range = z[1]:dzz:z[end] nzz = nz_int(zz_range) diff --git a/test/test_output.npz b/test/test_output.npz index 552a811a146456fecee84482e383ee25b19a108b..477115349975a9cd317c63e5f82b380cd61ddaf6 100644 GIT binary patch delta 3300 zcmZuz30PCt63z`#h_cI?03Oxoz`!ukZi(ed;!0Y9Ty1PNNN)%##g;_3Sr|FR~F&xXaR}C&PSCl(fR?$$M&vw zdIRAwzh8OHM^JI%XjK)BK!a9qw86zGkd?|CNcIsy4<6p4oYJSGXwRGB3bxYdh22t9 z+b%iuP19A)1Ktc|xh+n9GDsbLDD6(VJEMilWRlX$i*!+b;87E~LIdQvSXyz3oe?@< zVWl)dH%>&7PxtO`qq5MY$fY|A>UB_gs_%S_a#dvXylbo6R!P+3LJ_|gL_zZ&zqqjf zrzwyx)w7Q7>IY^2g8suN?!oevPuDV>et{#lCzIAGHH270<(JHj$1SI9P zjPBTX7}zrQ)o-38!^H8T4Jr~p0xv%%Z@ZOiBq$88dfzJH2i?CH|2KM<8+e%v%H46b z1D8#^7sj78he>O%$o-f$lwLNTzQ&USmjAN)&j)Gn(=Aqf_0CT=a|_R%nb$y;xzFPf z1|lWx@(`Iq9p z7OLDJm0aOm$(1%z!sn?UA)PVm8%gtZJG{F2=6HLnJ9nX$1F8M>^;0g6l9lm2Oz%X~ z+Jn6>_P>RiK8F)Zazo(pya;Q&HVP5vvn6MAXQ0QH^=$oi8md2ec3^St0u+(Q+>pzb zL4_e1qJq*9R=E_0be?~*7@B+rRSLerfmJsl;Z6F+-j6M?HDcygj$`c`nS9Y}o6a#4yfgI9NaTXyPcL9xjFtB#js!Mvt0b#)mAbGJZW zoi&f>F9seuBZ(J#xN_I%ItYu}m)r?lj=8XBxtDV*>_n?npZY~?*bYld(ZI&qP zZJGkLY`=#^ZX(D%^IgK7aWS;&LD5|~5osiGV9YAbUJgZeG;gxARz|l~Vmo@<)sP-o z#L}!a(F*RnEAwzBGJ6yH8YFEWeXQ_KE=ZsjTe$ZdX4T z&uojhoahdHgQM*-C%*$p|Hmln*)&ENXrs98YLe#y5W`nM7@>+(lK{=t5wj3t#$i!hh@=&s{3Zw zcG6ID@@3~Pt5J9xt+iPsVi^9@?0x36<};u>W@ON+I>9reIO-U)3>=gX*$oAkz(Ax| z-hgTrWZ>PVQq{gtJTWEy zGj;L6Z7X2t`c$Q3#-K8|#KY^N26%oK*K1@e4SHR>ULFghK*i>q^;SV2Y=-1AD^E%B z{)I+buz@=sbzOzq-a5;0Vb>dk3}ZO>G3T z-IBhQokB2MQOet9I>Bst>(*t6#Ot6-!>WeM31*W`ZM01#n5~Yk6}p6Aw!mzc$YO%o z3TsoVoA}JO?V1+bwcZrMeH8rEX5RX$uV6LEz2=W0W?}A^lt<+*<-r@yNqx?iW9aE##p57ym1@B6#oF_ zMA3;Uzi}OMDt+976mE>0sPpeKro26Qn~6IQe=H=%mq$#VlO&2Mljf7-3@im#oKP1y zmdzP*5jd7$D^e(-!Qnpk!O`p(j2$t0P90qkOqx zaZZb33~DQbQ!9$;2;A~dW`SEQ4bPg?7U(Vp`79EoM|Dx-%!pzJWIPE|jQPsy6bd^Y z*BFx%SeBS~#0tKqR;Y9Ki(yQGmVbPcK4b9AX>CEy{KGpNqfq$cn=c6yk|h%lVyfaF zrp3u+V1|KF@vSq4W%(MPX=5}O7PA0|_gv0qDMp$jOVihLM?Z!sx`eFIM7HO-i`_m{6I75y}gzIH!nzfy4J56>+l;Rp$i1 z9q=me#T->JOk7ZBKF9O-CDquL!=qwkq4^9y`}y2Le_GA)5yAAyJR%nm%huU*(idP_ zf@S@F)~3&XCqHM|#P2wF2rF%VWgQf>8zPBW{*azyWxG~T*)RK%gVBg3{#?sZr3rT^ z=UYN!9_K)S`Bf+>*`^z3olDJeF5%4Spkb0-Z^+1c9@s2V@L4C+12gHd8!d&sN z4mAM`f~sv)7uOP0ZN&zTn}iVTOj-6C&MF$Fz?gL7Z>6to5KfUK#xTgYT6A~93HllwLCjL*O(Dg_)A=zJDc6tFXqpxw&N9z#V;7dJU|0pD$z>cQJid{HA zhdr8GpU!f6727zq+FG2wg?*OcOe&i58~c1|Fe}KC2{eg#iy&vr*i;zLFkgHUGJa(#J71 z4OvA!fti?py_o^if8K|crI?WAFKj}QY~uf5(;Z@KI9WP6Mn__`EdLuWm1kAi-JHzI zoN5ok(A8I!@^3$;HSt7M%7j>$v@>2q@5+8XFHXOV{x|NyVm;fxFx=h$!uYSgo!O1Q zn6IabzLCFcWe7J#3+Br`s{|V8cjiUSfD%R&8CL#+Sbz9q6(5eok0(U+BxC z{>`SHou!R~yS=C7O_!Ta4mS^3b^34K2@q*!Ui!I?eve1yW<1+OJ#nlJ^_jzH{0xs7 zN%%+9O7Zy`Yeek7yzAHq`0_72N5_BT{Z}^W`Y*q;K*7jkG;5-XIvhO{+g+r~2}H4( zuZE17T+oD%0w+_$sej>(`}m~%2aoFS6|-}Urp0LfS2oRN$t~E!+bYo!8j6;sua#Fz zF!=X=?x-Z8qW(BV0CF06L(6KR_b-fSA?_G~zsStaZaMuQFowU${L49KPpALGrOkii zvaP~i)+@jt@ukS@v}XmQ%VqatKacq%2HG*-sL>nfOi9Is`-GbR!j*J##{UP7`fs?y zo$x`V2?4Te}&84 zbfaJY0s1HG4@o=>$qWC@8R5C575PIFUo~ruHy@ICrt;LNU9dZ9Rcd}rcSs`fc_VqY z{}+im_U1j&DeN54bEmEp?(%B?oqH*_t(u}G)mWlgEsjU0Vs-Ji8EQNxAI@WsA$4)a z6Khxp2HC97mD^a|%R31NLVvJ?#-uTN*NJi3_ctDKJ*U9E;z^`=ewG%e82&7?RG$$? zP5n9aXpq40~+b3*RFc2?Z& z(cMQ}1O$H-n!`|cvbpK*aMRB6_BA)R!&t|wsHprq@lzZQ`=2XEE@pLka3bFV5Y;>y-|2rWqA8r5QO{APGM5OodC#d;8piQ(Xl5Sdt+=;?@ zOT#PBwe1*lUW+Pp<}B?l?^-qDk80sJ3;2R;qiyj=zS4e0di&x!#wm3uNUu$gL$Cqi zOg48&a~n~~^Mig`sU{@o^wCcB-FNhDn^Zecyaf$hOmv;h`hhxZ8lps>wxXMl{^&bv zwW0CtkqM7>059oGLOPX zHpeJ+p=hc#jp4~IMCN;yE_W(+4_B+>AHGniOyvSdxZOPex^j0xBW@#CW?r5#VlE;5FM-xSVGI**~0I;9gm@5hi}=6>Vq#u%!l zPTW>f8b|Hoja9Gij3ZygG#*UpI7&Vev%j@7juequQ@_##!hd|j_iV~@0^J{Z!h51{ z0xk7^^tiq}fqs4q|IsfxiEh^Fv)s0sM2yBGEV9ohk+@i^^|78wM5LWoL`5@&-pqeE zB6Vg8T|E=*>U@6+oj==Y)le{n%HIn5pPimUPk9PTCRwKum4jYM%K2%ed^|VQ1wV}> z)W)=}^QZAB=jhWt_mOFIvsJ9=8TAb6NJ>Qgm>IOTwOy;~ID>L{mg5TJXOMn#9nZ-x zGw6Z0;qc_*47$4UVJnMi7LCSP)ZwwS$XIe-$;o~ekzVM(@Az~Uz0h~*d{{h-ihZl) zi~DBL(UWdXTL-hKE)b`x#XE;K?yv2<*Tm1E^Bu@p#&!-F7eC{y3YtSnpQX&rv*wW1 z_6RxE*Ew|K)qKp)u{rd!G*|o8pE)#<>-jF2Z63u&-0k;NoJX@t6|FwH^GNygi{GIR z^T>(jZ2N2fd1T#6zS5E~kLn*_A5p%WM+xQUnlIJPqtchkt1o-z(OZeY8^7@L^N3fY z+|1@+9?@mId)H6BfWG*x(;0Iwpdyk%fkCMSM9E!z;Di=X*OF^8rOpCU9@BdD&S(Ki zQ)Ij_vR*)O+Iyy?E(^$vlB>D?&H}0k3?~cqUqG79GjH@CFQ6nblGB_K3urNy&3`L) z0dce~_I4yKAceY6N_=tZ0-}-hDvwKBKqu&Ph4JYN=oT}!>{iAC8nU=oW_?)KSu%e! zJnT=cVvN&1949I3@{l;f9;Z!e(aFB+X@&I^dU?!{|lMr*Ya1LGS+5RdKJ%_xA=FhI*nnSfNZYsqW=g^t@ zZYpo_IrLttCr_Si4vnbfaL$a)qHCCv=;%+g=#-Z6X>9B)BGB+EopqW;^K1-uAsVwt zjk@xo0Lx*}N*)vaco@VN4~b#YtDZr#thmDQ=ou9Kg{jHKW(L{bxJ0<6FoQCJ*YDju zI)emr&NC2qOrz;GnbgN=)2P*}iu2gbX;krR`A(?HG`gI1j*X6Z8o9Rq7Ef)SLXLMl zip1ik&;;X~xpMO0#U4Xpa zEYEkLNA7&Q{1>~>UU*K97j}1?}jFu_oNAFME0v+Y(k60L8ccpzoE9u0ww{}h!k%xP=VxT&b=%CRF3j>#&sJmmLa>Lk{evh zA5mvdWazh~52*iGbmHerC5V;sPV6*eF%pnkxFb4Ui1M0JWo^F{pxYt2W)E`SBMKFB zQuF8eD3+s!C<^lqk>7JA;?#eOwr%th)_il2Zl@w%CBEtnvRufp$Y6SnR8+eTQUkKk zNwG}b6{1WeA($CQk&up#u`Rwk4=>SXB=eH|?h8b{YaMF)f z8R``c>9preL`In${*tQkh=>ZMceuwQ1D#4Awr|mh;@k(RT(u~aMcZn8u08@OY2kBQ zt-8Ze*SnkNTrP$pC##a1_>ITt>WQQ0xoU%vcE!<(U*&;_Al-JNt2Y4U38~pWWP5H^j1L+V)?Tr#zJu&ub@hC5@j~6tztRp}@IVHYI@E7X z-4Nr$p1t8HSLD-=vYW@8kv!ggV4`#82EyulnOe)*BMa`g=Mu_nP$OBiUCF)cXm^k4 z6^X|+M0z))FiP1H5pr>UpB^(sj4ac$Z-$M~a{{A8HDLp!HCs5kRjiMk)8o1GALt>s z%*x>Bv=@ST#f&^aW9Yg>;nK8s8R^qw6pLOjkSee^ zW6+b&TcR3$DyS&zO`mU`BI?9Al-vegTiZ)NUb*GQlP)+VG0%h zq{}3PW=>v7E@0$G2i+{sI^20soB&DT))`KeqaCj3tj3C5=<_pI^BK_>ia%YhiS&rE zM8j`+ng)*u{J7-XxTsL_=6pcUF>*walO0?ca1?C_)6AH-A32QbwlJ5(KbV6y71pMs zzc7#N6$=gMw=u_K(m9#mu48ap4Dg|91+%!XEJ1T`0b?Gd52}(gn7hIPv3d9j45N{z zx{S#vCW4`mFyhJ}hV6I-KXFYDMwzFF{=^`@15;~n8ekRPhM5oeG;`Iz8DlpcU^QFa zfT7aFePpSu#dLkSG6^Qt80lj1^8Ni!7=wDkyhjorF_A_?#u|J@n5(w}e)@jS#~5lI zc{$LUgQ3KI7u2SCg{dgi(^_Ck$K;qw9#M~cj(L24gOr&%3B#Ok5;Sx-7IWgWg8pax z+ei%A3oN6|@h}YNT}yB63dZb3MPF+V4!}HHw{#=k#ADLd*Jb6|?qcrWiu-u>oF_&p zQZ((c;Vq1rUL;?hoD=4o6Yt)ONn1=r@-xMN2rJCUTf55R)D{>oT+AnV9V5)kXa_CP zWPOYzrN+2*(nSmn&CD@&XDy8V;-@3{9%oHVB&Tps&2=oMV2$|1CQcP|zvd$Y8-W6b z`RCfoK&>=pLu81h_>(ARw_%bqNlO6p^ty+j>jF1M*vF5IWQY}mza<_@bB+N+L4lQ3 z>7~JRmN_Jvr&3^2Zy!tB_(P1jMR)i6?h^t`Z0!XHev&=)oQn-w>nb{*M51O8ht8chDo_w4SsdiHXp0sVDVdT{vM$I7}+Q!3;C++ zCTD-uSnyJ3R3W99u)U=&v2T9S#_EQ8x#RG|lAY`7-II?~2c#_2f3UVZjIJ|Kw|a2q zW1Wzm`inR;ccocN{rAw?ZK7)$>P7;j3$(Nt^?}t#BrRVQ)W^nyP8yt&Rwu*fhz!eU z39AnkaO=Ng;#Ft=Rlx1Oz^WcC^^yHJ;R$t|p9hgi9+f(2q1Uma`$yI5s*8N_I0E&h zqMEj$Kig_gAHUvSe7CB$d_TIUjC@XQX7{2nhHqT$3(LBD0^@+%C-X&8k(ds(iuv}H z9hMd~?|fY*iTN2*RSId;$lsT35TFou*cO-pp zni|iiHvZ7|L^YvDO=8b8qtsaa`NNO$W1&L^idFxBToWHDAUKw8M+;U#6RW;V^yWAs&r{Y-$PeFx$Nr$ zt~6C5cM{xoRJfaFJ~>t!w5W2{!U0ghxXNd)knOi>eJc0%s+R*?+Eg@or7n~5)~m=lF|FNuRjG1r zDw=)pY>|pUNcyenmTVP4I;A0e*Sj>8iaOcOAf0%XN`&k4O$%2E@e4_4XMU)1nb*bs zfx%rBhG#pFtA0ym*MiHjDA-m-@~!0^BRg{yBK1&~R3UwpILG`~+tX)N-p5@?e=LGk zshAVKVznl((xf!yki{XUl4qz$*&M{9Qo9xIeCIi%%KRTK9GMW6%EgIweDp}=5fyRX zv={}N9p#6#6p~|IOUfGGEcWt@naUzsJkkWhPrSviV~>%7n6apjvW`OobHh?Hm6d5s!#x0QTs>3dJ9 zFDluzotlv)8B^L)Q*c|8>QXAY`M|^**QCT%*XE3htCTz%Qpe@pi}sSmRL`Lj8E{_d)9YgwI-77xOb^1cJDnAj{37peYS)M;HQAKxeVgD^BHOOIAAg@t z>G6*zMNO{#38K%ZfG)iS>q+V~S=!d9G1ubt&?yXfZu>`=)4r zUTQX9_p>7H^Ty~Ae1T$1xtQ9=?o36OKN93?)rpEucv|nPSHcxP^dRx;d}yyY5^03L_u5R6q&juBtW;M~I&b2-ub+nE+>bX*y~T=nMa3^~ly~h# z6@4!qCD?hxsmLL~>(x?At7r!AyA(r+6+=%L2D~raRyZ&)dQ7Q5ukdDRDns9ESV5G% ztDl6aU7>6v)3)_%tpYJylH%Fc4+^5{!Xu|qjzUN8(dm24&lN7WcE-4sL@MNBO3i3R zA1dr_zH~EM^iq)Kta^^5bKqlh5(P(ejebbngeZ z1Larisinr$@5tA@3@E>qdqZCRYRdBy+UD{fWzMeC(do)}r54*5q~PSooZ@0#&dSNB zvoc2{oDh)j+<9TodxS~;IbpjXVG^{=9;`WE-ep2gB zGhb|z8!5OpkzQOQcd1KC1%F(&NbbZA_FUMPOu2*j%D^?dc)1RBQrf$0kL4;=K2Prs z-<7M_G`fLsfFB!`Q%t_8(p~w4PE(r-Z;6zPqBhW6J+Hw?!hFNikE9LvrMl`kumd3R+ks7xhWMxnv| z{!DI?jE*mE-IgOnW~xqibuidR#y8PQGOxE$I3mu6%GBYPnj@t~Y z%e24#bkmSpLZ(yU<(X({4w<9Jo{a`;Q^_=vxo2sb{g%#lMNe*bGYZ zwU!gSGijE-Vzh5bs!<{BnQj}e!u(d6jCN5xvNu`Ub+(u*Bsx?YGwgMK?esnAxB*j2 zfpte|$^`RUJMpH{tk#zm+vLwle|hwVh@(MOTIEJ^uKh((Y2%@3wb@!0X&0vK?=ceO z(ng)e+Y|{qQk`i9KXHGirSOry)b^LJcT2sOYdxT-t& zKw9d)%khQ%G%l%2gG{MeUl?HvezW`@-u3bRAlc&}Q5A;ClB8DWTg6huN?zGhwi+D_kbIJ%US_}X5nxu{x z?PLGDx{~A0`^QQrkmQ4@H_9-u zFe%|lU+v?3r9(oWnz-*pJHAGO9yy1o+;}gM5v4b2%KSpY05?8S{vk|aw)ipooUM<< zKtktfZZ!vq0Ku52(b9$zH^*l#Tfr#_n+}G`%Z_pq`=y28Rn8;9Vl(xnUYth4PVS<% z^xfa$k;|;}P91aNWs8P48n5?=U(?GAv*D`~UmryK#K()o?ad!3U%|7!6yKEgcsX?; zLcCvN5MOV0U)-{-%=Ek64e=MgItAe$4aKjTb#_$AiB>BMxGZbnRN*!hi^d8b(I2jb9gD>Y(1;wil%b@^h!%yaiu z-lm8Z){ASVC_NE-x-fJFugY~>jJ)!tt(k$fm^e8$D)*t# zUpJC7L|L1RIwzVUML#9p+j>{vEBZ5HUuDYehNxnDg@w%*1JREk3-KMnks6{~bKAH_ z-=#%=uz9rK`o=E0v>^8CW;~f_dE)@_>9q|JF)if={LfL5DngQK|BEdmgButOdu6G} zGxf>ggRV@G(8QP1()Iqq|O+KY@IzVWKwie%ptnMlb1Tu_fDw^ zSB9D0{NW)e{16`_ap5+duy~+Xk&E7ekZ30Vlf~84ijk;L!t8gie;&(lR~H6iC>XC zloYz?Q4(VHm|5t3z#X!ug+xNO62IscQkDfDerfqhZPY7RU3kMG@La9nQa0sVdi8w4 ziq&jR$JIo^$}sNbih_rNcrzO+*`B+ug0H8D#Kd?^1TT+t9^84WDOhsuevX+c_( z99hmT7Qvk@_WP^Xi3LgX&WD)LuL#^5smv|X>lH9fB$x1-{USiDyKgdn@|}S4@&2n- zZV3WRJH>A|X#)f(#Ke(bsj~pLEQ|#58wz~i34KINa8jW1bar_Cju>8mX3VQ)pqN2` z-ERMpdk2BQl3{5XC)FIkW$1X`BwGi61$D$Y^K2zQgX!S7am6eCh`nm!0mUeO{g)>< zm=}Eboqv5bzIV-rKPDtpx!pvUe>G^ea{0G1|NIl{>QGuf{w3Q>^>_*@el5=>z69be zK3baV^PX=;_+AkUUL29ZH}X-{+jqYqE#%8%tHiB-NanNr71lhS8OWz~LuQ2dyer>( zWxArChemw3TT&h54k!7B4}{iJBE|TQHV+U_g`MCdcgs|}_~L-qt0%hkLems)T~lc1 z=e1Vef%!GtScy{JeO8;q_Pm$8CKyD(T@b!tMlbF8T}FU3pU$CL+5ya&IeIaZ8b}VLv}~U$u?9mfJO0vGODLrbj2^?@wvm&oro0Q%xRo zH?!Nldr9fWeU0)+8gY&>ciyazYmNFz?#VJj)A4Z;?wAmQ%q3GgZj$MvyCX;TxYh|O zUozT{aj8*;#aI%5<5J@>^UpCU;3_|@i|cJlvnTuQYlw#rK&jJ zZnvlSSY~nN(Z!Z|_l0xPUWmi_<2*SdKYwtuJZa7;liqjKlv9(_K0$P~GhK{xra;22 zWSfq2)H02ij&P5|N(P^tqR>3bLALtqo1%UL$8nL)aP%af!#g`qugxuvL->uV(l=LM zj_D0{FZWbyj$9gnC*P*ean#u`2)hl+a&XkDTvPC2;aD~yqNuRV(te0AB+*hffF ztSPX2H?h3!7vo@OkbaeM%z}vB#rD4YK;A4HYPxuv>rX42*pmeMuInXi3qRH##yY34 z2?|($j=B`UR%{kanqut07GQYw+Dx)8+f6SOqOK|hw&&)=X=NVlY%Ye+*;3+A&Shb=3;$ET-{ z19ez!t*c4rEoE5=VhfbkvY1$N);32E91d8PPlP@qcr?z^`;gmJhO&`G)_XB0@6S7y z?KWSnQPUWf5!Zo?KV#l37?EX3##>ieuJYkFq6Rft;!hMbs(6UvS=4O1Pxn>PupE<= zaQz{>$-Mj3bpJ$CA9I<&tpb;sDrTv)w$vM1nasS&H-3<2Jz>tvW#7D>-raij z5+ktLm#0K1GZMC~HmVHqFs6HQeV}5aWaPO#ZjrgR^q0(qfI#(s?tkC@x4vL5kEr55 z`htHXG4B7C#Qawk2@^P(q<=R;(XVB)d5yG?`CIj#c1IEP@fl_J!!0tjP|!?lkv97; zrHXl%bIPG!8YW=XGMnF*|3xI{}`&>WUQQ@Y`~-73~+!kqbh!xq-! z`lFnC)xWV3Z0nKP`b4;!e2Y9k_Q-IaK7H{gh-h#VP1v<3r%vE5oBZ0>z_8$cnv-1~ z_{NF*y+R;QjN!%6W|ngwjpoO_m@8m)O%uei-{|~0aZU)2d-eXr3F=k>oF;ED!Obon zTuGq)=~FM5a1FP6>K*#Ya8fZ2yK|Ypu+aw7Zl6MTuu7XsBT>TR*xQ6go1;|i*!wEh zZ+1n}vE8zf9iDgIV3&^ls0*3!#6A|uQd~O{g#98W)YYQ!0DG~T=lxBNJJ@%BzS@|W zyI@beO-jq^(Y3VM3wW|2y4HP{Z+#JrzhST3Foa3aKi9vUot3mfsn zGm%w!V7orMN?iXAOoT4eCEMnM^Hh~-^K3rQDzS`kI==^1-NwR;EAOGhBSZhHdjUvh zGH1TtDuAJb==IxLG0`93s(E(9JI0StTiA5FBk?1!8pLggbCg2V2sB@R zQ3@GB#5vx4Ww4hM8=9L{1}6}20)c2b#8A^x$h|Fxt|z~jHKaem!K@EqK;b7~Z}GqM zQ?UYuH@h%L$|@k8j)R68Qwg@gA>D>mm4NRFPko%J@fq5wW+U3_KEq01Qtcvt6)@jB zogQmj1(gENzpA~i0{r0HcdfHkkVbN%__A0vM1@;4A2?OR602=!MqV}ePDqiftW?9N zydu+2GBt49|8B6Idku^%ami5?*TBcY&y8tYHNedxy1=je1^S5X!hU#v!Gj9-T%oy$0nc(R-rkHf7=B6zi(x~UTT70^j2<{rN4s)8OPD3 zTi;I8xFT-Qb za9~=%QRpXGjc*Iwpk^w%T-5>*r^pC934g#O;k%mq8b5%gni3-t_ybD1*=6xHbw5Ce zv(AB^xD{AWmOWnNZw1`a_~*J)t#H$Lq$$p}6}~_A$ubXag{WQYtJLQ3&5YdGHn0~PT6R4-949ht>;PIVof_T&c zPx7{$oGd#)*5!kyvR?=2DoGUYz370fJ2;8z>JIoLEtHTv)&UC+xR9GfonUg6lthWY z6WY&o&TMLSLWHiUU4u<0tcS5TzX;bwckxFN!9_T$nsw8661MJs*<&JvyfZ*Ve^gv<{ zsHTx#e^b^2mCKJeo%(wq!kU{)WWNVyPDzB4arD9kf-uD)OfT$|vtqI>df{5US+Enn z7ijK{eiVOk_}q==3}tmMZ2Y;(F*4c<6oe6cIY;{N&~kpx(uuDR7-b_8rTQcr!$0ev--e58{6+x-v=+3`txLF`asb5sRcE8Ka{A57>tYb1L{5# z^j5ncu=BA!&W`rYPI8i14gsT2881Hf6_@bOmR0A&68400U-SL;qok6gQ%D+IDI0(7UcPeH|2Z4ZE|4@iF z2qeRO-){dJgh#zEW2N!zL$Jj6B!w0;1Ud%6p0nmdz*Onz_`!DwysvnM`aB;3v`0~= zP&EVsE*8Nf5Oz;c^K1MFT=MAV6nZ@Z ziY%wM13n(wdt?N!G>pKYhCQ)G?+7S0N(9!;j)2FQ$`R4M5qR|d!rey{qi~ZX7w^8v zItsO>ovmhKqu}5wC{nLB3R*Q1x5afvp_@>+CDd#buAXp|nspil-mIcTll!9(x4^eu z7d8riekJNlrjA0#fkSCT-Y9s2E6-BpD1_yI=Q8{;3etVSrC*0fA=L1SiPXv{doj8kmUuqyOFU!n+;>&xi_91}ij$p|IynX-bH!KUE{y?EmA%Q%p)fhu zMA7()`xyAq9{yGtFa`r-gx@7%#^CZ@JMZ|6F&Gb?WZEwp14h1$Qq!+v@Y_qBy0mi) zrX9MUiB63H!3Q#y_^mO}{BiTi0m(R2y}3ttjcFXT@#}ovl|tjd9ZG2_pgazp$)3Y8 zXU8F8;-2Xrqj6BS%;>dt7zeR>^2=Yn$6={@s+Rl7INWp!c8^FJ2Wh=biG%EMSSucy zv@IWp0_G=i-0{l}|zb1_xK3Dz#dutQG#QA0WGSMW^ zZ8k*K(@#RA&TmEq{z*LiKBbqKB0mXNGY?n^PfdcfJH%TaDlRN29fK=vCSi_-Sx4Gy z60%htkEH}nf@JiMZ^ZGFaQG+n1G`s~pk?MN_4UIfIOlm;3pE^$|0Ix*&@%~NZXIbP zo}GlXxcD{4-AM?m45Voxp8~y+ubg76Q&6LnnfFu_KLs+O zH+06HwM;?R4WpYoLsKA5=kv~Dc?yQ@n(Tl4nSuj`aXK}cX*jDknUlsn4b!GBr-<+} z(_otxG-`>PhQ0w3tFL;~aN$Rxk)qW!D1=?YymXre{4^6QdB8M$a<<;P5j_og18d#&?e@)SNKkpMD@`&3mb90* zBN%7k5}O)1kcI0Hup44y*v4D7kRe2`vtII+X9>yf5IyTKU)v;G;lifKx3o^(@deH&A%A z&BDuAu|NHzv!HWnWD%{+!p!W2&g3I=5T8CLw@o((;0-G`_~zjA3jh z=EEGEqrU6mRzC;Kq3sesy5}HSgSK08W)2vs4eQf(=Ac3D@&{`2dC1F;&bY%e50}@r zLI*_V;iZs~ubS#SMB|^n#>Z;U!$;?$JR;M1SSBny?%_BOI_WQ*hkfRu|C40nsnB`o zX15f5l{^n&ow6xJZ|7mG&aU)&#XJxs3e|Kq&jaU3N0a8zJmknKcf4Jihq7AhA=*Fl zAgW!s=SsZ*FP$FZlJo*peUY$wja>lzqG_7$3B3it)~zF-YYT9 zFM&J({yQ0F3Es*-6|*~!UjnQx$HK%oR31%Kh_7(arK@3Uz5cAU|s5%*a z#QXUY1PxVdbi7-FJ}{S1s#=1Rvs0I|TbJOw^xPBX(IqIcnQ!)AU4o>|8IdW%WuUT& zeRYa%86HX*oO;2#3_2~jRwQ!EaHV&G!CPY)oT$C}rttd9FhRAHchPzo+L%L93q6*> zGrOjME^rxAtdGw*#V*4Q2aW!Y%w-t%YBxMnvJ7P;%QHoF%fKC}5XRBH3=Pd4I{wqk z06%5rN4A%N>G3`{2k{Ery-s*h@7M~gxKu?vIk5tH8ILTg*;gR+SN?}h{uQ`X#+!;4 zl30O1@o6lE3M-&oEGYUAvjTh^jvXIPt-v!~L*>=;D==oIAS!II0+%CWuD}3_D5|JhfjhmAwKp49V1|R+O|)$V zuJG21n)j~2NUyU;*ysw-7G)Wf&#r)n+WF;~)fF&UHFV|MTY=116M?2jR-wtkRw|NW z6;v>nooeY;VTDmdeUo(+s@kGo@bRreUyfPYMX^=rq{i!1Jd#@lDx&MI)oQCCA3I&T zt+5K+{050)I;+4oCVAZa@+$N+Z8C+Lt%6+`RoO@DRY)Ubez@qg3U>$$*91LRfiUpd zuKE2{peVLYj0!yLzlCY54O@k1o!2imVpl=aPmYu~brsO0<{Q1NRrnq-z8sRb3SxL! z)c&Pp6^iuxF78$yb}VD<#p@3T9-fBgKURVNRk~eB_u&L@znuRtyb5t*qFO7{t1xGN zpHF0E6-un}g!0(SQI`)GnSXl6cZ^ObB3{GXW>Z_uoPX-In-56c>e>#_Ff^Q?j5 zuvBB1$Ql@&T`>A8vj+G}J7*s6sjLC3c$cxn$u-akmiT6_y$05lGrb{~)_{sS{m}=L zH6SRONL;mAgUFmH3sJ{4IG%p;>s9wPsA0wYi1As27piRc>H^ka*weddGh_{dA8dc( zjadWgF;kPv$!j24#&9n@V-4E8-m28)uE7Pn=Lx@x*6@%bo0B40v33|#E@v&juECmN zSVCCyp@BelSZU{>0VdsLb#M*BXG*d~r`F*0&)A38mezo*QuQEqYYm#}eYhG9)*x%j zadw+z9f(Zdn+i~`!(fRNfj-kZaDF(`5yrI++TVk<>x9-pFeA-sUuqqMr?l9mmGSEk zeW~>BHS9V}^prh}IJ*w_0wh!`bk`vbQ*?LPXdNWEF|S3gtwXynlg~AWb$IY~XD#m5 zI@D<rZdMSnJsA#)S=_ zziH0LV|X~BKjFUF+HvoFuxRXjZAmqX_p6+|~ z1`N5&m;Y+q095C|BGSGA7>|kV%Y7Tr5LmAjF}49lp)Kbc=QiMyW7F?HYa2l98f7W} zYXj2mTz0o5+JvLb53zBSo1i!7SO0~66B>dT$(GqR!C%5HQ-psL*rYTUti(4VntNp` zQGOGm# z*#r}<5zm^D%|lPsTGnPZ!E%r(USf3<$dZFuZFe`J+GN!{nQ#m4w2NPBCEtSbPCEW0 zbXy=|ow1 zQsi%m$}I@Ayh3bOw}l6nV%~;b>-3c8+}m(#M`xl{ zcpFL;oLB!yZ-ZywO0BrcHY~r_yJw8sh7xY4`4}zyHmK;$bG7Jg!^O5TuA?T~5dNxR zPQ_{)w&H*JIo;TXf%EEP$!^==W9WIj>FzeT@kAfpez1L*0eksGL$-m2BEHEwdK*qt z4JxNR+lFNA=Z-(q5677q*^}mOgOO>oyh_nFU@aF~?LTdUy5V#6xY}({C2PN3jsLz4 z!M#s@?{;j%LV3=C+~79oaBMa@PHuyFtFdL;;x=SBwN!R&ZiBM6*38j^Z5Ue08ImX2 zfp_&s-&{Ms14sG79)>aQKqXJER3qmOJXEWAb3|wdlC1OBRHSyG_0nw$7o{CI>3wuH z4R*kbSU0Qv%nm4aofpLu>h3^_|MoozqaD~w5Kl9^x&yOynM3jRJ3tdPwcP5u18*M< zc944SK*RdoFm?YOko=j?;_`S0ZuYu)KabpjBlab^&51j3A;hh4H*E(ju0{)qW$!?z zow}fH!48C{2(G7=?Encfi|nk~0jIq|Qi`S>xNI?jH&Fk%107ec`8xIQ0GHQU-=y&! z$lnb&sGHvbvJ+28cGq`6R*!j6Zhr@oS=zFkh_QG{N1KYkF366km)S||f{#3@U$nw5NQF!3)nRs_v~xs_@boT}v{fCbTsW*3@ePTt z2D`8?xN$7gVi)$vz^L1H7v8^loJe zmZ!ULy5jRO(u7^OPw}l4^I{j&JLufpUhhIGZ_G+Y{w{bv)Y9qvxC>t>r2~nocj1Ee z>1R@nyYTkQ{B_H=UFgDdDSu4r-33F7wu(Ohe<@r6VC{J<~oV^T>XgSK9k5`@3L8y?UTqxV5NF+#3b4jm1$A1 zk5XYkwMpxtNu;g1OOo}(Bsxl}{DLwj(GK@hi953=(KL~RHS6Le${tpYd|Wb##%)QZ z(>Ev4`GD(5S~ZjCX-2DUNW&z0Xz^$E;)_Z2&1-pj+q+2;+R4qx#@aQBLRO4MU|%QE zE#-IjgNG>x{2+~ zZyj^+OP@lOMZDM2&rP9KD@NU#3#L%dw`UW~SErEO#o~i{H>S|c`=e8l)lF zqbYPCHuOf@vneE>n(${;%M{`?zbmax>X<^8PV*Ui_fDbxT)4Jp-Bzp zX|(J7xTlBqG}1`gc;*~&8VNe>scA8tM%tH?-muzEqp#kn?+sn1kU@u?q}M!UNUGQE#ZBNj73^Q`8_!ImSlch}_bi`5 zesv+=F0GkCcRweOe&Cuxn|&-MR|w6ZB9`luDpE6uL${wlKye1imH4pcYtEq8L3hU5 zuo-00I{Aj(WCo4x)hsr#oSl(pWbe1Du{+sgYhC{DRTP_=OeiFvqv416_%ZfzbbENGoUUVH4q zKlaQZ9+LJRea?Xyv`RjV-R##4GN(8H96dFIVxpKTuPmEIXMOG*Z(BW!@{59fnK);W zkA?zL6PQKG-;-tnC1z32xx2p$LOeHPh?Z|LA8%%Yvqsn5-fX3_KmtS-iK z7UhV3y?b@jEW)iheXreZ7U2b*o-q5(B9io&V!Qf|S=5%9rxzGGi{j-Kj}+{mrTjY} z-5^`Ec*1^G~E2uEb7>P>B+vFS>(Q=a{QztsiSlj6+Nx*TwOVf zjE}MyYSqpn`?~=@{2x;qC{8|^-#m-t#|(`>zMn-N)2Vb^-IU`6;~tuSn{{pXl}Fn?w6_ie=7=&7q#7 zbRq3>bI7#6u!>V{4$0psXWfXKL;0)B{743r@?R>RT;Di{zOzdcdhF-WtVVI|I@dWQ zI%4EN@R>vVyP3xQgGqC!SDU#cCwvZ_=iGh0HD(TRlYMGA4$Yx{o9sF z#?c6tR$d9u#qX)H7Piapd5TF-()J9k;O)vzKduEtv#)$2$|DkoB= zK!P1@byl)cKgog8afZi(*K#59!`zuLGB+xB50~F9%7Zvd8X1ai^P)c7pG<-)KXQ2L z=n%;)h}iZWA4#nfLLn#D8wcx)poAMo#G7PA(ditHyX@=5P>QdU$M#j?lwmT-W?_^- zF_vR2QjC(wU2oQ7Xcb8cMPI%kmdhiJer+8Qk5!OC;k%tawVKJIc1_k_U$)Cp{!zj1 z@t-prP7Q7ce$!S7X|r^UT(eO|q2&(OLuZuH6|(_X=Ta5) z{6mWW6+czPdTlT$UtA59Ki-sauvZP?nm)>$Dpg0~VL!taNF)t3AYjUJD7Y~PdB%uYaM zY5{a#6)`k9bKw!EBSKQ2u8qeU648^Hx;yfpiD+v7P0s^I^^o}OHn%n{ePr6Q)9dzk zeN=e0rG|8>*Z^$|J2X7K!w?DR-{wyv7$KgzCWnt)#%Mfzzf^Y^w_eGf2EiON-*WAU~;lV`fqG^G-X+#V_a1!vwy5n zRlAYK)p8r8ZB#Y#BGMN9xtd^e!q^UV=vt6|(Mi~&3YlYC$?F_Y-5HaQ7uIY-VGR+6 zy#kJCdP0<2is*#)+(^!M+2f45C`+Kc(FM7D=H0YSa|?3qGMkvV=ZXsO7ccwxZAE+~ zZ*W%YJ+7WYs;^9Xh2nQt%O>P7I}`m@4FjHH+9r(0ILK~(SKKJm&V0|C94D$U~3 zqwIPwlxSgRTyyvtj+ng)_k2ISA~7#dKw-D z#weXz5^C(9VY!AcxsUU+IAH6)JTMe0(!<^w@V-ye7r|tQEMgc_>9FMO9nbG)E$t?m z-|Eal@SxW$`s4o!9z?@q|8K#Aj^$JFn9zGF9{a;V#bb$mR6M4Bkc!9L_^5cyQkROy z`Z=k1%ut+)$3DhT@mT30DjpM2rQ$Ks48UXW%&2(mrXv-PiCIzcScWMTkMY`5@tDsA zDjvhhQ}GzLAQg{w0v<~PJjUNo#bX12#~hAO@z@u@V~+uk4bbq|>C;p^#xzRBW0wJs zT?IU52YBot;4#4nDjtgkJQfRhYzpw0`Bo|(>jpgb3h)@06BUmMv{Ug|5)F^p0UpZ- zJT?P(%nk5Z2jH=vfXA>aR6OPac&rETn8YzE9<%JH;xR|SV|-&&JT?k=tdfSuq?4(5 zY%MDlk8!q9@t8i~u^zxGr;;}-&W4(aKP5~a913dN#@R$+cu_+oJ za|S%73wW#q@K^!hF=@bKxq!!B10FL2Jf;PB>;T}glYqx20guruQSsP0z+1x&RkBn(W)FCbQ>1#(Lcn9=fXDm*kFBYq;;|6GV-kSJwgMhw zqVrt9V?s1M)&Y1-8t~Z4$5cFa74X;|z+?V^$I1bZsQ@0U2Rvp6c+3XySTx`~xd9$i06caS@R$wYF?GOW&49<~!l`)d0N^o2z+-y=k4XU@BLW`t z13cCWcx+`W;4vKFF?PUX)_})M0FSi-9%BJKwiobNI^eObfX6VvV|IYYz5*U20v;m) z9!mi{7DdBj&jF8x0UmPzJk|+#EC}$JBH*!Xz+ zfXAK!9VU_L z0FP~@;juBmV}}8cH31$w3wUf2@K_Pxu@ivD&Hx@e4S0+)Nf+=~5#X^zz+)MJ$Akfo z@d6&pqTw+{V=5jC06c~VJmv>@Yz**N8Q?KZz+*Om$Ls))MFJih1w7UWcx)2z*h#=+ zu7Jm$)9{!!;IRV0V+nxABmj@e10J&mJZ1@a%m(mS3*a#uz+)YN$20+t^#dMb;-=y; z2f$^R`DNx);FfX93QkBI>u%LP0Z2zYEW z;IV2N9t#0Hb_noTI^Z!p;IT-+V@y(1Jhm3_7(*Bpj|BoAO9wm_0C)@!cq|3*STNwR zKEPu$G(08&c&r2P*jK<~F@VQ#fX5h1sd(%W;IUl5V*-H3iU5z<10MSccq|+6m?_|~ zdcb4qG(469c#H^m>>S{+Y`|kNfXCbck4*v|+XQ$_1n^iM;4u-vW3K^^RRSK{3V7@i z;IS+k9;*R7HVJra6X3BPz+*jt$I=0h4FVop33yBb@YqejV-En2X#gI(1b9pe@K`+! zkNpBXW)67l0N}Afz+>M4kF^6H%K$vaR7J&O%kfk^rU-bf4e;20z+(k~$4&trI|6v@ zCk>DF0v^i-Ja!21SOwrQXTW3ifXDU#9_s)+rUZD51b8e5@K`nAF-gE<2LO-Z0FPDD z@R&K^F)hGjKLC$q0v>w|cq|R@*eAea=77iE10K5%cH&|*03Le-c+4E|m=NHxSiobWfXB)KkEH<~bEe_3 z<)^853M4k7)xQ`vQ2Z8t|An;IRzAV=jQlv;dD40vU`1U#k!c&r8Rm>ms|YR2k_Vvz+(=8$EpF3@c|yY2zZPe@K_fOj~xI!)(d!y7w}jk z;IS~kV^;x>eE>Wb1$ZnT@E9N9F(JTX2LO)+0UrAZcuXJgSTYTdbpswV06b<0c&rKV zm@DA1O@PN(0FPY+JjM)oY!l$ICctBd0gs6S9y<+qECcXZ91V|^03Le@cW73+X{G$3GkRN z;ITD;#~J{S#Q+|=2Y5^o@R%{+F(JTXk$}e{0goxs@E8x^F;~E2v4F?g0FPM%9xDJm zHVJsF1n^h^;4wD9V?}_+OaYIP0gq_`9@7Inc8rF{q5+SM0vqTw+Fc+3v)*h#=+(}2f%0FUhlJQfLf zj05l(C*Uzlz+>5f$07lbxdI+D1UzN{cnnX&W5)oG{RBMb0C>y?@YqGbV|c)0p8=2E z06fM5c&rxiSOVa&T)<;3fX9phk3|6<)289Ee86LKfX99S9^(W&Rt9*i8St1o;IRn6 zV~l{uYygk>03K@tJa!WB*ebwdy8w@Eq2aN1z+;kt$8G=~TMc+j6Y$u1z++s1#~uM5 zy9{`&1n`(P;IUZ1W2%71?f@RE0z6hl!($B6R6I5fc8Z0gqYG@YqGbV6c&r`pm@eS4B*0@~fX5j0 zsCXis$L;|h3k5tT0eDOi z@E9lHF-^c@k$}etfX7Y&9%G{6u{pqFj{%Q00v<~NJjMukEDrFPCg8Cez+)_c$G!p{ za|Ap_|AC6f%mI%j03K5TJjP4IW2Jz{5&(}q0X!xIcq{_&m>A%(A;4q9fXB`P9+Lq) z772Jv7x0)c;4uQ=v3|f~L>eCR06dlfcq|?8*bLyYeSpU}0go|uQ}LJp;IT7+$6^7G zNdO)r10Gumcq|O?STx`<0U92A0eEaH;IWl}#}WaL?EpMh4S38Q@Yq$rV?2PzlmU;u z1w3X0cx)2zSOVa&R={H<8XkKEc@nc6i-5@ED7*f3E(j^z+?q(d5x`?j^LR|1g2xg8kNquM4xBiBjRqAU!4)l87)3<*{44p?<8Z8%+x{@7 z(yiJ!D|-Tq{yAY6m466he6p|frt%K#h4Gqo9uoRXV$I};-@Q(22x8h$vYYPz5Lj_ZbW9q;r+oxYybyi@66FC_uIA%#ZRhO=%ldG(Dks=_C z)tkq1S@UVEF8-Xztfi!}dhwZDIiJSL^^5pBJ4za>5}o6RYk*evC_XL_qh+T>&554 z-*&{4LGOb_fNa@1dcFS5%uk{RyNUG2<4crUpAe4U1(pXh?w3T6w7l=}??l~z^(pW{8E9z~gq_IjVIn2>MpT=sUy7%>UN*XJi zqZ`){C5_cFqoS7WYW9?PU7W@1pBoZmZblvXu|th0TXEjuJSC0QsS3)TCQ2HsPt$Td zhW{sx)xnc&VYC5+3}qDgM9EMV_^(NFz&%7^UOowDvrnC*$S3otM{DQhlXvM6_IdgA zajfAV@(HzHs{bE^5hTf9gz*D`fB0~6E(+VhdF({#HRP_X`N~B1He#9NmAx}siv-V$ zY^ZgAghFl$$6oDxf_{E~HINYT963L@G-{Xn5_uH0aqEu0M&5VsI%xW|AdxE}n^=3^ zBm4vJ)U+LK=vIo;itJS#sE46`|FJDy=$IzU&)s*rNhqqUnlD}K6T0ptTIHAi8BO(D zNbO(Oi@vVbuMauZhc1aU-1{i|4O!;)9@+l5A0-`Ilm9C0J8FEF_|L*hVCsM2WGN%3E7h+0%<(S(?*ec4JXc^q=INp*j>^fXzN=;+em~kGf-Wv2I=CNt#gfTL+`nAB zoGx=Yk*Sq5GPwGX# zU?F;JZ}h6xV`{vX)q|O*5>rVjZ!tB)!Wtm5rEkZ}ZWzXm(DC;q%)U`8}+i`>@LV?mr9 zu-H6!$%gnyNNs)CE_TBh%X$$L_f;yA=PSZmZaz9Yd@>cYd^-`i;e{V|`T@7* z%V!Ep9xR&G>_tl&nh&!#pT?8Y&_7XK>1LPe<<_2<>FF-qWnWv2eRt=WIj17lF_*L| zPHT?fKDp1qEU0ctkI39y`>3w(HQ#+m>gOg@PxYr_<=A`LSpY1gq(}W<$u`x^75VU# z^r&NE#2=cJ^r)P9FSScV((dt`(HTD=?1?>e<)|6Furxhtd5s99F4f@DQl+&%Y!x;#RA^&B>#Dt+sF=6lQ~ zan2|3+>NEp{leE$Bks|=_ui*9N6UEopNM<^TpXER_Cj7OU9g_vQopyMW?0Pu`+MHG zQkYJP6}x5WI6*f2>_zn>RZGTpSKo*%tsPh#^Q-W$#qrOqq8IY&1_|As*faKHWv^Cy zVc$1Z`n()9!PrCBB)VElV7Vta#%y($E$QH{YpH{@=4e?(7ZTn6^M$LT#EBCdOu3g7 zx%|j_!hw4sdw%2MJ{e%k5og>VZw$M;mD!y~boG)pjs?nXBBz=NkZV1cmk$|b!i+s` z0!fF-HfEU%nSH%aHOtUWhDN8u=BMsIuWzoQ!OHopyUK4A-pyy-4YfbnGM{z#jLw+K zeAeBvzO5JjZ`R!*S^IQa8?@_7$+i2}_5J5d?b!Wt7uVFE7}=6vSwCnL+fk#&=#%^r z3&QW_P}^CF*+^BB)9Evo^j+kXkO!?f+Dq+!5+(n$oy#eQ1hjv7VwRlxeanbmm?MMF z=r&}EecaX@!*Edq3uHdGj{f(TCDZ1UdcN4wcK$Ox|2bbIwCZ>B?l_P2-Ll|quqePb zNjo()#T>^vnI?=&&iY{AmVYiW6FoOS zUvf{NsQLL4;HR-;e!hI))ob~G=L<{F`3PDYw0jENM;BB0{z?kS*4~alM3@nhIvhmg+s?yAHeH{JK&|et5}P2|mqUNy%l-N8ahZkWaV#j}2tn zbr4|@ClmGg%DwIK>4oR-a*yA*tI@6n?coKfq5s!2Wk>ZNN%FuiOyPtUU3J@UYmU#5G+X(55OLT$ZssZximR+`ku;lwC3PGP7&#; z%6k9&TLYb&eTOScIOfe`tl(YbiFuU#+5B3ch-n<_wbD-&!Q_4h9W>6KTXMxmPBpFj zAHhp$$!q@=ywrG|FT3KoR^5)t!q}JHpLF^CGz87_26dOso*fQg9@CBCs(0ydqr*9w z#qX3dVZzx}yUC4?v*8YgoND1q=EEiI+5bf8tT@i`W0B_T&+@q4>rYIu5H)bm_77hF zqk(X9Z!$lre>TRg%uyw^Gu*eqO>eleOl6-lu6Jc$X>X?|j(u`8r;=qGP9ylXhxE@- zT>cKrqwEKxacx_!h=+|H#BF^Le1Y^i2}f^N=6kbo%N;fm!%wZmE344W9_e3) z*GZkacBW-L-hAM`)S5>kc+V&W$A%AM7O94 zhvVs;nQr?(jK=TUs68Gbd=Sr6Y(prwJc>tY%UQmyOT(*`GK$DQI*k{v`EJMjD;K|^ zx~%k>Tp|9Z-N)To*`;{8fF7p&BRBBT>B*NB%&PIm8#vtpKRv*c?iZca-8tTf*M7bA z(+lD&{L%VPp-RQ?@#oreZ(>j?S_^}hX!MlV}{AlKjwO=Hr z@J0a>;wLvS5Z=C()ct9~L@3b9D~?rJM;NX1DtJTYB8V##uVbnZAgmoQer?PuPI%y> z^M$AJg}nM80n8!9`(qe<|$zPgnlO(5tra*wg98W8;JW4-Z2Gs1XEU1m$V z4dMI7kE>IvHW5Y+Hi?MsawWtkE-RDI@FHZ%Z9YKI3m`0)`w`N2YA2z!4EM7BNI1cH z?Y#>-F7G86_wE%BF^?l~`mN($?wLRkY&~9L^(~1IW&M`h{30ohptSy?n(4~Z1U;qj zp03_(0y>_fk@7a5khn_1ge$n1pxN)TeaPq`7K7BZwbuVH5xn@+6V_evR=Cx(oI;A z;2yE!&lf^N&c}x%q;mrVk#}={=p9E0@^2)|MXG)iq^wLmEYzk5R%?0=6KCkKM-7t- zLMqF#n)fSI3;kAM$saZz-m+~Sw!F2HYl{ab#+;sA)rfd8rM8P)WqQKc%1x)*A6tuK zHD4p6zrU2mI;#AL-|ZE!rZ*Qxt}v-$O#{E@Zc1okN5b^=N#DhBSm1Uub9^Qd8yvFm zw>LG!I*iLJi8ZDehsoX5AH}UO#-a9SI6FJcX<6@FJI)b{ZZh9`e{?f8G0sx;#M%Q3 zT{p;*U*U}v)}DRi@xUKTT~^pn$PUKp*#wJ?HFjZMnI7*Qio&t8?$4I9QBpe|7(a0A z@wt82xz@uZ@sro%u%g7Seq;Van0n40MFz1W*q1K1oew*YW2e{Xve?~C#mXwH>uMQK zVuxyr+1_=W#-fNXT1-5$u+6HLnRjgSFrB{H`bfV5?02gjDQmhAy*rQ;(6YdswFnwVXvwA@GJ#6&2Xz!Y*by$4R=0mnc z4>8pSZf5b32JA@a%V$#MPcgaAFK|sWFR+65$xd9MuP`^8(Tm>LTTGjP^Y;#ZkMS+D zwJ2WMj!k`@+84&!h3SpS7_+(eU=r3V)yJkkV-5AL+9^YQm{sm;{1h6%EK3xZ@s0nW zB(Ge}qRf^)f~h!it$ILtcsbrS`09LQ4AWe9$T%r%0=upJ@Y?mVX-vu-XLmt{4v`+e zn6!~&K;PX%<*U3H(TNK2Pe)r;ptt>;x<`_jQN?$?wO&VA(0%`|&Nr{uqVFgFJT$w- zhV~wltkJv1fu2s+)vj&jLIRsE$kS`_AX0w(m+cu(_)zG}H_|dE1ramdBi#Q2@OZo;#?AGH7p;SYZ zo4R&v;}0X`o$9pqw6h61R5PrW_resZ*4DlmwlhaQkxlBQJsVLnn~SVI$r5c>DDQEW zvPQ<`)^*82Ht70Vx1Ua?BwO^=);7dj)(<JpQSH)gCD{k8011*`pt!I9Km(2ecdw zso6Z+gr@IG>TJE^h~kd5yv!?dLc}|cj97* z3d1xc6SZ7X)`!^hliY5ov#fsQ{b4t>-yqy*#lx*AY%31;la%a^aH)G#LgYNqgHi&) zGtmPbO-VNL>-In!BhTxKXn3Ojsg4ZwgPzEvn9l9;Q%|J5|B~{@HC`w}Ca!gxnHO@B z6tw2t?}fDJ{#?qr?1jSS!mJryc_Hyyuc2SVUI?dfqg$KF8|8Vvw5wk4jqbnS5)&ru zjdX@zC6Zhyzi3ZMAMmC8L6o}6kDQg>Xj|%xnk1#*#kU`~tbgf+#tpAKWnJ<@-=a4x zFNyI&XRfh)E;sQ)b{W^d{a#J^-Wy+D_rw#ey25uYFwPTQDJ1oNR`o=edxt6>cX*)7 zhQm4JLmr5P@6@>_?SZ_iEN#~vb4Q(tJ*47;_qL+P4=)j2`rVKv^Zr{K*0~`QD*^f* zWmn{F^is^#ato>v(x?jb+l3Ep05kNS;!w@iuIq4=C4(GB{xXooc^v+1O- z4e~!~5EsB}jr{x`xci4%BGa!PzF%54qCG7Hmjn}Y)M@?3yX~PVn(lZLQ)6j@?vy&d z_4{OmQo3H|FF$OEuKEeZJ=QirnJYRjKOfgaa*5xDhUl2yGuQdhFj;KKg0hvIB$8Zyaz_IfC8)^2kJelc zRBFu3pkBm={7+aL<=$9}HpwxvY~I0wgnQqIYKAkT$7^~;HrA~`X*!-lVm6EjcXjzK zJst-1;8oFY>>9(?A929$y!(ZX2WK4XjUB-Z z)*g|bd;5cgX|_&;TM7?gPxZdiX}s#gGMq%rj6Z+IwzYdlaBS?s++;j!dd51j7a}|> z=6c#NL*BS4{Kog#zL|4}CB@!iu3It%3bwq$6ym!*!yY-2znDy3#(u5j{GI=<5VMY3V`*ekfZ>n+yl0`AhY9jD?k-cv!tQJTToqP! z8VgS!?Rndoj#bp|Q9O{HiVZOb88eg~$9PNzv(}C#V)w3nDY*SJ9@|>QESHlUhaF^7 z8fYQ;?ZbxhI<|VIM`2d>TxP>@;aH?p@TURxU07E`j=HrE0Yy61#eMWK&OqDaPim zFR)V35X-E$bHnW>VzP3Qx|=`iV!r%wGo;WT8W`V#!qMS=%0qtitBb zTqr(2!-)N$UuI_8OozSVXgWFfWs*R@Y*ylA%x^;8o6?d2@e#s`OT_`Q)&qn!yB$gL z!FRq8#A6iQqp)s5V&TO8XSQvGit^dG=jYxMI*+|Pz|HZJFns)_vdKszVgH`rdL_D# z2qT<92`$%a37>T4&MB(gA-r4Jze@K|1>u=Xt&k+^H9`{JYUmGlG2!VrQ|tgapCGq2 zD022-HX)n$W3FMvDZ<$Ux4gA{(g?O^Nr|`AVv`74G_Gm4tVtkTyoq&8tcoKXkchvq z#(FP7AW~%CE@C(#VD#7CvqC!wLsMEF%ewsuKJ8Pw(H35W>u1}Kb27LRa`n`G6aqF8 zPUBmRlQ!89zG`gKk7_j~_-m80ZuaOCu0$049w^2WhUjK-rnfZ+8oD~HM!S>=Q>0Ro zRcgBo;gu-yVsDig!Ox)8jYmy@V76Bh$!rco$ZJxecNzm|7n$q|{ zPq?J!;VpD^62B&<)!|p?FMPDFBiHrxLHsqvjmHJ2zu;ZC6i!6icHw(_x^McayvLj7 zdgOI{dx^g?8g!yAv=N_nH_*v)ALaeuNq#KqU1=5mxS$W)Gp-wW8RJ{R#iJ#7%?*@v z@z*cnmpjvGzRt?Ur(|$$%f(LP>&1EK#V@Dg-<>$v!drO+@4z}4Jbe8CzOdUl>a9dH zzV>17pX8D-JR$8u9OvF(d~V{tVWuP>e5HBb<;pBKJe%EIorvuwyu9#1`WF|hNO-AT zy)GSe#`qcIu9U=NgimGXz9^8aiGLiGW+ua{jDMhG_ka>15?^ES3-7&21pn|xrMjff zdc5wJAAbb)t;1KP+^RJ7S&1Lo)N=IeAv%1BUQ<_ zg}u14G2B{a=TMxL+ln9^oj}}1_ueTVJ5OARVB+Z&cb#yDWyWfBHEnR0FWM0#Ng~EL zm9x&XjWQVSvdYcJd56_;*LqIYIQ4J9?f1AhXz*JM*H9i^Z*`RycUWDeFwdC{moHpp z|2|+PuH&Y}se8xgbV(wwHaQV~>uy`VynlW2cinqM*X>3WKI!g}4tN(``ayU9;q6x< zWM1mZ>;1Zs_N89;#X~#8r1*QfSfB=xBxX~gyP|?SP@3Vg?x%cp>-eEOT@S&iWWl;q zy61dV3XMFIbrVm0P-Loz*NrXt@MDcxv@VJ4E3oDFF5TV3x5RzV`RfXOoiuEuch|jq zdQvsS+Cldgxz1ncu(|HwsDUZ6($mfVT3BRyTvK;##^B1`n-z4$;@w-0a*FG!)%lP( z_to?0@~o;Z5@B7ZTi3)Gv-d8e?#)4?vLE**bt?OU$`o-UI^B6`_WTOHI=rje4T|I1 zbsQTv=De|ZrNi%LtlzP_L8swqj_%>#)jB)GVrm`v%XQR5-&JKf73plel6WAMk-XF?+}(&W-Lmmm&`4 z6wTS|^c4Qkt`{}a3E(+&rs^0a(T1OrC!s-2hbgJC_1rBPo!JI3T`+Z$~`t8*Jv)CfzOBb?sURtG10}f^a*LLEG|fx?#n~ z-&#H)&ttyJ_G`to?B8K()~PiY>a~0(^OctEcUE&=fqJcE`WJV@usd4*x9ts51WUCx zF?y`M9(-Oa{_QHEy-!bS9n}9AonCuHYoq?#YSAtGwI1G?;rRV}mlnHH)>PntkCtDZ zzt5j|XOb4bLJ{BX3JWc{iv?@i9uu{!d3JUQomJDio&4+_p+`zf@7Ffz%P0A?p2pOf z7`L!$sa*<*u>VM}wMT1%fA^J9&6@QI7Pou*G>>kmn+rM9rrGM!n;ZJ}g=Vd(#|50@ zea*UzXYP`WH#C0^oN=LNFVfs`Z*%wYBWE=Yj*lsR+f7Q=Y-4sdZ2ldmxtU)lz<+PJ z=1t>xzDsugngjQ>#Q=PQ( z2+L!Q0qN@UI#S~ujYk1S9_O2{Y6QP{R*Fl@)rhauICv=_O~Z1G@w~iKyoT_}(6FM= zNR87BW^uFk0yGSEKMZ?^wrC9Ct_ZA|u++HsEy`xoZlXq#pSgG>yQ+rB(FlX1uf;Xa zcX)5>Ime~JY8Ds~m$Xu&>C#bM!@LP~a=Y;=;>Q8?gb%EbNg4Vb>LTXsS$6`P)i1{n z&8~b>r=HHXD|Jq(4t+*tRWKpW}Z|Qf|C7re;``|1I-RpW21oO?&p(d{B!%J5_q1vQh2y_&VW?6;*0M zEj#vI@Vu&)xJmWN?}yoHUW|m*Ts$di;g*Qr&(mAn9t;W-IlwjqhwmfFm z&xg}qEhfX<-0;4E8h@CnO0yg_kEHJ5#Z4?e%)QuDt~s!xhMv_h>TlHFWS zb?b$h1+`WS;%wkS^H@*7CPU2G`xoE6XJp7AU zC6wUkGvPgrDR8~5C-$8|NQ5mmS z;t_u6q_XPjK9x1k8-pBJJp(e)lPZh>(^nHDLrMz(dL&$ zrYg!$My45l#fvDbA3VYAePS(1Idt}gTQmQxQjl@>%H*=|O73FbVIy{JN>6MjW$-gk zl?JQmL~0&1t8QLM8K>H_UR<8A_Emw)3y(iC2oMI*tR*Js8aC@N}6os%b(DGoKB8xIr8Qw#~x9n39ER+Mpmq!@$W zr#N%s@aF!n+Z22BR#fSgI4hbwbd~EmXRMelcHZ2(OkL5lw{>;ocQM5}<~Y^&rfiBi z1-o}VKRTyyfK+-zzV^FUbP+;--QRcceU*Y2A z?8Z-fQWbny#hsE(+TU2+pQIj1|H^-c>h0rLLg7{#U%1 ztC+%;lS3l3(BZ^bPhWN*);;uJIc3m2X<^bn7C~ULJcB;St?UlwW)M zMOkc`g1pXw+vXl!yz=KAV$Wx0uaK8G=~t6g`)k8Q@m_63)$R>-DGXJgo1bk^TRqL3 z^6K`6amFGpXQzuBqSm`a;jGCUX8hwxo$If}ZqOc+-0=Q>@P_J4O`_9d=M9D1_MVpA zWVAuqrs&Hx3)KxOd+FG=7z=N(^Jl(#sch8-iDJDSNe{>5q>u9pCi;God-=XLndLyU zoWdajvlx4ooOt+;Yp327${kWsRBY2tm$TyV@4U$uCpYybKxOttupIk6NfkNFS!y5J{GWRhE{rL%lZRa&3m@Do8xhA;kT&2eyWO&yhJ5hq%_63$7_FveJcXx6T=#5Pxal7{t0c zQvCJszHiQ99^&tRSFG((HWeRjTqF49p^ErqMUxvRC!e^1&)NuOp=IKn_yE!ay6Ob=AQ775$>YP1S4?w_k}Kd&_@g zL&qJ_b@DobZSwh|+pgV_>To(H`p2y#yl+pK=!Vt(_z$q+v`LUzIKy>DqorA{A z2GPKr8s3j>9HIvv;9WGOr$u(1Dk0q%`Xn;wmMPO9)FjfY`;!!tR9h}$&let>Q+ZaT zpQmdJ{mz3TEYIOxt+ zxj$EG@e56#=^IOQXAsIcr_w)E{8jMH&914T1I>bO_?A~L};Eow@6MRhClI0#}D_Hi`{%4S-mf&5D;&{CoA;Ap$aP!HNjDqEN$nS>K z1_T=K?ypK6cqt%2Jm}z=aZBLpenn?D#vFk&??$jn_jmzB3W^P`-X<_uBL3&~VOxRO zsJHAl#xwa*G*jW2{IuytC6E#T9JRe0&wq~$Y@YFlmc6i=B&Evh?wljb4exA3c%Cp7~ zym+#oYnJeHn((lkz5PYpRi5YZ41dPa05+bm=ca^lm*4AoH9GHKkZoPBwyLE)dRguI z1hUAnCr9(w-@W!Fdi$-z>*X@tg_P@p*Y6ZPuts*9?fN_IB7dHd`E>!}07Z$%G#6za4scwo?^7<7JA8Aa8ZiTi|cxd>0++FI_HrO!HjA-}c~n zvg~Gk>bE;LqlIsKHh-Wam%4G2n>I^P|72{lHuV?PEy~nia<=f%Gm-CX5h2dklV|tb znuZeQ+0Nppfg(gG)QmzD=4HC#j*vYz8|fK1={Oe-nJ-K3pi~t^p4cL}AUIrxu7#%1 zQxd0Ae&&6{;zx;0aOIPpFG4SSiILXayuw=ieC4BrX$^L#2)+<+-i2dl2OIeReQXN< z!m(jt^T$$qmFLRGpe;a772_Z`yRwny+j=nEb&ghkUb`%QrUhw=?!V&!{;P{`Mh`yt zcO6AS3w6+T2!PjtI+XvqPM@NqvuKsfQ=-L>10}~SrPKaFWED?NYDbjZ_*VR(d~LdS z{=0>aT)On!m{xw?sVshyMgDUVi)!Ux{ZI_vPNG#YuWJ@Rx}vmUUbq%Z3d_*VU&+1s z&Qr9~e~(%xr8Knz{~3wPZhQ=f#V8}eL+*5AS(q8C|9Cy6%?!%W{nuwBdDT|Ah0#v= z^ga3CZE#9a+wh+s?S&AFJynEGT73flRT~{$l@wV>N{GDAn`Pm^LtR8i8lq?V>%asV za`skkiTU0Qns6)&q5R?e&wtYRE*vZ%Pe*6CFb>qQ8Q#h#TtWHHARutbby~Pg9*z@OyqmGFWLZOwKdz)XAB#yL;llqeL#A9J?+ZUGRSeEB;qY zOY>YWmVWF)9nQZT)8AEy(Jh?qeDZ4AMf`hE!ZwpxJ!Gj}CVB8N;9bdQJ@{Pbuhr6Y z%!>_=xzVQL|GHfaJ-@z()*xj8{?{jSE4iDGgZ$NlXJNiC&A_~Pq_4+(BY)*%{=09~ zt}PA!yjX!$EM@fm4hH>K1@rf<`FsuvMLZ=JX6Vv*w~MtV9a@;7lszDn zCA%>3=Hu4V`ug|9RZ1YQH54QZyRmrAH?%Yc>tdzjhiMDbOz*FAgi~A^S9P&uW+H9D zP)=t-y8W~-K0}-0f8Ux3M;4~@;{7Mti;vVo`O@@FsVL|7V1b+NT_%8|51GQS<>U+21XKi6XEz!KVA{ntHjVFEZW z-K4cxnzeMHG;G=WSIhGouV|zH_YJ$Olq}#azA&AZZX#N2AmTc0I#F(T3oByjW}wAm z17tGbuE@|d04FA-VJJRz_S0J?vN;xbm=~Z#nRq4|1;|U?+>7*dlMF`V7)~> z0f4AWPhg?{OE(-WmM*83rfol1T6%tosdh%5cczq>|M`3JJ6A4Des>Fd`+pTQ?E7V( zLiTXwV)#%=xk2%fbv&t$mCOF_9JKNM`}vb|hfMEAeV{Gfp0Ie%fmI9ROL^6ytj~p~ z#8u7byBD4kEY;*^B0S`?ZY-wrr@nO4!D1cF_x>8Tzh6K8dix-kxbrc@)Q}t9CG_U& zT)LfLu})$gwNA?Bf`45Ot|&6Y!bx`5k?9FMplNU{8y~*e|2=JW?!u1!n1`sXg7tw7l-U$t7D$1)5iMmrEL40?Cvc`eG5A7 z&BwssOfL48T$u4ot@>iiT`&JLO7pv+{&m_X*iDUT@1!dhOHzJ2^lJWf^Iv5ar8JQ~P9=HS0eFU_>eBrt{Q|{QKl+J!0pKtL00Y#KrKL7v#