From d35377a62c029be9ac167a2957ca0c9277cb518f Mon Sep 17 00:00:00 2001 From: danieldouglas92 Date: Wed, 5 Jun 2024 10:31:01 -0600 Subject: [PATCH] Add changelog, fix tests, address comments --- ...e-tian-parameterization-kinematic-slab.prm | 19 +++++++++++------ .../doc/schematic-diagram-overview.png | Bin 44263 -> 49692 bytes .../tian_parameterization_kinematic_slab.md | 20 ++++++++++++------ ...e-tian-parameterization-kinematic-slab.prm | 5 +++-- doc/modules/changes/20240605_douglas | 6 ++++++ doc/sphinx/references.bib | 11 ++++++++++ .../screen-output | 1 - .../screen-output | 1 - 8 files changed, 45 insertions(+), 18 deletions(-) create mode 100644 doc/modules/changes/20240605_douglas diff --git a/cookbooks/tian_parameterization_kinematic_slab/coupled-two-phase-tian-parameterization-kinematic-slab.prm b/cookbooks/tian_parameterization_kinematic_slab/coupled-two-phase-tian-parameterization-kinematic-slab.prm index b9b35311a99..a116d37354d 100644 --- a/cookbooks/tian_parameterization_kinematic_slab/coupled-two-phase-tian-parameterization-kinematic-slab.prm +++ b/cookbooks/tian_parameterization_kinematic_slab/coupled-two-phase-tian-parameterization-kinematic-slab.prm @@ -4,8 +4,8 @@ # are prescribed an initial hydration state after a period of equilibration, and # subsequently dehydrate as they subduct. This cookbook utlizes the existing melt framework # within ASPECT which governs the mechanical coupling of a solid and fluid phase (see -# https://doi.org/10.1093/gji/ggw329) but introduces solid-fluid reactions -# (i.e partitioning of fluid into/out of the solid) that are governed by an approximation +# https://doi.org/10.1093/gji/ggw329) but introduces solid-volatile reactions +# that are governed by an approximation # published by Tian et al., 2019 (https://doi.org/10.1029/2019GC008488). # In detail, the model simulates the upper 20 km of a subducting slab with a dip # of 45 degrees and a convergence rate of 5 cm/yr. The Cartesian box is rotated @@ -15,8 +15,8 @@ # that is already subducted. For stability, the slab first subducts for 100 kyr before # the hydration is added, and the slab surface temperature diffuses to the mantle over # 10 km perpendicular to the slab surface. From the slab surface to the base of the model, -# the composition is comprised of a 5km thick layer of sediment with 3 wt% water, -# below is a 7km thick layer of MORB with 2 wt% water, and below that is an +# the composition is comprised of a 5 km thick layer of sediment with 3 wt% water, +# below is a 7 km thick layer of MORB with 2 wt% water, and below that is an # 8 km thick layer of gabbro with 1 wt% water. set Adiabatic surface temperature = 1600 @@ -24,8 +24,10 @@ set Nonlinear solver scheme = iterated Advection and Stokes # We choose a low number of iterations and less restrictive nonlinear solver tolerance # for model efficiency. During the initial pulse of dehydration (for 3 time steps), # solvers do not quite converge to 1e-4, but achieve the convergence criteria for the -# rest of the model run time. For a production model, one would choose a tolerance -# of 1e-5, or even stricter if possible. +# rest of the model run time. For a production model, a tolerance of at least 1e-5 +# should be selected. Reaching these tolerances may require on the order of +# 50-100 nonlinear iterations, and the maximum number of nonlinear iterations +# will need to be increased accordingly. set Max nonlinear iterations = 25 set Nonlinear solver tolerance = 1e-4 set Dimension = 2 @@ -38,7 +40,7 @@ set Output directory = output_coupled_two_phase_tian_p # Checkpoint the model so that a restart can be done if needed. subsection Checkpointing - set Steps between checkpoint = 200 + set Steps between checkpoint = 20 end subsection Discretization @@ -229,6 +231,9 @@ subsection Material model subsection Visco Plastic # Use the composite flow law to simulate dislocation and diffusion creep # of dry olivine with the parameters from Hirth & Kohlstaedt 2004 (10.1029/138GM06). + # This will be inconsistent with the physical state of the model (the solid is hydrated), + # but ongoing work will introduce a more self-consistent coupling between the hydration + # state and the rheology. set Viscosity averaging scheme = harmonic set Viscous flow law = composite diff --git a/cookbooks/tian_parameterization_kinematic_slab/doc/schematic-diagram-overview.png b/cookbooks/tian_parameterization_kinematic_slab/doc/schematic-diagram-overview.png index 09a041e75a0c3a992b91ba50116b6379acf45af4..3ef4e7346744a57ae8a6ccf504713cedd612a5ec 100644 GIT binary patch literal 49692 zcmb4rbyStl_ce&3bV`RbDBX<`(jeU+Al==eq#z(6-Q6i5E#2KM-F*S+cdq*RzJLE* z%O&dEXXcqXvG+b_?&o)M;wXp&h%hiPD3TJQiZC#6q%bhBx6k0gzm)G`rh`9t9L3Zf zm28Y1UG(jZU}W_jZ7ppaEzJyGI2+kJnAupfF|sqV(Z4WtbhLHgWn!}W-vCA%dlROC zRMBDZ5CmHZ4F?z)G=1oA*nFWpGZ>iG&yu2VlwFhe7oF6VcW!>4Tp};hlGotn??dUY=)%v){CkYtOW^+egf5JY8s!Q83HrBFe@as5f1?Wr^ON-k{T)Cq zFLC;p=t}iS7Of8Qk&X%67Ir64Vkv$`0e*yNG$MOEwO|#m(Gnz&o zW_zKIO~4(JF7@rDx|Wv1!@XBH5l`Fo`QB`!TUFO+Nb)lzq;bftP18A!&%+%T43%{3 zhnaGHcIyQ^7Z(=;1A`x)SF8=r2dga)_iMY8MTuQ6$8ZS<%yH!JmR;9S+#%bM8wYA- zx`RQa-p^hquo`WTWQbS?Q2Y0EgN-K-61-ZEUhgpMi=lh`^l9E~u1b+=kX#7$)8O?)oiO zM9jkCB0h_OsIEWtz(A!@e^k<$=h+~a8D_**wfW>rv`kELhm(^N$nIq5Q#2BoHFVF@}YOp>&4g=IJ)MH(qTf zS8kP2;{4GV2;#OZalw!*5Zsd;*NC2&}%#u@~Nhi^XIA3c~ZHS&;2+NY*s!Reliij!hsU_>6)MMXu;d^f{7YIzWpO=QnduP{*8&{)~s zl}h1ta(%eF9JliRuuabgtPZAByUyUdlBJXe@yZUcSI_fFHFmQxf3VBIXhmi#jm*Zf z5%j#ypPLLN_uBdxg7Hze9&&w}eR6I12uMiMON_4u2Zxe*+AYeO-Uc8(FR!l7`Oq5$ z*`JNpX>c~8E7hnn2JR^zNiJsV=$QRBjKK4x6PMrnmV4Z*iq~D))YSCUqA%Y^-dWHu z(2(HT1LkeNLf@GY$C6=p_{p3_*~-R7KvR=|oSaewMX$ zy&B5I+N26Rr4O8+p6T*d3;EscE%xh^voq<->%n9mB|BT&p5u71g5hHUX5;mq$le&a>#uAXOx z?vkmCgs6MjgQzAvrlYk7$>cmZY?`wd<#|z24%WiL!U%X4zO~QYwNs8nBzbQc12?z2MuT(7pO0NbcA%~r zZFo3Y6gO&!*Jt$YrAmxJX9yh#>1y7;SDqT(TpaX)xcF_^`-Xuhqx}1KWsr%jA>4aY zB}(-4^ik2dq&~L_AP_1+_mN(&nSb(9z1Dis-m}D?R_0Fv^uF4NQ74UfJYoM>6 zfz6<)4B~o2)Yq?H&6E87A9r?li^F1badA-z3JMO|G3wNRV25S>0nFB(VaOL2?u&Kw zB%zvEFZb0CyR|?A8=JJ`9|CDV_(|2A#Y9DEsia~gnrmWXW2FEhh_VSGr=X~EBvPUw z^<7|Mi{pv2vCj0($b7eKM0WDHV2bvY^uQau+L}H$3bIksa{K}agx=14b}Ab=JP2n zvnlbnZ(o9(iY`rUOJu#G)u_xqh3qJ0b#r@N zJDxiG005P1XCmHr)B&b5uD=^h;)=E)nEP%Pedr90)umom%gtVi=XCuWU=zj!!dF_~ zF4jAZ7lUM({K))ZBZbd}^-BZTK*ap&AHu=Kk*Hh3>tA)tOitH(LGDacEz!)F&bwD% ztaAu*JKG$sb1)&`wCo8$!W|zwK9@)=SveT{QSv#u$UcZ!7`l&PNpVE)Tk3(#B!m^s}LN9#6uQjM%~7EeHj?} zRQ?qZ6QQVbtvkZu_V>y9iZRNbedvln@w<#iU>D6xle55*llk4E9`25OVqU7eoU8}l zad9x;4FI+rSl87~UYb2s;MR^K?;|+4aUD`3qU41-FA!-8rtt6es*DHME6oOewff1T zn{~;Y;}6iabJ#gK<#=2k?pWhuVn!L08Uk$m5dg$dC$Rp`(NT<^xvf`4!XhFdk08S6 z>gvj*^7D*)aoQ{;-|Db_==thcdj^5Xkn}Qu$g+4us-`xLyoJFM!(~`Df2^~;wH3GX zt}XDWtr)fVh5y$K=K$t}8msxS8Y?XTVmiW!`DR*tgr1@jy$31MKPU*U9OSJ;4vVgj za8JRxe>iOpLM48bu3cvz!=!f~U`VT4jAw0ay;Cah77c)7r1_7jyPHe6-Fz{P6gO^0 z#_%9C($&MmRAU_w6~`8tHr%S}b=DLt=U^k|T@7dlL4!@VfFI50$SW#Ii z<-!NN6#4qZa)j9i<^#X`{W@e@|^?Pu*kDqI3g1;*0KN~8QeEI^V8hjFJI zt{xsHKT2wAxr)@vsCju40g_wY*!X5l$i!5}j`m5B)f1P%V-f6kwZ*jYRh1`51@u-l ze>c7kJb$FZP#k2Y_S232I)^P(rCgZ|5CdQUiwqqWyrG{I4AeJKqNZdbuaewQyl3&g z|Gcu;;G(yk9GDSUg(cm`U0I(}L6Zz-&vifFy#;6{ftmz{g@q--`ivY3f%qZja?9<- zV>^x$g({)IUQdRF9&91o3IJ|L7fx-=+q4YYN~kLJ^iev1u;%NYf=nXho$NWt$``2r zu<+sEKB?+De8!|wJTR0|ZNXwC5acT&Cnw(?%@RM|vIJ?gz?_LK>XhTd#T-i6&DEV^ z-|n8Al$2C8KS&VL?R-K^O-FCYL^T!9$V3^dDk{jhx#IyQ7zI`bgTtgd00QV3e?tEt z&1M42f+v(afD~Zj*0PoKu>K7|xREkFq0Fo-J4eTLK&tWp=us|IM!?0z1&9_F%2F6K zzr*x6dSfe%oC*{x<{*M&=$f96Bj$4kIqPQ>m9%{|zJjW5S;bzeeM0Q?zS|X5EN)Y& zCL#jD3c&X;$VA-CD&x#beS^~}tO$D`og_XL6An&}UKa`GAIp^`?2BdS2S?}%*~UO4 z;g_&3{a{*Sf3_(*xa#2Gum)%pl!_8?n498=t1_sUzqmr!0AY7LW%LH@Q5}i&Xr|cc z-Z>=|mDx4Pr)yCckgE1)D&St!nJM+1!%Ip^?r_R}4hiWSy8;?Q^bTjBW_8TqJu=5f z4W75gd*dDN46LH8iNKs;aD_sFuQ}|?w_CXTePH*DDm&q0B32I?CnE_)G;*K!)KfO zq$wE0PHN0B>C|`lJw+SS@d03IJy(y;W*xXn^!)D8d@WPE!-txs<3v^q#^_Otz8#Yh zWe(~<3n`HY{_7wiVj zL#i4D)rP3B@FSRtii#b#U88jI7US)mSE&r>s1H*`*JUt-k25XG8FGl|)5F>k-~+YZHay zFi2hE@!TO3^kJ+Q)8m)+g7m@m#RX)t9o^|h$nlDZ`^CN_7zqOirA>X-f>FovGAR-g zk}tA&cfI@hOF%v4m%7KsASm3(%K-`s3b*IedelPA=m9BzM2Qa{C;=xZR<9s$rzQoK z^M#0q)$PMK@RoV*igO^ru4j;ThY@x!)H#d`OaWQSA!p4IV5Mjqf_o4jWP-Y(#JN9? ziOY>h!^qW~+_%V==V9ZOz`;yGx-9YfHXS^ToT(Tqw0h{=07N#Htuz2-o(!Eli z+z9~rk9M8By`QauKt$z(s_SvUku{hJB7G;&rC9k}($e%!kDyF3@jN8+_+sA1#Pz`b z8ee8sVR6QGPUAVqpV}78<`c9)J&QF?g3i9^vJCsK-qEbyap$ycs~Vj>k338O*u>3A+1 zk}yIpYVet_U(s>s)nB)G-{osozf@IKU1;)%10n!Oj)=~O3%&UYSz`+eo8N#|#NrMy zH83k5PActa%U$qn|aV{nk+ zD))!|MxRDVp6F-HTfo?m3D|>yL^xKiFA|7K3`ONoH6$`J66E1zqyAXC&4C2V<@C32 z-=6N3wVZ;Kfq{Goh!k(X{A_Rc1KEld=ti;h8p1$FsbqKg7#D{N3=Z&$pa{Y1XFfhYuT%U0 z?e>-+F?I*(rPmJzZ=zHu34`M8F$i~1-nBbjHVoLm=wKoT0vehLp#3uO%pZXwX%de- zS)zqYN=oVsOt|AS*5}Yr3T9?>a6rICx__Un0ix9p;n-MbGo=+zAHo8ZIPG2Khu(Gj z^_`uBoID74qg${oU+`Ey0U1NF#oJ5p{jsI7l_-jf z>+bhjMDJNv45LnBhDeZcyrWPVNaNtIyjKrE^5g?RcO4`Ge}DgSU|djx011J|>zW-v za3lbc;3tNCz_!a=yCWy3q$IGJpaM!s!O0m50NGfRhYOHQPh0@j_Sdx{QlCC3uw1fQmR?oW^Ptk~zIM9##frFM%k8 ztjPs@09m-z7#a%E8L_w5y#ac)fr&}nic^1|+snayhiFArp|UO@g&AKp z49`SVIFRtWzN1qqR1Um>vg6O6KeG@hsS0L=L{;POa9J$cp!X(gxK2N4yUK5_cnrLus53K3dHF#;@&>sssbJt$3PN>>Zd@Y z77NAc1h5L9&6wF{sZpj(OA3@i-UAZ@s=QLOryF!y-9~PJP)z{P?M;ac0~=!iOno4p z#ScXJmDSZxK!ydZ3o}bRViV|28oSEcC7K*S_vitkC6&mpp@9#)F;tLyn`DR|9X0$x#GRn^|w`Y}$|{WEm)0g-1j9ZAO(Y^I~F zV=fm4+r{IsiM;IdAOP}aZ(14@d6sB##WHAhWr~F=o>%lTz{0|6|FHci^l%#u*pMPP z9`Gs^K!vuPtpq$j(Wd2&Q?C(M7@LxpimO@RZ?Ot}l?MRYe}$|6V+HD8U8@6W*fupL zrUNi$JA3=S(fHs0%3Pmh6|e!q1_H5i;2BYq?Wr9An*fB~WHuoE{p`R1z@xN)Fx2nq zabMZm+dulX=u)O{xOSTvgpv7rzWm(d;^-FE*IguV)P8+W5sMaDJ`>*BZt^HLpJn9{ zml#H(Aq^y?CnPkH6co*v%$U5KnX51N8h_wszS6@^U)x5bg2e#q5zR1oPK@$(N8qa? zH!gXQCi@Lui)ssut?g;cE!QR?w%925OX!+{Z^fg~$Yg%DoB(N@80lvuzLEy}HBwAA zCA0(rt?udlsiGyiM`X(e`M?s>Yz}_756b#Gi@(Gh=%;(~?k&FW@jr2Nmt#sDPIeJ` zLC2!7SN><0PoGE%`ZHGqI{9r+7{PDTbJuzzR4s$2`9Cgc|1-l#G87_xQZql#lCI@R zuU&83-57AaW58Tnt45XnQ$bXM3cJ$g)7z2}BuFC~>JDye1i4KY@CX$Znn6_1c{1?9 zG(BQWEfWE$O|j`!-x$_p zkJH1{=Wte5dMR$_i%3T!NX&?1p>0%gASq5*6C;CcO^s>dwB*KlBK@b?QaUj{WnqD2k1+a+Bq|}m;Wl^u2-SD zShAu%-|bIv*!s7-ZIsX%ZP5>%Y;m zJ;|57T?>*q9FHLEj`-gY@J6D7Mm9%X*8=e&0Eg_1}wIZX@S>vB&rgE(69nY{jramA zF8GmL-uU?5q-5i5AM-oxDbH)6nxLVPe;dmV*3pK%?dV?bavP+H^7c{rNfFbtByO>_ z{6Dck5o*?Kpi;j*c8I2H1`&kOYJLcm^|>puQ~u+F!%%a?2h}Y{%X2K0MDIsGYxI0R zEFP_~hfDnl1kgnqzs(yTZSLe^KfS^9)mig4nT%-3qz&ncgl>9f`JXu^=N=pqpWMRy zZ1#CDxZ)C$ctj3efh;&Pe|#j<3k>fP?8(5_I`pj9p>9G;*|--lazFWQ zWQV8e1FXZ6TUvI9=wB}pC4-)0fP^2f-}ZtLH8c^xB#7 zV&rXl9ye9`;rcj|va9o7U_gY0&f;(0yx?~zfSXxWzS);s?T_k};yXeSumond)E>@!d6}xQ$Wwpw4_e5x4)3ZPWiA z0zGQsQv1$W$Wa49Q#G8+#euZ@CTzl_7AqdjOh@6aNBsa<^kF4Ez$ZLl!G5Z{ZOtLqN9n;~A;KBK*^R7Sc-)f(H zg}#tX-H#YU*JEd*5Y^3{MJ<`#$v6*OrvKK02VINC>b`d3O6$C>W6|6FJ3cTZh|gsB zzai*SBhlac+}DR83XF)=ffewX$@yElnJB^dssJ;=e}!2Q3aLH&u~O z-nl`ZIr7@t!iaa9?4ck?C}Xo*7jyQ$+;-#@5a5%~e%rK?PFs9;i&Z0^@^9})f{o4kKNRo5S#?e5#lx9%0g(L6w2<+~3P?4Jcrc7{Xf09FYXTcp` zII)tU>7wR=WwRAwM@RmAG7sG7=x&YMb2aPn9f;%JF&7RFC(-Q9Q^FP!5l>H|!9dhe zUUq0Gj1P|LJa7!j{NOaNoXYR^DKxY{bpPdPMXqd1wpi$JH7bSLw&O@0+eAUueQL~O zE}OP5_~?X5mnXS=;ra09I)3uBOBEBP7lLx@ne4MT&US`P$L-dctw5^1HGH z2Z?}GgIpe*n1TgAJgUE%u0k}LgB(8$=(!=*y!Xvn3KW=c*<7DlHZxfJIj`qg>}2`} z0cYgX$j7Gz9M~s&E;PEaweP>hCISK(o6!Als`*-%h-|88a;otHeeldeeRwE9yv>4l)+xx0wYT6>Urn^LAfZ~s*^4UdqEci#Ap>HO*| zOGu67;xz)6&prGW4XhsVdjPB^5w~AX^2ZRas>C>QOlDF${pz04#Lqj8Y z2k~+~x)1N@;GV0#>9w>*thG|&bG!MSty?-RQej~J2$P`D`S%czPpaQ`c0aXV{FEiG zRK_*o+0j+_S<nYY^@}P7G6t2HbpYo{qIXtXOEHQ$@JUdXDX=2CXCtVz z|M8p zrOHIE>U?ecQq$a*nr~_Q%R=A2y$XUs^twhi-W)&xJtOjec5OgM)d3DLw7^`sl+LeY ze(z!^3sv^yhf_1Q&o|j5G^43U3RUJdK7Te}jwnzHi%|cDFVceQo5+!&Kw;~kM1%kM ze1?$oXPpDjk)Z)2Dsf_OrP0UNn2Yr&@0S`2)HM_TO_M-L+V%mE^x^b3cM?)k9sNbP zO&-G?ZR8PoPlv?>1xISE{>c;hX@^oQ*R+jiy zm7|#uHdUltihGsmD01D0Y55{mZe_c!fl1B#Z3Fxu;X4JoFU}Ys}BnzyGSzKcf zCA+<@Cn{B-q1<#7lb*rqhpj#Wn~5_8S*w55shQ=2_wOSi4t$9Obai##loVdX*q|Wr zT#eO}+x)7Un)im$C!CwzZD?tQH6|&%c|WzXWR9pL{}ID<HHQkKv1hy=Tw1C1M~$xWdoQsx5F`-CDJ53fSP(3Nr6oEoK$(M2H9wQ}~K}o*{*> z7)}U{Wn_Tm*noUiSpphLr;0WFN=hiWtSzb8zY|wfU}y#zv35tae|pk|1qjJAWdBMd zhph#=n~#u#qA6QZ(R|sqT)BXjhrk53NlD}+PL&`e__G5(#R5ebIQZh$v9ZX-2CkK@ z>G0rZXLJ_8sqb6bgCm4t%7No)8*}V{Gk-+$C+j&D9s-c0z1^ z(A|?b6_6C{Ksvd!>FNsGtEhEZnJTf8NtsK`lP3juu{W9biNjVQ2RfaF%ei~`cmFqj z>7H)mY13atzgw3G9=n}B&Xx5u?EW=GTn%@-(O+&&tld-?KA}<&wX@bE&o7X~Vmjhp zT&&G)`<`$xv6@Uf5D|p8vkY41_rAPs{VvM|T*q>9W}p!7Utu^&Eq4{xB%oXr>vB+q z#0*`S#llQbd3hK&FAk{GEy|=W0HHJpx@9KAg=pH?R924+^kLzi3MSPp!~k2cH*$(! zb~~HU_fE4>b-HA^h89e;IkzAJEZcWQJVk=@dA-C;dJoB!^ORpIdv-!5kaY0;g;+;@dxeZwrd z=E)j{V!gDB-y9Z&-TYp2zpreTRntUc>CX1yTh-Ih#N}pf;dhgVLQ!`}vd776#n#V`iYV$)e zyl7d_{oSN!Ym8J3jEl<;dofwoF~X?JIgXk7=b%2;Fhq3Y3-t4}^mN&o8SJ(7-aB^l ziBC6|gLxp6iDp#0S6|yiknEC9<6>aJd3&1^ZtyxyG?R&ijB_-+QLl$$N6`f$J7Qe})UG`&z~uH^51#~NNMj+U>Zb=M?`RlQiVk|&D5K~_Uk&1@KW zEnr9wyECtLhsyT6JqD-Nx?ovs=X~oO_DrKjNe(>6iUbd0EYC}gO4QcYq>AeEU;DG| zzkjoD4H5sba0Oaop3M8C=^&*u`T@ZD6BM1Z$7Tk>91bYTJn2P%|tUk5lH(Wh{H59}9CrE^TVcTwI9_V0Gq%INido z!$mDS7ba=f-)7Wu5DNa=MCe&mTkKxG_Vk>4ESl$+DlkcA=kStTI7PNeE@dLiTwpn- zZ-huX&(C^sLNg*BkK2XuznO?0l8H5p`QXRImN@{?YLNe+Fn}pQI zb93(M9O=$EsiX_B40C8$6s(LnjFFK^@996YP!bP6oUl)r^XuK{|9%azchEf;>6-5Sv*P6B?;*pJ(@@)h?9k z8g|AZ*GHh_Geqj|A5>J`dU3Iv?-wX{K$oPLBg0?stP-7?O?lo>E%LOQr=a$->TZAq z*3O}5;47Jyts`W;E9}hG`*vx(nj&DGjW>1wHhEan~_<0Vfmr~JTC zMXd(NV9yiK-}!zxl>vC@=boN^nbiGmzo%Q(sQn+T7ZqDb4GhC*Ug*^Sd{wr-PZWgK zp2)H5AGlawfDn*DJz!%=BYgaJ@hDJ`gz zcwAcNXx9~rit8Dkeh@lbls4(9f3$t5bIwOjO??iB@hjLdk}5s2lH7NE1abE!!yU;V_$z`z>R>x5hmiJ(!o&V}jj z@#Bd(MA9@UPGc)j%9}V^{fsr)O#PMfB%bvNXbq|XwWKcR+aK2rDT|?v^`AO>3$V$6 z6*+k&JVmn}%a0`Oj!63}$7zfIPS&ao2Z!7X$ZRxC(XwP0=GPjl$+@6?DUBOK;=BDBUgE}U^W6>#@>LBb z`gbjc3=ChMbYrqYx@dYk3W`EL9LC+Due z$HR@@oFHC0HI2!z@zG3YXb}f70OIWj1eDZ;SOmt7{Cs@X1*)Qt8y}kITma^G`qWnZ zQbp?SQihu2R8Y|Ttg16~lI$a&GyDAP7}XVe-4M9c1Ud*$#-z`K!Xf^&;R`aJP zUj`rd4=l{e&jy{lC+odlpy6()+3VjsOcv!`l+fxDdH;TMH$5zWFteBw${VH98B%T$ z$J=O_!mGhonv;|Gfs!)6+hUOb^pTo8BN5UiHs_)Uq*h()&b#+ft+p@`(&!Zv|KA zZhgJIpG&l+e;q6)`}<264_ShAm&CSQ;T!zLzd&iySfLAk-*imWit(Bpc`zH9c$iI! zpVLNF{4p|`59j0K#fIJN;*iYn7p8}nmdVdhbwzf@20>*zT{fx4ase;*k31%-DWCc? z&8LufLeE2X&(!`KINV_#Ztilkal}MUD_H{r_A2A*t%W)#=$k7jn3wD>HUNlDN?J{- z;e4Z=ttM_pf5Po7W0h(&Gd$AWKB;myJ<8&-#h71QLi`L_`^{`+NpM(Lrq|6yu|{-z zA-E`k-lm)@r;;QCBpA+40+v+`8Yea&t&jMFLhs)6gZSAdv4oU+-$Gpqx0ibkNYOGf zU*?$BF09rNB%0jKs}z(RuS=E0t)k2gqd89`dPfxwE-{Eptd-Q&HD*B}mz7_D4>WIx zfR--x8q05g<*;(2&^b~nR?N5I{zom-yqVfGV8$;cB@7l*^tn!ZjFeQJEI!_9a;e&d zP>+NZTzlD{8ACyVQgmRcPTHUYBoh;BQsHr}m9)@Fd5!ARSz0}kp%lNHOEy8{(MIE*U<_)>e<+Z4j} zoVZ_OR2-5l1F_b@?LX%#S2MGJ4 zKLHLb43tn+dL7qYotJvV@-&PsPDC(p`sOXR+oZu_GSUq^v+l$)C?N#f0%x;|&Nh&}PoO$X(mO!+ho z<)W#y?Ci(z$lD^O+m?T|Q$!SebG3CsV2O|;p%}Nb+(*2;&dI_}?r)$b78`3GCfpff zX2Ba3_3PFxUPv$HX)kth$s_@q$ZCc8^M!`I((#>%?Di>~g~drp^jB)C7UnfJYC#%Y zoR913fDG~BWEEX5`I%VWIE(i!wyp=m-y&1+qgO?LVwawzloUDbHj44UZ(QFRtB20{ zdC$P)s`7UhJdKTCXkGw%n9Rc?|Ju@W0ztn>MLNk9FNPjrFtMQbo$x2#+1*kdb073} zZSrG!1s1N5sE)B58(T=e=2TbsU0fu^zfE?N+`K}65(RoM=Dz#jv3{V0v6wQ{X3*cl^tWC+9hTxBSF-FoL)UW+)*X#?JnxFds8A2GEx2oXnFFeaKKBn74R;m zCc(i0)PIOYv9@&X^X+dtZZ}N16e+^Apdemp(4?`u?XxKbieLrfq(Hi^@_8uxLU2wk zpVr0d4)M?o0w9>Y>|izv^wR4NChm^E(Z`JP4EB@+$6>u74_aByZjP2K+|F(&sk_D+ zuS`os*gj`dsyC`Pjh5-~Sv^wBMMfsbo|5S30$DfESTQiWZ3>5gL-EG@M_udrUb$hA z=FRD?APfpYzHDVU5lh+7<_@UYeYbgFi;d;nU%Des^JEQ~nNe@PJq3REXE{WCx%U<9 zs*8G={%$V*13(!M_LR;WKLL6pbdxTZeOjkV(m*4fcz>*3I$P(yCxH?+R0qN+(BcA1 zDvaUX*@5hIiC6Y*|44V{E>p+gb$myI*XV(Z!>Ak`K_ma=C~Ld`%?10(#Sxdz*5xfI zRF2l4bE#|2m3!Ubg8Pl2?{4k){mt-TCRA7L=~(ODsb%wq(8`OjxQuOytXuG=go6&K zWTxhKCbtW9XVOz8H~zRpG-S&P*y_#NO`|@kx@DX`h=__wGDdN_WrD?ZK6)uvo~xqO z36p1^#icg-3hrLZrK!@0)5`!&nL*PtJS|Og%KO$1IqR*LIp-e%tOpHH?I&&VETm_P zE-kR-9+#KJn&-{43)n9ai?s{F+0Xf#NJ5?iZr0jL=F#nSjSZ#7ahCJBQwAR(=*we@ zR*&G#(=ghTDR&MGOP@cwzbC~)Pmiqzy)xlk%eOx*|0t>%8~`VOJZX7o>Mv_4q|iH9 zejx0RYYE7tG^$isPYa;vYe zPxAK=Er9(#+NqxY1m4)V=j)qfZVV{{10n$rIgv812=XSYq(tBE)+z3F0Rfsf?o&chr!Us1u|z-H3c@5M9lpKF0NStkABp)i>N_d=6V-R`gp7xhLv?L>cgCv{_-aW( zvVlf%pgHVLn#B=PBNOqfr_@puAtMO)M(wr)cC^oG$O^*jPAFidVp2_blS;?tJA{E- zySyjmwQ4~kq|7jlO;tC2h^f3z*Q3>YM2@`DvGfUIag2At5nO*Lw($CMsC7-%+MsO? zY6S-bJl`5(cz>~R2njiMOh;E| zJ4(<<3%8P4#fTsM-d1pO>$64$tXb4QjIce7@c35XM_u?Ql`3O-cZhnhoJnWjH>B;E z3fzj>6j()IaXV7Eg`4*(NPn>L^8I z^4Q9{*pw8_@s+K-dm&cE#fuB8>Z{46CLWKJC#?x=Me#x{gp!gK>^4guRf=%HWxYT8 zHe!y|{Lks>-@Jr`PQN80Gm1CifoM}zVLa%r>vgUIuF8Ri!sl>L89V!SpBQiT=b2A_ zD%F`>PLqGVy83>Oy_(I{cwr%363c032~gz7qM=kT&LlVYlJ8s=$*-T)>YBV2fDU?k&7V?oQNM?V0%#vW+y3us zPW`p6^O^Yu@C6H~xiV&h<9QKKbwe&K5D2l4DUFtqn{-ps41c|L8bVbnuTrAHvz(H# z8mo3=OBQ;t^y)hk^2~RB%_P$oXVaGcX0M{se&tu|0lU>)Vs~z`uNbs?0d4Q9n#hmj z($^=nbENNEZ19NZucKmNnT^M$NaORsF&^ej1@+{7kzcCC@l_7sCytTkW^=^~Go{@Q zk^H{cb82Q=P%IbYakfMYBdlH5cmDelh(OrzpEVyzCjpnZ+GgwUE{{mxk7aL)iY<;q z#Q|Gb<3YuwcLHg?+vg)KLauI2zfVBj1$76Fv@Mn)%Jl=dL+SqbiB(%yE?%A-@y~7= zUysWHsPO8Ko(q-VETB?ocQ~`Qj)3d$W7=ryD=d6Bt-W%Y%`vL#Xg7BRi8R4?)w`jj zKEovOg2bQ!+^u-C^h595N?+GS;+qIKbd9G=8Gh7# zTRmBemPy|3j;7fh=ULfRaf)B`_y2z80G8XprCfHu+Kq-b?Di!fAgH?^nWXUvUxOm( z9e!kE<716o6V43IQ4}7e#sFb(`+TI(SP){?4~{i}`H1N~RQ)YVPJ_-LslhF;@ia4?cg z;(B*`P~3N@=CuT6a0&gu^I+(1)^sL^t&A{4i|;5PA7xIn}=fPyuqyjDd>A1>cq1` z>j@(TSnQVZmv+zJGBO$=ARx4N<`95y1ZnT?jsq7n{x;D-`yNFzltJ_B#XdroVh%w^ z$AZDGvd`fH59k1(c6bn!NktwV?ZWzcaWJXT4_0OWQcOYuPb_5O*GjtsC`KvC_@?4I z0pp+Gp;mS{27ChnO%3bou`S@HXPYCLVoTjBc(b$5jgEhzxVo~Fvz}ty-!shAIs9u_ zaRcQ+G8k~P09<5~ZgL+Cj_?5u2qkrnW<}6EA?iPea5qW;9#&;MQt$NVVRQB9M2)dg zR7n_EXZ|=0nix1Wvj5%^+OqYK25y#U)f9IDpM%N-gXY}7>8q+%XF9G&=5o+@Hfj1j_wQTz)Oo zRYUg-^k4A5a5*B7e@(!t9JIRp<7;AY540W z_=*Z=O2K{^qHV+3{vuFlZP4(J@a1-kxf~PBb?nnHSomAS>Q@_Il$BMVHX>uA4TB;EVMK4LZVU#p=M!874&jvwV3k! z%V=fftY7)}_=5KPP~72ckIRp=QsmOHkhPJF1nAD$*)dz-Xg6E(!xrcgyPrQo3Z2+{ zsruMtiyx`#9k?-&>2j3N2pM}F6@`hSvqteI^ zT(%8?T8bzCKq?9@BhU;Bl|Kr4{Tx%&naX^EOz`;bM!E%EX?%}iXVmZ|Qwu+yn4|%~ zh+jyBKo|oB>SxiQ&ncK5%=><-`151nf@Z65!{hHhw~lali1LY?g)5v|m&le>U& z;p4mOH}f^S$uD3?NRG#(-AIe-z@xX9yvVC9(8Lw6b$z^;4ZH7){`Oi)Lfe&acFx^% zIut~UoO;QLm%|MFRTXjgRQ8^F|#T3XIFlSpB}*E!6CMwL+L0~Ikrx$FPY^&QY$ z_u<=`8Z?kmAtez-$qp5gtc>hPR#ho9Q9Ty7>;tw*m9w!2vZz4T<)F2+pK@m2QQQ5h z7x`yi7*q%wN*O=SHt{k2^h9+_s0i`Ysv;?DS0Y*~EWR^XG+3=`FqRI(C{&#Jm5OB=)L##%S7 zle@N!x9^upo|u-$vPsy(xgR#L#;V8)tlR8X;g?uKW73mZTN3UkIg*YvNbc_(YeT&kFm#x8r>J!gO*9HKth$wVis( z`{UyO?Af!l)c9@X)`C>sTvyA_dIzk~lPJN#u}7h%VF%&DUnJB_r2+KV^<|ZCB{u&M z^_t~p_out#4#d9UY@_x$avfjk;Z7d|}t zLP_knw~!0(+vBzeEmJ3&eoXI_mb!g&BXww4*zurIuhVLMsZW; zXV=jkyyCTT`;m1C8mVWFHD>72lWTLE<=&=Uc(Z|3I0D8OI5=)Zh`Vrs__j|~1s3EQ zT5K<}wr*5o5cS?(8FWiGx3BepOF=fLZ0afriHq$5b|!pvhVEut&b_XaP0c<3?z-jl z&>??c-|M<482r;53)#4G)jaIM(P(P_bNVY_|~F5PseWbyq55hrU=itAXznjtS-J}ouZ;u zL%dOWhOXqPQ~oXZ$T!aG!Gz$^=$kZi(XHVuuUzZpMZ$-4wIm5i&x0qj^37Y)Z@#eW z4Q+3~ZbaTOf}NZWIW*arSKnWTnh10X%<5t~9?Z@%!`j^j`L2>CPk5kojf?JuL~;33 zpF>r(Y`#hId4W&t(+vvt zF3x9FmX=;a9pXcwJN528^-oE_$^8lkXLr!*)w9{QpY_n)`^w#Eqekk>y}~nK2S*~# zd2Y>e)Fau2^QQLA4OeJ(zG75V+&42uTh^9;Ks7GHh$Ogb6BOp<42WPZP-J?nUF=1_a+Sz=bR8wgNL>7wybZmbtrQ~qO-e(_A@+Db-}vZz&nOn86Oq%pduZ;u!zxkAt)9F( zh>g|cBd5>q-L~+>0y29~f4VTyCtz)(Xuq_e=H(UTG#pO(MR!&5#OcY}783GS1Ux`B zOaoDBW9U8+kU%36zE(UHG=0xFGOFnGhmx1_u-<93sbuifQ*&p~=+{BYmQP!8`q(s5g};BR zmq9V0F8GFw!-*5Z{ok}6PE2?r%Zs@|okBqdS;<9EEhs3+wrV?Rc%@*)f1&H&-$j&P zEGpW8LRtCeq7_@WZl3r}i=w(WrMyr|158R*N-O!)*2^P0s#w(q+eZ5oNPCW(n`&4uoYFYO+Zl`VZiZcMzBgM;yKtI(&cam5oSR^yv~`E9wF~8}NBCi38MdHL99j@Nje1`KI}<=U-w<~W9Z{t; zy803V7e6@GJg`zUFgOBH26!_>p>QyzgmEvs0!^g&EJaJ05CsL*NRjcgDzSRVIq<8x z9f$1u)c0!+$$BklpvSxJ{$9#62kw$*3P6A@LoW}-;y0mT3{}_Hd=|ap6+&4aOG}yM zrHyF@1rqDmJ4_EX-J6?YqU1UM2>Kg#LD(srKVI>BgZG)R)r|aO^&xwlC`T5A#J@AK zkRS3!{o#>auOYe@|53OZ&R{#$j_K;@L9|QBW187Fku%5H(q|J}4CQTq{0^TUlMH=m zxLsCVefPC%ixEfN{qQz0V^^^$iB!3W(?GhM9%fp}neFAiWqw8(cE;EbD5J+r%;|gE zIk->L$=q5y*(j|+A#TgAAV+Ua3`!4`cYW4x@K|+5CvQLEyJd7B8T)~({rUBx;%Fxb zt{dzM5BV<7OHRrA9w5T|s719qKpNIzESKb>Bj$pnM-wTc!SKVF2QZl!C5phCr zj(u*n-_`%loqjaWe;}0Y_ubU&Z==3{s}2Y_?K9u5^F&v0bkeTmW<#v<4t);NM_?|C z0MKXok0iDZZPwbpH@zXtz=ufavCGBMI`GdsUoNs?4=#`XjUpis0srT|X=KXKpYWgN z*-#g?Cg0Y~uB5-J`JjHW0aU2^L@vj zKX=Psq=!V4QFEoHKe1bE@M#!L!$Ws3ng+axiJG)cE%yzy1#2h&*B`juKUW)xKPN`zWA_Y@7-@ffY%-NYty`mpI-csVR+iJt>Qvh6+S18~#576^}oKCma zg;bpDS02(3n79SO$PVa=RDS#zwULQY&`N7F9o+@Wt2{IsJJ}d_=fiL&B=fe8cKz8U zOuWr*_9uERvvq?%vY$iO?2{h-uTt_s!9UTR58=+TnwoIMhv9-Il~2p7s!-?^kl`=- zbi1c$Sq3EE{>nl04&RJmCcEvmOJY%kjy|RO!FtO6v?j>x(8|!cmY2t>o%IZ@?1AMK zauhFgHvL74QOPzPo%U>(;q7`ahGaFvh=Z^6_~Vxm_+Go<@Tc?#Z+>+V-{0F24eEv@g?@{5NjZIBsMm)UY=@D~ysbXVGqSR$ur3$hH*3T?3*lpMcQ)iTP`Y}n?-B5Kd7rPTmM!ELbazPA*~ z=Wz2biY8Y-sK(bpZRd^AShIt_}^&Q}#LX zFzvbxoxXW_JNBuJB&j7ac9-xY)k(PaF{Op2w)yEIIM`g;l-AD^$&)p_E*-_*g3q2+ zKq^fjlxxkFZ7$lxi|s#~evlG)@iasV-n#gS`|2B4a-WX`t2u^*aVkGkA~}!jc1XIK z0^$TxRGlPW6>+3G-rrxJ@*=nCgjkM>(n!3yA0zEULqNyR_lZ3*`TS;w(X-~`pI@KJ zzR*)?Qi-E-g86LGUqx)d%?lMZ6SKc~KpX9^4jra^(66ujhukI*kk#7HP4a9Pp6>q( z1vY3LYZi(gv9tAR2wJgfQi@HGG*0b4cQ)M(XtmbSqcu|l)ncGL4*u$|rF>5I1w)`+ zYp8Ec8kAhwWM?N0fz?Fh@cyOAKFcqhu5fRd{+XhD)cjx9*?z5zw{epCXf59z67mpV zEQosZLCV>cyQ4)=MirH91@2S3X?Of}Tq?xLlVt;smhV9PqEXhZllLAev5a0gEQ7?T)eC=H%$P; zHPrpquQ{IiJ6P?R5^+?r;Y2y3$TJ`6J}!w#b^>Mr><|eB71h`;CRS0A+*8O7 z(~?grUr$-25BXW*36G8kv2|8lAYH2`Om!hq)t4l}!%s9oD$Tc*4C@%WzS#qVb^1rubpM!LhMdn)dTa2}rIAc$g-J~1{-(B-9Jqxt+Alzv^}jDDSf zh9eLcNta<^4ymO>&5!d2Y96qXO%Kk@IA8lXIIgn@cG~~C5J8ENO!}4VC$^t{Fu8||>w0Tu7kEO9 z|9GRyWPcsYo0>xEqqRT2a!^ERCPMO!gI4DDnrG*HC~E`% z{OnwnMry-4207z$|Ad95`CiI5{3=G-#|MZ5KPEeXugpd@+*(xq(x${}=K}$G{9Z$B zm0DKM8>+%TJ?kJ=y8MDRJ=ksWcN>Bx@ImaMd74>PEdmMx$`yLI5(*wY+J*l~K>x?B z`@B%DGil9g{6uzOKpascf|F9fO*C&GE01qwoVfvz%_$i@+8@4X7Z#z`PdEg?#hy55 zemj6VIJCQ01w;N~q~q{g*(*$pjD9Mqd$_qz4WHSt8X`b=&_jNsR&j=0N&~7YVi{;l zNk}bBY@~W=J;%O7Ae{J`c2ssRS(= zK8YqS^h4POIf~k_fb2m}Xm^Xb{A9r4Z4nP!Nz^A~$~*(yM|3V-pixhg;1t!;twx=-z}Yo8*w{N;}q5R2$48LzJSRrG}gr3{qDv zJBcd`ZV=4585_%rmZ?Gf5+@H2PyJlJO;=^KU9m!@ULHV=b$FQOQ>}YYDvhWLVo&W9 z=H~W?`sykH=t4P*pPOuqIAo%3<2g5J2H_#U)YN^*7GC|+lY;V56xUjfnxFIc)ipI4 zh^|W(Y7S6j6craI;>bPj-&u*oc;RPAud#^O?u1-%0)(*J+uQrTeLIPoD74_c(@hME zy;yK&p#wwp1dbS}5NNP7;PU7KVBVw`>ZmUYMR5VHl5t&gVo$Rly1w0G|rERVZnR>`S8l`>A z+n*!ghV1d{iB&>p= zPV{j|Dj@L#RJE=w%~_;tXCabjMG2Mz%rGG&)|`h9H4R-MO9Gjwa@0ehAElmjl9mcu zG84al4-!3ji0wI#w1;phM?HtE8BxcV7%vRjM4_?WYRnq7tLxHfgtjyu^WnmWs^p0K_Y4^c+hA&i~o=L`N z8!c{(Rf)T^uDL7syHRU~T5ac}<7JWJzxk&B00FRf~1Z5U4(4*3)CXsDnAxM$a{ ztGF0bRqQ?!t@?I*$NR^}8}W1EdZDH^J2#gSUn2^d#|_$#2+{TMc~dRW-H{TObAy{qeY928I# zQ#+R3F$e*^`yTor%sZE4*y6T`D z1C5I=^m==B#oOHCU6pW@MBVx)SN8E6Yl)6sax&-A+;{e*d2jru`Zo1dM7ThB1gU%$ z&Ac`xaLD*2D;y=$O^&)B_9%)Yl65Wj%#i}kLZWO1v&3cV^hM}YD?-MQNW7Dc!=@9_ zF&cKMqsm`HfN(Qa2C>J-xRFT06EC0=ceI*JQox^7=Sf1J7uln9OE#K{;`jfedY+yA zg-ozO`CQ1#1FfkbrQ{5EekpB(rB;Te^fh|=q8j#COsNbG+PQ#Th( zIUj2zaa#+ozwc@uEB?}8gXtg5)QXC=1B2kiSCp1shw3{$E31ujFU*TNR`5tIShgFV zsQTT*1yYSN!_$b&cZ{TGIAoqvbbE2`j?`}2DSAvvi3Y-v<+ZiE^&T=)*)ZqHmlD4Z z43t9w4dSt6uE>WEXNcv>Dp|U8q~|Uq+@#Mg(^H|vp}V`AD$jh#-7obsgyT==j$!{% zB#v|KaleI<1U&;oIrP;6^YVmn!LZtX#hcH)c0TmF*s}TYRXu&P)BKagm_M-MOf-*V zyjwC^C|J8g=5@?a$5c6UW(zJ*wSxTRtcQpE4ugX-e`H+nsON!7y7f5l@KViuS&-d5 zzp%)Eichc?I$4fJ2znu+#&YBdd#Y-q=~%a@$Dg>{H*cQ(Cb2Pu)@69p+tSjg6vg>$ zAgNG0#$jhki&F}yV-bI_Bd&6L7`NJWaz`O5kLwoIyI=a z1hj+puC5)DDko1rJZR%w0Kr*IL=k!2WAC)+mD&%b{(D5!8##)Z#3mD(`trh}BEjx4 z7OJ33QM!9dxZ;GVJa8T>4}kiv>(E7JLGAEFC_jT4|>(=SQusw`sF40e(JQ?_Q zG9VG9-W9Mf4=YaPUKz_MUweuAr&`iUNbi-R^g8k94^jO@C`5GzL5#IDkH`^u#Do@p zx6hhaSKsjn`Lp+F8!soPFYfmT05N241ddz)XelrZ>*=`!tlI@bH=PPGH1Hd&Y_<5E z$Dm1ajL8RR!oOga7lHpceJQ0gZy164>c;FBO3Q0c_%r>QtUZ_gT&E5QJ&ebCo36_i zE$M4mdNMYgtyzMu_sA*NH#l%8^~4S|U8hmsUrm+{X@E9qy3MH_;T zHW{1KRFydYxZm{q%a<>RfcyLUg1>ozZG8~GoCKX$1Z#n?K|weWJ8ZbZn4mMd@~ud( zSE%{w#}N_Cw!Ou%qbA`P+k&KE0?LAqf6~Nu84h^!K%eoLU086nDsRh)Sf)^gZNg{Z>r?!~Z`J+ z`z1iJ>spzHd$~^5=$(Yxaip}56cs{Vuq~3L7eVGH5nDG6G(f*$8!W>cBL%6?$TuA} z?Ml({3l$MKQS>(PTG73GsfG zI=L203WLXD(;5ch?3nkG{JFzscN2{2kc2eX<=vmmJAw|I+NvW=ls!+u?aXx?JZ)_~ zQgGRZwfV6N*v;Ub`Law-UQ#Z(!khj4Ig}t>R*s;r|?d2@R^z-cOX0uqQR#xLAPjkm`_JwX3*JmW2^Drl14HZ>g^7@!jQqXzP`S< zeeZ6d@#6xfq@?wX!}U+ri}X_ce4e!yYNK47oU0WT6|oSv(Sro-T!6w2Baor4U+wK) z@nJovfYpF&_wge(>R@Axn2GZu7%3!hKvBd50jMQi1kS=^)%b?r!)?mH z=B0Q6hFAFXpodTkAxvCZfJEQr_FIBHHN3ErUj=Moar%F0M%^*lt($=a?Ec zQ-b6C@gZRq4szk*OjlY@G_p>j#E#;5NT#FJK}rU#pQ*QoME%w_tOX7xf(Z$u{`czL z<&Sa1Y5(8C2o0RED)k{8kpg#>9>gJTK7%?Te8^h_IO>_h?b&qr5ZB|Tm=NijfNGYKlEWE7I~Sa&2MIucbQcDBSkDJVXhzcj4Efl@8<_^nfNMH0Qf z$(YSW2w~9ib;8Io5{_Yx8PyYSiR(9R6!$d&ic*X!A? z-?Zs53=W2x#31_lh{k6GSaO&lb+I)gvZ-MaumTaL$5Cy@!s63Vfi;F;?bXhLUIFS_ z>~k8LYKSkynaA80DhTx>Kgj+3WdRPYvZ|_T2!bUyu8tJ(uN5?}dxVzQIw;#i-?%A7 zk%>rcL;n+oWCL=5wQoOuyia7QF{NQ)Vd3KE*FqS(k$}HTTqVSXneTN8T+oB|@(>vM zDDFV7*YvLknxu^&Px}_0{@3JWqFrw>35%hWvK7kNs?v*Lsr8+c8%bh67ggdTC~RvG z6dNm;Z)<}i)Oyr=JEb(}?jZ~^oQsaGm{Z?dVKjKkxV;L?3xn1^xmzpx%X%K0+Sn+= zd0(|qR-%DFCuaYE{!3G&IlebUe#TJJ46n;b4_B~Q^f-j z44=WbQBae0e4XrqPDNxe(gqF zLnd{Q`6}+-Nb^L2gdln^ZDKz|YCZ%4fKq%<)T>wWnpuYNHeIdiFP?-6X}q^K zADw27<}p|oQu3bW`-yeAFbUoe4EH5s7oWQt2dWSF2MUiePo1Fpe-tcA2c}rS&N3#9 zi^R82^z%1iHc@3u3yk6>NC8NPgt$_;8u#}OgD1+O-}M!NA~SpZx9);H1YdQeVEEV6 zbIoZlkI^JuY`PDs!gYnvh1>q8BxEwSi2c)I$UkqPeSafUVE^c-owM)-5Yu$`dEKbCB?>3}?X>~Ia z)cX?=dHc^6wikKw?mg{Bx78sIPsA*Q9D*WiW<)R%`GlsX(BZa5Jh3W8@t=y_22OG* z)@Z@_wdY!FOi|PZgO&~~?8&}&?qrFd{C}JM8AOSL+3*5)Guw4stLDMe)6UaET&gqu zJw?qFAo-4>*ohM_`@X6_z<2r5-(UNqSvUTd4w&@>NESqN!$T<16#@M8^iH9q7U1;V zmxAF#FAuc<6cncV-YJuU!mI<2)f}2*cM~gyIVFH8wu9k!=FAztc>@YabW0_*46!{R zf#f%%`Pk0!4FwH&!;M zg*+Lz_K->CMNFI|NRct)|J?@vvoTX&o@_sdd=(vR4}>{vK!%nfWN0Zqvjfs|{Xx>R zfa}~b^#o>(U_#KAHg{i6PEMw=p!T^^dD-cEU*ECw=XYJZb`1%`Mmsw@wEDkC z?e-&(9YGz2XTv%&ViL@I2Uk~%wr9$uGXY}mTef$`9KX>DqJ<>t?dF&}7#kd{45n1n z{GZJqGIX7UZS+)ctW|HFyn#(#`P`;U7zkrTrnjV|B#~ikKm4OqoN9L{?dV7ucKG z_XgTe>4P!`-FqG@TfNJ(dPBkGCEg^W;C1b5MbYiMUcMzh8{y5sG+FuT%A)wOm-L&d zpbn0?Iqxu@Ow?hwP3UVVR#OR_e|t4%jpf=)ZO`fn^n>+q9qaxA*s4{WGZfw81hu{I zH9)Cgm{+y5JiwI)vZ{PLXtTxNX)S}|!gD0F$;op8AKWF3yD4mF$R1c$Mfd^rR;5uG z&1uX|0$u`zDrLlmEJrSGv1mvho3I;dwnx%`RnvRhLHNPJDHggr-=nozqD&y~d>h(K zW^Nc>WB%cVrAAC_Y<7NrDnJ7-xSgZke?JY)AThDPFim7cD?L3uNumo+BzVsS^9x}U z*FwY_lQxKp1GSX}V4AL1o*3%n?|=#3M5==~F(~p7su#wfwnPM83_R+q0mYC!5)e=w zprB@J2_$;|VEvt<_E)i*WY8mYyMMr1J6>RcpQC2?x(XN;`sK(%Zov6G7K=;dMrMj$ ztiPx=V>H7~h)^eDu)_l8ZXk$E>DX5!#m+NdXN`F!K>ffgLO%W&@mS1J_sNQb7bt{@ zQa)0iGON8*SQH%^yrsG>(&m-hAtlN}XFG6-x(THyitIyq3PDbeg9++Zk@qkvfvEBW zE>*_Y23{wDq#OJg79_;5(=YehbzvzHC16q^)v2-6EzHdRNl7N*jU5@fg7CIgxC0t4 zkI{&@5z`-2eU_J4#-x$ftfk#vPW%KASLm7(Oz|$TX1qEEn5x?>_q|&cb;OcnQDC+w zBBlgfzGdgmI)rgCWA4~SIV2z*uU&YBE`550uZiRljOE@w46fAF)C@OrQqg>75Q#B8 z3EUenMnoZR4D(r2@fV4Ehq(Z|>4dZM5zu(|K&TubeAUpQKO-Hu{>qrGwjHqmhBw?> zvyl;v2^#^Qk!TIg8!3oH2r3Ii_2W}*{#H1~pQ7k+M71CYDU})~gJO1tFJ=@0r9F>` znam~Y%5uf=S-pBSq>N+;T16ah-lc8D#l^_(CAzA3ttn--1M>3nxbw$DxiY~9GoKlm ziIFZbr0mYi%bWuCYs_k20{#vFa)VKWVRdp&z`BrJ6gHtOLFU#hE}H+{ zyUCc(g{f{Ux6olK>cxwWPLX7+O8m%?1e9LLE0a-1r)_D3I9T?wW5mP-tT!nR!>z7mu+-RQ3q@2gtO({$^F_3g|T;>_kGamQ1ifsftWW zwP=hZ*yzb15!5-za6d9r7ciI@QVD$U08SPedP_}BO+Fup1;$|tVvLTg8bk5?*c;8N z>=!Rgaj~$-(a0bKUB?W7SH8sPK6P8L+D`;GizD`L>{i7Ns0WbWyjzpiylkNAN?M#=h_ap*i97 zWF7EU$l&gx!jP)>XY9?mI?f%2x-`IQDP0-ETC5a6f=7jph#F|)Lq>lEOq7Ca2_TM@ z69noma2y;M#H2a$E(B%_Vu!@dPAZrku;f-D*aA#==iWW#NhwUQNWyZRY&3aQ9fye) zi1e##Gcz(UDnpu7jc}|FcV=*lFyDWJ@NO^vJb*>ybYEj#dlFzps%vemt!MFPD1Mn> z#ZKN5nn48vKUG*NL30L~HB1F0#27<0VbVC=y@668(l7xsL&@&)U%W^Lh)V!q@nXcc zq$)y1V)>VbRwTqn;5E67DQQTGj=*U|mJ;5FZ$~Z#2CpE@o#bSlZ}r2a)uD3^mV}e@-jKZS5 z38S#$u#fxNVvc6v#*(B516Kogs@<}Xa+}x%DL~8;(;(Q@d~#v6&@0bB2Z^gpd1#o3 zdYWZxS7#?NdL&^$B$G3hV-%6#9;K41Xc*o!ul;ZQud=#2orLR{Qf@Kqyk$7c+AR~9 z1^1x+S`-zgW00<2l4p&0Vt^W{Vx7FpeB9n%0MH(p0z=hNd=hhERdVZsgL`1b^sDyY zHGOgpVr%DZZ0ZBHP5tQ1lSYaC{Fl-==TY(KD4~1OXBocA>FVe-XZ3nvr5X#IHV%D$ zLybJ_u@QKP9WEGJa8*-vM&mG_h_%Hf3@Yuzu>taP7C9O{4htrsgpSQ%q890jrrUcO z{rgtVJ21OAaHrV)ow;9(pu^%(n`Kb2@`-Ycde6uXGI}~`eLxzE?JSH#RN)9f9vnlA zY&{rf6n?M^MQ5wR(+$ zuIJXs&*v>(F3BOO3L;4*ssq{H^Os3_3#>-)Ge7Sy=9h=GOoB<#km__epXD@<*JB|(JEgps&sg2fscpAv9P%*zt_M`Wj(AhxloQyDDFjX=r_ad95QS;OcGl0(j24{j z*+5A&Z(=s=-i?e(Nk-DV-x@a$mMpIYMN^)lx&ex&p`AiwT(eGLdjS2u6D$x_!YUNb z=mobXY-7g2jFqR@#Q(VaZNG|J|E*>YE!hgkj!|K<4b|cJuCn@iTI>Q)GGzeLsMA`$ zbeA$OvpT}(6PKM}MXTDQ7{+bvJXunRx|A3OclaA-RkgK=xD&YU6(*Dee?6o(f4^&- zJ-D0weel&}XDo+-MAoTzl1MO!Gcq#bi$-2sb!){2Uai%+7D2zZn)$w84}Rd-=Y1Ma zTLnO3Xn%gb@*a%VSyZG}#ku^w5@%ds{H7c3_AY2O7=N1p`W=g{m6_Sipdi)BZn{Pk zo<(tCP}#l>lw#Ex7&e4OKsGnSASmVnws&ggdJ8os_f8fA5e zcszAMBJRi=x)Tg~NMLT$FTbs^n%w5m_@KzEqHwC+R3rK?d~K8{0)In?mw4wzI6>9pC&T7Heusj|(| z&w*_K*1e9oQJ8W?=G|dX%T0KVy+Ji$A3@b55Co~k!W1W9naEgg6dCPN8&-*zQif|X zx=2^_{d0*A1SWSK%lv`9|L>>+R#%_JjMNB8&&yGqB$V|66NZe{L436igcJg$VEafk zhT_x0oW+>b^xQ<*W^hJnF~au{xTTbDCgB_^6R@Y9 zSXyKwViLKnyJv+YE`tf36jV8}_=j`@tkM4`vUOtJr__(|^X`+OL#z^&xFGvF*siPj z6IQw2%Y~f~W%m2Y?z6M(iCR1%%Q0*qD1+<MzJXYcH5mZOsCyCzxYTq^Md=)cXQ}17z|t1}|YE9O6Q)EJI~*Qt~z_NAOr) z1}#kt1MZ^lvd;74y2;$Cgm1EE*h`^@UFx1y20NbV!y)%cQDiGB>i5Hhetos=DXd0~ zo$3_+d8URQgAa%KIH-uF+I?=1d~xOMFAWH^@01k3CU?gVd@Bkq&fT-*WL2OB6)Wmc zXma4-aqi`U`VMd(yyUiqGEwOBo|x@_Fuo$TV;V&=I%RNA9?fU8VnNMev-4Ch$f&Z z6(%#(eTos<#TAG-Kuly>zdoD5LTBouB+)9%0X8^c1DaGBmRI6?>7qRX__4B(V^ktH zh$3gop7^?v=562iiENUFX1~6sK->p)9(7W;I#H@YL8m-^YybG{$>UK(ixFj@0ofk) z-fyOCI)<Y0W9Ev&(=1r3n7z*(ANq-aRD(>d%LLpXs5>s z8KQrdMyQNupLGW;Eoy^3*ciw&_g(W`piqD&aili6DFK)B?m>3vWrZgEZ5 zF(*MUjvQrjD+CgdbL3dHL6xnRj_ryc;NEqcOQ5BPzPvZuX77U!w0Aum{MjN%#Abh| zhvi*2#n@Nid30O1UPG=03e>xv9yThV2Y4#RfpMdw@-fDkBkmRbY_gY=(;iv0W-a2D zl-z>%uH-DNqo0~guLCn3gw2&(w=xlC`mN3Ul-;tmeZV5_%o`XOhiQ;2z;Hik7?z{8 zWu&`+9&Q-}}oX1RE%P4us9TIE-=|8w25PB0b!QD6+gaNDy7?!DVAhhAzZfSzQ zayTBpx#2$vj{VRti6lLQKGMUi@;Be_*Va}Ay8IOoIqHV}{ry$Qxv;gS5K+2|%EuOZ zdKrYf>o#qw1j%gZCBkGHmcy&z_Q<>Z`D1}Z+w!Ggco;4vE)i{L6~C+sX_P=GFd`M0 zwVRWJPtmv-3t;+{>-Y-_F#T0MJV7QOkoLX*Anh~_AgxJAFq$DSft>s_u)Os&G{8@| z%DfTX{#RU+m6vY@B=izhsY}raEsZTLBQNQK7535j7}E*pOTlCwMBf90=ZvpfnmR4C zVi!wh#*PZ4Po7zxbv4g9ajJB`#BOF;wXJfkbTi(ea>v$ja2Qs+JXn<2vwud*dbsuJ ztW{IEHK*^EmDi|@*RSV1c4}Qf;<}m(VcuIg)|*Dx#q_dk+7{XJRmYf4dM%kaek)Y1 z-xzDr;NC0a>?~l>urOJw7BbxO=(ytU-8mTLx3Tw%XFyNUjk#%!`n)Jp6m;%x7Ydl1 zbOez}5l}h8tF~!Fpx33PrG#K;W@d&Fw0G<3>YhqB*4L-xT@=r&Igc zcke#<0o*)CC&7ZL2I?nIL{+Uvj;X}Vep&TS)>f;?LpefHBJ{B``=eL0M;HIna!MOH z)hkt3SEEWFQd4sp)zi0tcdppq`II*~J#7k*tecd;H8lkQ#V*a{mw1n-SGm4ue(*__ zu}t9Z4tEnZwXL9`F)%WI9UnJpDS%YBy0J0GyLWrq--s}&s`{AMwV&#G`C8~ix}mNv zRLsl7oVYd{8TlMOER6A$s<|dDVyLZ>#5=WF0ayugRsrGY(Z`L$_M$qkop~UDuvhEp zw$$diqWJmQSs>Oa`&kqR^bD*-K0Mfr3%Oz4y7OpQ?1o0*^R7GoLF_f!Mm8c6;jcK@ z1*?`xk*B_x0lYR)}4xR)G8 z+8rP30|#4sy#JhI^MUSOSznLmxEauWy;bJM)_g!(mt8M6h zvt;RaR*U${_xH|_Rf)uzw`$`8+lzswtfA&X=-fc)pm`OYh-t3-Asaoqp4mVH(UYEU3!rl6@m1rv#^6yR|(^{WMPoMCGPlu#}hAq|4&qa z{*sgVfYfcFDSMZGwt03aY43lZkCH8&xa+2OKb*JF9ljUAT%WhYjCn1hYXaL+BQPzd zVHN+hw_;TvCigDPYtUzC@&4&{+AQg5*r7G|r({{uXSR&p;J(~|9_3kEpgDvF5CNA< zqJ5_IQITE~$6pX5pqt626va?>&h@4VvkEE))$`+35v)&J57FGPh;#Q7>2X_i73?ku z&hC16&qzzF5+#L_KT*v;l7ScC6~a*FpwhN|_cv?#FTd?XlkM)arw3(ZlhM(j)B{IhSJd9NpXtv4Y{q&olM z(Qyn-eRHHN=(FRR4>@Lve2g2Bq!mX&3DB_LwPZeagMWkyd4gZIUS@gaLC)P@6{v)5 zy<%b|mgB3_%2{>!v@q}CRKBhDw@wXJ$p(Ep%V&#pxZ9v1FV9vz$PsoPIs<;|@Z?UU zz6g_a-m1`1%CKOz8;Y<9o*R9@#Z#4~K(IV1l z6tRnpIo%ynalv?c$aSb^vBn@;;*Wx%_gSaWRMA)!aBLKuryuC=2J)hiWmu$C?T8D) zUf2q*2Rp-kqs5=v#4lbYWQC!mp!BQ$mc&?J;=WT+3XV;K&(e2AD(lbU4-e=snL6%w zkY~JebL$?%RP8LL6-Ye~f*ouTs?lMa@rl)NK)gANGWwWPNMv2}q` z0*>OQV<0?%&{h7!`88+i`w|Tt?6h}SRJ$Uf_}=BEZD{B4xm_Q<{K)P=N{ucbqu(a) z2gTl5b$N4BMgw!-m-U@Sjc|d`?n(zUt;5GipltvmO0`btLIC4;x~Q#K-+D=d6c&LW ztqBn0_>-5zDRdv{(r`_jZw0#7Xc zDFGPX_9TPP;?t9o)*qX1s!RU;cZP!mHprRlk4#_EO$q zXxEFoj*U$2WIAZj2RpK|?|sM~>~@CM@E_aQ$(fnQiHX;7`4>mV(~(-VTMMIJSpjMq z`4*@$mNqx%zFxNQKYz$}!sF8a_^m5;+H1CMO*6DBDGKgg9BZR;9{p=o8R*S)N6xkB zbW3ghxlF1aZD(3G7Y`jC);`4)C|E_HdAZ@ix3_xJaYcKa5& zZ+5>ZAj!PVLGB81 zZKn=Qm(y=un2OKV;nSKKy#pk~K1Suh&IjpjMV_nCT}q`zN5=>~!&hDW>cxOi0hFId z1c`+u2oK2YO7*7MI_uW0d-cr(Pyx~H_4(DV5-WZP?J{^VU0q!i^&;-yzaM>%cQTZ6 z#nSx!0X*})yOfuWxs^~lxZSL4d6~~a=1(ts&FnA5=Q`G$`}oa2CP!*&z5Z^%yMKQo zN^=9>zxRLna{Sb(_d?$Ot~XKCJlE{tXDon;?pB6LV7kh)xesi`BgMqTXqW~1s3k0K zYi@&w3en>`;=FzRdbnUV^dG6y(Wic~Z?A`E9qRE+J!E`>j`0pPa0RpOSC;LgV?u49 zXhS8?8?GP-MGDR%N5W8cj$|E92Ww*KYt)WA6(@kN6by#zKO64O|9p7|_*ei@gX7~> zQL~VflM|uW9*|EX{>2ehP}8}|2l%nW?8F8-IxaO`U0d>9&%IVMog8?e!xxufRD#Xj zr(rn8<}?b%9$`}>R0gxM0>=bofV__E;(Y=F6HM`MzTyu7v)K+U04Ov;XoliRjzuFI zX%>K6up2ViEiDJrTuzH<4j@1a9;|WZx!+zzwo+SuAw)z8XsbJHc6z_oQxO)HG(@fm zt^Ddz0P;bp`wDo!s;diG?{gzBLqC@7!q&5Q22GH?NWUp-;Rm6U76B+RA>gJG!W{a` zx^_s#zh6$7y;4e%}P7V%diELpKVSICLqkPIwm!C(jTQ;{E=g!@aEs8sS zqw-UyTn(WC0!B$Ofq+cfb2j4_w;_AuGn7o*u#g zh?DemL>V<&Dr+6_s;?fP3Yv?53yX%G7TnfmXD2;*@}X2M&4X29llk9i+IH*T9R-3n zR;)VN@>;G2b>>G1yBHW42ntWi^ttBsb%DP_$yW`P7Nog>oak(alMcue$8~g=f$iW0 z{ek}4qad_jEdP@Y@kh4ieQ(cyY!>-4+;uf{U#|ac78Umb#=z4pR=rXYprt^OZSuIQky?D>Q4}z=bmd zqJG@=LXZ>|;3r6}Cj#lAUK43eLnpV#h+08`={ zfRc*F#_^w-mk{Bb<=xvircBQ6J;{k(=eU zqJyb{N`O4#5K;TStBRG3e-dzG)Tz5&+F}R3nKvgZ@E$bmRy8&C2?%&$ubFO`Ywmj3 z^Rn>opLvOo&Pz*rOpADjOuBBa>mTdpzGT?O%6T}4or(X74Y5edQLDGI z#VIfTwL!n=uFX>IL0Kij7SCX4aFGA185k@SRW#?dxFW6KV+2laeA>^)DB0)NS2YC2}T*`W(p}997J^C4Mg<3qJPJ}&p0Qi%_&d$$w%?r|R!=G%$oXAK)3$!R=G+sUE zO;sa07p_~4u6<3}K!{vKEi72u7Uc=<-CK1#;Aa&`Y@++mHujbXZQa_jdhKO2Dh4*^ zT1d}~XY~|>$Hg@uN#x<>O+vxK*UvBAO%p<9>;y{y2NZD?@s=4Ge7PnHw2#RX<+H&n`zV4~421rSw6e6cBn6OwfDHgF zYEjPwGdMjwdtd_~3sUNWIpyZ%HJP3Gf(C$&5gm0Iil;IL9<_(xjyG2Pv%LG9*o6yg z?I)ATmOV8pVcEJ>|6=zm-kUeS4|TaxdCq8#8t8+PgLSU5&2KLD_5?x!SDEqg&_I>D zyZhGd+fO4HAlL*bS@3T?c+4uwV;perxkN-l-;DiKk*jdOU_uRbGFy_Cpk;J@u7zEK zW%C#VQ+%pX$+eNrvqDZbpeZOt9l1@AiOVDJu>t317B&_jkPAfyPtcK(%M!>>O@9Au zF|WG6{UmO*EIIvRV$q-lh+LkR`}%dQ)AW!Ec?fjYPPxsRM~RHr57shaeujX=%pRwj|K>(^~a@;J~>|29K+hh;|&drFENb~G$g-_H6OA_J^U01`Jqx2S`z zd3u)JuFf>KgFF#b2hp9Iz@t8W;>2}4!tX<`%GRv~>1dNPHOoZWfge9qIK|u$vL}-A z45*QO`0)G1C86NS$y)Tf#Y|5RM_Y|Ur~xLI+t2s|NalFs=JWd ziG=?Gv2AQ@JxcM#*8^Gisq5-uoE$x-HdKA@M6~7poB7+ z8@0<-eQBw2U(uXM)j?zYo;OYWy+Br2?xxja)dStBBrq$BxX=7J0C#TeN@sCUXZ`>hbLUpX*yTjF|6(;bQ``K(E z!%5JFRTD6ilh|RK|X_S&*h1(UMcEh=HHr4ElP}$#qaYY zO+>y$o+c|M{@87Pn;L(mF2}5M;%%%~{uOFM0ou2F_sY~}OUM0bD&$N#enz6ck+HBc zx6x3SFY3YXqxvr|e&1f;u9A9VV1NfDMTV<^tn!oNMt=F0&6UZf>gNddgsgN0I-Sw) zBfF;j9_5v2$t%r$<1L*1>V=FyUf$zx?Q5G8wh?8<`*gV z@=MFxj(P12boAToqb$+(-o-(5Kfpc=y0aU2^VMN;Kpy?WG@Cc`o7Y9Cpk!1Z$kT!f zpDhf8x`u{=ib|zzerQb0-CMV$p1r0&2O}2}Eo|3gVCXr@N;Nsvd;YWdsKK7LXE`~T z6cw^ls&MO_iqoko{bE-ZWDcGahZ=&2poQNPw@+Kq<;G81<1b`C_xODYG3Tw%wg6-` zIS6E!Mh6s)0m=53tCF6K!BKy$Xk`A_)HddJDc%lZKYjf8w~~M=(LlEPq~w9shK9v; zQQJ;O3uoeO4}Kz6O*QUr%Ox`e!gs;N0u~Nuxs6!lko{BB%}zS61n2l|Jos5w#CIv@FW|dUNz?we~9x~;+Irs#Ez^47n&Hws! z>(h#jef|CNZnJ!d9xFy7E+ir_18Kpp-h2;81SoQ8sH?*&6(~?Xc+k9YlQF6CW(sfG z>A#1^a_y#DV_o6AH8mS<-y5x}4rM=|=kc*%O15~xk=M-Z`xEUf(ovqMoBPOw*JUCjY>( zvLDUtCV3u?Rjl-)UT%VegWuDPq_0j)=pM2B*gg`WegpEZasR8iD~-l-{nu}a22v=) zHkCxlRLVT0l9`MlgwRHjA(>LC%tI71SEfiJwUIe-4F$+WUE*;U0e1@Vjp8v*-tC>FmU&9Jp!o<|JNi*J}hE7^mzs61tX}8Wa-3p&p%u z)`hTVT4Fr+?`dgOZxzYsK!9Oz$EQZ#b2`H>Zbum#Th?|;Tnh~~hsf80i_<-J>5fL< zFcM1UfweK#H=v-F#Jy8!_ioB26lwb=-xw72jvZ#&0gFE4-nhDp_;F=zdwYrJ8OtkE zI#Tq~gwA|^;Ws{aVDUHi^yf_5RSXQtyhhGNkV&QB7zN0B2i>0(!dUK6B8+8`)#Z)r zx6gy<$7Gg}I3CT{uVLB8DIy|{#nv*b?B}M~Sf#QvgOYrF%3fc-ZCt&fu-;WdV9Tww zEG#AHQ7g&XERkBrB+I@bZiaZwwHU%wToA9>+Z{nJh?K0Xme3kt2<=p?Jf%Zb358>vXDfVoIE z;HRXL|J*P#B5?aIwA&Y#ve~bKsd?qW0~KL?*KIpPdbl6vfjsx@O+}ESsMG zz5Msw?tp1E5$0wGn8x%D9oi8X7>IFWaTm@YNSY~FtqhJF*&oBFBKYZc)SID>7yMd#_sfTwyW`FWj0V>=%7rss$^+va_tadB?gMNL{9 zXBYm>_%U_->TB~RvC>CEii)cMgRjShDePt%avm(|IVF@E>q4NDK+fm46Uqc6+So{8 zB!srKT(Yb#xSBb6=X9^~m5>pel6eY4uz~999I-WwDVJhmVtS?<)}v)90Tv=HEv>%? zs+Ce?+?wo)q2mNzBE$LSr*!{GAz|T&;PokJ6wW<#crLfhj8VNGh1<^H=^B29hzN-z zIO{^$W0v4ZIN~Axk;Ib9N<|O`c*JrB@`!f%3RB}cuCCI$$&aYiM!G|X8H|j*u&sJr zXXbs6gS;TZXpi=IQV4m3OF$(>T)k7tkfxY75rggQc1TE!L@P{gk&uW%&5j9I7)!gu zLZjR*Gc*6Vw~0;c5Rq@VqgH2G`4jxmFaxh_QG!ET6V)pLDi`O&hth^?Yp+_MY z7JM=(Zlh#pxC*b-^1lY8Srm!nqMSK%&E+K6!Looii{nP(`dnbIjPfM9T`ek&;12 zoha=^3BC)29Boc&7hWe1(ZvZ!xlQE)_v!bm+~;Qu?qpt;J2%N#a-wzL!tixL)4~AG z`nyQ&$OHd_z_jK5@A+7ahlht>ARPt8J{+SwN`74p4Ti_S7SZOn6TPT!XRk!x>ou_(5m-$vTJ;WG1jr4}h?>tyil#`kaxe zcQVm+m4Aof`ofTQ_PTTjuhhe`59f1xr?Mx`Xhk9K@DfBL5A$lsoj7Ih*!XZr^-n4z zyauSkPzGI3O3Di{dh0>k7`7qICbF{fexJL#%(J#m7B(qxO+PSNmtra_e`so_tmnKE z{M@!R?t=B}`hUN4KhQ@_hD{N)KPeO7e>MAkDf=x52MRI%jC$<>>49FAtdk98D`*RB zaiOT;^DOPe)p>V+CcJy6%A=C*hH09ZntHD1krOjzd757zhe4bxct>DH7ok?m|KKz0ASb-!w=y*sEExFC=Tz{{rob+2;vVcu&(+GumVXVFo-qO z*Vp4QuF@XEl#pk;T$HN(H#7~F=^S{vhME7lq=IRh9Vf%+s6@&6Uy78#3$ z=KBOR(+c@bLVBs8ww$UTEY1@d*0Owz845+0r-bYGiQs-I^p?R}utjHM8TEFswupKH|y zcjS@m6B->iQc~K$F|Js{$hb{edFjd94Vn|T2l6rLL_lT2ei>eOxp~klp(J45hndF= zD6rI#lflWUW}C|LS<2E(B0?`ej^sOtqMU#Q8%LSYZxIuV!XE0(rR^Myu9Iu&>CCn0 zlRy{$jqnie*KVfx1Zle>Zd~V*^I+9GN#EMa)>G=1-Ba8bXDh3z`+Gf{oy-)XtZ-QS zJh5ng1kDgFIPXrxtSEANErcc9;QonlxrZest?lir2o)K%+R3njY;+p;KMoGoJ&T{#N$0@@fw>s?1Iylx zo4G>wy_U)Y02B*p1jZWJ6H2U->+FhgH2oBYROiZ~vty=kwm=Uv{9bs<@><|`=eHzF z2i@O=R0H)tmv=1bG5bdpf%W|KsW?|0M}x;YHs`-deXf-tv0o-d>8O^kGKnq3>`?Zu5##$6HNNT8t&Kge}nr}<+LK~~} z&E7p&aryfg|8|)lQ}~5=p6GQ@I%XO5cU$++aiH1H-m2Om(5dLMpY=AOAu zLkvI?WAX$f4NB#76E{u>(H6>OJxXLsmQtI?7GdEF7Zm5-Nl~(LmxCvI2c+dpCV-;# z4m=fP=C=+*?bq766B7K(akHbyKZ$`Gcm2i1#ei3@UK#IqURy5-M?uxXcmmqPB+*oK z-6)KU+<3!jOY^M7qv?KM=t*QAueteyl+=9hA{qwJSFN*9RID_AEE2z#t?4H3*y{8P zOTPaG{?f`@gsWLoSJz}k-e9WU$eY}^$9u|b2H!-j`a_Vjl-d=W;dBp`**?z;lr0Ls z==VsS)t~Py?$r*B;BI@F$&DpoJI3_O%nQNmk1)R;D~KP7<<&AjQz$ zWR>L~?HS8Ll=AWmcJ;v{#`~9aVVv2$yI$6J4>!#?^NLGW5~H;H6qpvYQrSMXki$B! zsi}!4xKclerkz-5_De46PaM$wL>4kST3vU{e>D4~?@$yy1t(atkG)m8zjEb`1nVvG zers~|r?LWn@5hMN8K|!6h(6m7l&Ff{W5@)!R@}{-p@5xO^Gwan;f2+f^%!0QOGANK zaZu}ZVQjn3YsgCG?(aFG+!{nbMOSpPvmz@;`I2o{zgvs8)^k$nq2L=D>4s{l<2M)Z zV%$t-98naG;P@a^ok8;tD=X{Wf`YyH75mcOu$Y(;h`faB>*_Lw>W$qgPn9<6^?YVp zn|-FLK3@VSlnmK8YZ42~pSylv{zuL^vI9Nn9eIKj?zaCB5t zRdoe2+F-519wHFVW=T zVRPBEvm@LyRH_k*!;z;~^_|`op}aBg^X%-$q}zk*ID$HYM228@1BA07`i8u+Bi7b> zD374z)HXMVfCJrQ-?EW#ytHEetd`k3Ch8DdvBl|trFqu4IL8lX%9dq@)zq59BAGs+ zf&?53@9R4Sb`68aq#R9)BqEGp1N4?ea7g4-R-(LgMGspfa%3RCUIpLA#l?lT-vNAk z8pW%TCUt(HUi4T(WXGpggp>#h>D3LM8E|WnuwG17iwN2&B7X>i4l0_%xq7O28fhZ- zqQ-%F(&LXu5xh0=28oD?ZI_o%CKAuch`OmMJF+TveG@ z4`u{``1r{mE57%-)eyab6qaoXrLoWln0Z~L6igqF!N~{~A{*QiX}$RT`EM<4Z8)4v zQKyH7hf`q|0$8BcW*;s;J~y<|eLZc@EWwZfZS6As?c3nJ2XDc^GFSAUSH2ZE+o5{r z&H^be-#=Pw*IjrDX5W+EqKTtW)4A9hXFj#A%ZvniLN%Cz}A81|QyvIr2yD|ku;iJD_J{>p9DIZDqZ)t74jcY*!xYW|3 z5gy?(HPx+i_3Bkoi$|ePpB_k-K8@BcaTwjC@reY;Clc7dyi|jG;>A*;m-b#yr(F_k z1DWx9>3s)v+9@bE-J5xT#3q);)8I#&7xFHIaaIJz;S~%wRxU101B2L}V(0h^60W6v z9*ki9SJ9VGHPIU?k@Js~Z}URJ|wL ziJpbgwj0!YOIzCzo8?}lQ{bceKrIl7K~6%mnev{k*#H`z5aJ6k>>{bQhPe||bC z^IK|s!QLW6gmGy}SN7Z_9RwZ&!r}Xqi(I4Q`|=z5?9+psf~VNWuzn;#g@wOZ9O{V+mYeK9j}!B2a_kv$d)G3ZYMP(rXXBE) zh~spyF$r(pWRIHj#(W;HEo*+h!ITo&-3@Y`yWAnnLZR&f#{$Hwh zPqMCJh>d^Dr=p^Qm*X?gFW-eOmwTP}(2WO;)cFQ12P`lFo11e%Q$W5qv7wOFS&qz^ z*o+tDv%G&KN0)VsDhM~TS;rVd+znv;A9{`c%@b|7@SR%^r6?eUn|Falrtj=0av4nD zG$_HyWtBU>0Oy+X6?6}81eo;omBZ+Qs*?z~Lxckw{A%f3v+hZnRmDF@kG;Q%gI{RR zo|}vu)q51!L3M85>BIva0CBvpBV4|16fa5q@Vv~8kefL0+vLU zZ9;Jg36Z6xrB1%i6#%@~Q&KRPtgtq3k(N#XBu?RN+l2I(&48BgCg!G;ZW(0d)3!^? z%eyKT?cs42*3ybbs16!jkC!PngZ0Tka`4DuMF9caxNlm7s>7w6cVBW2_r8Nn+xm$osbH0XfveH z2#z!C9Z;Em7+EjjoIW}bO7+g%5;Etuj&-Y37)l(KR09@OiY&fg7<`j$UFCJ{#|EnavwGV6CvlJrS25tcU17MuZpk;`jx0ec;4ndU%JN zyLO;TCGSb$_th}f?Bc0C!FSoa!q@d~;@tvw# z1`IIn0aEhy_rIQy&;+_4${O5mB_uTDgh?1l5592&p4AVxL-_0WJ_tj{Kc$^pxH36CA9d17A+w!Ib? z2U`10O;fJ?dBD}IoScf4@N z!K!`>ig-_47tz8QFZ1uu%p)>v0Znlv2$aWDiQz2N9Oix;t>%n=4{~c803X%qf{YfbIi8PjAVI?Kr zi^AzXLD&$jblHq(bJ_Sd4W;gL-VH(J=IW0dANDWjvK0I6qFn2bZ9%%eVDSqn==$H3AYhS}~YGK)1f&%!GI{Y$pu zw-;VanV?7V0v4xV@hiOrij14#AMgLg2b0fL#!Psso2^pK3Pa#%McB+xUT)#0hW2(H z;9p(6?+5XbV`JU-^)gCS)-Er`c*PI#vMp$_GiRYx(}TfmZUB$4)N=D2+r`H{vrzbBnnhhAU{F?;b3fP zii&{^PD`~;8s6MG$65Q)R}PAcUssi6K3nhiG5#p1#ZO2V5WtQC*_~KTo|!OU2q4=P z{Shn_60op5pRqhT!}G`Pb|@ecA8~yT zP#ub-5z;jUsEbwda17tXtCchTQX0zh!)z4drVsqpJMiRG04V0(LKHEUkgvX$)>MJm z8Y--2=iq_w6sVn1RG}}Acu!IIy%!ku>GHr<$*~e&te@Q+LoO=!${LC zkUQ`bn5MiL_e1&Nt`%57S@O`BXp0bv6aNDkWqG1Sc>kL-F6;oN<}30F3haRa1Li*_;x&_X3E_JA!+ip8G|8@v?}5bOB+`Tn zb_FfmX#c#0#+aPvPC#`yl<7aH4h*Lx2p|zHQNt?QKs=QH5@vcx6YIgq{$Y1x`dw>o zdzu&<7x_%vVY~(7Bf^h-620Oksb*OVw7DH>T!N6iYY7SDU*WYw+hij;RxG>ap0%}w zg@p~GgPLSIfpB8jz-Or&rj{Q<<}FT3#{Piw*kuraSr3z3XjD{oW#wY>zsF2X4PJWh%tF7%{efIC52^6oAg($#K3)yo zXI3}NvYcttIJE2toTV+#o;~gCi~)Yyh}L!|r_mTbe^{STkSY#wrs0h+>2M-LK-F6;!u_?KP+ZM2kg9$zbWgW?9 zL~b_fDnX!E8^oQU&*v3Fj66I&%ZH*B+M`7bxBqV=R>QzRO}=LsUT+r{Nt#Bl1n}16 z7Z>nD!kE>FyM0q>BhdZXqT*5^p*%uwiXvtoG79Wp%=>|%T01+BPH!xw836H55T zPop(01RN5Hz<+_W5=*v)8O!y1W7ozX8MRU?@^hgYQFmzfk0vOZOoV^?fA^O`_e{6F zOmr8KYAwP@!Y^Ts?>v|Gp(6ULk!x81)PFiWU<$X=Rq{lkWwBJM9Sm-PGb8LniCv&^WpB~`t1j1I+amy%NPR7~xZvl|ygrdEKdAV}Zqtv@euTs zz`ui}5xS6|W3UrYf-(PopAjF8^6wZ`M(`j!s;5ukKTx&L&T1Vf;|>1$g$n-X*Dsj^ z<#!r$RIIERuY=J#@?;V->V<@a z$O+l>yL%!CE4?mlWmCAGeE$5o-evcTV30BeH8pk8OB)i)>9X|A$MTZV-E$%yTe63T zTf5Ngt8Sm;*z6#HCOGvxO<6R(!N$(FO%b?Ti&CK-(wpNJ0;3H4Gh#1yy z)_~A=?MGs&vYSKcuj_0Vqmz=5AzSlxY_PDfBNaxk&1Ty!TKtfbhG^Wm=$`pOc=TZ+M8cD)EQmj!P%NQBj&mN4KZG5wp6Xp2f&AWF` zU%h^PYCT@6!_$0!F<luiXlM^vxUU;^20%L)4o^(uYi{l20@6{?2#9uO}W1+!Gsljn0mN_c4 z|2g!vBfICC4UTk>EWIYzYs9DpGBUEGIjiynNsC|uAOC7Pn6fofN!xmNB9kWU(>T(%xe30&CqkHCxw_QS^( z_s+J*-LE$imGb4ulTuQ$wQDSixXiJDlhou8BF2sGRyP{{`m2u znQ(1E)fXv<*M)u?4;*Rexm%D&{xw^WJY&#bylpkxI}I>Rl=aWk4Yie4!{AOhlfYCOqZBH9YbDTervj%l7LNLH`fPiR7IS| z)9CVG5f>&D(*1=+hrr7Z36lbEccL(^YWTy0ZnOJwSRVosQUs)XuEr{YULgZ7Jv|-O zd(W!8Mco|T!NI}!9N>(QZwLl?w@#hy#_bh!9He_?#n9pAY&%=Mc+ia!tQgO3shJ@> z`}xD!xI#7+^0n>sH#t%cV4=$gM&Q)blOHtG9Ns`Ar+{CU^YyET9_SY8{kAvZR?GpW# z#`et;<-%?Ncn$6A?@Ud@;gB$xk7|lE$}`+oJA=~v{J={+E@GPQ37;`C7PVh0)^7S> z{#i0pTCS2Pyqvdd!&hTe3%DKMc3!Kt7(emZ15 zF6AX)){J&yHZwCjez?C4vh`bOhZ+6(vv6d3YSCLw=xQ@>JE356vVvT@I1@U+hC=&Dw&mvO!0{G9?FZcU;55PvZy1IH`cG@2X zuG{1EPng%`0lwSe5{rE5OF1$QMoNZP_V)Hzq+I4y+FDv+lzw-WM7X%L2lI6b zX@VZ#HFG2*(|pSHnqDrGd?w?Ib;8lglg2-;8WzUov0ZphA?U8CTW6bo1pp_EBvKNd zs6}{sDRr6-sRyk8>?QU zs`urho0GnRipsbky=I-=GC>rD5aY$}v;sU8KObM%n<$FGQ~ElypeW!)8xz)+mW-dt z1=Mw0JZs#+HapDM+N^;?W^g@R(klV>2UD752bY`k|gnk5D|IyI$U zU0I3Uc7LJc09vFH&XLYzPVt4hP7O3;mXSrm&v%?DzNbgMmnxLvEC_``?Ya zf{i!&6M6s!hJu*iQg70S0UjZLw~=Vz2pc4sgB?XK@a52Vs@#CBXmTfQIm~z2 zo#I=@8&2lTZ{G*V&h#xeSI|M{;wKseZ1T%xg95YR4175`ISGg{u!hqGJ9eUO0v7F1 ze`kCU=moAR`~aTvE&q^DenB^U-_VuR)YNpfHJT$hk4C~J1y5BpX*oJHBu~a`pOuu5 zaNhk=OG^uCU&x}>b#;Uwteoy`0;9V|%(B)f*B@*M=QQpP^{!*s2`he1Kp|Xa(8ijQ znz{*Jmmm0R+(oE7H}`3!B2Qj&2F>N0yrCk~-8BCi^l2TG;j(<)|OY3IRDEmG82gEx* zm|VcqJeJ@Zx{x3a69r1y9fm)uG66kIv%w(aJsN1td?PI_J+%1*avUY{;Qq%1iruxt zhgEbeEOKWT7elkA075-#e=f3GO_#CgH@UJ5znrhL8*Ui8=_R%M-5|opr{(PI4DhhE zwe{LJI1HfZ<+2}N-(0=JI%g&&_FKR-!{3%qM;I9yA?SZ`17kzbh-p|2S_PMJAxd~X zHCEGqtVMveZI)5HT>*+>959bg&aVnj1+}hWA13vZAG3X{T_LDKODi3L;N zuZunGV%5S}JArX?E~wM=K*QlOZkaV5wlHuG#T8#bgF@2w0iAm6_Gu_5*qg5_Z}&L^ zCeYpW#$h?I7lpFHpu7`y1ZsFIyD0+RssKG?Wn_N2QRd61_s<@H)%%$)4%1$mt!0S? zW;EPKlJYo~ae=tk(!S5C-<;%hUIVNU1cOD^L^cEc72oe|GvMA(PjEe ztv*d7jRMYFRIp=ex-^U=7`*h1&l04Afa@Cr1KpY`B~>k=Usea0)Z6TFI#Oq6;Jh^g zvpJkOTona~-^QJw@7)z)nFJO*8tKcD_FDi&sG=m|*Jq+j?|*%PfdJp+W=lwR8=wAQ_iyL)+?IJK>z~|jUZ%OyEI4@DixCPU_ju>K z(|y3qG-=iD>gX3pP8Qr%R7j>%gxQD|7Nuaz5=nB%FZqmOpA#LxV=0C zDFOAc51=2+M-B~9IurXSDd#(Y4koc1!K#;N9(R9Iyl(KiWZeOzdg~<{4JRiquzWa( zGk`jA4cpqfx;{X{B0!)(1d&L>+L{@Iod09GfZH0dR)~LfHCHnTa#C;Kc5G()=^NKT zlxG;NTM&_ukaz_I)&WZKl`FYM0@#CtRNF1Lj+g)w>fr4m=+r#ySnKzV)ULBd$Hk3< zt#fjA{s?$9h=^z1yEoOnk)(v4&47MV^YJx|Xl-Qlfp~MgQ(O)U0eqs;umctlOQBPT zhRnc-2$`?*AYdB;qlHc+2$c{VPoUsDm-z_QFriq(dq09dm(5iosMDo78ow`3N_FbO ztI4q)iu7o5fTQAgULHgai#(D65NO|@ufrxL4g)WSf!HoJhbahO{|@VWry!6j5ZxSV z#=@9*n|M6FOw`-g*EYgYl~3ZEq3t-WG_ikj8|(xo;J%%gcm%dWnA061YH7K{S?7Tb zBK>f(XHETMhy)!ntai0|09ea92-+Pnv~TPIN2D)n^gbQr36lS5i6)`-ZU|%<{%8E} zWu&BjcnoILMD6t@J~e4UieS3sX6$z=)2`XUjryseB!Y{QD%?~M8yky-g%y1&e7hbc zl2Pxem9#Ltsx>Th>mhPqzYu=YjE>laTW4l!N?_r+nU8fO61W-FVWt%urXi^cbC1lLs;auhSnih~e)aF$&J;oZn)-ml|Hr06;$V zRuUKQ=>xg4bfeqh=M-+MwVlagsONy}FL7Ltp=rJ!xc9r_;iV-CC}05#9VyY`oZWd| zgn&x$+xvBc`|%pA#+mQYidfjRv6&fCf;2*42i5K)HI+{eDFHz+_y(+~{PeelOEXR) zGz3>UkFzbMn9ptiz|CyqeAHOn>e}|00wiYpst(w;boTP`FVgxuy|AM7a{!?)5W6^@;O<2Qb|yN39KWKwvVkf=RwcolHp0)UwT^M7nHdwmbSef!4F!vn%13JwlV ze~%fpl1G=VhIrodOSKOR^gcFo)i4l{yXQ$QhDg5xvf3yG_tVX%dC!?f601MV5&P5MEE5t8gg2tx0>Uu(Cr@>^PA&B%U`3v+v01fQX zv1>tqR8&;7e?UOZnBCvDXvlnhVY!7677+nQpc?;c;>d^dGnk3vJDm0)alOvgQH$hOD_;N z_-fqt=gz}!zdu65Xd2A&_0J5xLaNzekY5$H`vXE8AQ3@;UmvyBic~d^t=L=}f1mR7 z=~F!I!-+x_RBt_*aR9R%05IOYdk29$lS@%uVRj2TJ~{cIU5e=c%B-&us24!ih(0BX z?7brA{Vsw^DdhRF?cwIv)-8R}7L%S!XozsVpUG1>tE2%5%)I(m z4Y&=TipaxpA7zq&n{B{tqu--%7+M+y2s%S=yP!;YV4DN5)L_6yiCR8q_diVC0ie*V zw>QeTJ@*s>>RQHZni3Nm{{f5iO9zuhaMYm<2 zKc{@&SL}dqKtvmNVs3LuqI=KYWmeqZ?3P~#ryRCDJ|t^=a=ZnQdLpzf;&Jk{_4*Gf zWNi2D46i%qP1KMFW%J$1uWQ!`;A{ly<@$PC$3Pi#@_hg@*69LUI1PYkjV&vUg-gSO zfSqf(d_C2#U%&Ecpf7puvVio2d-x?(xB}4w6;hzz!hdsl_}JOUgrS6ojfWTSgzyJ2 z^WLM?E(3D|z&eT84IvPVFNJroW|ESKK!R~BLvDGv*=5u8Snd5p95;@B7^fiI0p^0T z=nIf!rq)I$Y~od}IN2%C&p;NwHJ-09^wvK$l@v&G=uAvZ;B;32a)f^RG?0uq+1ra~ z;B&5)aX*m60Y}d7Eb8d^0u~Mq8c5K{`SAdy57IC1YjE1EFzP~N)~&|^(JrW}s!Ay{ zM7tD9e^w3-KL7sx4k(+39Ul;gURo0YdIOJ(i`(7Tw+`L|6|bNIPu-Fy)9_Cl3u2&5 z3Ev&P&dtxK1rjJw2VQ}!%gDwia=kAK0QEyVNR>x{PMCI#dF|llhD3$X4TM6&y_w32 zy1HRL20-I|V|EIGMgn%65*`78T3DC@%)A|_A0s($FjY&n@qhmO38hJ2*bF{{<%?WC zu>{KkCu9arj%{l!4;5I-NRb*7FE6hso`i0zH!;AHE|AVUwt7$sxLRjeNRg9~f#~56 zR&EUFAyjy}m8WdFN{KBjS%fnMgQ>q6gJz%w6lau&|EteF{&CL)gD|T;? zsq>A?J_5Vui-Wwgn*B}(eZn0u6DWfLxo`+z2$&#e?&}J{gl_m9V6Wz2Grk%XfPJ@j zap?gd&f#$aT<;wfFE7!{moIs}E|?$gFAci8y9H1A@Hd5$7C%EQ(xan`jMFl7efGQ)WhX9=LQ(xM2g2WPN-vJ;f7;lYy1thf#&|-F@ zPIxkYXJr4ghr`jF5T#6omsVC*1Ie67LY`-$o}TDkUe zujVS>JjhCV5Yb8Swy(X9dk8%!GFv)5fY|Z1=G`G}8-8)shJGO*BO~MRJwUl^={&}Y zR_fQmp`<*v%+7hXE&7>Na}K_Kj!f1;r`7-|vy+pxzLx7KsjK+SeSCuJsP6AOMVTN? z_nA^ee`5Fi`T#1g20c&$;p{58T<2gKnpS5wG9*03fLnmXB;V6nGL!WP{6*VpnK5xwGy8B ziHUnu&f_a_0I_lnh`KWkPV}u;Yp4M5;h_E-6ofEg$z4)fI$G~wf=2uT9*7LRz`;=T zel`HHm)GZp>zRC7C}5abfZnllaEJj*8U;A=iIfMUva(XLZ6`)2f`A1Yv1C(UqJxur z0~no=r6o$KdIU(IJ3#dUC`%j{hs#kQ)2Qy=c^*D8qG-L`N)E<`1O~#%%E~I1YI9#* zUmGnodq81cUgYFMfhuE(K^q11tJ<~JPe2-Hj81`pfzf+?x_LbMhVmy|o!?{L`6aeWL4$pIh_;2}gZCqORdM7j`QHb4>X7#hN$prAq2ee#}x4~2TUQvS$Fvu0_ z0feg;So={wQ~B?Sm6a8Uh=@oDz^mm0&1;4sFZ(|O$yWdsM!&vB0{jrK@yBA7H%NSW zK^h|&NfHD^gN}g#Gy(zw^Re6~Aew}zgHw$I3}^;$xcu(;*^8kQXjR2Gnn z;x`VBqJ~5P@5=%@Jll=^-KCjaJiXpA@4BF(_y0HTl!IyO7j0TN-_kR8&a! zX6|4v6w=|bMGrBT-$STD9S^1FDIEeN!vmumW>dW9_xQDsyK5>0!prW0^Tp?ZTU;r= zfe}n8`+LZ0;1A^Hh2;-X%i@ssRbnO0vsE)JbVO(Yg@%Sv2N0T1xpBWVoe#O6`1`M) zni@YeQNHkp?%zh|k=Q%&vVDEYsMZ+&tj@xoQx3B`0rstP$3&1ki>p2Zg+pSkkSXe| zaX&<=0teR~oB`i;Iim01f3rpaT}a9mlnCjvd@E$n>t(OP>dgv=)&O*J&Yq zK8S(#bEJuGHVQ3`Dlrsmbs>BeOU>da{dh}{VKTO2B;|1=`rimLQ0zo2 z5~z{Tv`3R2V4__ecp^(8h!%CQl+eTsJtiEmlBaQ&YC2ScyolNnj8cMrB##&v;>zC5 z|IL#KlwFDW)Gfk8qisgC)L)1o(a^kzqhBd%@kt#Zv>c! zYOki2T2~Sn;c6cK>DGXI)m#hTZZvBJac=vXAr&Q}~dT&I6%;=6kjPY`-5SeFz1ccB-SY zmJqFwAqk%(xN6?*$G=B7K_5Y1t&RE-(*QlcanVppp(jl#$DYx&_}{5OtJq3>#>U3M zBklMFT5W{gpPUE$uHTaxwEbrs;U8B63(`i-oslUTiE2pCLy{80g34RMyNaj!|Gk70 zdMd9iYY4;w3A5zm#QcR(W%?p}I^Hw=XMl>pM>gW6^EyKP`OHVGc@0NQu65Qc-8yI* ze*N`fEa)V@DHmdIiI7L?1j8c2skc3{UJPW=|J%)L=&5ePiv@Muf6KB+9rb^1yx&J@ z_~s1l-vjysR|BPBz|R9(9>pRqOkg0Xt#@b~akTHt!hwBny#*6M0mSO&@`~JMVF?jp zeSctO{a_V4payj^A#fDAR3E6a-0CPtTgTr({P#s3{m#HwYUZM#!f% zICYZ{;y)we+p|OI8~i-{g@z;;fKLwG^>3Jt#IP^Riv?jMnrSK>TSrb1JvG6sSqiUT8q!2_GBm z3&KccSMkxV5bdZh zrsP% zZNpbQ{RTtCYfcScLy*W*{`#yyBl>@_E!Q48Ky|#=hr|7cGBSyZ?tCkPk`Ws%CyU%) z+>q)32Fh5`i|%Qtxq?z4J77=iD)0?4gqD_0;@^O;&;j)~E=|VnDz&LgpJ8H3Y8z+&HQ<;|XzK^i8d=;vewsn98A81L zOwsKZKpS&&^!49}4`77SXA*9HDrR?7E>l-``;+I-&?$oBZ!iDt6dtriBx6Ojl)E#~ z#L)P_Kg3MI7{H&q^JUq;N5DcKL7_&TN(%)ynKw`Mc?12B5NH^L_5aSYAPpQwwl|-; z#}n+1`K&6(S0@{%S3HKuIMB`yC_Xawfprv!VkD3ntNw6J(w-LYulaT5Hlp}cx9wk` zP?ZO6r1W)ff5%vM<}bLN_o~T~@DjPP>HPyp#t`s~tz>I!OC#AW49%IzH`;e{vczPe z0)(vpPAq{F9Ej@Esi{ehXU+o5`=qq@2R5?Af_QbqX)*ui@di4NAIvPp#ujp>V6mf4l@<2wThT-_{I#$CblrxM@uN*gk^RGZMkMGw8vcUY z4}uvrtB{i1kWdm191Tp0|7N-l9TAk2GPky|B3Gchh!4O`mXrgB@bcilXLS7?vF$uX zgkKG#X6X_ve2#637_#EMXZzQuvtIze2@RJht`YAJttr@9KnSi@9l${Cnza%52P^$l zz%7%tbbjfH9VaIx(RyDos+Lca2cK~b5aS2LrPlqoCmdi;UhFIohlWDZ(l^R-ueeB5 zKBDObr=WYj9vEoeVO99wPJ!56{FJ>})j{uJZ=z71N(y;{>KRl)u+NLBN zmFauq)!p*dqayBq82la{8X0wUA(WJ&`EafW@COSplAjbMAxgz^?2bp?6J8d`S;L)U z4t6W;vG4!;4g?x+om>zaoncfN4M%AV*Ch>4O;;X#82l#+Y2AW)@|7*vtyh@2tc7vd zw<_1t>#%tI{xjzRuyN8c#Y`*PPYNd6Pd3_r#*OXg%fF@my(DD#TPn-f$77cZF+n== zD_r(yt{RacyWt?xq0!|4-@li?Fuuw{DYR5;@6L6ODx8$|qY#9^!&Fp2m{};7Gs#6M zgkTpI^H{Yu(y0pA3l86LSqUIwNtE24(e8z&{4>?3xN;6|(Xl7Uo-K1kHPN3~_2xc$ zUCMt|$#T}MUJQ8o8ypBKfb#R4+ zV=wa+;20RjYV*`I-r~aJ-+A4m*Jy`yIkcIKYh==La>Av*3Hqf+0#KAVOtX6*UI{Ov zaUT~47Z=d{QyR>?7bPnBY%A-hg$wWFk;Nc47YN!lvA!>f{7HHIi8=kx&c@2kz?OKr z82-dza4t|)mXoFnr%|`@ro45#M~1j>M;@DBynw=V42_`Y$7@zr_p2l75>&13k{JzQ zDJ|E%D>TSR-Y@#_CfK!&b6$5KXtu`1aSiUDZHsEB9DmLmmD1p*5d7hlNd?J%u6Int zhhsi2i9#ujK>?vchV^v0ryxnd@N)GI!o~HIc$1FM>vPe%kOvl|#33Xt+kHo4t+99H zu-+T_Y8*X*5(q!IpyGzQH*-w<{rlMJ?rlvrtd#oH`bvykUTm-5$6(m;Z*TQ7!ey}h zJ$kOIk&du)81!(k0G+`}}aRpo$2(X^^{)ERa zR18$cxT!J(klT4RX`|1~!WbDT42ho8eE;r`O6Z^aw0&Nu5jIvX&(nLOYZAYrG6pe~ z@c0+3#PN|k>SE(toBd%Czs{}Eb^OUG5cr?71!MNF@7MJiuXVqL_`@MRGu|H4pbGCN zj;?T^t#GkgnxC|6pCUlUOl%ax@7wsLoXDi3!{Gm}TsK3gEax?nMQ`buGu5+acGD%{ zor^xqcXzq9FFj5Uk)L?rBj#cHf<&gMJzQvhvP!aA<@)@vP09Uy(mpgh^qbu4)k6_t z#stT2Y8OY*>Lh+C`fPeNq2aI;!put?x_SWzi+RP=eV^Z&2=49q2czeesSCJ592`yt zYSq5%%@YQx6?ejf67h4Pahg3K)oAnCM1pi=J3IxQNVq^zNvR5pf&yNNE6ZVhm5A#( z33$=7x4SdgTT?8D`<^Xa7qpCwkHp0pbY8kh5H356`UaSbi!qu&|LVg?(ts+n2$5cK z1`LX{u#I1VFhnbdBCw}_B*b!*nXUzM_M zY;?#(1MJyMymW2RCyBzwM_6dOyt^GF7iCn7?Tu9}XVp_{a4o&kP>bLC5p%8C=u9N8 z6ea8vY&e=t%+9c%w-pBXSObT5v=WPE1ztcv;MZxhZ2cA(WVAHY8g0P9CCAIiN1EK0 z=FdF~i&dSu>^DC@aFLJ?Wg2`)=Dd$dc;Bsd{jS!IvcftU577|`wuXeAGtRWb(V6RjBNJaRBg18C z0&ip+Wm)oxzGtBbBBnlUf^-Dvp==+vqZ<3v0Q(e5;(}KL!+|uSui5I&2ecE1?J-8# z%q7g^WJ%40@=v_(SelKhH+*jFiiNpx4+Sd`c4bCeo~a zeu$5iAboR(`Xr{;CjP0I`&|PG7FAHB=^*HJ+pe7|4G4gPp`_eimPV&SC%H17sZdia z9t0wkr}Oq0X%sJiW~fwyu}J_TxgtL&KR2 zKfm^YbU#`!rl&2CWjlb13UsKj-RU|i{zjy+@wnftzsy0MmVo((h<{QNq@zQ0Fon)( zZ(&^&7A}v8?Co3S#EOFFD~D|X$((r-XaX*=-%t_b1zdCWTUd!Wj9IuyWcHR?h%+*T zqT_SZY60DTJRG66c(_gF3AkKFT3gSPX!tV#^MyIam(QRX`4c(QJr1)KzcWNW1HScE z9?VS;%#B)KKfTgo95E?%e&x^LerWiKgtqpFh4;jRKjYdL8piNL7c|_|7)5hnK78<> ztBLem*uz@gW%5~#OG1`SMI;kKS!iWJ_iuB4qS;1K?OP890Jq<9(}eE;-kfOs2~P-%cvbixu88vB?MN#hdm z728_51`nG(#<M@$lN^pE*aTJk9-|@ZSRWMdiCmtk4v>t9a(#qa8@P+TpI(_beg3NDOBmGDsb9y{ zdc$yUSN{A7Y;sE&2T;#g!fln?1_zsKP{u~Xjem3&fAhS~k}eXxQ1xx2-KDJs*-yZe zl44=t{3RoEZQ&Q2;9_DF4y>;1w-2?dPDG{_8k99Y8g-$R^H4#|XO(XnOx_eOyIy1# zspW}gkH6(Md%xjx*Wjc)UDX?z*E)p)UDQ;OnhdCdAk@}=uJ?On|0a9<%JIae;n$0l z2*Lnj&fMYTJ6I{htTSR$=2;o-g=PjH0(`EF*LVReI=RIY{&>8O!yod*CyQ!cnQsv! zrSPu}B)_(IIq%({?}tW7{IS)pvE?3ODGEY*7o>W7F%sJE=GaHGW@A2)BrpPQ()S5G zw(*@^WqP1O;~$C6?kgCbcIR%;75uhM)ogg^QFFG=EF}H|=yH?Ocml!39?Lr26Nw~+ zOJI_L+p$}MS1;|qy<$+wHQ36L#J_iW4tV{K=P5jP!Fo-`hGrV*m=t7v@eKM>(v8#r zmz}O(u=_MZF%(NI8umvtwP@V4UOx2;`{4ZPvmXi<_!wf*!fmjy5I48=2_4eoH0GxP0329=CvAto7;GeBbD`-iGkZ#x-WAV40jz4R zic(V^>PMk=e z=-D`4k0%n{+ux-X^k!maNq4xv;j6EA2R)!NfAfH0rrY0*Sq|Q_PRkM5XQ-?qHJ13z zk5<#UVmpgXGH~tma-L+yuKTgI9-btZVu88E7^i>Uw&_+4Cez+4qoXm;IWdJN8hrU2 zZ8R_9j41Cj zOl4IuxXMvl@x3S{m4~gfrw0xQ{)y1E){|m(d5)a-XziAro8_uC6xl*ItgUTAkBM50 zf0j^!Ar#`_r!+VFwoPh%JSa>UDz(91R{ER`tBYZiN!`a5KREq z>DT+*CPCv}x&Q(amT05)>dp<|ZHsa{zu#_Gxu3kVD0vdi<1Q`P%x$ zno=#LVXLJzR_BedJ^n`9+iF{xp?e^ezLfHXx#DP4rcls(UzP1YUxbdl$DVEzHc6iT zPN7Btk|`auCe-9+C)O3r67$6sG~X=Mdz@B4=iTDfK2bPkrtAyTJ*?1pbJXl;a4+}P zquCDRedk^;+@I139fTuDWmb|5t}{_6_)VD}zOi6tA;-q9L~wOwiHWu+B^U5q9iiT+ z=BARGi*0ROacYZ=qShr`Ehfe467%Jfv_YRrm@WKTTiZvzj-hlsJ}}AE=p^m2sy{y5 z{%(utpVlCWI@v8(zWjbJJu8=&5oA8}3%)E{m7HugFk(ajdtv}~mnvS+TV73$8F7T0MWVQ*n zPwgFHzp%vlavaW3r!&i@GY;yXR__^%a9c_(6EVh9g|oS$#Dc3!q@<*R-d9*tQ&YWV z5qOLY3{Uj+g&bDeG0eHXfJ-sxY3Tw#^}dX1YD?>JF{$;pHH{CYk7<5iW9NwJ_4`dI zn8F>`I>W&kC%pcd(y!UvD1mW|H(t9Ym$mH8XT0*y-&uu(o{WwPIqXhhBqeJZ4W>*U zovHp!t%9yEbPDHZ%`)Ha#>8_<3?eq4h~@QzGTVG@X;4H2&-d@K98RWgy7d(svzW~E zv2g(f&kJ=_eQ{l#_nx~D`e4s4=6*-YAS4j1wT>}`_r5wh#%l*~lwFeje2}WZnLi!> z)$8eOPZzt~Z_mDZxrtcM;nDup^*rdF!Gi94bVZ= z5k{yoJ@1~a0`l=Ma|)r@vZCI;O-fN+y*!6^!vcS@v#jE;5u$*AbD%XF^^0QN$cpLS z#^Tz$kF^}s7vPb~#RpRsqca|mg}g4n1(crVGmz)MV^&g9P7G5o%`MX`Z@F-fZhQq& zKW1=$q|&JCg}c|~hmPpr2B+JVi9$8U#P?i46RQ8DQ+I+%1&n)9GLmDfNsXycrbs@K z$Njgd1w-Vm)Z3Wa!soKG1K4=gWJ5zO(J?XM!ag^K*C&&I;~Kp6B4gy6xAKU+pAn$+ zKU^A=u737{GON`zYjC-p9}yaHl87HfK24JmkeawJ&wMf#>SK8Pkw4&(AfmZ94~<5Q z?dEC0*V(c0x|`RGJteKNf-@V2bvq%G0*2q? z189FrhRVfC5`^R|E()PhNYCWsTi>ek24lu1s2BGD$}?I#Z5V3&jDVHwW6m4GAHc+N z^WnzxWvb*i2t`wP62D=w$wT*3BuUgQ3CMIqlNN&%QxkifqCYN_EaFR;CtMxswYkyu@qdecDM)|S~B%Yt!Oe3q= zjW4L2SB^Df!68cM+~4OB*Sxfbm=5-fjKphL);tNuNDu0b4Cylb(HV;U8k<~5=dUEr zh$den^FA1~Djlunu-ob5L10?E-elX%V*w-tH&+_773*x5n8*>~><`G;{aULrVmZKW z<45hSRC;oQq~AG{O6qIZ1fIQpG=W%&fkC5kl}QNKpr@}&G*eJVR^M6=Hm|QB$p@;H zWk1^ElV8+-k6+ampPLi}?_^_kA}pK}M$>^Lc*xOdDm)?KhAdycx{L-N2^Zi|y3NP~ z$trgcK|RG$HMCH93`5)0rN$Z8y`OG2$!c zW3U^4bUra4YHsxqyHr1X{*rGwuXrbX=SO-#mv_wCTZvjlQRSveRU9*ez@79sn6k_xAM=kty2K8Wb;g>qs^ zS4aZ)Of$4VXEvOfHM&*dby)}14x=+NKC`u6i5`rBJHXXsLX7v{2M!(mCX4zhkhS_j z{2IW4hkB#{^b zuHM1m?gZBMHZgdS0{$pSNw zPZ_p2q+6&~HqM0`onz+apr;CMdc5KID@NLa5Mt|n)Zu=2SH*QY+3M3~*b~k}6&`FR zSab)dG|CeNmEydPWg{d1gvrx=P8M{jw;tYf4HZu2>GC_@TsBbbw?@|It0-0S!nq7v ztftKli<@^BJOn|7p!VW&S;Sye|NGe=?ToOo;b zkrMgC3&*o>a${7^Y~?_oE{1^$xZ)D{OD;PS4UMobUxMT~%8rvWt~a8-5f+_PS_E#) z^XcX=X(^nlC?wvPD2J=CeC3gX^o?c=mfep%LqA#IIi?=GB@FFQ5u<$fr+e0~taHmNhT_ z!ncLVRhHxSM&j$NNlR=BijV7kIY1~i9Lcg!EQ}w`t+Pvn4J|*6t1F%yn+(~5I6nH0 zUS2zkDi%^i?V68?r!2WW`exbaSlq_|0mFt>7Xq5H5sT^RgnXNWf~p6vPy5j*X6I_r z1GGzfPBz4nL_|nHFE#WJlF0vNVEhQfSoi8xu+g7^2E{@yUKh)ke}>2rTZC1BHf1`^ z3T`np*O{>q5fMBqomWzwh54$@y?wa1Vr0V-gb6X9ltvL$E0$gt2FlzN$!yi;&lMwE zR=E&iWB8S|ljh-)JBvw6B1E(#II$7KpRrn7%YH;JIRA<{q0*=TO{@X}V3K(3Kn(#m zO{fV{Wrq1meafkvU8z|5etS%9ZjBIR%TZ-E|!SeJ$7?hcPsC7$i zY~-9Tlv7}@eqe83EOe0#4CVzPFD?e`P6y-wms-OM|#KNk;`e670BzjqHhaFOW zOZZekpJwD0^X_2K{a0rS1iitXb8~YuF4#he8^1J6Rz&^Yz0b>W-Om8>BsM4}La9B_ zQ(JxVGMzFPNYTJ)kB)--8N*8ofEYGXo0^WTJM(W}$Ry^9G`Yzbbw*p^R4y#m3=Du& zBO_z1EiNw426VwR8uawwA9{R2%2z1wy}N#3CmLJbUqmg|VEb5D2#xAWLDfD2E2}Y& zKK$uiK2PzB3TAomCE=2o2k5#ukc``PjW;P_KKiHhrv2JRrtWi9ao@cYZmK%gDC!5Q}07$c5 z&r9p)7b+~LL+SIN@%C1%TtOmVzE!i%_Q`<3J31*41xSg6!y{J%8FS@QjluPNyWD({ z4{>pX9|AQ*i)%pbJI7)UkBC;lo15`i0Fl=M8AS zyA05yDV zeLV==q5{o%C_J{^RE+(Y+uL*%&x2(@QzHEt9wsGIW!y#)SHw~*o~pXcnQqfO0-d)r z^~VlJKPXTgZ)>?;zQ4T)#}L@no zwe@dkK(}3-3;Q8ukef#brEh8c2q3LRh+Z(^RAsdA`5Z%=$1U)wV?b8hu)@i#P-glI zOM8P;)Zwtt-R)%_%ohIJeb9!5R}J&**-Ow(kHqhCBJqKP&=Ay-ltF)#@>3oj!tQSC z4|#dyYwMSWI}>G6az39WIQr{3&HQJpfOcY8!-W8Z_8x&?Z5F)%P+EG`^Jh>25i@q8 zWc$O)irk3hzWv!2MSYW7&;6#BKftj9SkPtCW4YYQQfs%2X0_X257GjlRrcn*af9kI zLzzgXgRzZ@Lgb>L7=2|AZX^*dx9PK?OSSn?&NL=;W%kz{Lq6!P-l z9*35n42s9pil=bCe&t~2b^%#O$m-5cq)uHtRLc$h6acKcT&wCENTX4_yy%}kO>6)C zGeb&Nb!g=ylSZzEYt{FYvf;1FFt39lY%OPNw|-X@SKR7&pFhiC6V&3Hn#$MPukz~Z zHruykqkjaO%LuCUYRI4s^p*9VIxp{pxbYgXZ1P%I2xd9-7ilD>e zg9H={K=HZYxy8mwJ{&US>K3k4`=_#YkcpiEWsamxVBULj-8Nqe=;Q3g{yuI?g2P)e zXl-xrlaN;}XuU*~l5hP8#3uX8GDk$Ng+@Ar=xQztGO5}j6K0 zG{ZdA$rKh1!N@bvka%aJf{TsKQBi5}Nx-eFb8-?6{N8|nO&qg;MhAi3eeS(k+(6Is z6l(7CXW7}!_O}-t&{a0M1r(|T#K!{#U^ru}XhZd+q$CaWGzLm1q^v&Nr`o$+c8})# zvyx%ZpT)@h=FJ;W$w5FTt$n4fEl5M#Vt@LF8u}>U8jzDMr3aW?FSjg}{IT+4n(x)l z;^N%|GxHz_{}`pf16=>Fo|z5(^)+XYB_$;lF1sYa-^A3_UziT=1e|X2{DMA=UIym( zQVabw%+Egu7~uBTjfp&O>Goa-XGwhmTF&M&2MY|17pMfaw-M%4aDdZK8Oe0!T=YqFSdszi(J zuZ}nQrq@iayp-RP1Qt!cO2f>EVk|WKsGeC?U zGX(vGF|>)Dvkg2iILE>LTuNr0TIs)eqBI(qoG|pJr0YHnXh#QmWIKp`yuRLH06^rl zN!We(?e9H-ZKyG|FH{C@ZDRLZx99l941D8d@YP8D(T7A|xSMp^V7N%KRPA_jP?f z-{1H9xj*;gasP2&_vQUs&*%9(&htFZ}&mR!?*Fq{T}|UUbdDl+ulfrulfA-uey7d1nXK+m70*r(PT{dgD!VV$fJB zOLITyKaQ1qwmQqNWC#)z(_;{o+)dr?Opqw@-4LN0boNi8oYyg43hhl<_R^vY(IhXw zb*p%&1ywNQ<-scg9%|7bQLd?AWMW2dtS^tl;Xv|MCZeJz|jOg_pBY1?7z??VUGSKlu%e6s9UBW>^sT5uV* zX;6ymIYDZlIAN*;%1{8U_@g7$n>lni(+gGxUs)wiWuitKff2?Pz59u#$(O^t+m7>7 zY}?`SST2;zz!vyB}(@YYu=lkMK!jRH(p(hQP+i;Bqs)^6~4VXOG}Y+ zYvbJ$KN>I()g)rn3jFb$|9z~bd?Haw$!z*kl+?$p~3f8Fkt1OIw zY!Wcwc%E=GM>;ueoThv5)OEX%3E!QeB^eN<_~xb&+dsYTj%ihA zka2!k`Tf&F1iZJd`CeV&7V9wqKG8o$+t&6tXX=Ca$oB@rx3Pp~p2$;ifFY3wl${NG zW|%QKu%2Q&Hx-2<)lrW7-Z_W%AwfI2{@OKR?+QuKrA?mGy*g%e{E4X|G@_cPrM!Nv z%g#Z(`X)AMJ!it1ZN%KS_Z53kyymJUhfo?C$H6*kad*jawMbAg!>dUcK_#jMbD z*pQi#^_tU5$tv5E8WfNOnoVCiVif~Cc! z(P~*g8p_nj5?6Ix)a8w9?u{kAG)RfnEN>XRxcYtd;FzVb$w{Ttg6xM6C#UYcJs~A3 zisM~n2?MR#V2ugOUmxo^Tz=SZZ?G_vqF2AH^dXfuZ4B35p9-L7@jo~H*iUV+&5VP> z!UDyj>*)~`g zpfaAsiv1)H8oRKzPR*s=PtX@`kn>yA=Ve8qYf@92Q_-7wX8N#gll1rAy#rJ8nxVoJ zCns|HSvlx5J>APOJAHq-#(syAKxZZ3nkm9CU|}aQm_Of_#S?)e|bcBexXl;V=JgH+h-?P zw4I#(y5SceTVp}2ckKTZfa)W%hE25zR+kQDExQ&Qi$+Sh&zL9&7d^T!aWcDdf3k@v z34H%0Y2eIdmgPi?VNcN{6C#zvp?Z!@)7ST}JNZEs@`J~|7Z)KY@qzIobI3Xp175;g zwwNM8>#pKLRcYXOB}(e9zJ;%J^(!L%lamW3O}FCHmR}R8eqcO+Z71(FG8~2Bcq;$V z7p)_s<%8|5?3nfnG&>Po9oK9-J5f~maS!ODk6vFCh*z~*!*jAX-0RLh#=N{rD|4MM z84rb)(VAseWu-nS@Q>bj zvmkKH49Y?q=6faQYpLh)znVwFy2_f>Q*(;bY}S(8=fj7pRs<6Qo-TS^7RD`n6EuSnGhj^7*CoTdp0KbT!|KI)_9buX>HA_qB4IiKB<6ggrDqq zIEnGm{!mOQVQ!N<*X_5qqocy(r{mHI;pslVBbnGbm{zT_;#bLAC7|r1rKt&m2B(haehv&mX`f?J&7u z;%73tk%gHb(d|CyCnaBg_ijtPd^zp%<#8Y-txh`mGnk2`m6_$LWuEUOsIDrWKdZQ? zeM)9D-IG|8Ffs-tT=F4V=K6IR+p5*<$O1ITCDDEs5w%|{y=n0K?BpNlDbV)SW#ODQ z8o;KN6vO;$#J`yL%tSvkdd0=<`A5k$E1Pn*liuiwK?d|K<}91sB*;98@{e(| zp;ud(kM?aqXEqXGnM26Wm|XEi>Iq>{dnp&!JKu-fPk>BOd}S4NL%&4asfx$Ym0jIE z-Y@B${Lap~9DD;;ta2mCTfG)*x($Hyf$@`G%<6NuwX1smV4_}-pIk5rb9f*Yuwarl zVq!uWOAmupjMO_VO=K1o%_;6E+9tTOAZTAgf#-J$^$$UbN|dj{lz(RColmcv&pA4_ zX?c0^Sm$$Ns^PJEZ|~GUpmgeAQ2I6NXgXPM(aXqUM^DMLv^HVwr{}TC@>^EPxmucA zzkX>zp7G;HNnlXg=Vsm&kMx$0P7Q|%D1SxXRGCMfQ4!58m)0Ea-`~y(OGzzncs7Tg zzbiljv*YpOLD#Njdz>7H-h=2H@4?7TR@A>cL;j8ilgRJhy$HT@haE|G{I+WLzB+eg z0>;9aFDWTWF?tlt{`vFuuD)j4!DjLiqdiHCgfTWLZxk9@wVC{(eoCW#dD!dz*3%ocH?v&jVK3<_ zJ4L>serl7boQKSs?flF^A}rt}V!T;}y1-8Y>?Ym)d~B;l!I=%qmIVrj`AXX;@iVgR z0Zum-j@=yuUS#5~7-3#kUphwivM%R5OQ;f{m53~GH-F(qp4q+LeYM4NFbm+e+l8& zSXg6!NvGdR%bfG{i5t5c&CAAG-KlHa3nu~d1LI^4zQ5wSQHw%UdUJCU&C}f3zQFE#K}kimGd6##F|P*4 zr<7*d{VCoAd0-G<4Rg=mV+=c@(Evz1dah5v_VY4Khb&q&vHP;16Vhb`QZHQ+0Jm=1@fDheeCC4(I>= zjKVn2*K}kbb1D*#SJpqNt{9!V`3)1thwa)F$zPWRv(mhL8J5egtX$mCpg3tFM_xA- zKF*E!(Nf-(FNDTQRs4B;ZTF>+3gywUzbh!FcC~^Q4v$}fu9(fF z-228D>h-%i&beFHdL=Qk%|{-oJ`}WZjbmEajidw5Jo`Irj-LQD4$jDEwv)?9Hs|oFD5h?2xfq8C>I(pqg{R(XU3}KcqhnS`0-1Lq0o0P zy|A#-ta~{Bj@vYsLaXkHXi*ke76rY$=-S#NtPH8Yp;0itoA>FSTWlKpyS`>7Tja^z zpt-bLJKx8yzv1W_3YSyey+9Ey&4+&VhsU+G*FQYC<$+yUpWj_wu@x5PW2Uh!R_p>@*Y_=5hZ()@y+@Rk^K52`|k4yNDB z9W*t&>D74Y5HF5qdWeH@U_@OAO~u-~x32Cu6|9+FcMS{w0r#e0+waAHIp|hV^SU_( zQL!NEkMIomnr3bIRCM9PmjfDDai;1)jB~(p<*2MH`d!@YRKu)pX7+}JqiwFgqBMfI zHJ*I7|NTv;zhw($U;rDzgkkIP8%tF`>5WnXIK>{PY?vrTG?YxDq4jkb`1}*Y*HAi? zl5@~ib98!2J%3|1TtH?gpJeis#J%$u=ESzelAH1G<9qwCIjcIBliY}Hp|5uZgb5xt z-tMZ3^^V`3svs{fhr%w813i8iSfaulH)w+t>Fxd7#FWz966_%}1HqE<+p6oSoD*UP*qY zvgM$q*XTi;n^-BZ6Uv~_*7j{l-1 zu?!LiD>N{`Xm5X2E5zraLxtO!QGdJ+NBLtM!eN))8(v9#dc14k%Of%?cB8EPiE55( zE31^$xwh9?JO`fnXzZ@Cq5h81BUMe!?OG|>6=OZ;x~R-favKbdDF#wEdAugUD|us0 zO-1WhVLU%VKsl$lDTjHuUHOv4fj8glADtmTbMoZzh)l$`lrM9+lQewzs5>=pV6q%sx@Gr-){2z`ZBJg8%)0FLc|FRa_UOytN#hN>0l~ggX zeO&U$DzdG+=j!40BI4bF!A{uBln|%JWf#-}zU&~{3Syb{m)Q`ku+@yYeIlEj!qL~wrYi4E@c}0${HUg-KUQYdO zkauocsjW@Lb6%0$Q(&>Wvqp&v4b6CCA{iYuafBWbWGaydJ)r;m*-cFTGfS?U0QE>T zsD6%uaX^9%3sf9f2*GRAj=W~O3vXiGT0k=&s8A4y(AX92y(|V^CdK`}kR> zhRF?RMk-<~_^~=SHyd7acvUPmK$&v?mk(>_rZN$Dm%q$m@SVDGQE>a74$IX51EKfs z6c$3KS>fTqEVjzZyr}+TLjO0SyCi9oa_jiCC0)B4p?E?s5d?r80$A$_Qq=<%(WF5E$I@f zhIKDdDrOV z@5q`vxrkhta}_6FgxkQt0GgD{!^1WYo7GB(T2wuU+jPn1`Fwm6=gOCBww>_N@tqYT zY88V`cnh4d;^Zs1^!D~rk)G-~^>mC8?>BcS@|!p7?f=C5C9R#7k1x5uJzv4> zg)K4*Nm+!_%8jJB=hnW`Yp{u(K})H(ZX_y3JPs4x64!Jc!Z^ilK5%QlSXs3@C`M5*ke=3aqF}OFt7`TzCOi6 zEmf(~Vju9uqGoq<)19Ejg%FL^Ss9KW_VX)k1uaEC0KG&yhFXnBx+pj}+F{lUx-pT| zBu|9gYX8Bz{ z*5Q4@a*D;`rJr_P+wcm3V!t_Hm>vcf)~g&4ZKJ?`+uPzbwHD4`Jx z3@4fjs$(=IIXgOpN?WatpLjbzzM|tyk1DpeC?B>!wDiu|{!Q}}4N!&XC|s3s^5@bS zWBv7|Wj&j`1~+79gZQx=@=zqjy}lfG^`RGoVBOQ5OY{_eMLEvun}y6qTPv-(B%<@@ zEN1@ha&;a5(a8zU5Gw;Nk@y+28QL#>K{lA_yRQ$>Sr)QXo~ZN4R<55K2Fs6Z^_>i( z69(+@+kAdzWB(-iv`?4?J_SA8zaVG)QD7RzrM<6S7lmbQ%%mSC5u< z`q9a6!PywPo2>R&K8O{)Z+p%1!0gWVp_ae#TaGeK%#y))`s4NXatv(N6b-qOFs?tJ zp4D+UIs#5|#%?8__5d_j{2@W8^{Bpk{_pzu9SuSU4?2^`{N>dEg;lVCS2C7+uxV4k z{Z94APpEN#L@&ux(8r2EBfR+CyVY?rsihMG3rN#IrhBZ8!KU~B?_d-4Mw!Wt?o#!g zv4@h)r-6^2lfQfs_E`A%qNvB52JKs;lg~Xy8Zov}qekEAL!OxDeuUk*bP-bVm$*bh9L z6<+@nGvv|BO#j2+3VAc6)frk90VN~5ipWhF+1@y@L*z)cC&rz6!R+Iyq)PPksVP)y zMDpJFqGL|pOg8~CC!d`C(^8dkEDPQc0;>yM?;nKXHYlP58&rtBPQ)4^|9sz)Y}gL z0ljNXAuC&!oPJeL2vnIKc1^NuUe@XbPg{J2UgOw+>&)aa_5uM`^km02pq7i@1)pd63kwzK$7xKV`V zBbxt_K*tBdR5 z(!z5J7$v&mP(i98VdObim;Apx^h5vf(D8ScW(cXQLqnFLB3+QcpD2B(Mc$7NT6WrD zUEjYWovcAiWz*+|5mBXYYnzdgk`si5g;87EMfYcwYjpe>o&@D{4kL%ihdheCFRGIX z-7eHT>k}_aAx*ITwpGZw@D26I`uK5(L$HX`^m!go1r1OI0H+>lomG)tY|L+EMFy+b zT)%l~I`wK{C~-xc%x{u4{evY^^%|ciU{;%B=g#&FZz&ROH1JOKqai<)Nt%m!XT+qW z__iGNsQ~4ItenRU2gSB}t$1!et$Y0V$~UK{%bxBf{|D4^NYu4IXlZHxr#BYL{c9&? zSAV?MzN)X*m3qO{rj;cnQiOv77KmouJ=IK~O}Jov1G!*5A^m#xoFt^*L-)gWoSFqf zC!AeoM|=Co2M>OF_zYhT4>y6JnipR6poa)_EUK#y}CKOb# z)#s+>sntOs2U^B`P=b6U5+!iCyD#}Keq3F87a8UqlDsGBxS#ajpq6MCR z{)*s$lusg|Hcr1GykXg%+i7^h0>aHAv9v`TmOwjDP^}96B2~Cx)qVeR&(Rk z>|={_*f}Poj>xH@Bry->hDa)K{=h*`SaI}-`_AI#+Ghr`>hTlaS3UJ{V-~!!$!G)X^8gDOd0bKCA zxRqfAvn~}bNT8#~S-<7dCDyW!w@y7YcE;S03w&l+eJ?yE%RBq{HrB<(QzA@jqz}ro zEX=WMIgti|1S9^;h!cEqra!%#E?K$;W*z0Rh5j0_H;GUWdTSU(gl$FWO6XF|HmQ$v z1V!Jg#o1OxwuOl{4b1{JgYmJgN*GR9z2)dE!Uj}c^aC(8EfmapkU4gD$mC#LXQk7 zorji+sPDdu7u$(3O}`g6hMQ>$Q4a;(-eRf7$e>@pPObw)UQ?)Ttxd?@91*q=Rzb zyV``!1ZynJ>HUi}sD-I)lFTcUScsZ(3SMp7&MPR2604iH&XT?VvYZFnMl5 zFm-T5#0to6(<&$^z_h}6c;%IF49JI8>O=K346q_2_q3*F>2d-`-n|>7ZkqnDtEq6L z_Yw)V4&*WUm3{UM8(#@0 zDm|5~UMjgpQ7nYeLo(rcRvP+zeZ`WQpfwW@M;q zYtu=Zr-i?JhhWpj2a)2~WP8yb*_>J2#pq3R4~Fl$>Zcn&11(%di|%p3HpRxU0ohtUfQ?= z2%dt`lxCb3UrBm9p}X}V4+*OpGoc z|FJ74`{gm{VUPa#qh)Ei27oOR9x?r1Z#0FVbxf2biKHt8i?P#609_5p^N>9xsR7-? zuUX;@lfSm1K)6?1*SC8L)(Y!;vRG5Cy?$Y-5ReDh`Q)Ab|`$i-EP$6OR zLpfdqnyQwmIr>qJRKaf@m@P!#ycvvz;Lz}N?|SRa(PJd2#LE3IR5IYf|AjY30pW-h zy7f9vP8;qfBy2wULm%!!D$P4#!W9@41kAD*`nQkZL;y0l>vGdQCi2u5rvXh-&b9f= zWedNs@$pS_2(;4!&{09_e_)k4Sj99zcUsQlySn3%BM+fn37!(&)-^ROzpG)Aq@khV z;p3wV*Fsa8C2;~*lGxTA2B%|;)nUr6VFGSm#> za;Vjqf9PmX@c%z_w2ai)r%6d0(Lex;iDvKZttHn)w88{bUM4pl$3eonV$D{6jQD(d zsI?gp_0iCsq;=E}B#&mIOZyv!v@{zK9Ks>Q;juXHM1t1mffpPrU!tmeTEo&(9&+v( zcC@XQ&hcBFSC*EQ0qoq~?i%LrPf0R)uln0KrIA<3Ya{kNYfyK%DcshRoh4# z4n&*rogOhpxuD+`#IORas=c^Am;s=sEX7u4{*#>PZ}?xhZq&M$*8WfOiptZPPEKoP zr+x<=sWt=IPjqj*91S#m;3RrgO+oYF^gaeK-vAB7syFPDTez}#7NOxn#l=a zUv2`*Kv0VZ&V}i>c}M+vlW~*BJw!?G}Up8 zkcibMtIPO5^5aa0QAiu>spzQ?@}r}iIeWIGqC&6t`D04=;NVqEGB2YL3o=ZyFjD)7 zq2yIntn@~KnhM758EgG%Q$(HRjYWJ&ZCx)#^^n;qco; ziB$W@HMCqu-gb`84u1Z)W9{~W0N2(H)FKz^ZtZgluiRy`vrl)Nr9JFdlViq#UG?CB zu(s}Q8&N3nF2{L(D@D?}uQs6`W_lD%x-3P&RAjUNH5ckT*{H!G2Mwn)gQ-z86sw$q zpqe{xqhCuJqu%%FQ3bolsW~Rf=zsI-p1;&fZw6KCVLWu%Da1$ZSZ7h;M^mnd{Cs7q zkppCOb@uP)U}ZqX(MX5X=uLu`2l_8xRa6=xba6QaDT;&3*HON=+3Oy~Fi7}C>`OqB z(r7N!RXDx#sMVYwA{frnG%Ae z_RnI4&hvq=&=92HfQ+-q9oD+H+V~23wsiOIAt@QFkcEmyt0hhe=r)JIWl@uZZf0%{ ztjL>$PWCOw^-dNB=l}lUgcSLU7t%P92?-`N=`I%N^yu&0w-n0M&|g79)9XTh;dth3 zVTg zyZjNKdXjeD>%y>ms^YX9pcn822EI3mr@+W#8P+)>JQosJAO$V_7X{Pf|iQ0jE3cxu~KD#r`TSySuv| z;LrplfWrbxE4xh)vNyYYl}l8&PxW!zH@%=Z9`XFXkq>s5 zR6Zb6;Bq(}gx*tSlyJH>UU+UANk_0lMIV25`AC$tBwiU)jAQ3{Q z+heR-C;kpJJyrPOlNp`=LLOB7iqS5q2c%78bln0;EqL zJQ4;mbbt(7C#PT-*+9kwi1IKJaLGREMJmf4iNTzPiG{_g?M?Qo=;-!qw`U>445DH4 zaOn5%5$M#Hlmscmks()JP5&=q5y|mj!oY^c9x|xnjt)Sbs&)UJ?Wl-n*uY(N_+=75 zNv(}Hnm4nBUTL8-Fi?g`f(WT@8xK(Koy>nh2t!nQK5s25Bxl)E& z{{LCEf=o1Vb7y-@{^IlM{~k~s_tW_V4r?qkM#SNJyefNTWu}D}E?uJY8v4qKgc@n8 zq&QhilGkvEiO9J9o5r&Ilg6Gx7=Z+sDyq?B8VP0;G?+$GgdSMyH-Mq z@Er~BjmM6pU+qRw8=~>BMuPqDKC+HT^7`DC@Hg}FT#>xe2&%pbCjvRs!v_ftdzqbz zbwWr2G;WDnX&VP;CRa|*UXSVFZ6tYpcjdb4iT`u#}Vl${gaoA?gFBx^ARwbrx@JOdC~cYtR=2B&h} zNYuUmAK0P#A6{jHm+%n&o|ToA>Zq!4tX@r*YGx%`A$kSz5hlANW`SBwOcroH;v_#p zXc|WJX0R|IlTZwY#}6z2 zg#+M+vFS$gFNf>FBf+5YXP6G==MN{TGrXYT(_4ygL*5(F#!Y0d{SWNZ9CNZqKGsM< z7-|bY|2Hv-0YMd&H@0;J(< zWRDPaUF}Pt$J+JPt~l0MC2;H3^7x&FQfbyL5KDio?&K89!JR{6n)0Tn{t?}obEnHA z1QBI-C8jK= zzfdsSOU3>T2SpUmY||YY6S2J^aWHXX;=wE@04&?ljbh*g1heu5 z!l0=b61;eK=q!OwD@h9#TbQF@Budo9#btD~OUuH_YHdE~IoFL~nF?GC!B$JK%W${L z*KD64Pxs*CN0X0t2co3Rn}8x#8{CxENk7zn|8E{9{8zxGhmrEirDf$miKcF9Dwbt` z=UX)jB5yvin_y}4f&DAkGc*)vnLgEicBdw-nz!O_=h7l0+DGSc|GLS)osnq}XEh{B z#DfNo-Qj#&o@z(xj7uKh{=JzNs;!`pp_PA_S?K;w3@1a?L+!50YfZAAAjkp$W&9<% zf$`QYXIx$OrnE#s>H=s0J4zm*xVe6%$~`i)2j_{O=kCJ+h4}a-R>-D<&A`CM9|}Wb zO=N=ORZ{T01gSqdoJ6pVhv9=v2B6V{MpA5%QB0Cz0>y+QbVItzcoHuz1Q1CFb$>Ms z{Ju@S0OMM89;Rt1#5xI;rQIrL+JGr2pOAN-A8~Jb!G{k+%v=+_ZAPj11Z&$JlQlRtZO^nfj@T;fKY_fc%j(S=xVqteS_38*7$ca=b z-xE%Elz{uX>6Ga@#(n;KwY4w*OF9l_-CUBWvjrJ}NB`ulV+~G^+nz+1CYh`O%wEO9 zbNl;`4}45)f|y)ny>vOEIDe$*&L7vaB!ffpVx0OL=7*zGQ;@0g_Jq^h-PYejL;FEZ z6&{L$pQ!`9*kIRWQ559~v%{z;b{q!u>dYGOZkBkp z4~a3@VGmKqfGU+doklTEO|(>Y@~dndS(b=*isOU8+K7m9@%!hoQ@n_s@XIEhCs&D? zN&hWEm#Mi55U)s~HG9DX-eL5=nKiPi%T6_=B4G^J><JskTT>4Ic*{Uc@*G8b;`h^TEf93qDDD>gtsF`SE?dA0hG*^F0FR0FOxU zhug_t!3sP!%5iu*jIe|!0}ZIUFNbM28Kg&Iw2)}nd5qAB(`A*sf6-B9NUa6V?cMZ! zYKbR9$#W#><2gt-2lvwObRh_Wy+jGfyxIJROvGy4C&p!+OKMP;KqCn5Q2GZ(XJ=-9 z+TG9M73i zJ7ROwrwHm6F#;nNSFqpecuQTJ9S|qHZnFGh-MrSx9{MXRnkaf8 z*(5YyJ$~gabTs$)7`t?9AsmFZ8&3xM4u!s(Yo1fO@voKgE_L}3&`t+nTrCo=i zzWG0Ovb@SuHk9J0dzUbWq*jIA>y?db(6Q2o66y*{Mot z&nsJUcbpV5O&o8SP#J{Y;>CM&SweU6rA()IqgILU`_a;2EBE8vIpRHJmI9*Oj{eD6 zm1#H=4KuWW-*tb^uLTD2xE*%>tFx{#+?F~P{{m&KkP~`rE{$M zBM5E|hV(Mqo*3;o3T=0fuKkFi9F+ZUFB-NLU6Aiv41BMAgBVxdRX)uIaf6W9ST0f% z#kTB$HzWPd_fQL?j2BwcP@J8eG3pzHYz6dv=xd)6{DVCCyiiQ%Z!iog{`lAXK5mQ^ zig-GdFbiK`4IHj0sl#Ya8(|1kqmss7T+fHH)yZ5Q_OYSbh2;tRh=nR$Yp!1)W^26$ zr|9Qg8(8WGw`}&i2oaHkZ0hP|!5BH}HCngcCJ-e(m?J*Wx!U`(>CBHx&aKD4Frms$ z9Ig-|OzkUt@5$)5i2qV>?~vZg8f%IaA_Uod7?%{xhh1RyZjWGyb$BvPSC@4AwadpR zn}zVUx-j&p<1%FaX|SQIZyHV-syu}|X;1Y{mAyW7LWDkQC;XI!ug9coKB|YMbQh9p z=nZs(JAfCWg<`I)94}-L>NRr>ckf;{I#czo0j~%dCGl5;VZEp-e7`r;xnkE*{0-t{ z><}My&PpDEQM)UKXfs_}Xj<;PzrV1hMtu9|P ztQnVZRF-yj3Iz(=4t#z>e1>a3+COF`_{eefNTC878s7>ObSgXX7$-P%aD$~soSxdh zXFY6PkT{@?h`?Y4C)_H}_U7QxiTxE$XSiR{ z!}#{@+|cjeXfe@H&=@nU{T(M`Q;H%J2rm$Pofe*|Ug;fh9{K#b?6rl!z(2!R{h zJy#+}Ky2YAovWS#AMtKI#(|$Dmm}8m)?A6>QO%hcUplrp{nq?JgY($86OTM*M%BQ6 zgNFL@z~YfN8`HzX_2M@_TAOBb^|EkMecl+w`pqkW$4k{MsFI=cEID3j5&?o?qR!=4Zz$@3* z?RB1*cMk%peTa(LF(@EIas$k&bug1Hg-02L1J6d9nr(y2m{u?=D$##3#+JV2Py25i z)L!Zt>53G6kaSi-s(z1+i^3-FShE^w1Ig^GHxskX`dqwglYRW2eao`?{luqEV!GM2 zxhqZd;YP_YCT1qKW0tpr6NAkQghhofQGCifk264E(AhQq$6) zdid#Tn4Xal_qJ`~V`F0uwRc-w6VFh-G&ir};wnd_FgY#Fso0-p;&*FYcenjb8JpcO z1zumluj+dPsVhQ>9tuB9Y-~@TKM(8gKNgAexMgtb)~y^oJbGw3RzwLsL$w>_Jm{M$ zpTvob%HLfjjiT4|jEsMN#29XO)m6MWf982gN=@YY<`4GPJrCN#!kA%Ry-`-SqOr+=k^`&@t;xGcg#Y1t3j!~^-Vjwdw|J>qez~|n^rkhB5s;;a$jNyO}#Q{-h-Y|H3 zKY?3m!06~fuVH>ipFPeZVPRpwYMeg#32-Bf!EV2Q|9(s64Nu3blatLrZ0hyYHACIH2qJBWYxHlwQkQKUMPE?n1?PcVwBW{164W z?iqM?|9x3m5GoH)_>C93+b_E*t|}xnRI#QfN)7LOrJ&#;Jpa2yodeD;U*6TPHT9l! z78MwhlB4Q7OvZ=eCNRf1M^EO2ViB1dSYoqLku*;uA{zU zWM<~%BmJ5|Qb<&^+XPmU9My5+Jp%(;M)i{2-KEElC2dx|nP*|zDDk4{+{R7qSM&1~ zXc<(hiY}Z$S+1h-g1a@~_BPn;^Qr#v2q zF3LWIv*QaGEM6lkdzejV-}>!SBEJSxK1QyWF+O;fKgFBw$dS3R(NVi%f4z;;)_b+I z=xJyY&Q7(o+*Xq`e+--lh}HA^cMUdj-rR8h4)0$#@UkV7&ZC8I(00005t%+#TTg#I ztG&J8gS52ILl)hRTwHpYmLaF6hBQ=ETamxM?7ZJoW3b8_17`|EUc*4|QE z;+B+*>>g;jx40l=KhohiG$k1r_-1SCl1L0sz;Wl3qlJ zF-{e)o$yJtb##nOP0Nb|Xki2Bo?#*FRPjc{28))h4{PA0m~NYn&lH2EmiT_-m$9Lt zo*!?;=;r$*)I@HmdDq@Pn3(?duhWiTqq@BclNiF}o0*x3zcDmN(jkh4GBK4l*j^uiMzDnvde z+}SrbWrIyygCw@*-MjEIw{-dk264ln9I?2&yT`=X43#EC`4!ME`;_URhfq=hm)Yd> zba8q4>wr*s+aIh#M{Ferhq|llCS)dClxw$PrK~nGGWycg^f~>o+DTXvbQI>dWV%^= zbKI}3ZCL;|#|E&@xmKQinsKsKDXk|w{aw(4z=b+t~V(4P&|%hcMdnpflPZ7zymbe z?%*uAB{25<1xAN~RvQY>oklyq4BtWkJBTs6bQJLAw_mIueCuAThs}?y=!HdQAVN4kgdX~|`2A;j zdGW{(5wbs@@BDgVzqK_^L)OX3NfNE^-sKMr4BS(*&V1kLofv`kqGh0qiHobjLVA*D zFT}zf{&T2B3p?Y4$^MM6^#+3PyONHy<%Lv!v_>`0ir680dwY}RYi;e~^D)KP+`Iyq z(EpEn+uIYUOrv-lLPN{2-_mloh6YVUL`0i+v(3jl3Q@HQvTAV>a&p)D?wyN7HgF1t z=Sb5GmKGOosh-<@@#4ij_je}1SDJidyy&s%>An(1p&#gY7yPh3@N|zI-p7~5MsjAe z3GLj2%nvt1iz+Y*LvU%0lVc-a6XituRbj-y-RQqzdS!Lts6PdS6<#}pB)KS5v_(s?rnj+X{E#MyQHnA z6ocz47pB#ApFF98Rd?pxIfty4nE^w;Dm=!O$Vgn}-b3^Brc)t)h(>TXjPM zM2qmF_s%ku$<4r*H#sSZLaG)JrnrhZtF_r^CbID>s%vR!H8|jwRHnWMEwvLaM8Daoq zWkJEBdt(X!z;&q8W3{hXvEtqP_lM-qC$CRQOJm48_9P8oZ`C`k_T5y4CPd6`Swzw4ky&w&g(`uqZ!LbhXpcoJkfR64VuQZ4ZTjL-9 z8XP=`Q7{TE18d%MwLv_tjHDn1`>puz+h=KJCMMeG^E3!N*?VVu@|jo03mSTW6S4%zWQZvZWs5u)Q9a8(sfs1T2c*+B)sk zFAfe4)c1V;kYj(Zpr8(YczQ?*`$C;83M7KA04^ko;$mV!Fv7k+G~~i9d*~w|);!A9}77~OkFZ%lR#O``sf|&HLa`W=)yAAtl zYU+Oe)@%O2{>s@s|O)Udg*un^hQ3Rsk@;iteZ z^cm@BZ}UuV#qJ&(AMXa90eF8kFvjvnKA5f0OaB4-PN&;E8T}84H_gq?F0=Nbtzc<( z-|I8| z=I-6|^?wdu{JS2pfmKi?9?(uf)h*M4SHX>cad9zj_03*0OH1bBD8KohtAZm9&)CKE zBE7zOC)i=K21nbnGC5-Qj#0}g<= z+DB}cyG7!<7ccGv(qcu_g{$~2G!(l>M!5Py%KBUKYm`g>wstN8oS$>X-5|ssdvoX2#RcDIhw@*4b;)_! z4Pn}}?)u<^;nk@zt<2{d`iTCvOmVZgU(spsK%&Fpd z7Hb)%qB>i0err({M?=_W8q*JXew#$(K;=PX!v+RLwa#voT7AAvzO&x3V+U%Jlu;D; z{@@W`RW-hL^D_Uo!bF0+r5$R!c0@A@#)p*|8QOO{Aeem|vT(g_r>MVgvH5BDSPA2O z#q2Ybyw|Tek>nyxr$D+?0Kd7f@M%2w>nX`+3g*YFB?T2IvB(XbV^vgDckPQjT6IHw zbGYvH>w-rsBF|(wWJ?5h2}?K&D*N2^)0+7pZg%MGPx?6PR{P9Rq=e(24P9mdiCnn)BJ;| z(`^7{{9aha3QcO9-3IQ(P3bF+@5B2jj49?JS>7p)Ll&vG%jzrp?yK$9(2#a>8`jp7 z#c1KPy&vo|EoK_-FTPvrUOu*cHi(P?+x4BndTN_3F zT?y{1m)UpDsN#zs22W=pESxn$V(I9bHP5d^T!=3VIe>7ye*OB(A7Vxnl$)k~f5w{M zzUkPW=|JL-r0}L2WnIFXH%FtI=8$QDt@Hgdk9C&)?N@uKZ>1OcTgSN(1;}5IP%1zC z-RA7<{%ggrBVU^Diz&_0L<;V$Xw6qe?d50D-1}PsgN%ivr7Wa+iS>hnSAWl*cU>Ow z0tQ88Zn-hnt`#F!U)(7xb#rvoBg^5+2NQt)$2Gov2}?f`PA}5=*GwrlrCqS0|Dk;R zL^X`{^gTVfi;5P(jvmXpAF)ljiR+H?X_?H-OnfY({rlCSlEtjg0-SOsD(YEA#^{LD z=S`4|zJ-i|Tb(S+Z?Eu?W2Zasvllm|`Z-PhGB;88zkOR?R7~8X`<>IbTz1E@;7Mdn zs12CA;x9o;=@}m0XJ;pZ(!tl?K3NX!pMXmb=6Q&K#Z-(63{rA=*Kb>qjua;MMYSPh=w0NU(2%v@{~h=Rk5iSFY%K64yI6@GRQz{^NPR>YJ)? zD($(TdC1h8Ic4(DZiBb$gdecHetm0xrsDK%fx$Pbv;B{x4;S$R?E=1n8T9Dk!;&Ej zg&W8^KwVQjKO+Zf84ALq#kB1BK$g5D2)(vbSvs!>*BZovp$>sv_T4W0k!$-_jaOCm zSfoC@5qn`Tvf|2`8feLhAe?T;*@Bnai0%mK_t-l)Aoa#FTk$Y0Es_-K1eAl(!Na`3 zo_IN2{_`@jlPg!OS2@4I({l`85y)F}b|ng2H7y&PgqF;2ufBb=OiWB9;2iW+tddo0)*Sbdc790K)sNsU z4)48gO7B96BBr#IuDeI5!p%MZ{I*PgCMlaq{-Lk-|8n zDLsRO2A-aq>+9=LA-FA&viEs%vIxov+XZ%ZDuXDfmtt1IJc>ZtC|A z&u%)z>@Q!w#N$FJgKiKN4ZMDRrIJ@af&FD6xgD2T-fjEmdAPW@-$ISPYha+PBq-%d zM1rQ4R_oay-&rpaPELoU=Sq8!CH?3-TA%BFILqPn$52kO%G%ne$gXbOxB+!s4R7yk zytPZn8|<pIa`~HSdXyjm7_fzYKW8 ztL`YcJ5^rG$?^zLq?Rq)MllW7(kU)6=Q;FMH}R3Uk-53E&x;eFV(k+Z(@IWNEpg3V zAtLK8hRCF#eA79YO<=d1n+!p_GR#t94FH+2`3lfr&&LyWv*?ljz`|;EDF=Cig zg?#%1jUfC8aMvxYO#5t0WV~|9W+~(C-98O_6{6$f^gKLTZSvIbYLVanF`^ssCIsUz z^Ycn4#~(*jR*rr7a(6T0l(MfEOMQJ$Qo35%_mZ-ff*nDC`#8;c!NUCU^QR|%3?--% zkcW}?ylBa^y!`pIF*P~gpSSa_fnvVGO6E%$aK@836!LQQDH-upWm ze3_r}$?R!LfBqc2$elm>`cS75LL8HNLgLssd+-Vt5%Hs90Rc7=n-)k|L3a_~WW?0u zhK_aHg{|Ao#$)b>+ciA(1}hmSZ}~42kI|MVG$a#K$To0-le+1x?d`4~A7b`aubv)k zj+&o8?>jZ*jOA5MmY24SL*=y+myWxTa&MQAYQ5aJ1|KrBB~!)|6?L4mrEcln_s+eK zzV#Y31~t5}Gd#hh6PfP(HVlmITLLLMD8wrMQ@&ft&%Xy)7yYv>*p4U3b`l`9iA7a)6n57xJ3Gm%JcjSx((OLR&wZ%0(82zziEKPhO1-#EOfDFr zMu^>bSJAU{$YN6|RLnvnBcCFaVbVb0PK_>B7*cgtf9eNOw; zs)K}-@-ZK2W@fYKYjX4(`eLh?CTuEOg%A%56aJ;ze%Na~)N$PwC(qD7d>@-b2SRI(}v& zo1SXS236hNyG0SppZ@&+n!EOYDDyo&uASYKbxqjjxMg=K*HKe0C2HKJZI#fhkn2&3 z*=(a>Y|3SG$)#M=p6uDt#W3iQu_jDylLm!w3PnjSnukL(w0t+3YpiZdSP%+PpX<@soyu>0`@ILCKESKARtB_$g@wV% zlZ3ER@MY^i8b8eL(_GJ>9XMbBA{$fx+?KDRReJ)!sa$u&lWgeVkQ4twahIp3d(cyF zddZ{7s`0~vYg$o(nk#MgtjE4Opmz5A;_GG88fceFM#9e|`SA(2}W zd+y{(qJ{8$e7{aQOyQt2h?XK83$C1A;b4PajFdsZJb(VfNtiUUOG@Z)eiL^(S9+)S zjKS)`V(HX)F&^!YVs>(#wB1K@ilwC` z;FBqop!xY@pD$l7$#|xYsc%k%&p0d=_)lwxqOQ z8XME{^P_aB*@F{3V)5HEH8mU=|AB_RDQOF|n=J}#Hk-f_YHQtZJUIea^_Hz$f9m5G zcEd#>0U>Rnl(>w{M$@7CI>uk=T+Z)Ucf4H)&s1tSkZ_o`_8laT8%>E+NzCHat8?Pp z^BYKT&*xzf-jYAxyKr)0x}A0WI7|@Da>}*|A1g$3eqr_b5Ox|vBclK<#Z24jznFMk z#az113GZmRkGEqbV02J$d9WzuI52zAljiqtXCK}jw?Dl7*5;s)Mb(5JzxNc(I9{fq zOc3+g#l`M$Cj(gK=87_~aq{r6JD|1!mAV?k`nfsz`DMs|!P3yd%fa2#vAG)3RrXo$ zjm9cKaIdT0O^^{K z8Am-p4G=~}nS3g9Eowa9x-AS>B9{NU8fsp~>(^&bCfO(S`=Nby8?=*=%We3o5j~5l zMU&r?2IAkx4aCi@2o;^Z@)SHzwvk)D#Q{YUCe{XyKPX?zt>0u*lzuIuN7g*FO&3d! z!l5zf!YITg2`bCvt}a8|`x7m@b~rj}4YG6C&VOhpHH~#I-9mLL`xuQSPs?llm?4bo z$aRrmFD-Ex5+MgmLaeCWXZvcx)nL3JHOvp4MsDr*Q1g_^UDl;DIk{@AsittP?Q(8; z^k?H0YD+IHfE=lB_hTl!JqefmYe0l5lc76%pm?-kub0CEwj}7Jx5eR&&Wgo)*rf)Z zMAFn9lZvT(JII%grG!f_dtm|k15CV3N3wQ=_-mzOtHY+f``0y2M0~k&vx|!ZQ0h*U zi|lC^d4QXS(`v2V5RmXBae;2%p2!Cv_3}E0@vOmGieOtk`r_YQ+5gI6sZ^s=#E?d8 zI_AD@4GdK3S`*ZxqBh`}6CNT3UfiEev=JBOnVr-8aT7 zge!}M0$|^O=*vq?pdxF0FPZr_6&wk9@fRGO{qg8SDCgO1dF*qeYwrE~rs9PKp^Sg+ zhBIfTtQ8cZP=c!I2=b?<3Pu#qT>&g58lll+R135t)Y_oh!&}BJHK00GBY#NH*nCW99@&1PEpqgWH-@q03H=0!dlfOjRvUMboD5Q9^V-1{1A-$i|wy)v<3a48M%#;)AGz=^+E z1qKvc*+i9|H`qP2AMu&b>3GA1xY>fAlWS=8GpHAa^npBcbDyEI{|VFDEhEm29IK&V zM+{YD4B`O==_7|IYmjL>)FoK}#R9_%+$YOGXrVHLn}aSWx&1~OHlpuK|dY+EXdjD1A+5pRR#uvW;Y`JHa&~qG)VtTwam3^KWr&qXr>F$r zeuq~x0_~{ZSs80at@(78MQBe2zW}9&NnotNz_N-8;bkV&CxTc)kuvq>O)zSpG_u2h zORDq7%G61u@ZT(hF7`B_F_a#>(cVv(3BUtHi-Dcpc^qLAl~0?I=4_RyEN~L=bHT>< z4DvC-)9u^$erF~y9Bw+%-~R~&;e>zK6__<5EgtoU90}>r%1RAY8@zGNQ?{(Ma4JtCDO9qNhA+XyMeKBGR)hqZf<8Ja7PC7mB)~$c zrS6=Z$$wZyT@Ghyq|FFt7v1uesc8`0U4T^9b|6tVU@3!RsAe@e@jA}N$Puy?6kZ=# zd{`_NYReS~XvOAZq_&1~I8(Z!oqSQ(T|b3IX2n@#H_i_&hw1ZBjA2#?H()P51E2D} zbd@4+`>KC@3RI##`@j7qGs3RJT!~PO3)iKALIT@mziLy{sf?*O-gWkFFE7~+!H+K^ z?1+8M@R}Zgp}f32`{j)0aPpXV4*anCs#TrX{=!tVm9k8oSejW6B<4!rn3%kPM>aBY u(4l&`uoWl>KexL(PBMvq`JY^o)Fp*mv6lq<1V)x5{J1#Nop_G^C;kmqjU;aX diff --git a/cookbooks/tian_parameterization_kinematic_slab/doc/tian_parameterization_kinematic_slab.md b/cookbooks/tian_parameterization_kinematic_slab/doc/tian_parameterization_kinematic_slab.md index 5f20ff7e5ee..93e81355817 100644 --- a/cookbooks/tian_parameterization_kinematic_slab/doc/tian_parameterization_kinematic_slab.md +++ b/cookbooks/tian_parameterization_kinematic_slab/doc/tian_parameterization_kinematic_slab.md @@ -22,8 +22,8 @@ schematic diagram showing where the model is in the context of a generic subduct :::{note} The low fluid viscosity and high contrast in density between the solid and fluid phases leads to large fluid velocities on the order of m/yr. Combined with the need for high resolution around the fluid to resolve compaction -viscosity gradients means that this cookbook requires more computational power than some other cookbooks within -ASPECT. +viscosity gradients means that this cookbook requires more computational resources than some other cookbooks within +ASPECT. In detail, this simulation was run on 128 processors for approximately 2 hours. ::: ```{figure-md} fig:schematic-diagram-overview @@ -40,19 +40,20 @@ Schematic diagram showing the model design. ## The input file This cookbook includes two input files, one which advects the porosity (free fluid) compositional field through the -darcy velocity, and one which advects the porosity compositional field according to the coupled two-phase reactions -solved for in the Reactive Fluid Transport model. Both models feature dehydration and rehydration reactions using -the Tian parameterization and both input files can be found at [cookbooks/tian_parameterization_kinematic_slab](https://www.github.com/geodynamics/aspect/blob/main/cookbooks/tian_parameterization_kinematic_slab/). One +darcy velocity, and one which advects the porosity compositional field according to the fluid velocities obtained by solving the fully coupled two-phase flow system of equations. Both models feature dehydration and rehydration reactions using the Tian parameterization and both input files can be found at [cookbooks/tian_parameterization_kinematic_slab](https://www.github.com/geodynamics/aspect/blob/main/cookbooks/tian_parameterization_kinematic_slab/). One important problem in models that track the partitioning of fluid into/out of a solid phase is that these interactions can be much faster than the time step of the model. To model these type of processes, ASPECT uses operator splitting (see also {ref}`sec:benchmark:operator-splitting`): Reactions are solved on a different time scale than advection. For this model, this means that at the beginning of each time step, all hydration/dehydration reactions are solved using several shorter sub-time steps. In the input file, we have to choose both the size of -these sub-time steps and the rate (or characteristic time scale) of the solid-fluid reactions, and they have to be +these sub-time steps and the rate (or characteristic time scale) of the solid-fluid reactions, which must be consistent in the sense that the operator splitting time step can not be larger than the reaction time scale. We choose a conservative value for the fluid-solid reaction time scale of 50 kyr for model stability and computational efficiency. For a production model one would likely want to reduce this value to be as low as possible without -affecting model convergence. +affecting model convergence, since in reality the timescale for the fluid-solid reactions is likely much lower than +the advection time step. For example, at 300$^o$C it has been estimated that to fully serpentinize a 1km line of +mantle would take ~ 10 kyr, and this timescale seems to decrease with increasing temperature {cite:t}`macdonald_fyfe_1985`. However, smaller time scales would lead to very large changes in porosity that make the +nonlinear solver converge less reliably unless much smaller advection time steps/higher resolution are enforced. ## Model Evolution As the slab 'subducts' across the base of the model, the hot ambient mantle heats the slab leading to progressive @@ -109,3 +110,8 @@ can partition into the peridotite mantle. this dip will impact the magnitude of the slab-orthogonal component of gravity which has a large control on the fluid pathways. Varying the convergence rate will influence how much heat is able to conduct from the mantle into the slab and affect the timing of dehydration. + +It is worth reiterating that to extend this model towards a production simulation, the nonlinear solver tolerance +should be stricter (at least 1e-5), and reaction rates should be reduced. This cookbook serves as a base model +that showcases the use of tian approximation in the reactive two-phase fluid material model, and these +simplifications were made in this cookbook to make the model more user-friendly. diff --git a/cookbooks/tian_parameterization_kinematic_slab/uncoupled-two-phase-tian-parameterization-kinematic-slab.prm b/cookbooks/tian_parameterization_kinematic_slab/uncoupled-two-phase-tian-parameterization-kinematic-slab.prm index ed230596188..9b6315f6776 100644 --- a/cookbooks/tian_parameterization_kinematic_slab/uncoupled-two-phase-tian-parameterization-kinematic-slab.prm +++ b/cookbooks/tian_parameterization_kinematic_slab/uncoupled-two-phase-tian-parameterization-kinematic-slab.prm @@ -10,8 +10,9 @@ include $ASPECT_SOURCE_DIR/cookbooks/tian_parameterization_kinematic_slab/couple set Output directory = output_uncoupled_two_phase_tian_parameterization_kinematic_slab subsection Stokes solver parameters - # Lower solver tolerance to improve efficiency. Since the problem is simpler (i.e - # no compaction), we can get away with this. + # As the system of equations being solved is simpler, the linear solver tolerance can be reduced + # to improve efficiency. Note that in production models this value should be tested and decreased + # until it no longer has an effect on the model results when using the "iterated Advection and Stokes" set Linear solver tolerance = 1e-8 set Number of cheap Stokes solver steps = 1000 end diff --git a/doc/modules/changes/20240605_douglas b/doc/modules/changes/20240605_douglas new file mode 100644 index 00000000000..503e31c131f --- /dev/null +++ b/doc/modules/changes/20240605_douglas @@ -0,0 +1,6 @@ +Added: Cookbook which showcases the use of the tian approximation +in the Reactive Fluid Transport material model. The cookbook features +a kinematically subducting slab with an initial hydration state +dehydrating as it advects through a hot mantle wedge. +
+(Daniel Douglas, 2024/06/05) diff --git a/doc/sphinx/references.bib b/doc/sphinx/references.bib index 5c2df4a3745..023f58d7473 100644 --- a/doc/sphinx/references.bib +++ b/doc/sphinx/references.bib @@ -12214,3 +12214,14 @@ @article{tian_et_al_2019 doi = {https://doi.org/10.1029/2019GC008488}, year = {2019} } + +@article{macdonald_fyfe_1985, +author = {MacDonald, AH and Fyfe, WS}, +title = {Rate of Serpentinization in Seafloor Experiments}, +journal = {Tectonophysics}, +volume = {116}, +number = {1-2}, +pages = {123--135}, +doi = {10.1016/0040-1951(85)90225-2}, +year = {1985} +} \ No newline at end of file diff --git a/tests/coupled_two_phase_tian_parameterization_kinematic_slab/screen-output b/tests/coupled_two_phase_tian_parameterization_kinematic_slab/screen-output index 41a87386243..c3f8efdb280 100644 --- a/tests/coupled_two_phase_tian_parameterization_kinematic_slab/screen-output +++ b/tests/coupled_two_phase_tian_parameterization_kinematic_slab/screen-output @@ -70,4 +70,3 @@ Termination requested by criterion: end time - diff --git a/tests/uncoupled_two_phase_tian_parameterization_kinematic_slab/screen-output b/tests/uncoupled_two_phase_tian_parameterization_kinematic_slab/screen-output index d9cfc8219a1..82ba15cf00c 100644 --- a/tests/uncoupled_two_phase_tian_parameterization_kinematic_slab/screen-output +++ b/tests/uncoupled_two_phase_tian_parameterization_kinematic_slab/screen-output @@ -66,4 +66,3 @@ Termination requested by criterion: end time -