From 84f9135e737147fee779962a3b3f33ea40479657 Mon Sep 17 00:00:00 2001
From: Alessandro Maccarini
+This model simulates aquifer thermal energy storage.
+
+To calculate aquifer temperature at different locations over time, the model applies
+theoretical principles of water flow and heat transfer phenomena. The model is based on
+the partial differential equation (PDE) for 1D conductive-convective transient
+radial heat transport in porous media
+
+ρ c (∂ T(r,t) ⁄ ∂t) =
+k (∂² T(r,t) ⁄ ∂r²) - ρw cw u(∂ T(r,t) ⁄ ∂t),
+
+where
+ρ
+is the mass density,
+c
+is the specific heat capacity per unit mass,
+T
+is the temperature at location r and time t,
+u is water velocity and
+k is the heat conductivity. The subscript w indicates water.
+The first term on the right hand side of the equation describes the effect of conduction, while
+the second term describes the fluid flow.
+Spatial discretization
+To discretize the conductive-convective equation, the domain is divided into a series
+of thermal capacitances and thermal resistances along the radial direction. The
+implementation uses an array of
+Modelica.Thermal.HeatTransfer.Components.HeatCapacitor
+and Modelica.Thermal.HeatTransfer.Components.ThermalResistor
+Fluid flow was modelled by adding a series of fluid volumes, which are connected
+to the thermal capacitances via heat ports. The fluid stream was developed using
+the model IBPSA.Fluid.MixingVolumes.MixingVolume.
+The geometric representation of the model is illustrated in the figure below.
+
+
+ +
+", revisions=" +B2UkBoI`rmTta`uD46bIa zzbPn?6s(7@NzThu`IjMu014*SE@{OD=bw;}5E2s7v(=92avk7Hi|hK!V=Ej+uZ8sh z#OE|KET0Nh)D>`(T-j_`^zIr^0Nf|vxn(O1gF zer0}eDC@0sJ~ |A z1SIkF 3~b!Trl!)8liA Z`|rf+b$3MOoe%uMVa}@c_n#y`Q-%8mDKo0 z__+}u#ZLHPmyMi5HQ7P3g5FDgl;fL6zjH|L=tdDxPJM@y!Clar18Q-j5~o2_p*l0~ zN P|Ti#6u%inbF%& zd*PcmKVv`rvqSjj@t>1}KfCmQcA@dGc&gGZ5Pw~&xf5E`%F4=xfPidi9&AEF2OzV4 zZ`N7jGX&R@f4 c$d$eA#^hY_2c4TMp;f%UVmpBhs#k?9`8q&HQej`8}L3 zv3IiVPfsf<;_0#n;G%M!RRV`N@Qk}EdqH1NHvu>+HkRRmd+oR7Ws_DvBOP31*!mDZ zLPbf*u*qZSuvHEcFyJUCap(dUekoVQ-O|(Z&*_PIPL6_b?r;MpooFp=ruAg9_10GC zGyB%2CK~A{G **Ub29>be8xk{R)O&7=#t_B z3VLeIB}9DVz4WWEhsE8h>->8Ze2r*Djy3aD!cXZj&WB>)U8bP=#vs*1i-|OGTbrn9 zYeNnl+W|+Al9KX58#ocFYV2Srk`(8nXgu~e#zT*e8pic4rzHy~zsVyNrXG3dQM`C^ z*vhkV%zN?mD{;;{7xC~Z7 %r-eQK)vK>xhpQVa%|sP5hc#QQAdSa2b^|PWnC^F^XmX?>=lz!L{o^N6 z=cR)vn0fT#6B2e4WsUJL7|Et5z804h;*4Lxy>aEk@o`fzS1@9EE I8%NiO_i{Is6Y4Q}q+UOx-0_AmkD*x1-qP9u=}5jGf*H6AD{M;|&q_VDoN z&sBtWwYOfSaG>S%XdUuSwZ>g`*C5eX9HeHF@5`A~D9ke7KMs 5lCG52YH zKID^z4e>D +Qgt;o>=?qid9jdO_KE=}Sp~U^ zYO~ooINWw#ilT7j`Fa1`3uz)=Rt3Dw;ZewypI#6TP1_X~6x2iME0zw9otizC3>)>> zotGUFy2;0fiG>;eJXOOv{1+P}#x;EtpWTj=-v 7xf^hq8Lh;D1?>P89xqyUtk9W@%CG7L!jkG;vgocKn$ z>i6F19KjS`pFV#=-3?~M!HA>D4)avYJB^?2PhMRO47@Njjc;aVcJ%wVr6p7FSit8C zrta=_^*Q&QoW4u?xAGx1t*!qHdNd9qCF}kF`0>l<8ah6{e{TK$-L`kZbMgdQG7cUI z>WzYB$|_5Xi~8-e+&z7LuToRvTq6;c%qrNy0ECwZ(xfLvJ->e&Ddn=%${Q9%17!Ep zMylLl@d(WjbS~=sURim!;!PC^21zI@o< 6@#ry) zvckk*bbUFZtv$+fs5n1AkBNn4sm&aG%;Aj}3{ZRi&JS5x+1`4^!nz3lu zN2`GnS0zJU#m6(CXyOU(c^|cM5gwEcd 3HQP7;JX&A|&F1s#j zHiQumgEI7{a5<59SVRPp#n#G7ZRuych?SzVvz}g3RFFUZO=mKc3{qMVTNgLKBZe$TT?4A{6$w==jF?IJy?#DbJfSR55UCU zRaI4%4}fVD8;Ys;$w)chzevYk1GTje6P`r~MHu?T%;5#4&`gJ^uvK{~(7G=y^+$Yn zqQr85-Xj@Moc<>EuBjUl_6z`F$nRMv_gRYYO5=}WYFTmKyeUeD%1mzJGo)8k-&D=d zCk6>OSje4=Q}L?!a>QU5d6&U#^p0CErrKK2gmV!By>~ooiivLwa@A7Fy_##+ALiu) ze_RYF6aVsM3+PQMzZ1>JZ7?COkdbALSJ%{Xv~%RPupnn|-;?4)HEH4mQ#kQ_u`MW@ zNk3P`_c?OFz`Ce+BGtzGZA$3x$+BqUwAN{VbJK52-w&78P}J*RrdvEurVMT(2zN+B z^n%ICsdPFsW(jwfJ&2k}jnzDHo7ktRsGszgPjM*t_QveHReu&u=a5nJIobQNd#;@l z^A?_+Q!MkueWe+8*ASr4!*Os%dDw7pao0CD6ZW1VX4 g6WgbGud2htE0Jk3X-xJ34mG z+`Qbh)?2DjHT(gm3K^8|1$Meg0R!~0@Vx8Eit&g)n~u)I4AOF5r7+Ua&`_&I%4=h+ zVt?OFpBZPDa^rF^=$fpq0@Ub4C>EB5Kf1OQZ^lu$)9>K-1No@%u1SCN2whLWWCt*~ zC>_%u0`}8d`(YY0yzIP0$|HA!0@sQ-uD8Jk?H@goC-6Wyi|wLF4h W^a?d6`Ec6PljYHOM?(J{SOyi*}B(=VcTW||a-LAOhT7d)(4T2?|vcI6% hK%vp8hhPKsbr2niifKHhs~fE z-)i80KeB57U@!c`DmCziO3KPTmleVe%ybz5jKfxRF?AJubsWoX@@8Gk`15r+wtmf~ zH{ohtI6Rn6Yd=t9X2`6fd+F%S#!RcfrROgob}A|=oZc~~!#I3xc7KvgM@tJDW~|B- z;`!(?GHeK1Is_s!3Jg)-iGxLiZTewRQ82M9L2OB?0tP4*cCg8ZQ*&7FW;@7DQ0JbS zs-SHQNnM$So?1~sDc*5VK}(C2TTZEo;_U^~$B%iaFk1N(3e8-#M*s~(OS0J(ze!48 zQLju+?SQF_ZQW%KYSqGRs|r1yW^V`#qCYK=kI*TB3IFoxZ~s5r6ieTp$cce$Y^1lb zQRF!$K1@BPz1DB&o8as+fxnoVnp$Bt72xGg8tEJRowlCpqDn$+Hun-0qJFw^-K=WK z^Qmv0veQj@Pd&G}ljmV$ezRPDjAvJG+)2ULkbiousOV xquLA;H3B2mf?tdul?44eSIhMHg?cbI z(ZU?0#5PJ#YTQ<#u%6dfO$`w$EIdnuua>tE@AWld(K52b3_s{aaBy%#gRs6O4zKV& zW?LMqYWDuJ!S!A%1@;8Tb={juC8>(yn#}#aX}50RlfU7*l3r(#pPw7vI+42isPHTi zBDf yg0L_#EyP(+!BAjY9M?xB{cqnHgn4gtbPPujXRZV%24n6LF}l>_ zJw%&;7h!$DK~zfl+u&?w{ZjV(j!5#`$2-#csW`l3Z*ham%j(DzBzSmYlVdd);YV4B zgUCa9@GK|>Q1zXX1vn6ghW1~ZJg2>q%`kpmy?OIyC6dI9&W#q^)ZZlIP}$B#F_Ldk zV=agep<#W6*Vc5W4G|j9M=?x+eG8-|DMNtX7e@bdSqu&jUtQ;CWHIFqK*Zn?Mv^R| z(a&SZy&%^tWJf&-F~vo{dXgwQK#c_OmpmX>;FD@XrFuZR0W&I#eij+;eG2TkiD+tJ zVc~+GK-ro%_2Icj!*B^^*a}(pd >Yzq-+-3Kl!!=R6% E!w}17V;04FnL`WVs{5cGPZd`>uf0&5(s_Ff pLsFv8 KHN1q9Wl8{qLCkFiXB80dKJC&`16X-1JnfLDFt%#s&aCe+1R>Mu ?@`d5DbtUg?=lu}fTES?Xx3cZPu_1Mo|#jpqeb$&h} z0Ule 4;GoqS5Fr43 <)tk%j>r_B_;S6(SO@6jnLO=mJGxAmI}&+g(YQW z m8;q=- zr6ZA}gQOL!A-u}JU9Dr G(=&>Gsq_W>;(Nl48}b( zP}qRnPp0fz!J9X6+8GGdg&0YzFh(&^Q3L)USYvNQ1e%9P%9kiudW7IZcqu$UT62ct z;nB>;GjvP16SPP04Q@z5f`K-`xXf`^R?6NU<})pv)aMw}+e-q}>Z|MZ3^biNSnlwL zLRyAE=?J-GIh!w$eR}plwFk(E^B3BK9h+N{k|^+)C=}w~TYmg!eU}Uj3}kg#gk0B# z`}(46F@im?2`?hs-@ktk1d8JUGWI`=Mz+$Wf-k=h7oPp;(WTcld!hqNRlA7Hw7t1^ ziI?WS+Pgu9PHem@%O5kk-euuki}~<38OcHYk8JQYmKhuwxt4mKG5uX!z+id)+exzu z3ill)Vei_kmh>5xJ@#84Ax6_*JQ~E3U4F)R4kLhqh=>3FlXpB>cr+^Ue^bnMeM{|w z;|K*cb#+FA3yk)1k&;$pBjd?&X}+Ld8WK=R6VcXJ=amr1Qz9?WEpleDVE7rP6eMx2 z(aYs Htj-;+N{!!NEc3lz>yFe)Fc6CeSNx;sFv9dGLXM1vgt!t}AAXHmzx8b7+&C zjLZ~U4f9fBFFFifCu?x;Tk9t7k_pGg#)4{@3s$~`{48n3cADXMyeUe6Ay&%9O;1VZ zcaY{D-*9Q~LV{8^^=HuzZ%lAk@oXn^uRJLSG9y8e$1V%Cw7-GG92`M+v(itzEzr@> zsZ{7LDlQ()b$37X^mJ}+?n{|TVsrqaF4LQV-|(KHt%FUK-u2E)UC%&@m7F}TX5w$> zq>=nsP17b7hv;8%{rZ36`d4Zz?*+!snYs}%R8zM^>|Ca@3FRKfa5d*NPlbQMol3 0jWZ&}wtrj0(p&7Q7?d` ;MNm*ML7y+T#eA>K z_yMMPzAasFJ+U(u(t!8l*yX1Q2}y@kbBI^M&Uvq2zy8si5c12@atnEjJoz!mhyPP7 zKKchd8jQ-ek;RI%zmit6;pUKmBl9%M-qYD!dbQLw4qk>vil&OrRU7gpq@<+XSCC;J zfL<>U5)g@y0s#jr42lHU*Hv|ZQMTqXqN1Sz3?iL8s)jEMF70Bdc@dc>RTVZsV#Fq* z;ga9^tWnBB=3=9vOVf&90)f6W osmj!XH l)|{Md3*Gf?dRp4%@^XgM%Rr2>UdSp*8jpGw zeaooDP)+U0$&1ZC1<-lAxw(P1P}2Vo&_FO1`3WdMgPVkwp;5{ubcJ|`x0mZv?`pk} zOZp0O==s_oBQ;S-`WG>g3U;!(J_5q}ZY`&$FaXWT$*HUXshbZMf@9jygr1V}k*!s^ z-gW3?D*N>TYsf=#DK6enkn8&S^QWjoYe`=+8yUPyGjagmh(<)j#eL=H;9if);AN3J zF3(o34ryn%Yc;BEqrHIwn|%sY8>JgDB`&V6J($vAKXc3&K_hV89189cCUH8lM1` zSp2SXS{p9>`s$TpBA2VX`^Wx)?*oTg-I8~euZnQFm^6!vuk;RPWf{cLz9Mt+g5gky ztn3OVsFa#+o#*(+fbP}B#Y+1_I&16emZ^*=c^<@c`*?CPC@6qWtE#MwFH@lv_u}W` z8hDYYYoPggaxzkQbkd V#Dbr-{xmcj5W&L1iC=SA4i z07KJMfqi3bBi^VX0cy$keQZYxOr%_MLyLrj#Psy2j=uWyv@}`}c~VfYQNKyZ%*;du z2L~P>gCVnc&53J48mqKWP-w~UdW(DtTZF`geMN7uIWtbZWj-bbdx}Wpnv3#E=AiPg z?t|M04CpuhxUW$7*un$25ySoOr>7tltL#V2gCK;<(cZ1LwXCcnD+?wtB)ibB^wxS) zD@_f&5{^ra)6~}wZ __L>UvgMmjXcTFiPfcszuY$q6TLr;b51 zY%+8G9n*hQ0a6FqE zD3Hyir60-5>pWl~a0LOa?R?$)D8!GplFYwqiG9#C&&|Qy2>89;;1^ypm`LkZt6Hlj zt2V1{s{yNFs|gqyP6p98UU0!)xiTo687V0juoc*Pv2U6;rq;TxvZ)W=GMpwn!Ye+X7P}?+BG>W)x-UD+0p`e-=|7KAw$1gsPZ70b*CGs)jMc{+T&PbN1(5R)TJ0<+iYkq$d@F$}N#~%+s zM^E$M_yk5zpj)?aX0i~3yL{m=yMU|-^AW(PJX2yOYO&CP5W*l {)+-u|6*cH^6 z$ZdCLb8~w)H(M8%b*;|D4-oif)~U(K;e|IuVlp!s#v^SPWYpBOw6tKD0COMk4sAhO)nN~NW^WakpJ4RNxOs?#R?4>)m?pWSQxYKcG@6I(jZMl${EZO@+ zDM~&*VksX-|J-TtAi}@Q9oSjJbNi#oL9*^Y%X3&-`qmxLJ9GQ`bz46_V*w8229Q(s zQwolL_Vzxu#UKr=$<3|JTWSlzy%Kp%{02V#NlzSWoyk@k!(T0(rbRJ~CJYR b2FgN$c$f#dgc?6K~wQJB3Loq0zgM)T&x)<7pcEn=m zt@x^kknmwh*UM#CP+0irk^R=j1`Q2G5(Rx%pSR~Dd8L*F&sJzxM09 hvPpO zLGC)o_4xMm);}=L2VtPaXQpq)K-=Ye+BlPttT&9011^JN`v!Rb>gwT _d=Wn}_yVQ`+B>L^Fjr3+WG0S(8Ly?IH9QKmgsrWOO}jS; zyPieIfFvJ<8Fn{(KwX}S>qWqz4vibp@|nc+Y2Y(VdV~avkKSST&ilEA_NQUbi0RoR zbbRU$!*->Zr{<&+>HkbgElQCgWfj7n$>k#dIXvlvcJHjCEMhYZHAeLt2@wi~Oi&dl zNDLKLORv9uEEN;46Z-V2We$^t9PkxgGPZh285t=lRFrgEo11nHHa50a1*JtCoSf~) z(`FSG8KQ|VQzG6y7U3Ge#;c60B7BW|jC0Dv!2$JN6Pqp^2Z#OPM^NEDc<_KJ_TL)i z3&gUnuD6I4|E*94-9fos5C+Cr_ihF6)mJd}YiMZ5?T`(QhDuRR4o4Oj_rm_hUQ@-V z-YtIpyUzP@<)&j|wpaiE<;iA=6$@j&G#6iRLz)$MInWiBy+#pla|%tdP)GPv%VD6U{Lau}bGKwcl6z*^^YTqshz4NH<%^ z)^%AApEI(x?M%+^A=Px ERv@HmA3j;KMQpW4J*8A?+PHi zM2|1Q;G=X{zVm!8sJ*MpP;T#xTfyPVb OdgaU+qM8ZPnLSX+bReBd^20*N{} zKjevm1f=3Uzpe4(S5+@xRe6s0G%UYxu6$8>lceesN#&>cn;IWqXNaezSME8ijn`Mn zSMabivaru-wI#0{)^{I%J``MYv9+*O3T8~t9&BVIeq(BhcE4M~1l&x |u-`u8#aRaD!1 z<7@5WQJCtltfzDuK(#$P>j&-FbG-)zNa#e}qQ1~x6NS`wOZ)ZvnWJxhn8>bMvJA&x zms@wz9Q@^gN5e-$Um#oxooi0Etbk?SR=v9g1K1Ju=a{YEYr0O1Pl2UcNMS iwX$#EA#G}{m-?Zam0|6Ri^HNfX_s^A^jIOPJ zX}?0w>Sl6T3yW`XKgQN6_-iOkIdbi&WW*%zddv-GMOs?Y>5i)y^wVJ{i0_4cy#1#M zJGd`+_#X}@m+RAh`X00ub@o#)N@5V4c)?oC3H;d7Oe8zQ3uZCdwnga1c}x`Dt)uC% z!qk5n7NdtP9ZV!yt^Phf$LNbs)O3vSv-EZOQocWrBnvjm2>3`e>VPeSRmV3De|<%h zjD2zfnE*Yu?x~gYp~FNLh5QfAREX~PFDhve(us`$>gqFN&mG{q$t-;l3E<%reHCg+ z&IhXan&&OmpxA&xO?bLsH;@qM?2LzpM?gS;iyJBA$p Q@Qb&C$9waMwm`o=nENEhVZ==BZ7SiE6wfoqO? m=tXJ z4Vst})^)}E*0WMj0CQ=$hT;>MTbZ3*Q{!s&7!fp30>ixsrNZrcz(k-_wPpPa(C3Kg z;W)geyXd?(R9lYfi7~3790b9+rv}6}Y5>fB_cTy=J$hthF*g@9eZLwSAbG+qC;(J& ztFB#I8sbqLbaH8FH~RZ)Zry?+2vb536@~djfnt3B-l@pjor#HwhqYGkx^_Z>T}g?@ z2)c2fRw}G(bX&ON87B`SAamulJiP;k-Fq(g9=8Q3G?m zDf PqS`iC0 zrLpz^P|15z!90_sjP2Fc!oothtp^u5pT}8oN_9QUfBjlZM`wR;k5Qh(|88Neo2y%| z!p1f|1w}n@4wF}RDPiR?RLKa&x$fu8y&6;e-VDQ1KD|r`u*(6Okj%plj)04bpmxB$ zd>PonYpla{Qn@fuTWfM!>|87Qt{cC`GLcqXe|id3F5FuRwS9>ag!Gca@n0tCH64@= zl0FJ-I~Xi(t!}mG>AZ|jqCg-RSik6~X&4!a0|9xAlY6A6qr 5$$hGS9!7_4*Id&eEW)Q@pc(h_B}@L84&3ItennBvTl7=TcUBm^!N2FB zuX%C;3%B@iRpK0&^PUCpkLbi|LG?QQ1*{1OETW9B&XYLqjSUq;pQoamBcY9Vw&Mh_ z2A&U$N$*j-*I-poPi5fIZ@7li( ju;b$pKmcJ!_x8p&&VZuLDwVnvu9|m0kZ*!p849!s0kImIU}WlEyb(^=%=ni{ zgM0U^y>F#JXAf;mdODq@rNY!|O^>{mtUmk6$!`!NeYiS}3sbC^UYIFCrw?qV*V-Nd z9v Oqm0vRUD~P3h0`ebyovu8doiNx$=#(iH7Ek zqnFpwr(DH?XeX+zZF*v24Lv=f!I`pR`NB7LkL+L849?sZ6Z7%)HL6vekA~8=Rsfkd z mrx @e0r)|W}I;2r3VpgQ-lM4F7M?T zyS3>Kga}CZfY?^~baLWcS@qjvrNHJ@dU|a$Q3rYwP>q61>aT~@Va$1}=B$j?V|h2e z1tHX;*woBz2pO>RaR{Spqwd>Ai)&TsiGBxD?v`(_8K3{(%WD&V3j+#GO@L3#OY5L- z (8@S-S9q;OD`^iQ_gNLDq*B9(ww6oX?F8#_@G@ zPG2?ZPktT(gi_f-`~&_Whr%4-^sKGV25~{3ro8N&zWYc8ysG{EJ40DzYi&(qbO?A2 z=LQ%Zyb8kv%C~ZeIPU-RbIlD+)mr-ML^L$lm$s6+1O)`L+c@8Nj#dv3v-&4s(fg>O zZi~gr*1hlT&Com$l@MUDCIy|Dp}zj}&6i0azyEgtT#Of)2MT7fvlFdPrSt*SyR2-- zr}^y46bY(jJWmNVn7SCqx?X^|r1v|_vTa2Fa~1z}#woQ28N=+rld@d97A}7b5-+$^ zz3O=Y;jXFYHWyf+m|0M8SHQ^Z-v9ET9gZA(kVXnE@4HG}xNw1klf&h$+$D@GeI-Rj z5;B~#U{&hF`QOKbtP>Ux@dEjj(&95LOdV})7_4!wdjLyN1^OydIyx?EA_^Ed-bhYD z{Ck%GS?t|_Yva7B*|!ABcR;}QAkAK{2aHA{JUGL>Mkpx+Nh>}&ISC(9ObzJap!xl0 z==+>mE}mFniUlJJ{N4*Sm&Lg7);}vN3y!t-vi-||LBr86Um|-kK7IadY0mNhElc72 zG@B*Oy+Th$flW)E@m-ow{59#sYQ>Xtyxs*x#adUgWUjY~^fD>npLMo_LEu18Fv1@B z{vXEPI}pqMeIG9h5pG!}TVzX=Ju^x+k;u&6Wv@G ee!jo&^Uw2-?(4c<*Xuf8=XspRahyOF1!7WO0OYrV357!CS?TEy=bpt8 z!cAo`6YqDP-4KAd$AL34A&I@W;Ikv9^icw~6#lOyE-E}Q>BS3fi*+~ JVateKr@X z0`b$*(qccdx#q)n1IpLhACf#g_C9vXdR3|CQO`95`v7OrumtY`$T*qwFYv=cYl}q) zdboNVW lQZqoABE#(JMTwKCmp?_v}558QK1_~RNb~zODykjpi@ 3kIXO#JAGn*MZXHMG5P-8aB4ox zxXLP$-x@ZxGuB*AvCGds9`vp~QMvdMqt8_8(;8+x*^)pcy!NTaRYPMD7|*69=;a?@ zcrDR^tBL-Lk&?>|ESuzgF_$H*AdyIC0tKUVp;IzHyr<98Nq9{6PnqZ2n=`s(Y%~Tb zFxC*<(wBrj`w|iPzrQSP#SAT_h|@h2KgW4_csaRddnArjKbx}7Y7q$94#^l;`XsF( z0gvT4P{MjzqzxQ3rac+FXLF^trpJ_yw<^KnTtpT}lZ>FU+>G%-WLWQWvDv )8*H5NLz`*B0=~$H$_TFXyf3Z^JKnh0IwPjaA&WFl}sHN;PjCZ$hf%` z+WnBzg^Z2LCh-{mH}=zZ0JhGQS<(c>5!*G i#SDP+!wTGZ4&r2ZZF=2dNl75G+#yG_;8`hln zTcRL45antO*^QqV`@{nW=-}<1Lp|>Ts?q!R=lIuhGBPxtl)#- 7siUKX1!hyHWT1gvaqK&OHIsWDbFOlst9Z{c-APif5=_`TZ4phL{~&(A+q zV*3#>l`txS3!9x@<0JpxyKjuCkx `n@j77Xf0@Vj_VTyiSt_# z0!voT_YlHxDIjk;D2@*ry}~X8fnoL6E`CG4QFAqk(N~v(A4oR{ULN2^(kF5y8_~A6 z;fEf^Jse)u??%t`SdtOw155dKW_I0C>#D`Y_jVI<&|mn!=Gdn(c8{L7L0bVd0MUn` zqZNw~bO;j5$jd%FNCXvC_}jJs;5{JP;)|}y#1KKW4!W9QKvsbw13FdE+lYS+fhPTO z&D}eCdd(%euLy$NZ|QF2rOkf;qlk_pcW*GkGJ!}uGp8YFa_Ao(&7CAH4ltKP4!{W3 zjEuRm+pK5K9N}ja6#Ve^?NvP`U&ZiRW*nwrYDHJ1XyQR~dL$Kx@<0u#!sFPsx3*MP zlA{eLzKsID{p2!EDl02SbD9O@5cv4`;Wj8r$~C0WexfFL7lP7!Egj!vGJL<0Z_%@9 zhOmrmi-qKa_P#e02ym?&Qdpo0nI QoVQ1gSVjIPZ35a?Ynp+<^>g#OgZ=$k+^K{1nx@H{ zlU3ioJ(|S{nfUK}4tx?56QT`0fY;o~hO?X9ahv@zUZBGJuRn#ZJza*ot%2vV-R %yD 0*xJnO=y8W3*LF-asS_?GgfFm%`*VNWFDXVc^&3Fo@ zAr<69I3|E&A_N@o!=r|O{yR2Ez_fwn 8$ZiXLW%dzwD#10&G1|88~ zX?Vs6pdw>u?}7Fn`P=ayKSEk)K>LlM_b(> gf%z+Of?Z!k^6`2Y?U8D zpaVDb>IkQlqCE<>B94`Ng|pO%Q`KKZz~zPU)?s7d? @c>Q#Q6XG4f(Nz1U5LA?*IfrGP+K9Z2xEy}S`cw`;nkijpz~g9Fgk z;vFkEb_hn#L!5S_v0j{B#F528!W@{WHdfXuUmoE56{TE?fTJ1Kt$s}}gS)o}X!P;v zc7;F)$KFZW=_f~Tr?p6I(2gb7)6vjUBykrzjRefHgnWhV;AQ;}AyC+5(p2wR{k6jO ze-D(E7AZa)rr*HphOw>f+Yn-4oXLw7pOXUi_C`v@UKjO#-%#atVBnDCYC%0V2fK)m zb;kPdWG$)Rb8g$a!A-UFB>>pPJu0X!1Vq-n5_s;+`NL`wym+|gTJ{&b8EXbPc4o*T zq#NLIkf00cL3vMG8 z1#4}$PHcX>B>qA8{D xVwy^<5Rb&$fDKPO$GrWU(DfMpsTW$>AQC0199 zT}|45oF}2p&v{<+NHx&?dbCi~u!a&jKrT$~(kblG!ISO!0?US^wg5TU* @HdWW1qaqi&8#w7k;+*Rsop|s*xc4tA>fNs@QSIeZ2bCbxR^%&Cqa$tp-I9 zHCmgSrJL&O#}fp-%+1wsqpInFPm}a5nkMKZXG|0FSgG(>x(NFApS)!cis0x&vGS5L zaC|fbLH-8&W3wQ~mpGVToK5f4T?0$omJ4NLL*=`9d8{s_;Izez0ycO*H1*f1Yck*p zEj!(2QF{m#aa$={E+_jSvF4@pU+jatI2Lk Ru}-Dk3E%Pa4uPFa`W2+5 zrJ-RAWDKCsfUX;`uth4&gZ2mKQS*$`p_{8*a}|~noX1GUllF&()bUt!Pf}0?Z*YpE z4SVoF+`e0zb`0c|?B0Z<;; f-}kIAo`j-;0eBONf!WxNgrJc~PLjv3)80{6nV`C2~x!28A|mu{6W_>$JI z)MD33F%N&ZV5bVYmK4&{F?g+|&yte(BvsxQ7P5&Ejh=za|F1YU&bTMoFM+xoBiHII z zmlZm+u*+DpNmhZFj^hPe7LMq1JapzjmUpwe=-;N&L|6 zxRCQ@h$I%rue@8_qBg27RlZsJrcA6n-l}>umVd0mXEAeXgCvuBh0qtSHv-$qj~ z3Eg>MSu(ojLns+;?o$O)LsV1*7&<|U1&tY6DL;c|!EE640=c52qvJR!DT^{E`Ym+Z z*1P*V&PeKrt5XfYK93rG=(jtDWfK`II~i**lC)PVgr4j2Q}MJd$%nQE+zs%^^2QjO zk(~8_ff)&>CHHp~r~~L~(zp+90_l|v_v{=
`U~faEWj|A2}31DdtBx2;(@IWxiREv vH0D4{3Li; z0vG{oljH#1-@GX=%-ztU%d7hM(IZXpw38A_%US?chdVHYj!G7T05`Wy!(0P@PKK ztT*kKr~5k_u|LL{rMBCmB=%O@F&3@fCFi90Xc*bq*&q7iBlNb+0t7%oL<3M(Nw6Xp zI|Yo(lMITtFRQ9PA=l&rMKSbk=s)%`;f^>^O38+Ygv8}vS>8NSYa&iXn^JReWDH&6 zI*wxhxrp9Bqt1=tDS6G^Th*CpyZN}%d&}JHbo!nlI}fD_DJb^7&Ibx3P1Aw!KNy>8 z<<|jhYL!CKfR=i0dOGjNPdYI-2Pt04ct76)c;~#s}8^yZXI0W2cj-|LTspt ze@m)z`}(Dw_m6t)T)h~Tm87DiMCho~>`oW7gS51*^Aljvv%(>`5c2 HlL z#%1GlysNv9r)&~?fvzi2+A531m@x?4Y!I`wHjg<)#n^N=#RPhSO6=#PQMDU4Zp09$ ztn2IPah#5;#CCMV%!)M_LJIgFTph}!UFB&Y>p`^j@F;%qVS8ud=TBl{63(|G$rAqW zD{~g&VxJm$Sx;Q>t!aITq79Bh?+(?}SaOQYIfzQ2B(bBNmi(@B&-3g%{9thbH)Og2 znB9~A(}fxxbakH=KMZA*l&mPvR?)L`oSzm7GBsx)4e~BsC~L?49Xfh0mgw8^-h32+ z8A5}Vj2pu(wLf2h%DJ{{3LC?_{RY7Ql9LOHs~nZFS83YreV+GxW_3y0!`H^9Wg}sT z3AO)u!;hV6juYj#G{T48Zm&$UV;Vp#3GN@k566p5Oyu1mLr9Ci8#!zq?u>$Vw)PQG zQPla}pFVxcBcisjw1{yvc*Tb=d&^OAUaayG<8I>~X_$N|VNlX`5q7>!k>T3KNAog1 zdC469VA%NeUC=Bv&wp>wW0EHb?m?%xaEJw%ei4xm*k^=wlZb=hb3_uE_Y3Z&3%e2* za7=E_ytEgQeaVgwKh3Z=S$~k?h)vO#6?{D9vvbD8z)usbT)lcldl*3#2C^{2J=& zjBoN8#l*Z#eGK&K^|mb-KTK4tY^+v_?`3B<7C8>RqBgon$fsR2H-wgs!QH_%fix!z z>&{o1Wyd)X9}9rZ-sTSC=`k2I3((P@w||o0UN{jE_?6XV&+&MA0(mgjmA0!3X9T=u z5;hVl(=^iB7exc4c1KyTD`)-dL@pt}f5n`r6ARyOrro#HBTP Jx!Y_Hwf8FR5+LjU@8_8U!b91Z#H*!R48H~@L)e=7V7BjDJCPf8CRCLIfJ z$)DA`KGs)4;|6dVT_!9nV_nzZEFiG!?$G|+P+fDc=a8q^UXs)?Ih=k%bWsbr8%ZOZ zrCq9+F`Dr>Vfg3mdlM~O@hTg$;9HniwYdZB8Vt$3S;*CYFBt$)E%@>Tf~sD3! C6ycA86jCX9TJofbj9xup&={Ouk2tB2t@$>Tyik$yihU>eg~#Q*_8 z_G _@K7CEe?)^npc`qb!zfRHcv7)&yIVo{}k z?l6cFE9_t!9Q7aswQ@7Ob~KBD$}1Si4hZWv8gIi5d8*pr%H1h7Y+@c!hjW{P9PVKt zw&n>jt^WzMTzJ(@NYB0^&M9{@`(5TEazivW6Qg8)bf{wqc`~)VKJ*5d`*V8BB$#Sz z$5_z;`#V{e$?dQ0GUAU<|LMbEea?9=HF$YERN--J$%EUQLYJWy#ZHo9QQyZGvAdS6 z?w%Ni-5jhhRA8D`Gm28%xoL@BN zGVVlhHUSLskyMSA~e zAWIGUml(4yCM+DH#Kl$tq)7-HoLrlcV6>LBS#x%MK}x;H=~}4a{6^10* SNC{bR?pTYEA zm%(rg+CCuK8ga 4=Huep45`P#QTMFaYAlL9G)~`sTL~LAv_)}o7)GKT zyEu+$tj7?{0Aq8G!*@nu-8nb6@U&XOB-qASk-W^zRC41zpmEYsS}$@LdG{~;TtIsJ z95_sM&6fMqIY8LiIcRTHOX!2a5DInoloR}?Z~Hc{LveMduTp-`b1iBA>7y!xxi)eQ zkD(~vSu@}$Mnr72V7H`F&Ow!?C&JE6@)JJ9)3%yvwz~}%&z)N=-xz$nN7Mt*Rc6B9 zGQ^YI|GB2DVofV-D~|ChVAf(Z8J1?Eq#`*fXr39zR~byBcS;vV=&+yBo1|;5u7N z#RXS0U|&Z0{9xHiT`QbHV}qq?w*loZ q>T0v#QFQf4 )!cv!Y8ST~@HTF+o zO$~~j4qs*PSj^YmwtP_OwO{*G-2>ECde204SI;Cf$tfsc_RXK3^mhWmCU{SfN6CO6 zvccz4$T?^-Udx<3@$7CIJ@w7Ie=8voHd?{~FkX~W3RN68FE5GK2qgZ6+k;(-;EceZ zzBgz|Z#9!}zK$YIU%PFUVSiX^G{VsHC#5;AT|{JLov0DEa=&`>v(!$wqu;C -__|$AtKSPloY5^UhW>*PlCRw&|%x&%00{$ecoyPP>&g*nc zgCVo&MT1#yBL?~Ky5jZjt|aaIl$NJ_+4j0%{LB*7h=~bC-*E%+;(9+Gj-6wu8TwBS zkag+iE9b}Cz-j5ppT%QA1wanpF~fHqcpwGp6M_V7bievR;h+JXKVWKx(hbC`6je%s zU+jYIkj8G+$mf&2hMVYp)l)jO8~E56^jIu 0Ryr7X5sWnF| z-M4#rIIi~qt4ZDuM6}(_wGI!YJz#XtR`nbg6p9ny#ffEC34S%YL}9NYA11N6sp7_L z&8F~lVI28RZ gs8VSxKUw=&FD(4t%W#3fR20Nyr06b zcXIg^a!I0#8%f1JTY0kvvvmYPgWD-GN1?xFvZfY@k_gD_0Va2ewuS68@D~n_NGM9n z(jT<*aC1|%n?iZ#nG0fR^lC6V o1*$wm^*R04yuomWm3=pDGI$;BD`? zTjw}Cgg%tBh@=zTsl4a8xUTwWrl}@8Dt|Bl8?L2`;{r8+BP){Y-Dc?J@Aa|Gx?F86 z?0rbz?JBS{FG8C-^;avhxEanneu2oSsHC_!^(S6$!88$xxRRiWrd+#;X~^S>*lEUY z)=N?4_@PxsoiWO8FjbdJMS%0CEWv-+lr-KmxBT=iGFONGp7mOeq6QkK#G-j+3 z8}-dSJj!t}+6WmJl-xQ6$?VMz5%Od@PJLT={&0XMkbMTiuYs3KYi%nNEh)=fSIfVd z9X_0{JLEIohe4S|L5?U+OHWUwcpK^0e+UFGBU6y_C^5MvT%H{@!#TjrZD*yE1cI0u zDwZ)^qZ#1`A2t90TH4SX2lrZ}-S&(2P@*_ClajvfUY1g_F||>I4$oYhmfGVpf7v55 zR pPbK6s*X z?ek2I$_o_Gla2)iH6~)}8QuY{7y70}NW!vNXXfFt`A*g8Syzr!#f2C%O_Sq-Vpjd` z`waUd3^84~l;w&a)-2pz` kf|a;2V#d|Krmq>7AS> z)`l>PK?=fYO2AL5n Ap?*ABsLm$fj^N eH4Ir#0p|=_Z;`p(f=p>Z&Ii9pt><)UjS?9~#PV+()bi!!M;Vq}WurXCq3QHde z4Gqmp>EW80nHdd#84_^W2+SeW G6upZ3RvT?YhLQ*;EbfS#P|`||ReKGMc+ZaV81yk%dD z?Q$gzNpi6-chj|?!VX!VCXj+K8&-el7&)u7_yBA>P}RVsyU{76*}3+hnJ|c%c+_a@ z+xb>v-BVzuf+U27<{e{TEOZnP&$RAEY;3{D+*jOdY+#NnxfH6l`+UtNxe~imS-DrS z;{CoPlQ5XKC)fDOO-s-gULf5DM-GgtrCNy5aT;n{4L@dfoYTvR&hlCsrdfw-w(Ft} zMTED8q>16nzE)hCDAWxn)0Lo%0)FAay;yzflVQ`=>4Ro^%s EX@p>i~)nT!b?a zqv^QZ|M5<`qgcPYAE|=8!e}c1rJO&t@vmN0%JOuY(czN9Znc HbRse8Jhxg6-l(suhdf@OrrZ`NeMj)Eb%E@8;Fj)ri3ql}3~p vYdN3RO2b01yr`PvI%4ahtI#det3%*?GQgYe~RF6C)rs}e}Wfv z;R1kdo%+qQOpk)c1r14@5KG^Kh)B+>QzT*GPqky7a{cTeZftFOmyLWR_?zMDky z5D@Pjx}x_Wr3WDJFTOxFNO~I@8$NZ205}HUv!2qzI8d@fXn~l*6`p=oo M+ L8 zBnxTjPdYz_pIv*pwP_X=I%~Oew5`QwT!6huCA1H0-
+%!;k z)G-#e5DB_SuMJnl1S=I4;>hbj EjDY zUd%{Ktj0A{@&2&a31fV6W?c}Gh7#R(SAVDXD*>8$La=QAsDquf@u48DrKACIMU4Dv zB=FToo$ejM&&kQ=2%l^QS_twUS*Dtrn!^MMoWvu(t2>Af#)68rQ+yKP2#e7r0rgE& z(N%--#c9XQ#)Z9in#I2_VBx4uZfR*zBQdkWLuF1@o^onxO4tC1-bx|^GIBD}dv=vX z+rYLTt@m#TQy594G@7I=L;UWBmh2AzH|nF+EIE4;qH?=$V-|qc_-o^bWe@wG_fFXc zHv*cICqt}VdhQWnm`etLq)E#~oBKVM(wvME)}z-Dk4;k2@QAE2#r)Q-VMHoYc)yCL zq_k8XJ2rNJ7QQ10x-&++BekVJNxYPkXGIEi#c#SBK*IcwE3|KBb0_yTVKCLXvf~;U z4!K4FKUymvYHbS1{cf|%_nVvi3z0(3#s?C>p2J%gVH1dEa^buZmSA_ |1zmNG;Ho~y zN}e4?M1TT|?8PXkN=is8bp;glQ^17>C;%j^6cmrA0B4So*X9=!JB>dLW^ IltUvejY}69F%^y7=C;{cz(c{3)AqAk|dr1vZf0-Dj;MV z191Ki-*RDatgJ|ZnIEIlJkSLUFVkNd`{ZrN-J!+C$`;nb3X@URWhB8EHZ^j(UpoB- zJ%um~0|V!(e|T{1dN*!_ZMi_hA9|>HA-#}QhvQz;G{}na&aI0Wv#(^_hA&EF$@CRT zZ)>YDZ2hvePW`;o%hwKG{RCe19Nbyti$Cw|0rV!#B K)gFcBFa%;Dg8^$G3GY}Tv{uwrV-v`Kou=zuQr}+1;fh)T(@LIu`{aKMlQxk-t z7@CO79BGB)`rN`w(3Dd8NF7q YnGSR~JfhL#Z`e}R zg%VB@=-b?hEL-1@Z7ae_IdV02ZB832?eQBZ0;81(n|aW@0D9)ABG!C+d%b+2Y?#*L z>K%cAp{OWYekDnGRUxOs;iVNc`QHcKUTrYWAYzPhc 89_BE QS zN%bKGZkXh1>m`^W@T&;Da{Bua^r339ZUhhAx8}kTAknC(O^%7bDksA)BH}cBmuw3~ z0&@!ZY?jBQcW&M=pje~Y+$>N4Pk}dKf}8qD{Vmdv$;$xcqs4chXlQHQJ?6F&Qtv>H zbkT+|`bbTU7r^a$a`*BD1V&Le0f(is=3+AwGvF2k*M#fsuMtW(X)W0Iku;oz5c?OG zL>sg55_rrTev4y&pgiCZqtKKrjO%t^imj|vzVni{f%fZ|z6DVH*7UzDDeX%4p4D9ZAEJl`}EN{h#q0p^$gzJXN0?jFIgFa _vP4zGwS7o-~JN#c@X`#z;B^-lYaytXbdx9 zaPWO6C)rDv+P2ukXMy^BNF%!-RbJ13z7q`VG|XTE@-53B-bmHipNO|cHt;U5pwLT( zQF~n%AHW8nrX$?e9F5($R-$*pqs-vN>+JJ9{&~!8?+?U509# lewTkFVYck{QT(%@MB^w@WI>foDh5hUg73 ~2sp(MB$dL1mlHY=R*_6VF{&i4A lOmDtY4C0bk`QSg_dxXrgq$?Ru-R;#KkT(Bs*f zpd+a-27CME)<0F9Iz_KTH8(4P3ods z^d2 zv#0nrmmEySplE8EJR9nB3@TZKhP?Xv9~$yM ?CExaU!ud)w>1lO;R^s_8jNdvNE$1`uogow zw`s*q6|Oly7G|5v%*pjrtT 4GndEbDtw66};;-`93PglIE@=P!-K?A%U-W9}b0Ym2p0#T|)Zy*t zHOH7$F+;nv(l7 yc_q%*>`^<0|ISe0|^r{W&jq zN3CD$%qK^Q9qwL=_Lm7^Uw7}qgt;EBXT)-nGmT3yW6%l6zDn Sc!2!mK zkej2B{Ky+P)GiLQnnsFOu2l6*HmLKL_o>3DM>DgrX-wF8me;(Lx%l{KU8!d)USv_A zlXAZPCdWvZO2kyzMG)r@O$U5d5!S7$iuZ%F2R>85&dM`8kK5o6{%H-~09$rGOCjX^ z+-$@EIX|LvRBHV5sQbn2wh71wQz)ScUL$)B?#8Q _&GvdqhERi`u5sA{8Mr;aKqVo2`ZWBr5IGmBbbNd zAay0An$Azu3R<168>5lbE^QDx>%U8_Yw;oB&w!4Ac&NOs^SeRZb6OS)zd{iB0g {;leOn39laf#3)mJieReNOEO>J-pAgyNrNOMTcWCy$NpCR4~-gY z>6)^+w8Kl<20VCN-3pRt6h78f_;!=&dGQ_ObdZW6LP$1ai#v =w+qrC+kGb &(CY?WOXjxR0X88w;qBoc4 z{l})n!ejwwGE-42AgBKr7EIuTVmZw0>)Wp #kLU6}L;Zfm z`V;IRoNzZ`;z^f37tzh`tF?pM%X~o*re4vryR1<4df;9 jSYoz_LCYDH-r>O}NtZ(k) znN#yFR!pmssqmdQow;QPj&OMmmJ#9MpASz}cwOU2n|E q| znf|#)=dQmtI8mo{P-6U3605Km=B-)jUtD?l7?)DdT3@$T5Sc<^YpcZ;?Ls 0 zg5~5H@vDb%y?MdM?(rj9QSJ)3!Oq#dIVT`Mj( 8z&UaU~HzI+`xE$Kdoc{5KR;T ztnPQO+!+u{#G`l`#NTE8i4M}L{}6?q8z)P%&pPlBl8TumlskfL5Wn!bmyO&Rh^Tn< z(D=fX FwwUx)l8^8dh01q(s-y zaDR+NksdwE|5I?dUTWXjMG2Ysm4Rh7p*=t?{>6TN`Ofur*4XOLf9sE)&ue|iNM O~Q!<>rkX&C<|&w3>J z2jDM+$4Df-*Jq>rR#|)qK~VfyBF4gJ?DP<-??F8zRS4?ufe18px+xK?;9uFZSRpI# z-hhEsUirrH}XDq}nIAIoR7@ z!%PiQq#Y((DNMqv(V}3cH9eYTv3qR#Rr*5&dv#M0H91lZRs{g8sEOl&hV^#DU(A;R zGyVSa;jXb0KPbK_5m!o &us3eACa>IsGM?jv@&f1~w6--gk#boR; z55N#*fQO=yfWL}>PRz8HeW>)46lc_SSe|`-cjfE@=}in*&alxfS@A>Smo8^m0>*Hl zMz1e?Y0j+Z&3pge{<4w)B%5UhuM _E;TDuc47jR*Y&AVcJv^qf86zXT-3UeS=s|IOYkn|$ROHD82AE^cwKW#90nWXmD z%YV4^iQN6Zx+*SC%uCv<2`647bRPR7U`m*m^bPvb^K*)Ugh*1+wj2-DnCt-9Z}zh& zMzXWUl&Ap95`_&Mvc+MVwdQ)#eg%zI)(j2~jy-jDz(GdlFo`YzXG&WG6u7MFn4)m@ zknfEr2qGLKz@9sUd-+1GyUYT7hGT*+Z(li?R$P+DEQy1(;&(Xg&*7MWT7>)dAxPzN ztor^JrOJ^7q9CW&rQD3e=`XamcQ#zyxgM)#0uRj}ml7;7+{W%hlc9-{zy#!6EuGJ( z#O9CZPAmlUI>NdZjBFc0H)LI_7cv)L*=RiU%;TMZ)27Xj|HV(qIUeL>@#E3lO~O78 zK7Vh-EF1TlzE(<%So##sNb0l=ZE$ijb80^qCGwrjRtf~2ZhWsJe}GMM*h3T>1ig*X zR_~Y}9~!^`j-fy>;uzYgiUR6PrNuvpLD>aVb!^7AR@<^Ru2p)D?s!5+Cm@8hsUaZ* zsr@o4)CBd}Mb$8GNE1`8puA37CJ^+t{JM^XecnDOrIFVKCF0NDikJTZNt6krs`ZJQ zt}JPo-aZD}nV5{T1s5bG=@q|xwA6y12e$|I$iOaxXaaI^XPQ<$_xnx-pdad7FR29; zF=HvG1o9+V|8KrQ0q*#+oBK_XSs#oBBa^`k>zS|S=k>yiw9L0~Q~Sx9P7?4?5k-Pz zJ5bYs8TS6 FX ^7Lsg$$-pE&K{X5AXu%*x~c}E=2Uu&^-mGpJOgqW zAgx1|Cgb^RVQBXKo&?i^nup@V#CXlp$J~-k?qv@0Y@wt{Fwe(GIE2$9OAo4-WP0eI zHjrRW!`Ypvy8`$^DaQd_VFqX+h^hiYI{U&(r)<9SfJS3{`FJK%Go_XjZIXfjHlW!V z04pbyQdoHG?F}wf+~g7)S4#q2eSIx6`#1Ygh35q>E5PFif0rVbb#Vi$4JasL5B!0; z?N-S=Y;sRfQpUzTyxSj}bmt _$08ei#pAr@rmFKfNT(d-23z z+}tNo`v}%ijvB4L6K!8`%Cgw4i?ta$$@l`$x#tr{V{1#LmlRgWn@q1V4idz>Pj_(S z7D%r;1HaPL{m`Ink5M)}6M~tD71o~3aWT&vES55HaakQ5O`d#_VWO)`O+xZKO>-WZ zGUD(ZMLOVv^H{KH;rNKnj+b|~y-PSEs9o1E1iWmmEJ|X%JIQCwnHA59RkY$uXejKw zp!8R-m~wk2_tgnnUPvaqhzT;jDAdMnkkCWut#B-nhIKU8{>uOb9j`j(9wvOaif-se zjPc6L?i=id45aVHP>!!=b98>&GO3{2Rro3}lWu9jZ*g&cuyTQy0q~DwzN9DIs1Ng1 zliWH(pY&InX(1@?Gb277;)s5$)9+*I<#mI@{a##O+)Oj-otuf@V6SvhvD@79+3YeB z_>%eMg{S7{6O-cG0jgopzBs5uRwQ@)PI5B#sL~dFjDE5G*Lw&ex3(1D;xhw2-rbkA zAAyqlRh5iY=~9d+xQDTA0ympwV7a-A5x6kV9`t6oW_1$=T@D=01=#b3Ub_qPi{S=f z1# @XUt+C6#MKYA{ zI%qV2oH05Dm@^o>19lD8Y928Bh?|4Axupd_&Ql*$mq13pgk5 8gj;2)SOqepLY;gvNAB1%@Gnp@%}z zwc+YsiqJr)62e7}8>)5{^*d9K&Yvn!oVlJ$NdWUrWyC*Kf^e$ PGcrm z_`!g$QI1{(FpjktNgEqw%VrHhKhrs oWHdOJt&c*T7k<(OWDXs_(X-wr VnaOC{lpmfn0Tf~FYi2ZY zk(xM7+9`=Xx5+^SPzG9m5!Lh0vT7lo6Aox)c08itqeEl^^T5&<=DsTF!_=ufLnjQa zvY6ZQ&NPdO?)N`grS398PW(!@D&n-zt>B-a6@IQywlxlcTFd!1X*48;+)a&*6rL?{ zY0pioZEnd11%YlP8fY}&p67~6$|DGR{<> 5o6~ldZVnDJf<~_qz@UAp7Y(K1Y z1iTD5exjIK>ydyLfn`063}^dwdkRO8QRLa=WJP>&{_lflgVEqwAQ!sXl<58<1qHTf zAf@qc<@RrVZS925e<*xC2)F0s}|G6Y4|_`OY(W3E|29>M!kAP0%mr z2^QyDEx_FnudR7-R)l)NtVY4z`}gs3Np46`L}kptdPUmve@SapD%~@KaXSc@8#+Z7 zn{sahuR0=LXvt8R4JNc Ha*oAYrZai?8suU%^UvsBrbyNoAIYid}|EYRzHdQ?1t z1ei$n>RoB+>Glfn(Mz`kvI)elQ5Oap0Ugvxt61}z2s2`h02vf;RD3y!=d&h{7?-ap z(12s@oW1?!2a_mpNo{jH$0`bTy>s19qJ1z;_ufEV$G_9q++3ItJE$UZ!^4p==%a{; zFyrK**xQbjxQB7*7mmLQ1sZqq73gApalioU)g?Ht{1wq9^Al+}#sDH-apt%1u4{}M zZ2D8tloJZaBr`2DRoE$McXxMWD@|wQu@!*T(80x3`bxjMySnzMDleqd@|209Mums1 z&akp^{v13by1raGbZRtee==%+RSgSnE-}Yi_e|?}d(EKbGyA5D2v0W6`bE=edYbd| z{dd!fvfNu0C| Fvd` `TqsENH5EL z4R5bGBFM=}SfKqZkrUHa{T60b)+rbpGm;Q&3z|n~*yrvvqrjxFdoSOzg)>Kb;HB_p z0L+;f>NWF(8r8ca1Py`?`Sq(xN+<4FeSoyE(XWC95O9@Msk*2L7zo zbQ*x1&lxQ~SNOzg=%Wh<$>`Y@WHmGl-zlq6)*2s 4;Ngt)L3_!bp(&9ssOxM(T?VEKzplofW <_4o~ejNJbC#B{1ncd78U(xmHQU zH85^Yzi*TZYT~?>W#*Fvb0@B-1xa)qW7%KpueivFo}CzjhQ1bi3CC_c)MCJ%2(j|M zdfe;+!>@rIl%@LNKaO_{85uYaJA_r?u>S-jRd!^sSi}kO3bSLKHD+Yg8Ue3^U$RS* zv5=4ETH>1MTAn4jcWOaB3N;xccIkzvE1{w>$Z0sz-cemjD`0Zqmz3n);pPkQUnCjT z5~{dmo_igIRjsc*Rpe&8n*wXo@jC3?ao^`;Wb(&ErN@q>fI&W@sP1cmMD|WAIH{T6 zw@E`*fc9FW$>K>Yq&-*AMhx`fS|QA163~-mrlkRPn|Q4m98G4Lek%wH??3}cp0{NP zOwye7OKC9xKQd@q$*;{!PCoQ~LR}?YoJzJ>l4D5e?)7Q4;rcD^)YV8>8a{Ff6AKHr zWhO7#2LbF$^FIIF%gwlS$~)? V z^1F2{SjI7i0r0))n92Sycx`AsIArh)WJEAoplv!&nV*=Llys=STb(UqJp`Ald6}s> zRL*ke(A8QSj|XE|=$bK&b1+^qJ<3gqJ}N;;07|gxSL~c{MG59)I-EqYh2u;PYA}Z= zCq_%igO}@TFKTZETsM+_%5Um0!i32NOrso*UH)K)t*iOHA?#}y2Xon1vCSlTg<=aO zk-2o#dj(=HycYM>z^3N(5OscX9&l)W_PRG^(|mZHJg{(x`hpHQSHA}Lng=88>C>N` z!D3SrGw^1UcXD9 U UIz~c zM@ClqCup2_&6G359m==^pp@j>wiGvacW{RTKSM}H3c2FISWV?}J_4Y!FFWt#G{l<* zf~`MBd@V7|F(S+{DC~}}0!2;b=uWcwE%Dbh;waYS*2})vr j~rkQS}LE{7imV P& z=0y(r(@d@Fx%jvyGJcYgxLzky`$NOp?|K&_0<{OtI?zP5fNewYfGRz90m zcC;N|mSuL$#T~k^<`x!SKh+G x_Ujj0s+~>j>1RBYg7-uIZKgnFadQ}TUN1^zg1F7IZ z&(^&c7Dkx*qNEp=vva+uc6m7uxd(FmY77RG|3yDU-&$ju11M7kiY~p7Xtp}ebXp7i zx~_@N^RfmLFmyP7ZH+@$SEt6g+i7;3yw3X5+|wW#9Jw!p^N mlp;mL$AU;b+QDS^?Bz7+ETrD8{wD(s9%;s zmBcQA<^(qnZ?Eb#-zI_1(}zH9&GmD`WupDA>PLVYC?_(+J@n{|E(&TJ4rV}m&`oDB zVBTN7@$}*;FgOp%3G~{Sd!K#^HkgZJSWc|YHJwrZZBMnEtKq-t5N{vM8~65TnM#60 z(?VI(xAohj?6 TV%|j{ee~`ZR^)M&@;J)0ONWTa zK&mHC;16duMpLQa5ipn1y48&!_QGhCk2}-X%%iU&L{8_r`~FV&iX~rbd?=!gU0ts| z%ifLz|GHsD(2}$NlwG>n%HTZ-S&RA$MmoAv8n9cNzBCsWN?{QZti2}@%~9E)Pf$$l zhO c|tVuc_W!+iEvfLIr9 T7iCfO$AXWcw3FYK5)Vd7K)k%Hp{FqtcIQ4Kbda`6vm9lv z$?8aCX!Cp=$3^orXbVIRO!qGz=8SrpfY)E;?|;~58@QOr?t)R}$y5IBgPw}(28RBY zI3tk1JJ%BkVby~)b%8UM(MJKc7RO_6EB<10h|O7ZozVIQJVNv8X>;b@fw^W>B}`V@ z%24y2uH1c!O&&oopKqhgl6&C20AI T?Ji_kj!MOoX!OA_X_NM$LlTniZQg=WIF z>O*Q!inzGDn{A}{uFVsHpFadmtB$5HcHX#W&!Q!s>JGn=?rv=btNVA!$vATcdU8_( za`vA4TbH&KL9S9?o3`6}iIn3WT<}ZiG3Ucxiu-3_Yyx~aQrfv$?AIEV!ke&{&M%?; zHk>R_EeSa~%gC@q?Mujvq9&InsF7I>+A>^w{$w04G<@rviWV-6zzfY-i7UDm|Q6 z`ghtMRQ3ndarAC=QRSX}>WgnTVM|Q!_P@~dP&-8cQ2PF>FlalO>cC7SJ28_B;U`Dw zDh1?}ln5g~!F)Z ?wW-D4KZF>d?#1d3XG({#$HhK-_ToiW3=GKw z<3}(L+3?;fJs$MWRSygss}?{BrbtUelXso_q^}z#(LY!1i- AH(|DRI z1mr^3`VNC%rP4>f>jodz#J<72ZY0vBH(a*O?<<)Lp3N!4d3ppf^@9m3kIBiEL{itS z{^L@PT=XjAi#T-nvxlO*3P}KW!AAvf@j ou Xoa!cNFv0p;)$ sN z_X}2K@T}e6zu)uiTE`ys$j!(oLRaV;a3e2N#DB4F+Qu3gi*&+~MTSn}z_Vur5g(43 zp3(Z{3t4vWmoH>vP)De$sxpY-V5|`5Ph|LBKm%zRCVlS>o8wKrBOZ%y1Am42D*UiV zx<rac+a?i_+bWynv3F=Z zM#PR zydd(0Y-~6&oEGcjbxCyN3!rBYSJe4`fzFa*VPRm6(?FeA0sO@NTC|k7G>Hh6DoZZi z#-IBTg9N#R=tt`fZmO6alHj1Avh{zZoze6Cd4@(EXAU!DJv6|;G*1#u&$WcK+7^_S z8npu_|8w{bf$oR{5>VxbfNBZ672(EAv;s{lq-yp2J3IdHrMpnFDDl8UkFsl?`a3Ir z00>=ezL&GY!a;;K>QQDL4z3juMyjsSsY7KYCA*$^#uZRp)S4G1Eto~ +8+rUiF(GoRX3P>?fR(P}?O!e|H1aa9-7V0xIP+JT##c-HvO#7!)R2 zeluK@wEyLOosMt<2a8khy2Xc+ |?;9ylNC v0N zK+HiwOI SxsELS|C zhz0ooPz4|#0I1o2AAqEj?Uh9RpYHz`)m1(5ECW9DrDhnYL>h;a;VYAA^-pS-Am^I9 z{%GCj$~fx`Bqe`+h8xT&WBb}PkosxV4hg=m4}Y(T_$^{$_XVjTDFs;uk&Z*qzMY{( z2cz>!Kp%pemqXTlF^3e*ez$OhD0^Hv)@>1La8HOR9aU6_iLxynJLG=HICo1+Ci3qe z@2>kF?;IZ?DP|c5V{siK??dF`hl}dw76_8&jW1g~%WR9uJ *U^e3YPIyo6<79H$Oj<&imHh7w z=Hsakf4@JS4PgxVaC}%<>ecm7aTPcy8E y2J8e)F6KkO4asb*>!22MmE`5=> xN`})<%4_VTYw~I22`8i5jB4z=ADgFvUcrOKju6>I;TUl zc3d2+D<7!ZC@f>FU8Eq5Q}5LKh{q>>DIAI#!QjTFn`9zl&K$g-@-Z^rF8a}(UQs1= zj0Cw(24myuDvo<$l$gE7ocgsPaAQ nH*Oj3|M0r%;Po84R}wJaUpMDF(5A$DfwuyfU0e-rvm z78VvF1;WVQZb>NM5Uu%Y5$l2PL-|i~6YHA)i+v?wJs-TU$Xz0Iv&LPsJ($OjY?5AR z9VtjVv%*5oc>B}OF;Qmw+k6Z2X{}P$npzA{O?CvN4*2{BkP?nP-i_Uy`of=SV3_sl zzPlZ9U|B)IiSYfdDHnthrq0Ito5)brRGiM|Y322BbT>$C@@@q5!tf?&e8a|O{B8jV zX$UcBc2ye}X?Gn%@zwu!7WBeaaM!LqmWa{U$admdxhJeC@3a4{b-W c^>r|4ym)MrfU|)>8BoZ$Lty06U)u(7-mSJPfZK8V zsSnI3=?d{r?0(q$f%H41dj!D%x!NeQYGCrXor}9Hzw~Ge%V&yTscO3!6zLXImXAfZ zBwPEG{i7Zc0d#u^$Wr;}r(o$Y=NVP576#n@zK1;>xR?DPN-kw$WhHw3%8X2d1a0Jc z5Mq(maq28WcnTl2fJmwV8QZ6(5P$2mYMh*#vs>8>U!QpKzUDRgmC5n`mx(u1OilG| zo$mbwj0iG;AeRobdO$Db8>Pysjmi6Ls%E}dixJ-bk)*Z|?$$nq ^WTfqlv%V z8rZEbt}dUFl!tT+DBe~)b>V7W+XikR+ct>kA+siGUSka%dv;_}FXHlH5=#zr>+z?> zo%Bf!G#;5F{DF(0ETlpIKu_k;eaN1}!70zTO@vO3k}5UJ-^$)(SISF!uIjru<0b)I zYHn#f&&z~Y@047os@Y0MmEnJXC7Gv5U!1Wke)bY-wFI_|N_Tt}32+?3rkEXX;E67q zU2FjtPG7LSqgcC30m_MRCJVQj12>L#zPAl~0r4wT9oZY_^{1!;?{tjuvznq7F|wqJ z(=#`4;o!8sVj^O__6{nPt?wiRbM%He-My|%NB9s)^d9`}6108FMoXK8%m8To3~!D; zXM>&H_^mhaT=m>MOlBUqM u2akQ@Et2;@>zx?1t%i?5KiWY_gQ}X)X6}o$QI+%g zMHQ=KD}Pab=9WH?1 5wQ`3^=P-_ai=jpyc~8aE0?@VakT-W~ z`C=`u3niQXyegXn#<6pM>Okic240HuuxSJG6J!nk{r6tcYNrCQOCuwH7)T3ohg-C* zda}S4Ls_1<{|bt_#Kn7SO(HaM9>DYaoabSA+u9#_FH)t1$Y#x#%|L7jvJ)~@bj8C= z9&+>G1GHYS8TVqlfBT|-^_m75P=qF7iXM&)jz%_iR~sWvM}9upm5dW0SIVgM(JK*P zI*CKZ@4-b3!Nx#C&?HO@OxEXoS6hI=h;O<({+@sCd!gTJ+$`!`E6Gy6Z4Mhd8?p_Y zkEy1qZ^S(}dS^&!38b;`h|vnJL$KTm3JN`z_U>VQ0jDm+jDsU;DWr3DtABsr+n!>1 z0!en2R!Ce-aw~oKVZXHzLF1&=_46&Sjwm8gdGN3-2~n`gFFWMb)EKimW413rDP0uB znfTMWx!fS5#ZPvw5u|;6Zx!U$R|;L&bJx=92dKQ`%eN!_hGW>z=*8lkbZ#m@)`C3S zHEIVHemRe^(<)bXkiHBWIZY%{?~dB@1DLsd2(lVro3A)dY$Ex)v^s-MT3C{kLKir% zK9nJV^O6E-v0%uS();VH-33m#FY4;~v`yT)3qSIHze_?r#h|*VX0)04M5_CNtcuT8 zb7sc2+@}52w^g`0pfk2{c)d^pK@=(*DD9|ksazaqCHIkcg?&0^>O!|=V0Cb5VtAP4 zkE}RK5q3!1wd>V5Iq@6dtAhf1*J$@)ET@79;^ppO(V&4WgIi&_HuCN1Q&5jf6o00l z>4tN7RrQJwBTiQiAfjLH=I%>-PXx6usNGqvaDpotGu?%mvEQ8WjEV>xUD(HOc@=nm zzbSshwLJe#;hQq`_G{m+0JC}rdNIk?y6*loAt@%{?oGHsU1VMhi%^BkNM_wea&G{r zZ d0}`)plClEx3ui0wxB(7f~SBOd2gY2`ek>i~)+T`4LEn zzt$>nk-c)U@B7fhrEhs4Tu@Up{OlQ!TW%~ahQACMPU~#h^12HUy8(Nd?SbTlx-DTU z)zy@!Qge|dN0C)$;g3z8Pv{eGVFsI@8*}=HHfifTX DwBQ1|eVp6~X`~ 7cS3fFnbB2%PeMvc_OI6XNQJtI z32So0I3Ro_r7uNC$IOr #EM; zB1=^hJQOx}^tNYzc>W*YaV#IA 7*>Kz>9P5<2tAZ}E3;%vcK^QNx5v14P zkXU9;xc108H@+#xg0XM=Cs41fw$|k4&EV&eb}V#st_5}8&NxHX&(Qcc`EUY~j~YN_ z@;MOP1qg71l{yT1Eo~hw-8hNBclQ_cqn^FoAI0QC{ !S^=T2=`?q`4LU{fzFJBx(F7mZ%MuSFHzQB?0O5^Gc47!|SV-!!tDmMoTs- zwqsn3_ZFm^Sm`;0IwAa=-a$=8{V7u>H(Y7VqrIbJ( vz=`2AmRIj8P*PH2t(C81l!a6L#{LN!G7PFbC zudy`I5*HY#oJLRF+}w<0vR79dDzj2hH?yd|AAY)3^@XfgO{7dQ!Z_U0EfBIoSd)81 zO6jW<3*Na_<3AJgq~xGdwdw)jXKFu?+Cc{`J*P?F$1%n%p~_|xmzOH 2)r5ibS>Vu4L*SE`Ff!A2K4ail>%b!ctOB2)5 z(0Ef^4E5TEq<-KaGuQ@v5DRQpn)|T0c&C$S(s!$?xZ?)B3vuQOZtm@vmW#3pneL6@ zgTl}D)pVg$Qt0I=RTZytDo>NzVT!*hB4`nv`e<_;yUD!^?G%Q|9Ys@@RBN%Fmy)8J zGsE4?pb3si(#abl6+uJ4$4GH{B2_^OGBs4SD^j|S$|Ecdd~YLz-iq%gkX4~lCKf(O zB}O;f#ACja&PSKU!b-S}IN23Y`SN^_rg!w`F*)ai%{F%=Db^jl%MMLvsZ>Rep^GP- zUWa{1hnG1Z_|oOP^m#ufCMJT*i3t4(U@>|9WUcoD6bY1rKL?^vC>GI;#3p*)kd@@( z jBYeX;Kp`F();!^h |H~0sS%QrM>U!Qc5U>%)e(#($1!=TxN z72=QoZ6<}$OI#cYScP`2XI0nUl7y)t^Ab (w&hP^W1DWPQnQE9 zDuw@`JvLKRghj{?-=u&!)eCy(*zw0>JeVm^Z(}-)6Xd3#Ooo_@jE4ghJb_`XB&WgG z)6+OrL{9^=EzvuW`%&T1{yTDH{=n4(8xDGMoq 4v1cZV7FwVH zeoJ86uCD45M&~8o6Ty X75)y9e>*K$Z(JPUooUiR*?`$l*a4`fW zio875V;&hadePj_(86hTc=5HwbXKDNo213py%xEBV11cnXpe!@8gBO-pOka*DCog> z%w584VNFd P-N8MFvP{=UYW4Rk8^Ub#Zm#r5j5@K zX4>n(LLXGfv{=chka`rL9F>-oOq2I9)VdT$Bfj;yv$Mg#LjJN?S>4+lEj<`aotk^? zqG#_>&ROcyr|G^gT=qM>oB`1p(1+(CYhQao`|H?c-7)Rk zu1GFLJq-6XEbR@wU_!-UPBdL*y2DS&=s#C~(SQ0jIX0%Az~be%^`t`R%)#24% dhRPgZ_Co^jjoh zltf|QQcGyo?ip+H=xev+&NnqSIwjtPF^0%p&4@ucPD6$Pb_nbTfh0Ur&J{X-8g0It z_g%#kz8L7}nvJp{O}T{5M9@bn-K!X3pa{Jm0PJe#_%@)xECt;Zt)vBQ#a-`7xD-ND z+f*>o9n-h8vb?N3M0gu5%aW?)P-(1`O)DyrhmPELp=S=tG0>6afGWf;E25k WDUy0IUzvrpcQ$L4(&kLl(I2x~a~)p# zjWDqgJZH|5A`x#&NAJGp_*$iBgE89N!lLrG7K7>yT`bxe?TU7P&n${?kofojc~8jE z^2$p8y-LHvOSI&`Y@Tf`rNNJEvXL2)a{R40B=Q{blM_J Um(^fwaN|8PSlp=fEf4;|6rN_?pR-6?$Gj$MIf8f%%eZY1#SMC z%SIiti-?LosO31Eo}1&9<)5Y@FhxhG }Bq}Mo0 z!LX9qw*) WuIxw~1!mm=v-nm-PCym4l8xwI)G zGxLXAt+CL;*4}ULM$bi(PMeCtk18PO@WsyUG6zy_J$-*|ps+6?A$ct>=QTfu_rOvl zMT x7y0eODw9_j|Mk@(RQY;4;TQ+&*I9PnmFYz$ZEPrZv|q#Z_PsesJGwb zx5xSn+pQtQDr{?xE>V4tU|BFi316o6Tk gL zbRnF-zrSBL3?|R>E+$@0g*(2>eSdfHHO{s78^KM37v#;c44>E6Z%d>e95>WVe0~Lx zlb2sMHtxeJU2X~s!K+*>77MJNI$d5iw)%-I{)w78FW L;wlBzr_6afV!$gznV>Ng8m77$~Am<>%e;A|?C@9n)fL^A#HN zKYzZ-OL+DAx(Knb63B$Y{PR4s5E%;Ta`?Ke&6m)7N9CHN{N#~+y8mw<`UpEvj4K^7 zD>gK~`jVcW;%?~N{$?BJo|Xb)-hHFK(mp-7j=_tAB?Yq$W)sT7T3c0r*Drd*Q|MrS z|0|QU&y0gVr5Go&6PjRM`4X~IGoC;Hb#Nq>_RJXmn3}f#*m9w!>FMb%Svp_P!xO?p zSQOWj-?s$*SPEc8?=Uz!`TBA!INiUmS2)1^RJ?a&1f#@kGcM_~o%O!y<=O0#zNHNe zu|f`WM1VW?6wjqNF5mDIo0P#T88UY~v#{{WurR4(-drR_Bab^?jf?hKMut?X)}2a5 ziq? v?8J6N+CB7Cl_gwl;b4a+jqWsJnooM+> z{Q{LTlsVlyatmTHyr+u`4J!J>0?i(BcE44@MEXaA@RsBt0?*aRGZ1K&ebu*QEW|_+ z3N{Hk-sNzSEqX?xf1#F*z XS=B>{QsE_ZL9r1nPd?S+kvw}$_^vh_26HOUJ~k5$(5K9X%}qA85!2_>3y zt)E^Fii+CsxL4_Y=Z;$8fa#|;AX)@Wf9t9=fTddKnBz~eWpB`lpFP79GPOxbN)$wq z9ePUn8ETl;ja)8-Gl(`jJV3=*iZBs5t09M{*_3__cpm>f(p5MzC^qbvce|R|V<$0F zBGqh7d(=3#ijz$@EcCFmQyFh4FD<2coV`TxoV@VsOYxmqUmG!NDLaW&!`$BGZT$iU z>X#6EG_?CR{lRo&zN|pDcyM|9FD)hMp > zcu1OdTk1vVG^WYxfgjFLUcui&=zO9Ee_LB+R8(PyW=(wcz#F|N`+Y5XK|&T90~z^N zYmv$Kt-~?x_LS`S09{=jF$YRNc*JkdcNk07RWzoLjN+Y}{M{ouLYppAXlBP1`*n%o zI$BzYC!x6Q@x?ahmXhF q}V*&fdp zZiGTjd)UX>O|Da7=Yg#$@B`M9ON~~5oW-(E{luh;Dk-bIm$yDW+352Sx 6rKM|aiFr{qpv9VlJI|tMLzE{EMf?E5wr`&y&?Tbz& zATJH65nzf9`)E@==HxN;92i4J9N}ogMu1!LT5hi%S;;^7ea%*Pe_9iYir;LAYAqAH z8bmWQyG!qzw5#JQV;IP2<6(?3j3%&VVWE5+5@_Se!uvRq;MOpUoD>0jP*?Lco4&HZ zVJ$}?sXY)fG#3GyDO)4!FzagTAJ)t^*C!??rEnZ@jFrf}=%H!>+4;0eYA;~%PJSXs yIZ;rdPJWLHd54h Rqf(+|9rQo`?tqg3 literal 0 HcmV?d00001 diff --git a/IBPSA/Resources/Images/Fluid/Geothermal/Aquifer/Geometry.svg b/IBPSA/Resources/Images/Fluid/Geothermal/Aquifer/Geometry.svg new file mode 100644 index 0000000000..4829707cac --- /dev/null +++ b/IBPSA/Resources/Images/Fluid/Geothermal/Aquifer/Geometry.svg @@ -0,0 +1 @@ + \ No newline at end of file From 0b83b17eef3455040af2ba12954ca6d7d6329500 Mon Sep 17 00:00:00 2001 From: Alessandro Maccarini Date: Wed, 24 May 2023 13:12:10 +0200 Subject: [PATCH 03/85] Added aquifer properties as record --- IBPSA/Fluid/Geothermal/Aquifer/Data/Rock.mo | 37 +++++++++++++++++++ .../Fluid/Geothermal/Aquifer/Data/Template.mo | 13 +++++++ .../Fluid/Geothermal/Aquifer/Data/package.mo | 5 +++ .../Geothermal/Aquifer/Data/package.order | 2 + IBPSA/Fluid/Geothermal/Aquifer/SingleWell.mo | 12 +++--- .../Aquifer/Validation/SimulationTest.mo | 6 ++- IBPSA/Fluid/Geothermal/Aquifer/package.order | 1 + 7 files changed, 68 insertions(+), 8 deletions(-) create mode 100644 IBPSA/Fluid/Geothermal/Aquifer/Data/Rock.mo create mode 100644 IBPSA/Fluid/Geothermal/Aquifer/Data/Template.mo create mode 100644 IBPSA/Fluid/Geothermal/Aquifer/Data/package.mo create mode 100644 IBPSA/Fluid/Geothermal/Aquifer/Data/package.order diff --git a/IBPSA/Fluid/Geothermal/Aquifer/Data/Rock.mo b/IBPSA/Fluid/Geothermal/Aquifer/Data/Rock.mo new file mode 100644 index 0000000000..5694cd33a2 --- /dev/null +++ b/IBPSA/Fluid/Geothermal/Aquifer/Data/Rock.mo @@ -0,0 +1,37 @@ +within IBPSA.Fluid.Geothermal.Aquifer.Data; +record Rock "Soil data record of rock heat transfer properties" + extends IBPSA.Fluid.Geothermal.Aquifer.Data.Template( + kSoi=2.8, + dSoi=2680, + cSoi=833); + annotation ( + defaultComponentPrefixes="parameter", + defaultComponentName="soiDat", +Documentation( +info=" + +This soil data record contains the heat transfer properties of sandstone. +
+", +revisions=" ++
+")); +end Rock; diff --git a/IBPSA/Fluid/Geothermal/Aquifer/Data/Template.mo b/IBPSA/Fluid/Geothermal/Aquifer/Data/Template.mo new file mode 100644 index 0000000000..208030b1ca --- /dev/null +++ b/IBPSA/Fluid/Geothermal/Aquifer/Data/Template.mo @@ -0,0 +1,13 @@ +within IBPSA.Fluid.Geothermal.Aquifer.Data; +record Template + "Template for soil data records" + extends Modelica.Icons.Record; + parameter Modelica.Units.SI.ThermalConductivity kSoi + "Thermal conductivity of the soil material"; + parameter Modelica.Units.SI.SpecificHeatCapacity cSoi + "Specific heat capacity of the soil material"; + parameter Modelica.Units.SI.Density dSoi(displayUnit="kg/m3") + "Density of the soil material"; + annotation (Icon(coordinateSystem(preserveAspectRatio=false)), Diagram( + coordinateSystem(preserveAspectRatio=false))); +end Template; diff --git a/IBPSA/Fluid/Geothermal/Aquifer/Data/package.mo b/IBPSA/Fluid/Geothermal/Aquifer/Data/package.mo new file mode 100644 index 0000000000..42fe1fac8a --- /dev/null +++ b/IBPSA/Fluid/Geothermal/Aquifer/Data/package.mo @@ -0,0 +1,5 @@ +within IBPSA.Fluid.Geothermal.Aquifer; +package Data "Collection of data records for aquifer thermal energy storage" + extends Modelica.Icons.MaterialPropertiesPackage; + +end Data; diff --git a/IBPSA/Fluid/Geothermal/Aquifer/Data/package.order b/IBPSA/Fluid/Geothermal/Aquifer/Data/package.order new file mode 100644 index 0000000000..97b1637927 --- /dev/null +++ b/IBPSA/Fluid/Geothermal/Aquifer/Data/package.order @@ -0,0 +1,2 @@ +Template +Rock diff --git a/IBPSA/Fluid/Geothermal/Aquifer/SingleWell.mo b/IBPSA/Fluid/Geothermal/Aquifer/SingleWell.mo index 1f8e516921..48f8cd082a 100644 --- a/IBPSA/Fluid/Geothermal/Aquifer/SingleWell.mo +++ b/IBPSA/Fluid/Geothermal/Aquifer/SingleWell.mo @@ -30,6 +30,7 @@ model SingleWell Dialog(group="Properties of aquifer")); parameter Modelica.Units.SI.Temperature TGro=273.15+12 "Undirsturbed ground temperature" annotation ( Dialog(group="Properties of ground")); + parameter IBPSA.Fluid.Geothermal.Aquifer.Data.Template aquDat "Aquifer thermal properties" annotation (choicesAllMatching=true); IBPSA.Fluid.MixingVolumes.MixingVolume vol[nVol]( redeclare final package Medium=Medium, @@ -124,7 +125,7 @@ equation if nVol > 1 then for i in 1:(nVol - 1) loop - connect(vol[i].ports[1], vol[i + 1].ports[1]); + connect(vol[i].ports[2], vol[i + 1].ports[1]); end for; end if; if nVol > 1 then @@ -135,8 +136,6 @@ equation connect(groTem.port, theRes[nVol].port_b) annotation (Line(points={{70,-40},{50,-40}}, color={191,0,0})); - connect(bou.ports[1], vol[nVol].ports[2]) annotation (Line(points={{40,0},{1,0}, - {1,20}}, color={0,127,255})); connect(vol.heatPort, heaCap.port) annotation (Line(points={{-10,30},{-20,30}, {-20,-12},{0,-12},{0,-40},{-42,-40}}, color={191,0,0})); @@ -145,9 +144,10 @@ equation connect(port_a, res.port_a) annotation (Line(points={{0,100},{0,60},{-80,60},{-80,0},{-60,0}}, color={0,127,255})); - connect(res.port_b, vol[1].ports[1]) - annotation (Line(points={{-40,0},{-1,0},{-1,20}}, - color={0,127,255})); + connect(res.port_b, vol[1].ports[1]) annotation (Line(points={{-40,0},{-1,0},{ + -1,20}}, color={0,127,255})); + connect(bou.ports[1], vol[nVol].ports[2]) annotation (Line(points={{40,0},{1,0}, + {1,20}}, color={0,127,255})); annotation (Icon(coordinateSystem(preserveAspectRatio=false), graphics={ Rectangle( extent={{-20,100},{20,-2}}, diff --git a/IBPSA/Fluid/Geothermal/Aquifer/Validation/SimulationTest.mo b/IBPSA/Fluid/Geothermal/Aquifer/Validation/SimulationTest.mo index f59efebe9b..301a553e86 100644 --- a/IBPSA/Fluid/Geothermal/Aquifer/Validation/SimulationTest.mo +++ b/IBPSA/Fluid/Geothermal/Aquifer/Validation/SimulationTest.mo @@ -15,7 +15,8 @@ model SimulationTest nVol=232, griFac=1.1, T_ini=307.15, - TGro=307.15) + TGro=307.15, + aquDat=IBPSA.Fluid.Geothermal.Aquifer.Data.Rock()) annotation (Placement(transformation(extent={{40,-20},{60,0}}))); equation connect(combiTimeTable.y[1], boundary.m_flow_in) annotation (Line(points={{-59, @@ -23,5 +24,6 @@ equation connect(boundary.ports[1], aquWel.port_a) annotation (Line(points={{0,10},{50,10},{50,0}}, color={0,127,255})); annotation (Icon(coordinateSystem(preserveAspectRatio=false)), Diagram( - coordinateSystem(preserveAspectRatio=false))); + coordinateSystem(preserveAspectRatio=false)), + experiment(StopTime=31536000, __Dymola_Algorithm="Dassl")); end SimulationTest; diff --git a/IBPSA/Fluid/Geothermal/Aquifer/package.order b/IBPSA/Fluid/Geothermal/Aquifer/package.order index fda59b4ae9..084ead276e 100644 --- a/IBPSA/Fluid/Geothermal/Aquifer/package.order +++ b/IBPSA/Fluid/Geothermal/Aquifer/package.order @@ -1,2 +1,3 @@ SingleWell Validation +Data From b073350fe7952dc70ff3f60a014f4431ad3f5d43 Mon Sep 17 00:00:00 2001 From: Alessandro Maccarini- +December 21, 2018, by Jianjun Hu:
+ +
+Changed the default properties based on Table 4 in the + +McQuay geothermal heat pump design manual. +See #1062. +- +July 15, 2018, by Michael Wetter:
+
+Revised implementation, addeddefaultComponentPrefixes
and +defaultComponentName
. +- +June 28, 2018, by Damien Picard:
+
+First implementation. +Date: Wed, 24 May 2023 15:27:42 +0200 Subject: [PATCH 04/85] Added documentation of validation case --- IBPSA/Fluid/Geothermal/Aquifer/SingleWell.mo | 10 +--- .../Aquifer/Validation/SimulationTest.mo | 50 +++++++++++++++++- IBPSA/Fluid/Geothermal/Aquifer/package.order | 2 +- IBPSA/Fluid/Geothermal/package.order | 3 +- .../Fluid/Geothermal/Aquifer/results.png | Bin 0 -> 85904 bytes 5 files changed, 53 insertions(+), 12 deletions(-) create mode 100644 IBPSA/Resources/Images/Fluid/Geothermal/Aquifer/results.png diff --git a/IBPSA/Fluid/Geothermal/Aquifer/SingleWell.mo b/IBPSA/Fluid/Geothermal/Aquifer/SingleWell.mo index 48f8cd082a..6270bfd501 100644 --- a/IBPSA/Fluid/Geothermal/Aquifer/SingleWell.mo +++ b/IBPSA/Fluid/Geothermal/Aquifer/SingleWell.mo @@ -22,12 +22,6 @@ model SingleWell Dialog(group="Domain discretization")); parameter Modelica.Units.SI.Temperature T_ini=273.15+10 "Initial temperature of domain" annotation ( Dialog(group="Domain discretization")); - parameter Modelica.Units.SI.Density rhoAqu=2680 "Density of aquifer" annotation ( - Dialog(group="Properties of aquifer")); - parameter Modelica.Units.SI.SpecificHeatCapacity cpAqu=833 "Specific heat capacity of aquifer" annotation ( - Dialog(group="Properties of aquifer")); - parameter Modelica.Units.SI.ThermalConductivity kAqu=2.8 "Thermal conductivity of aquifer" annotation ( - Dialog(group="Properties of aquifer")); parameter Modelica.Units.SI.Temperature TGro=273.15+12 "Undirsturbed ground temperature" annotation ( Dialog(group="Properties of ground")); parameter IBPSA.Fluid.Geothermal.Aquifer.Data.Template aquDat "Aquifer thermal properties" annotation (choicesAllMatching=true); @@ -86,8 +80,8 @@ initial equation assert(r_wb < r_max, "Error: Model requires r_wb < r_max"); assert(0 < r_wb, "Error: Model requires 0 < r_wb"); - cAqu=rhoAqu*cpAqu*(1-phi); - kVol=kWat*phi+kAqu*(1-phi); + cAqu=aquDat.dSoi*aquDat.cSoi*(1-phi); + kVol=kWat*phi+aquDat.kSoi*(1-phi); cpWat=Medium.specificHeatCapacityCp(Medium.setState_pTX( Medium.p_default, diff --git a/IBPSA/Fluid/Geothermal/Aquifer/Validation/SimulationTest.mo b/IBPSA/Fluid/Geothermal/Aquifer/Validation/SimulationTest.mo index 301a553e86..d7313ddbaa 100644 --- a/IBPSA/Fluid/Geothermal/Aquifer/Validation/SimulationTest.mo +++ b/IBPSA/Fluid/Geothermal/Aquifer/Validation/SimulationTest.mo @@ -1,4 +1,4 @@ -within IBPSA.Fluid.Geothermal.Aquifer.Validation; +within IBPSA.Fluid.Geothermal.Aquifer.Validation; model SimulationTest extends Modelica.Icons.Example; Modelica.Blocks.Sources.CombiTimeTable @@ -25,5 +25,51 @@ equation annotation (Line(points={{0,10},{50,10},{50,0}}, color={0,127,255})); annotation (Icon(coordinateSystem(preserveAspectRatio=false)), Diagram( coordinateSystem(preserveAspectRatio=false)), - experiment(StopTime=31536000, __Dymola_Algorithm="Dassl")); + experiment(StopTime=31536000, __Dymola_Algorithm="Dassl"), + Documentation(info=" + +This validation case simulates an idealized operation of an aquifer thermal storage system. Results +are compared with simulations from a test suite developed by Mindel et al. (2021). +This test suite comprises a set of cases to assess the thermo-hydraulic modelling capabilities of various geothermal simulators. +
++The comparison was carried out with respect to the test case called TC2 - well-test comparison. +The main goal of TC2 is to compare aquifer temperatures under a typical operation of an ATES system +consisting of injection, falloff, drawdown, and build-up. The injection phase represents the charging period, +while the drawdown phase represents the discharge period. +Intermediate phases of falloff and build-up represent periods of storage or inactivity. +The overall operational period is one year, and the sequence of the different phases is the following: +
+
+The figure below shows the temperature vs. time comparison for a virtual sensor located at r = 1 m and r= 10 m from the well. +- Injection: Water is pumped at Q = 0.001 m3 s-1 and Tinj = 120°C for 120 days.
+- Falloff: Well is shut-in, Q = 0 m3 s-1 for 60 days.
+- Drawdown: Water is pumped at Q = -0.001 m3 s-1 for 120 days.
+- Build-up phase: Well is shut-in, at Q = 0 m3 s-1 for 65.25.
++ +
+References
++Julian E. Mindel, Peter Alt-Epping, Antoine Armandine Les Landes, Stijn Beernink, Daniel T. Birdsell, +Martin Bloemendal, Virginie Hamm, Simon Lopez, Charles Maragna, Carsten M. Nielsen, Sebastia Olivella, +Marc Perreaux, Maarten W. Saaltink, Martin O. Saar, Daniela Van den Heuvel, Rubén Vidal, Thomas Driesner, +Benchmark study of simulators for thermo-hydraulic modelling of low enthalpy geothermal processes, +Geothermics, +Volume 96, +2021, +102130, +ISSN 0375-6505, +https://doi.org/10.1016/j.geothermics.2021.102130. +
+ +", revisions=" ++
+")); end SimulationTest; diff --git a/IBPSA/Fluid/Geothermal/Aquifer/package.order b/IBPSA/Fluid/Geothermal/Aquifer/package.order index 084ead276e..238fd56b03 100644 --- a/IBPSA/Fluid/Geothermal/Aquifer/package.order +++ b/IBPSA/Fluid/Geothermal/Aquifer/package.order @@ -1,3 +1,3 @@ SingleWell -Validation Data +Validation diff --git a/IBPSA/Fluid/Geothermal/package.order b/IBPSA/Fluid/Geothermal/package.order index 9465d41136..8ec14f19ff 100644 --- a/IBPSA/Fluid/Geothermal/package.order +++ b/IBPSA/Fluid/Geothermal/package.order @@ -1,2 +1,3 @@ -Borefields Aquifer +Borefields + diff --git a/IBPSA/Resources/Images/Fluid/Geothermal/Aquifer/results.png b/IBPSA/Resources/Images/Fluid/Geothermal/Aquifer/results.png new file mode 100644 index 0000000000000000000000000000000000000000..ebeba1b72610e5e06fad129a8e418a7705e6fbfd GIT binary patch literal 85904 zcmd?RbyQVtyEi&%kp=+~6jVZ5x- +May 2023, by Alessandro Maccarini:
+
+First Implementation. +p! zZ!Vta-QV8l-S63ZkFozb dQ!LBrn#Fsy+E`|@%rij$Y NnW0*_C!va0ap0-^z$5KX+#1o2`)s-%M~S9HLP%~uQK~-2 zf#0UP%QH+%;X cw5(HR~?J?XP;neCMB e%sJI&1+5ePKApG%f zu8Bk{ZJI>~z-;hB`wqm0s_c=-grJb%QyX@hrYAqW^)Fywj$Psm{JgiRUeKPjA8C3e z5&$LHrj?Utfi{-pfc_semHNv-xDh81Qr6^5==tG6ori&1$8R=s*8G_nJq5HhD`w8M zowGDT@8^adgKAB)J1fQ?k vKPAB2@<+s1< z yd1rps+@ZC zvyea4k1*rxQmq2vF`7Ri1f32gZ>_niZo3n^le%MxeD?Hx{M*UCO3+YGR4ItH ML7_(=58G0$GP^0ki z89nJ?kWz>JS7=RumYnDdQDqv2CI(emGkShzANq$j54V_79?d>#q>Ju(`HsjzRQuJO zlr3Z0!+EA=)pQmu7R)ERcQjakK5=5vQ~j77q{*v(psrVtP^^{l@x=$(z3-$3J+i@t z9jZE6>X}aJ>j5Jb+;1jd3aJY{hm<>n?qjXU++phs`u%)SHRQ?p6Z0n~>VA33rmt&% zI{$pRp}iqU@SG y9tHRF0w}P2GBqT+?LOxYKy!Qh+^Zg9!VT`r)1MV~9Gv_mr zH|%ee-!yG7{^I||npgY_kqgP)`z8N#>etJkbKiM>dj1^!S^eGh`~0Y+WrC%j)$;qz zVZr=8CqkRgeKgz2UxWfzXU6>a;JlA{H7&MBX@-)=M#lC>dq?s{gfkD-VMSXRDS3_0 z&;5?R5x(hClu0jtzWQyzB;a;?bksoZy92IlCqCy(XNJ1NI^Q}Tw;yh*ZvJjtZpR0> zXKa`Bm-o*a&t78-V?V^k!Y;(FCfEyLYroU3A8;8^93B?_HjMwSGnE*(1`pxW0RdaS zCMzla4+7d|BSsGJJ0`u{iafJ6Zj<=@ )>&EM8lTnjrHF6WPbB zejH8d?d#3g=C8U8dfJ-Xc2!TCxQ?d2r+wSp+BAOU>znD@odGu<8E)YkE*;y$4*hZ4 z#KtVeRZ0XRROlj%I2fC3)!y`&IHjb_R!C9^?aPGxK%OA)t%lcb)H>+>(a*BQuYOV8 zXshO`@A}!@_Z)-BlPFaty2~%;yLsN~vq5u%F@s0NgvI05=11~J%yu2qgY_NmLxI1^ zR$gy8&T9?|OjiGL?H!!6^Pgw-uJDe$=D(&!(?OeoBB0+TyINv>JfD3_$9neb8PO;B zC+SwPhxW2Y-A9r2NP< !udqzRu8& @5< zw>h@SQgQ5X?7r;ua;_@V*D}U3iOIVt&ZUYb{!Uek+E;I(P{&KnIm;Q#YEr^eoK8Dp z7GrD_!-$Vc_@FLJl8t-LUCD#w4lM644=MlFmmI&P$Ioz^l+X0}WN2A^C-3CoPiV)B zG6U|UrMTD8NJatSW7$-P)TeevwUhgbeN6kF$Yui{omTx}-E`fODkG)!wr=LwQ;G zZkm6JIdsIDe%Qf%D*Y-gt-PO9aIkTD)N+{VD@SPzCw(nF^!lS4N**HJDZpQ2ojV;g z^<+Y@#&SD5PmxyPBLk&Ep#rP?ydA=}=_D4e8Efb{^`*kJq(8%YP@`XCaqQ3ZO!fUI znjC|iw9^mD8Xa*=WO_cwi^Oj>s2FRiR)>}-lqu*?=}qdJs#I2Pcsk1;%7$wyx)+*l zWlxRQL+Xn+_a;U)a^7g_9NQKf?ZDld&jzqF1429#JcZ65@V^u|@`j())yidluUK#Q zTJ&1?3BvG%e}==^lh{oQ?eo7m?;EwX_RW9v9-aMRoo|;@9sdA5kUsE%H|%8Vu>Dz} zdqm0S3QzMBeH)z0%x`2HdUHa*Vn-OZKL2?je(z&Nbyrr-l6&9IG3}_BhNIrf$4T6_ z;ICw9N{F;%@mS}L_4VaZlST!%o#bFK$3o5qyIt#{N#v@>V9NUOrUYc}@S+rNfbz4! zJN?bNW&7`rEnewMIOnv&*!QtT0%-#4I7RfgsZs?zeutdLj%qd*r`A*0Z;f SVVS0Vw)v)&P)VR}c;05yWv`iPHCljZH!O6!9xv}WPoK6YWYrXUS W zPVXM 4L zZZPX52BwLS2x0cIg-do$PF;Q8pRJ8mG>lY7_dKgJ_3>BNh?ET`A|fIYD67Hco^oY9 zI5U8He^X?>XclXtt~qZOjgBK&2)xPG-zlih#EWwViNXe##iIA+xxTWME1*SXNQMF4 z&|6+29wM_!!iaY>Ebm+a-%$C~XMd}6pG}gK^B3weDAe|rv<2Y|qLUiDgoYC*f$ss5 zrW?QL2s2Pa0DzMp;^udV ;{q_@|)jA^QdmW@ZrPcFlU} zOhA<9!P5Nrc3bVS7t(i19{{+)etkT4#+#?u9~f&+GPh*K6kIr%tb-#+8pN5SeH9EC zkH4OuQ$icK*e|EdBt_~)U}^E0&H|DTr5x$}`dZw4Zq@fHp((9|G)S)kzi{`^(y72F zMPEr=?3Oev(Bfjm5OFXJ-N!lDfc(e(-CdP>B1ysGtCO5j+BDDAX}rw?ofig=y5kZs zoXHid3mpipdsnt3DP>3UH}6b*Vxe4WVTeUY-Wliqvw4UhbxbcBIb4W)kH){kUD))h z0?9Hf4yeqIPTdDiLcPz{B@P&;Py}%BpCUyb&1@Yr!)9p-d0_Xa+w83xNnKUfn9K1$ z94ST*HkIQDW3k95W{+G;BAFWJD9l_IrrH&Ur+-Q&Lgq&K8z)t57*h>RU@1dhMPh_^ ze zY6m7!l1Eo6 3#8oHalV z?P6W-KQJ@eXt&7L=XIAdL6xb0x8CN0&xej}5Z&*!MG8f*kPJA)QEwEMONCaoWpJKL z94{nk!}BM4^jz+XVhe#+J66?m+{lKUtz@&%)LYX+ePt4(Qh~Ng< NI>PU18>*~`4VgRi}x;jufdS9uylH#5Yn&8mj?=~lu)ND~W5t9)O7d0LeT ztg&h(_ps(cHBDt(zx?>ItSB^kT5mHnfjLEP@ GAaCewCJ^zovk+l?Z$ zHBvD9`izFc;VLr<0Ykfy;VM S(-9MLOLT+gHKnSBl2{h_L8~eHYm g{uDCNIj2x6%^~EKUa1HxYrAF^YOJT6CuSyIa|wHp>doAR!Rwv&qww4YqGqd zrBTg{#DVgTbag$42wjDKch09T8*7v2rQ~541V0MKYb0ovoeK1MIuQ&p31<$8+LHoD zy1TV5T>P)`M7jdD?qc+k!>aYo#IGw-Jep!sF^dlQyG5JtYa7@+Zhi-Bfhx`V^r67K z(%aN4vB`sqy5qDT0kar}!D3eIo_&dGbLR0MAEzGY_47TTj%Qt$O%yHv`D}j1YXR8? zGZd?>`g2grA?5koPwHnFi<$XITBTKXzv#K8NY!MY _9`V{o7mxPLdiy!aL+BV__&ie36F)@DdwK+d6zRbJPMM=@^^ixgMvSDO)2;k` z^Wovecg{wI(YV!Alb#?*rQtutkMq+!jEK`KP()1}jPQ3HC_ngq{-5jz8cc2mM zEor77&VZ<=Iqspdc(|lHr^-sL0m)C B6uH@!_boWTPMUj=ACr6Am z o8APkm(NHSx|a87zp+z`HqcMK+ N~ zN5lfRy0g~Y&7fN 4jvC}gyKDomUjgBzg63Q=sBc`xPO-;jpkQc zSYhiuEUHBeJXYzoADJdwkVs$tXGd~7`?@laaWlDfs1e@hF%iBunojrjT+#RyvV Q4eiqJK+bf0RNpi5jeGb i%06+5HVemLkoQBjsJD=TfH(kVark)`GC$=mCHY7^O(nR_u<-E zXdHZT8UmR^%a$Pt_C^Bd_d0ZVUz>R;W}w{fi0J|~Gmqh{UR@m_h#z=yG>p{WieD5H zEF-ix#9!*OSg+pST4v2Hf9_O0A0cefKJEa~^)(B-KQS4Kz p;Uhg`IC**-HDBfBC4 zF?a!1)Vwq=-+9k_PhKTIHmu;I&}f`!aKmYn$9EzfEU^6J$)uDAsZScNpjkH4q)a{V z$PBVdV&ma-no_dh7xSjIO0V8ZSy&-7F>r7+(oY@R23jI}f=`#?8re~xOWuw?{k%P+ zn%B;C&&;WYc)v{C%~f?tw7Hu4L-o92y&yz1Y~4F$@xG}wfkvS#H%-0g&T=g^?E$@4 zTwKkhrEm03j;5>h)SPX y(HJ2_9Nc}N9nALg8q?Z<+XeDlEUTEb@ zO9}V$6;26vRjctgN7S$5dU&|FXkqJv>m$n)Wv@uu2YEF&j7Ll$igAgw i8M=qYjUn z6}|BBSX0aHtr-qWV@PwD-eo?#zhe0j``a@ZI)`j*)*Z1nWbZGW=m!kne#Y9-^Gdw2 z*yyhOBvi;fnD*VBi%^s_2}6vfNhw~TA7^T?9uI2Vw+sR>sFiI1g8Egz51k)^*l$G} zDk;x|z$h&AOrbE6rU)u$KmZu;egiE}L5S3e9(-yK10{(|`ngT8Pz%i?XF$?|(62=6 zt@upY42{;S3=7@2gppF*gaAGW?4QLNylbfTWQAcUoPq9!$7}|oi!s4yEDV%@mKwDG z?h(lv5Kan%8{jYGiW3`)@yZb*qmm&M4cclXb*%&qLZB_LmgEU$#t2Q&Z?19!?uMeW z*HDGgjgB|*81Lx0aDdP1;KDp=2#8xvPiuo{=!!?*I%Wtt=JGP(Q7#d8ppZF`o)$!7 z#6{`|EJO{_pu!oq9RZ0Ny-6BI6?(!U(wl&0I_n1wAh`50{XYgZG8Tk4H}3zdsOLXz zka77ty4cDU4FHUuqr#-&9yoYT0U8k?NI?(mu>>S2&~jUL5WF3?SK)#KAl!BWy#+7# zGUq*;HI#vU(ijZ 5FZ%ifvY!CTwVvzf_Ko;8Bigi zkudp+`S#H{Y >1tR!8(MZdr8cy7dWYwT$z46AW>h7R zHhe ;*cB}RHc(qi z i za~Nv@Kdcq0O%Wuo2rF^o?%+Usz$LI~!?HAcV}DG9t((Cb55N4ff@IGcy|+%&bhE#u zka)^%zWrrn@VJEM`NZe&M>5m5IuShVDy=GK?`Q~zB9u!KiLTpR@Y4yvcsX)Y^`o2b zAD!;3V>gqZ=Ljp~?TV{++U7=DhYgSSTQiK!{g+|v{Dx6 s4BJM0PeK2g Qvqa+W=jfI8h zGIzkJC*mtaZPP(M* HKx~+qMriV~fxa5_lHeeYPVy;#n9*Fi@XjNH@xu z)XiP?&YtsPu*S0TM_87oLpIqrxvTVKlb_?P>x3ZUzZ(UeGT#UEKz;tIHC`zAonfzc zLuc#DFFV}a#cuZkjlmf`)!{FF>1nvmyJYS48}n*kzs8YnX|OxW2PD-4TJbrRB0xoY zc o1O`&TsMUYnQxoowFgYo z+OAsbz}T{)%7u6|mugd`zAjRJ!lgI8Q@N+obU>YP$^!uXVD5(w=Knm-cs&I&Q?N#^ zlk1CNEG^ 37zppMp{>|>*(z2KFw(T zqPc| `6s zYQmzB!zp9^fQY-LNIFeZ-ynDPIy&2n#zz-{fRA6-%W@h7#`e~t&Egalm0^O*duarf`p{gDVb8-BA#jh+|Y z0OE&6d-~k!inEV@@T`@kvl#}RILJCyy3HobwsFnZ54cy)ui#Ag2nk1n%EzttYm>{~ zowo2(QSx#g-{DRw&uy#7F;SZP v>aggfEDyA4`*og-<1H-?uoa47 A@Io<9E1h={Cvw<`0 z)Hc=Jfs#jMJ9P*T`FPi6DF~pJ%Zp>tyJXt3PTqE5<=bagUQDl|zGFOTTKAsC)&8Y3 z=Ylbl^pCf@{2_4VSV2eMWQ|&u%GH8*+^>h !Uc?OG1Cq|5L8|MR z+Kuo!zB~w@e(E93W&UhJ2#y_LyLu;fM@K?1CcEFGEYaK+ck E^(GqQu6Oc4`l0bL;Ex! NCV4C?o zw0jwR#wX 2Z z;!!vauufV70VV?o0Rf|WdIZoy^*@4?U=RP6X-8pB;;BHxhLhb;%1U<-Lza-ErdQK} z%H`v^>(pZa^%B=dg;V+`p>~&RNHEvFwzo*hIp>Wu-7)~%R4Ajmh~ajGI)E;x&v1^G zu @Ll N5f zQ*B6EI)n6&@jJKPepEjuL>y>B7*-=g1(HF?>j;e#Ya-BJq8{oXX7Yk-&OI1t zRslECflkD`n5zuwi&F}h=->N3@S@794a(7&1c#qpiR*!1hh&^q#7|*{O;2TZkV%k= zaS>ZI3IvE|U?oM=T8e8n7fg4X12Ic0xg9mS)PKl>8GW+fZ4AvrRR)vu%QV!l^q!jk zWb)SDzAaLY1)$P#B@gpI#!KWx@1-r|K8SM{-Xppr&@=YEn?_1&9D6l+W6Gyuw$?uH zg=dQ-48|^V_fY>8z*sJdCyAF`Sb2X)a&@psGt*Pwtw*9YubFh !fIzUyRkK4643(tVc9n+;{|gnw~G=X zDxqM3qunyzhn +}Ht^t)Orx73a^R_3woA@k z42*$`{5f!E$n{I=>r4!5pCOv`4}ZdiMmVKE88#`Mkg**x;W?-CD N&%SMba47arrloK}BFhVz_o{5v_g)+|c3u0%{aI`YUEfRlQ6xTJ ze>vd-#DQ=sz3J30N%!`o% p*D zH-)k)k3K)V2`p_wkQ^Rn%6g}c(fyt+Etm5f1Ed(o^{`3nIqV3(FZMr}?O6t3MQkB) zoc$lhe5}Z4Fo>ukkIPqo1YnjZf+GC`GB{OY>8;m^7!GP!t*<)RShCJeL59mWJNqXQm^39Vr2)NVHQraIzb=JDXH_lda74hq<=XOncMx5FQN{>lw z?|Im`@(lKNth&GceFX0xX`!hcOWjDBQSvUIpdr5q*x>2w=iKoXgY3nXr-r>JKxNz{ zVQpTe70Gqa?9Ym1E^%$ajP85(lQ*>Zci$}6?7k@NZ;BIzDBgMTdN0fSb($olTN;nI zs{pNdl|)|*=PQ>npYmfns`I$L@>Quc3XB+(nESXUTX9{>tCb2BjPK%?gp}`3JRcUK zGNnDrKf!FsPb$cfZCiV`;#x8*Oi3ytMej{!+}!ubveF0cf((#l)8*@byrzwdy4 zVeVmkgS7wNEdFnWO{Xm!f+z3VzNnaAXGyrUS*PW18hJft5Po_}kuLoP5{2F^l_k`o zLEprBh~WF&(cHOAn)4>6Qp)A!6?CrW+p#^LiXLu^%3FuKTJt6X+}~-p=k6KaFQeVX zO?&xK2uPv_CqMskuKP#xK#y~Z2FJeRB!-?t7r)U@>VyV6rkCuADX!Tp`b{Nj;_le| zQWk_mpMK1JFW#GaTq^4L^R@26rw&u1$ZjXJ0~P_HS04ouNYf47c_=bGjCF8cvda@n zKFR=d$JJg3r^$iWr`2{VHueK^@iTmI-OVJHr7nD{i;l0k$MYWAUnhy0KdOG61oQa; zJlS{)Yx@Mszb{mr2otZ_snYJCTPK3hogXXMft2$xt3m9W;#5Eq5CeSqESh zvoG%G%t^rb4mtuFdlik`x_yH^_iinRQzwT&Co^n~j4?^RMC%zBv+ga?*Ed#3Qf1xH zN&oyflj(R~y6;|axdPtd{ZpQHPsf+-A%e!=i(FR{R}9hGn-J*z-$rPl#BTHUIeJk= zzw?~cnU9yPB7_3lu$0}aq6Z@GHIgYs`vlHdvtQ0cNW@A9nL$f_E5F5XTdj_{&O4~8 zXz=b`9Mcf9h9W>u7xY%jcJ93RPhCKilU=XNr_H(%-}SlZ+yhN9slqJx>Uja=9zCq` z;*rTf5Dy-0ZixHrKWYIYJA(f>otMw!4 S5N$#bI^ =jMpw7pvZ1yygN|w7E^R$F^{19KXU-Kg5!UN+ p3vT=a9o&%c z?vwZ@3{4WL$J$*I0b!LX-@G>0*bineuyjJx8_i!t%VRnc6yQ}ub#RtztoxR{$N6;~ zt_ 4rFWShx%0X0uEIF>%?E)uqL7atb7yhrk~O_AM%N<-n~+;y>(7yP^TOwF zP`V+*VuYKXJ^InzU+Z|&Z7rIjV&TgAFagF__Y^4N40Lrp&~&w~L)f3VwwWcr$x}#t z&-kq;xKB62Q&wIiSgi>ezhe4mJh$7S<|W^kk4BEOT;~qj34Ez4lk6T`gs5=d4U9y% z^9; ;x71NcEJ1xS)ey=Nrud1lb%`2I>*c0f_>o%V-5 zBa9$P@L`Ka_nwuUoTR*ZC%TN|0Ou+rU>zy%nBMd^O6(pC7Z74_^4)l#+zwl5QPKh? z1aDv}^&>JIAOHmg0RXWN!I$I#4JvpJ0ROYhDEPvEfU*B+50D!GYc~hUccHxcU8yWl zJd1^Dn5X(b9uu|c0lWqaG2)vWBxvFwCzTY8)m6(Jw*0V>k8dB7>0pEXZh*=0r%G=P zJZ1_kio{$fGrw#LU4l3r7z1++L!av}#Tk}U^cchzeWV??4)pt(@lg-?JChIqePXmo zL>>)i9UgcGY7By(8XLMV;;eg3RV#3o6Nmbms>|GPkZLzIG`yxt{9;EB8%2b#95xLo zCOY9}tGu@BU~LhYtpzhJcek$W?|u$mrTnlu)wyD)hQ7TtC%p(hEqMk)0fHVaGq2af zFz?Sd qzISy&IrU0+C_`pTK=CH$Yt8Fb4)A?%PKE (3jltw_NnC!l580^%~^88StpSG*oSdCIK`4E-UorH+!n$0*x7hp zM(ks}P1NA&@u2T~*hu9hpPTZ)gL1ZW)=T33B?U98g3NjbwGojX6sEd!xoA(tFu!$- zQNITNj#I_i7}^ycNEEs J}GeltqoldNH_LpO5Q)w>w|_+w~0~Soa6_ho6KRgp+14e zL!lqV09dN1sOZma4B~-&iL3R4F7x?Qm*H7Y p$o!k%f9~*q%!t3`WeVKxm$ J4zw~04Yy6AX&^TruTF~C2I+XNc z)t}V+>wJjK6k(!0dc#Jj e7qC2sM zR^J*?yT~_veE&2aJP2k0P#4h8eDIvL%OB5v5NUBtE8BUTOcTSsZ{HB_-SAu;Yp@?e zE{=xQ@yHoj+C}hU)q0Y8Fzwei`T-Bf-)?#~-{=y2rMKMNxP1#@mgOH#qgJ%1qQqka z#Me;r12q8CWxFaC5Bio#ET7z%LC5a95p08`vU435P6Tm%1lvdMp%pjtYIl3xQMKfp zaeAlAhZOQInP6p(un#_&&Drece3qFt8npDEa^Qnrb_$Gu$R8m`ahI0r@d^C18LB(F zP _kqe`x8m!qEGjGsWNWCbzeWRzoEiJI1%}(I=<#Z-eNmv zUzCjm(%(oSdnt+z _X1yDiWkoM zD3(9H@gKykpifK;9#rlbzw2`wp3l5j+RtrbD>c{*3e+w(C{S@tdn?bg#aQ+pjy}k8 zcYIL!?$~VGVk$F9$gNHQA|pWiYg#FJd)A}7LS}s~S%!RTF41BZ&+EX+ZXk1{-t|lw ztzN)U=C618n%&nUoDUdx!1tYsHCwi^B>vG!OV*If0pljTL{#BU{z=#k=G2NmAlCuK zCnl!sl>T)V v91bt~6)fiNxWxc)+El{l$bx6tQCLkb4oo*r&|+J7BqK1qQqf@;_T z#3Eo-8$=#P_o~B~3gT?3loejV{7|z42X4qFi?fV|`XLk)`oA-AzTIxBM!f_?Yp0qX zzFb+CcE2G`iHj5x`!HZBzF9r;W*~d8#`C*voaM1?$*ZTH-h# hD=Ud`K*`w!CI%X)K2d$E8|JS|(i9*5@hbx-JO5Pqzr{y|g zK0G nPHoP$cU3H-Pisw#EN2QU85&|M%`m#ky^?_r*0>!1Y!c zHFP?!r4|oZ41+rVXYlmjw$PYhItA6m|EA~se`1`0=)l0hDbI98P)(Mu!E^@sUs99x z|Hh1e03;N00JMI_8#W4G?uD-ozCXkqNG4Qz94+`yQu;5c0R~8@<9~r$|KDfEzadiP z2M->+Klb4VYp)`%y}bq1H&L-85W`4d74WZM+5ZkRS|~36T%YfNHBOiKO5l*F!aD2E zDqlcfiMmn(uVS43akU11W>SJX-+&xT4s1SJU#)eM!Z0Ju@xEL*?7F@tx?Yi)@Bb6l z+VefkceDbvf=X~tx9ZPQbQq}=mmp{%P|0b9$SY#M3|TE~6okWk$cw|Ts3!N 2RITl=0SO}P`nCfWTrMg6iy8b0-8995d?%I{F}ODFwJTd)FZE=jfL<1D$(H;Fr@P?ES- z6o*iO?LZIj-(&!j#M7iU{uJv({i2wSW#q(q0P5R~y#wr>Gd6Ale$l`l#2WQLHGM zlBG+Hx$@Zi>LH^QkwyC4H+r~Z5;aFJDP(#!1EQ~>sVs_@maZT|;^Sq##HpK+F9?wi zPp&~+F?N1#O!mO}EYdP!q){k>WvY>9xQ_B-4Fw{yTn8HjE^eMAM#y^ P~a#Tz)(4AfBU=HyV$p9G1ax0o3OcGEUmPU~$U=xOa_`IO4u?Ka)K& zfbNc=AlBU#&}wC~N1*}W{ 73Y zvz7gqx2LVAdpO00)F{5S@Jb Lx-0E>BBDD?>0rqjL_Pu{X?n?;#aiqaj{@Vq$NCmx3AK6cH%o8v4n)V&YdAr zs^7b!psj|S#oXtPCmzzlwH9B{0&tAq&?ZORk3_Mi1vAcEE|$xRgHo4Px2cMyFXXqw zt%3b#n8$)IP?;q#$JF2g)&zhBmVO9=_sleuiX-`M44iso=Ju%wME$Ii6AIh}gWCa5 zrIW4i6baqW0>2{z?aK=mc)M^gEm8Z?!Q}02=Dfq0qmR?(*V^%6Ek_95b +V8ck%`%bNnYJ#T0sL)u}SkVzkuwti8Oa2AhaS>ryojJGxJQcs`Kiv3L6) PKXcp2xWI7;Iz!Cxgs4Cn{LadB&8G3O7B?bcVg5G-b7i%E@$~sUKm`3%*zi8H2 zdovXyFl2*n-5r3t-cT8pS+nj1X;N|17*(I#f_T*p2R|D8J1taYVQt&=+5t5%KPuGM zG7 o@9jB~uXz1g4%D9W}pdw~LHc(o6v$~ebB7zrhh4&XBhmPfKdsy|=@ zld}Z^)#|9wMi=RHs*eb1-Z^rY69)2KPHXqGMmSx%W<1iDFF=+zvk8G;bymlq7C-cj zuh5d?5l%V*IP6~YCbhEJ4TA((9wepdFXJ_Vx-Ycny1wS2W6L27N^3e-ryXXJ7;k)A zvA-@Y@#w!pM1F-mFxPmPCezG4r`7FrMIIZ{F;unzxAQN=`fH?unZWAD*S_+*tbq0F zquYB9469@tf72@cA-OH)ZD2a;i*57K#m@Da-8Et}{(3*+`U0+b{l~12K|Ipy8Xn(% zj5T{0n@sEIcq3T_vq`p?Zr}9oLo4qp^@B umIp!^4T^xm zO?&InGmza{{yu`|W_QHjpF uf`n~sh#?6~g m1o7NREJ6n43s{Elci3ZoJbtzAeKmgVv|@Pm@d)XC<*G|~{m1a) zV5i6hF Fyp@8& z`2Wrp>G1M`b#NPE?85O07(>v$a~FeYGw0#j88Cc-Vzb>(q80lw&%9n+Km%QvJ71(N zx_mAFLVS1P$N_mhICmMK<1z;M>iy-(-N R3 z*ZC``$c_fRM2(YWi^*M${@gfwn+FOEnJAZW7Y6Q8++lz8X#PFVv#f}$exu+DF#+$v z+Q%#?)7%dR(k=aCQ4!-*^a4a5Ms>I+SCJP25;ND&WFFnD!W5hjV`W*;MeC?Q^ZKkp z{E<@f-QPhhFEj?q2bZtJth(X8s_bA_ueR&j+uLg+e4Kl`UNRF+n95kgxl*!Bw6w(^ z$df=?^OaNZA6{G;HfvlB(U~F}y2nnVmdw>Y4^TTGUsZWnIIf6xgJ=%QE8DK|++sjW zk1;{Nx~+J$viwkb*;rSefkNO2q|eR&1u%&8M;L#Z5%ubT{H}#eP^9#s-pq;&DE#4+ z0TF6+*eHPYVe)K9(6-A>TgRS$w4jQLw~<0Vb+DZ?RvqfdL?+CLqO2jbLatrI WAa^n#j?fiM9oBf< zY=IH TgU;KK?&D+ul@%v4+V(6BA+!L?5|wZ zS3i*8ACHoJcD3we_T>%}ljVN0QrVc(WDMB#&x})q)8x$BP3Mdh(FwFdUW|VB(`lY` z?#Jie=Uc>-IK}%k1|R%I!)T)7LD3^9+^rghv={BCm1fDTMk)AW$0Bi0v0(2C*D0+v ze)PyCiRMHX@agks8RzfCG1K @LmL9N9gs?6f_XmlJ z(-t3eYc?BDy4cy;DGgQ=nlO71mb+HJTrF<6n3xOhA`KfrsR3f ijUJJ<{ z4BM|$wp^}=(=_n1zD|uavVCxmb+ym|nT}F0NrmVD1gzWgTWxy_;rq$~`&zTq$>!BW zMcWDcM0p+oLsyUrWKj_@5#y6dFxTjmaKn#_?-!5CuRX@wcZ93Adm-WiME{4dH;;$1 zf8WM+w@?yow3 JH`Z^+ 9lApM)SVdB3#`XC*XzN?tUytZ+-)7hWdL~NG@S|_8xr*d`` z-mUruCXq0&SB18%9RJeYBw^S#cTeJtwfp0DLX%9&6+uv#Whc8Qh!Ua}De*2M2g|rJ zkl1Vs8|`qOfA?M;@%SX%;%2@P(ACq^b3xqtMur@F$nA}Dt#$jdd`x>uG~1?~x=%4K z{J#F>M~2eTw?3}_R96+yQ^5;TG-J62*@y2cLb-NR9X({8suI$5WbJoyXTLn;vL7FO zU#vbrWY^>GkV~gOH>7Y!yRXc**_-~goOiF8c;r-C+M`c6-A}g7u9jO$w%pUEQxvOD zd~<8}^=9)z*XZZCn(pm8_uvMGP5=gP&Vc;pWV*5=2R`bRy5DD3$poW~u*W}%j=bg{ z)z=Zc6FtMtn15{H9qY73a}!bk4`AVv?Ry;UTvn6ki`?I}S-*cz(JyYjSl9Kb?!ClX z{10BR+r7I|H8JcS+jhcg@&I+Dq3Ky*Wo|vb5-WY05tEf TcgE2)*r!l3g S+GaUPRty4 zb;P)Dc|(1jM(C-DO-Zr5x#MqU-i-a2d$%-EkLMUX5SsIDTBS4n&D>*ARu_rvzV$+a zBQM-XZ-6jb6KlMe<2%RBRM%zd%OP;p*{C_~Rj_aT_y%is+5 HB1}c!WUMcH zxfOY?kBt7mYpf8bw^9;flGpMc!WO%~2Y=AW!;zR8BMrH=&|g||Bc-&S9)Z;9P4SvO zT9NFy53bh}Co|~#T#lxnHm
3`@50> z%+RCo!FNHWFV=2(C?z6>6a3xg<=H6lE>TKVo5QvihC6eP$#17K+ALcQ*59~(FMd~QJ^rl25Y?9jp;U(->~G)&VfTmmypX+i zeg4xRO>WPB4*(N-eW^*GVndxq!_iM%T?Mww-=SWr#%VGcdcXO@4bw z6QJsy^)Yt* zST{H482rRuJ0YJKBO$GK1()d~ZHW99c^z>@Qp;V}AB(7>vF>#yPc!`PcX8C!gFckW zeLm+Mr4sYCP)sG2U8Bhf%u^cq4Rj;>eoU^SE|yG=@y)V5jr0?jW3tgN*YD+ySaTag zp2dBHXYDR||E_p93j6qkHr7gm^@K(@%JEb_XBUt@Q $riIU_g>ktf|BEe9!)04)Yg|RLXHZs ziR7}{w73N0^r*XgjQd d@LnMhTK89(cU?N(@fRqt%;y>YzkM03tm!3Er3v8H`{`g#jnX+!cweIJZDQSfYs zi;~%$m_K>z6W|P)zP~or3uaZB(X)+&DNpv&`U1v9 N9xlok+vhN^gyGfexpHBmp9!d7soa6k`s@FZ!g4yEv_&`_g zp)WKpFPI&XK)}Q1c_fvnx*KKD2>>ik;I!HfEuQ(|0URm-m3k$(K@`U+tZiMZH>0dw zC%VniDCY3@6aKuJwr^~<%nN_%725;!L;I|H z_?Y)f80pqD=43($$A@yon&A4w@}6y8BJ|C-8kWg^sw)sH;|0g^!sl*n)yAG#k*K#W zaWq)paZi6EMk6?Y|5!(yg}u#i@#m`jvNS2}#26bdZlkt?pI?5JAZbXZtTZXm(iE4$ zOnl9s7N?iWL-|=&HhbS*UwY3%BSu!2c@If8Mv0JY8q;3AKBjip@Ty%3Ipu=_J5VuP z+P?gB@T392Bjb@m@}Q+ZFBoOElCCkH(YL1A_T58Bg#O$ekhZTn?1WN~iFUa}be4-c zJq;k>)-l$3LvHzN3=ZHulg}l;D?Vq$jg^IG@;!MYDP~-elgZC@H20If{J@|8|18;9 z;+&A=9)Z})CZ1lVKb8`tH{8ABl6cjSy+{wBQ>-})(jxs5erk8#ZRJ|cJ@}nH!Dr9t zJ0>58f*Ja(G;FX?Y-g29j>RuuI^7`X=Ryy?i6@NW`S;`b-v`MkYQP?OU1oY~Z5Y3N zZa-&>KX=$Du!;94iRO2P`41wIR(-rlGvRZZv3vP~J_;=36!D7EbZSZM;R&>}B*}fJ z-Sf;#j~kChy~a3cW4~%>4)_zy?bVAcb=MX1KegvzNSpcDUKvzlOv#A8pX4d zQc&LAoMp|nCAYmiKhhVs{n;g%JWR>)L#o)9K|hQgOw@+Eo!9R5X500a1S7v_dtv+c zgJlHFT1Aoi2HQg%K$ ;UXI}U%6quxkCFYox?%~MrJe5P)MUflhy~cAo<-E *@7>qmSI<1l%po zZpemmtx3;{3r!ZU2xr~Pb$m|3mkhgS`;Rl{SKxVnkBX|=`v%>>%A+iry-rPBue;qn zUtQwsN1gt0wb}NEBIn#mSd%HbkxPIdLg-Z+dXj=#%RIZEjlcPQGP|(02_z2C!KX;(f(N7e%HOCrd#n(wvT+hm~w7f@~P3Ag-=*M zhb=a7OI8ZTh|0`LLk@L6?s&7!F(X6^)vZ?*j^%f;1VgR$A-h{?NcUnxF8&;npB9UL zX&vU@VaW?x&l?t-InXuo_19SX0xvW%X2a+RSz4LZ#wk3*pmVH 6(K_vU>m zb2r+fE3OQw;$XP{%o%J*5u6G3O-Ed2ettB=i$TZE_J)+Q5iNbSAG@*hu%GkABrUGQ zJ6ICl+oqdm^DO9?hV^ aVk}m1{zQ+u#l-WyWT4 z{rjo=bi*g#hHo`|Ju2#celO?!8r=*iplcB%A+$o%-J@GZ9XO;Tj5QZ(Ul0X7O#Y4# z*^I^yk 0fC&d2QnQ6aDqVU2*|YJci_yk9V9a=Dma>e$8Wh zZC3b%>^a(RvmxnMJNQC74Ty<_jMA_KABP>9mgoU0^4uW<%L8WkwEFkt7pkVm$N0ao zjZ|HaRIF#?JCs*2`es0%#c_A}BbWT|$} 29ARo&d&5= zM;VVoi)%+DhmR$yo~gyvY*JNqyfb_?X7+Tz%xE3zKncAi=Jj~@E%h>9i>ES_%%~d! zRDWBusV>`@Zj68~j@(RaJw0X-c8fLPfvQtlkG;Vt6sK=QviI=qA&&ml*6sG?8#agM zEA2vE{ubA+HT*hHi#x_%$!_=Qzg~Wk`Xci9;~16BVxG u#$rRxEArXDXMGx 2+^0Hh;eI>-4Ob46=>F#@txc+qNKF#pq zf#jPW$Bbu?X12}u!iLlh@e_NBpGBIgEYtO(Q Tt!9g`~xP5Y_nZ3)9<_mb$|_tp%qyIbnTB-=Kx s6dl!SHhnF&pQMxCd#B6=^N&9txMN-k0}SMMrhD zZk2ZbJW;le8K;})`UI`~_6caHR}&w9)?{8d;@9uJuQ1>-cLkPI`pumt2xxorn^UZ{ z;tSOLN=wx#LQVPie2J-D9hxb3b7U(wNn7TpzIJm|a~{GJ_)Gt+YFLq<7UP}A{g^UP zUd#YWH8ueb4S(h1U>tDeGhszhw)Gap-aa{V68k?1NRkxOJ#^|gfgKNqhY@XlxbATG z6@Q}hfk4h>8BNq)b=Gd-@I5i?>l2z{TXD>-B=?r>6bp?!*B>X9fwUJ&jRtteUQIVj z=p?x2ICOC`JHl~oZvxr--k$j3&1Q*;$6pfSeOE1Rv_DZPVzjGk4;T_|-9q^I6YD9= zvuKpIjs+Yvw5Y%4)zUKY1NX^#9Lcw?*E5Y ;U0MJ${q2F6(V4Vpy)e2|T@QY$K*nOE>75GBkM+hch8fFz9un^yU4E z&D%&phiMIZ#IB!NTpNGx%kn>KV0{dDsFhLd+_;8i3y&=|SiAe!@1tSK0ztRXi(v=q zr8V5lDUP#6KCHN>DaR1EViI~P9V`oI&^Vw$fwQ@9pb>HJ-u>J~-Yy<*$TbZ7f-X9$ zvVJn e;l$^|6fgHaj$H4sj~(kVWa68~<%hj(Gg55&Y>m z8!fMEapOy5iXWuWhBNW)ZR)eOS=@1`>(IEr8L; Y)n1G$L)^b`{SUjQe>acl ze*Q1-gXkp4pDy9}4z|LFLEM!3`AOpz5r_Js 7~tP5GdA>{#6jbzNG}8 zJkg(yvAjL?sIR7N%XQRTLx ^s ios7eE~po-eL%U0$Y9V*Z;nnChIcs7q{dk zHZ1E>dg!Jva*)S(W>%yVix=r>>8%4wE>AjO^>`gRB%tz!r}D!8@eb@fr2E+3oP7s% z6kia{KCS%S^P7JUwP*g(-##AF_)RVz`ooV`>#|7X82bSI;C`w*5*&H|!y{U^imt)P z`FfG~AmaB|+VcJlH8i)CGWPuEqf+|1&gUIHyElbS4G^9z?vK%X5F|FwwHB`bY0Mha z{^6J7PuCU(^!Aw%8 gY$b+P)7FrbB_~$*<9Aat;ot(+sNU7@s{v|;Fdr`7= z*?)WeM-neZg)kkO {0}rjd1a`;`1;={Atbw``TswYY3+b^#s|ur!|k(Ijt3IRQQu zq1)VvP7M||;uh?2nVfg`&xE$$Dxf=!93#;^l_dUQ7=Ls~|2gq9NmA^OGGX}O-@HzB z9a$=L*+E~sUa^(-hr^ydyZc3~@U45|Vq2q&a`o+$1j2Q YOm(l_Xwruv zn}D#{FWoJj^B^>!c0S}HpjhG^mqn1(HSnK5cFd0t_1eadA0_o7T0$BVinCT=Z9zb# z1RO>}Qz~%U+|Ek7OG48xYp)v38vC@l_$Hr%4c{s7Gm5T^Wmsf<6J49cof()1j2i-= zrMr(U^W*LeT(rlwAA_Q__Oq*IX{C4flJx{2oKC`=7P>f(m_YbjE!8>}@7iu(vMD?7 zq13*~6E~amj5%F QZvf@f4|C-FRH#xdUKNx&~$xz%P9&u z$5NYBF&`}AC tgUlnOTEC!Zy7`Ecug)imdf%T)Nq+_T-#b)1cCH{{Nr zFQC{&NlNA!(J|XskSyzVbIUHUlEkCVf4xH8v|mun|CNC9ZQBP`yWk*aA-i~fKx$(d zKBj=j&vz1)iA-21n$L_e5|K#@9}*J_QdMB9Z{80=)+$Z