From 793067f5823057be4ffa1fbab56a691562d2f502 Mon Sep 17 00:00:00 2001 From: Max Luebke Date: Fri, 5 Apr 2024 08:33:26 +0000 Subject: [PATCH] Move old benchmark files into `old` subdir Move new benchmark files from `het` folder --- bench/dolo/dol.pqi | 67 +++++++++--------- bench/dolo/{het => }/dolo_200.R | 0 bench/dolo/{het => }/dolo_200_rt.R | 0 bench/dolo/dolo_inner.rds | Bin 0 -> 23156 bytes bench/dolo/het/dol.pqi | 38 ---------- bench/dolo/{ => old}/Eval.R | 0 bench/dolo/old/dol.pqi | 35 +++++++++ bench/dolo/{ => old}/dolo_diffu_inner.R | 0 bench/dolo/{ => old}/dolo_diffu_inner_large.R | 0 bench/dolo/{ => old}/dolo_inner.pqi | 0 bench/dolo/{ => old}/dolo_interp_long.R | 0 bench/dolo/{het => old}/phreeqc_kin.dat | 0 12 files changed, 70 insertions(+), 70 deletions(-) rename bench/dolo/{het => }/dolo_200.R (100%) rename bench/dolo/{het => }/dolo_200_rt.R (100%) create mode 100644 bench/dolo/dolo_inner.rds delete mode 100644 bench/dolo/het/dol.pqi rename bench/dolo/{ => old}/Eval.R (100%) create mode 100644 bench/dolo/old/dol.pqi rename bench/dolo/{ => old}/dolo_diffu_inner.R (100%) rename bench/dolo/{ => old}/dolo_diffu_inner_large.R (100%) rename bench/dolo/{ => old}/dolo_inner.pqi (100%) rename bench/dolo/{ => old}/dolo_interp_long.R (100%) rename bench/dolo/{het => old}/phreeqc_kin.dat (100%) diff --git a/bench/dolo/dol.pqi b/bench/dolo/dol.pqi index 408339482..f05425129 100644 --- a/bench/dolo/dol.pqi +++ b/bench/dolo/dol.pqi @@ -1,35 +1,38 @@ -SELECTED_OUTPUT - -high_precision true - -reset false - -time - -soln - -temperature true - -water true - -pH - -pe - -totals C Ca Cl Mg - -kinetic_reactants Calcite Dolomite - -equilibrium O2g - SOLUTION 1 - units mol/kgw - temp 25.0 - water 1 - pH 9.91 charge - pe 4.0 - C 1.2279E-04 - Ca 1.2279E-04 - Cl 1E-12 - Mg 1E-12 + units mol/kgw + water 1 + temperature 25 + pH 7 + pe 4 PURE 1 - O2g -0.1675 10 -KINETICS 1 - Calcite - -m 0.000207 - -parms 0.0032 - -tol 1e-10 - Dolomite - -m 0.0 - -parms 0.00032 - -tol 1e-10 + Calcite 0.0 1 END + +RUN_CELLS + -cells 1 + +COPY solution 1 2 + +PURE 2 + O2g -0.1675 10 +KINETICS 2 + Calcite + -m 0.000207 + -parms 0.05 + -tol 1e-10 + Dolomite + -m 0.0 + -parms 0.005 + -tol 1e-10 +END + +SOLUTION 3 + units mol/kgw + water 1 + temp 25 + Mg 0.001 + Cl 0.002 +END + +RUN_CELLS + -cells 2-3 diff --git a/bench/dolo/het/dolo_200.R b/bench/dolo/dolo_200.R similarity index 100% rename from bench/dolo/het/dolo_200.R rename to bench/dolo/dolo_200.R diff --git a/bench/dolo/het/dolo_200_rt.R b/bench/dolo/dolo_200_rt.R similarity index 100% rename from bench/dolo/het/dolo_200_rt.R rename to bench/dolo/dolo_200_rt.R diff --git a/bench/dolo/dolo_inner.rds b/bench/dolo/dolo_inner.rds new file mode 100644 index 0000000000000000000000000000000000000000..18740ebfc14e12a4569669dc6719965168b5e6e8 GIT binary patch literal 23156 zcmYg$WmH^Eur2NyG`IwJ2<}dBcXubayE_C6Zo%E%H8=!!ch_N%dGp=3?t1V1I;*<6 ztGjB~uA1pX8jAq&KM(X(58`kAF=u=M-=JZXADa*ZZ&-SH>4=TPj(gl{W}WCA^7|IF!a}_v7`s1hn*+9h*Xx{aYbNU1rUZD z{wD-{J)Y&N;x?cu|9##}SdX+HlBnMhf9-E{qSM=BbfTbO5HwFvPtw()mB+{YT^mT9 zn&W>ZoXh@IYE4FWpWrEwm*J(k-7jWfW$0DpHAVK&zJZ-xr{ccgT$E3^K&{Vz&MSu> zef6Kx1Hu1SJvrA8#C@>ShO!HRX3R2$vR@0HaAW)V=(pg%$DtJPb43sCf40|{E5P== z?F<8 z97^EtVsvhW8uj%_P9PT_=d}Xsk~$S|aaXBn^6~7rLe{Kd+hUN&-$1*G&gx-av$%yI zqdmNDOj7EYsP>u~&^}i|mXb20vx8BC4B=RNdkh?+~Z>SkGh2wt|nPfMHzo-6mO^A zVOhf~>AOnA%JeA10q$c-{+RG%Ea+iF_`c)i|Av}q^)Ce4ZolGx(e+>lhp#>bj0i2{ z|4;4|my6<&4@RdE@lRu#|1ib>6%5S(H_?dc|C4vKdEDP$3;o-eX8-N$?PhO3#n+{J zqkUV~rP{=|CEKKHqgCT?$8e>4)iEO0vakieb_HD5t+8!MuQqe~oyF^S{?`mQb&CAs zUCpW#U8&+`dBn$4Qp4@|cFpWzS+&YxzH1#jYiGjJU2oZ^tTd*IZ)UNwspkCalC|Pd zPEFnNegDk)I^}Y2N}A$|u7k#foV)QL9iLlV{BskCuM+f4cD^W3CrR$WP`iU z@l=CZT{bzc0J3SaJ%$8_cLvhB2OO^5c7VvAkch-1&vCDFG$ zM&xxnb%N!Gno*!_kIsN<+R3tA%~QgXP8HL#9q|&?IGkIJfel2>7sMpy){8B^?{@BX zRdzuaVzHe(Z9AfaN3!rjL1dMhJvOT^?L%^7(xIt``R{_gbR|x(#S8_b7u*dAT_ZbAs$?kiBn5 z@}efOUB?>NbV*nx0dPhq^}9rVk?0=EK`4>}k;AUHR5$WJqfH$ve@!cDP-8}t3YHh= zWE;#_T4iWZ2e@^I*ztF$1qo|phRMGl@`TqDrMuEJu3`O#pF>6VnE$J=(zte(RWzE| zYeefyHBOvA96=d3;q0D}J7pmX=3e1UIbH~tLrVAAi z;O2(qROD112mF0P&O^^k;yyXfWxN-HCN1;dX;S! zv!AXl=wqPyH0W^OodH6StQSL3fDUhC&zVVvfZ(2LWOul^_&ejlevt%{?bq%~J0VXb zIMm8f73gPyxzd7Eeomn?rD25k>&hFnK_IR1X0S1;1Oa2&N>6xMElu=#)L<4 zxk2$AX6P)UgOXKSdR-esIBif&`_bh@L}kWix<4`4=%-vq1wBG%B@Us`IxD;3^`p1- zXZOT#qGE%3{m|wjAGX1Yeork@CspxuSb<(p>pWkTY;elTY(~@JNa?T-=xFTJ{TnfnZ~m|p ztHS_OTz532&xCGa`%2`6!gx7#G-uQ=4TD*dmOu@6fkCDf6%8XR0#uwX`n+n` za5!1bZkgu1so6hhvV7IA@6n4!%7PWU>Nr(+gq_3Qj99A)t@n3|1-efetNvgQ-WV=C z!xN_=Rwn>s_%i-UzLRy@!&oe=!?C+(A*-Ibwx7{s(`VBKIb4jH-?C@;SWws#O45ZQ zJ*sN7$ofr1X^8BaZjLX3WBVYsDN_JVTV|1vmCO;#ZF|yJjn*m3_Cy60OV}BzrBYTMZ)kp|^(4W8h%E z8mHDj6It0)wBV)IMT|p5T+!Tx^JG&B$6Vt-*ErZ zHGd zx}ipG_eR0;)1Abn83Cncsf%TCy}g%AQ8|Q-#&7`H59X*S9->bb&fmL9Ty;L|Fs2@fT`2aG>HvjP({jG z{)aYSx<>;xx7O&q%Uto1V75JZ@lnO&GlVm1O&ZPuym zH)gF04#hBCT{4IaSRIKI8Ug~fJ8S}>!A*Sb;{vCVy2UK^Y&MeC#;U2!k$@m?24UnO z&C%Y-x$>$H=9wWXi?ZEsvOiK(*a$1>#kN?wjHS2eO%BhBASQ**p{@1T9@Jm#XWME7 zrUoU;voyHM_65t8M9tQ*1F)=K$0G;6U(z|n8zLWITyZxj1v9}(DO^k_nvtewk3wc~ zF;S3^C{0Q5K%6o-SS+%>Q!Thm4T55mvoO)kIm-NPD}qh2q6@^OQ#WFr_I%%hrqhrzIZbM&vUO~CMQd02!9nWae6|EuERhlTVWHJ6xadAbZR7d8oLM4HGT9C#f@UW6z zCH&*NnWsy;sz{F@3q6|PUtz`|KfRG=jLV1#YWu9p{_#{>9693v60hfRlAof%@3k+9 zhH>MpN_&-gq_TVx7IcMzFxKBn8WTw7$OLJ4SqT}{=N9$?^>Z+wg_`kN(C;#lWSdLZ z6*tkfJ0X}p6Js-E>hy=o#Baop5H*I6OfQX=Yv-Ee)?3A9;~By^(WSJvEPD2%rz$Ry z;njH?nu$nZ)N`PYW;0?&WIb_K#7`?VG|F=HWf^0havZgo^;&8*Y_R9@{7wqM zB|r7Zdl$(;8*U#-uEzJcCGNO50qKoypuPAt%qQeNq+Q(XZ!Uwb%=Gl33A>hEoov%W z`(_P`X+826Bw}f$cSmK$E}z4XFV%r9xIE0Kkyj}iVjlQjJc+;x5_D5+6GVIn=-7I8 zm?OGSZYasY7PN8hu=Ahk5Xz}d!Dg6!X!RG6Lr9BrRo`_dk?|T|U4p8VtZSD~HE}ij zyS$CjmSMHSH_1wVx-}(~EEmZtRXF@q?X@P5FF7PsFe+ibu!bFjCj7ndugh_0bV;0+ z5b_nyBb{JxAk?J*_UKo<94VgfVfb3N(_8-+zn_z=D`c6Z?NqlQ6I&UWbi;J=7%dm; zq9$Y_#8=V7N&zatF015;9%r|}(#6%i3BhwU2nHOzXGFLWrkBUR$YeRF&@DO02D8vN z>C^m>)$j#N_9W?Vv`#xSd zoTnCor`Fdg7=7c0UZ)1}hA=!g1PA$#MA!`XG6!wi^cadiBqL@H#jfG4t#e5DTkgmO0?HG;RdP2WKkmeYfp_wL={OS{^q%JPbQm zzhIDuwVJpovpr<--JDMoUHe-j;OVK1GKF0XRt~V%K)_8f=t>b<+H1I!B;z@tcnfI5 zg%@U}$_5)shFK^%ba}N?Vd)ecizddnG-jbbI^|WAJF@Z8rbTdJ>b*`E9qhM;Ule2+ z6t`JA3HRkF$q4+F5f+My_Zyb(6{cKa{MPbL)|K#1e!Wh8Nl5HhK|69hP>uN*m)G!< z4gX;>IHqk<)-IUsC}vKQzcr!%CmXH6_XbNWr~Vd8v=%1TOosF+M9Q0@=h6M{TJNk4 z?k%`$%@+V_5-(4}FE*#tU>YMmUoI==)~S|dR0CUw_6de={R@G0**Q8&YAIUcGXTCYwh;TFIu*Z|;OdRW-(WqMRO9?httjU?JXgE}zdoKENINq-$ zWsq!fYFt%xzQMIpH$QJL-LpIj%e}wz5a=%_`~(J)4H4y(!tuxJs! z_A!maV$d=j-$Q(A#ysg)-OebQ9s(|&nf$n3OpEu&GWsg1b2GySE6BFmy3Uk@f&)d< z!9SGKk;V4%-7(@<78Wi7}N&)R-;I9@hOTGkuq}H$S>pg{R*A-p0;9lGqGX06@xB4Yf zzOxqe;Q@hO( z^I!;gk5pg0X$#9#p%jMW-b=P@>4y$XCg{Z0EVW2osn$G1@eQxDIJV@}8U+DjYUL#e zSrbmqDRH&J=gC>fP38=!5b#S(8_GbVK2DV(487WLTtIKiXIN<3l4^jyz%g?r6{-A) zp)#7F_E|VNW8DM^^QZBE3S$WqaQ}k@1xn+Fa;ZW5JC$;2QB1#Vo{%AJsw}y=$g0k> zaDPbIXhczY&?s}tFez`cX-G#(q| zHcS_QU2J$8Q@Z4g{&-LEnp-DKVz+8wZ$ zert`rOcz|00o|K;Le9z=WD9WNP{|E%5}ai&Tji)F{{ywl8*9c41BG*M3JET zNt^L^ZnO6H48z!`2@GVOW88SG-O;9t&~jrbh7kTM&J{pC6l$qGL7YM36Lk$eqpTz~ zJE_jFB%Mm#%3sb;N0{Gy6eFJ0qL#Hnyhc|m#wJNpY|C)#{HCT6suswHZ7*K4dSPU+ z@RIBE@+yDpMZ-$VR-)Kld)3J@wZ=z7aQ`r|blue5F6P*)50Z9hN78c56}t9#{Sv(H zdfHRRUg7y2N0?)yc@!LA$|1e7yewE##Q|1*UEFDxxNwPv?7j(gmcTc@R>6Dx{-Sc_?At7sDZhex~xRC%79 z9;9O(51mJ-=83|KejJbMqty2(Tgw0n!_5%*2Q!-(isRQG+RT#ixe0v6@)YF-^S5~k z`0O+&(%v`~iV{V$?{=e`%GfKzKS$^itL+~P#EiXI7z^ph?PBa~^<8MT2MGU&o%G8+ ztMF#DQB&l)QrpCnXq{wd<9Nj~OgpQ(*QyZTnB&Lt8r6~#*NHIbRQ7Yw-=5-&c+4`7 zBlUk96neMjOzu(pI#|%3AC#TGHRONZuO(%n!#j4K{A>Gb)}zniVZrg&QpGBQmGiT# zDj?U{i<{yV$*_A69l;LF3hdluo$p5x!&L!t|{hSsKxUwPbgDpM)e)TZWn{IZf3IJt`MvLwl} z5db7*SdDIE0EOi`Jkx2Lifnrj7P56JvZqEl=~XD;>lw_$UHr3?SXq9V^}?goytL~SRKR;G zB3xXQra8CBuZQq#Etv`tTeT9!aqmWPtqP!KB<&=>%^P8u@H)vPaF3T#3yCZ+zC;VF zC%eL|gRX9;Qff2HAg+U{WkbI}6=rl6&zo1)P1?hsl$`y+t2YbuqU(}JUGsowpz(vS z!X1ZadDnI>GHZ%pQ|k;LhN5p)$aXCqBOQ3Den9xj9B?WKWSr!IB^8%-W|G{j zbro=*-eiF@;kA1cg{X?gpn%&nxepPsv5V^1FN_-BNku9w8?>0apo< zpTs}-cL#?$B=n@G4r(cJTrW1-&EY86X=+h^6=ZA_T6uIWHK)c(X>%v>!#4mC&7} zl2s^^9;W9871x9qsE&QX5L~6=^a3WN^XaI!;bVnL$4Sf)?iuR1g= zm%DdHN0#5}gtZ@~JPjw9B%#bKqCWoU)b4P+ti#XnM+$b*E$zjC{dQ`5x>xKTwR5n% z*pcr1fReu=YoC1tHaY=plE z5SP^<6YYbY82&eyi84T6PKtdo(I2jv7K#$xr;;;5N1>ufSho5LkbG9Dsl_aT(A%FlM$_bo3*rX^5+s(hTf z`AM7<)}NWA{DDbYP=YZpY51S;Kf)Z2Hq9BTR!xD9+W~2jDXb7vQNhczgC$j^e5&Ho z!=;Dnzg(N`h*~q8jnT?p*^6j++jj(=b+9F6EMc7_*?kSO{WEv)WCry#Bo_;_3BV^hXj$z+s>EzUgp-;-5bCS}5(!R;kxmDH$#aX~sOy z5zc=hN`-2*HulaQNVp&yOufE2Pr>`*z^Q*qf=a?lmyQS zRZ1FtxJWM*u964Xw{EVoXtFp+YFQHt@BShci{7%eCSKU+*H;>I_0s@|dzM?5VquNF zi?`d%TE`QRb+Kz}!ZJEdouw^&a zT4Bq2@~iPm)|ZxP)$^^r7CDSyY-ppK)r`Tp zX+8_Qos=(@kVlc<-8%{{m-j2*Zf&a|k%ZowIUVtqHoDW<>;`>SZR?t=oi}$Q7-@Nv z`8R7B+p@GqsxIKDn*utl#)5w{rRo2**B&OeSTWb!GTzHbHN9w>T(7Pc_uI&GgyNr< z3ped^y5y+xZ5+L)403mH>@}=u(KtFs$}8C+8y$6oIkITt+iYi!InTk+%&^kbU0U9k z?OIZCuSA9dj`Py*uzdR${`277Bdsyt>?yb#5 z1yS?pq5{`u(A)x{KQ`cda!w+~rrZW`i)f}b66x`G;Kffct_f-i^&n{8es=&m`wRXL z*CUq}1zWzlncZ*szl8ru^%t%QQo=N`Fam7!###(WcjOBl^V5@Pe)rw)ro}S(wIK#n zcW9P$#jtc!F<0C6gqA@9zHEgqCrjD&tS!~O*RU9#8s`zGZ-y@KEroDa?PTcHzu$}g z#@&L7?F(Y7zS_6@ssS2cE(26Ua;kW$k+_Tp(Z~oWYP6$9#R1lN4U{y_cIrZnguG@`!S^#As?b z+JL5W->4|?iC<4RbhBGrvPipBD78m_A-Cx?e1Mc+#jj#h;fc*~o%7nsaZ+VGr^@0H zwY1bw)CGaoSNCg;i7LO`?e60%X5Wd%RT+IG;e)$pt^M8`VRQA<0nYw$(gtO~$7*^U zmjCZQn`hxo`<5;SbBN*IYGHQx`LKT*Fx4OS2XuUXPv0hK`qb^z`R&p?n0<@d;Av`C zA>A3)y)clZLeX^;54A30aWtA@H8vy--_ErRl%Z?!wVQmpL2BPOo6rkQqYpV+7L-3t zHlNW#JExf6WMwwdJciNKIfSH*HkH4x2dh)v9@csf4UPGFlL&`r=7&RA;jG(#FXB<_ z;lW=I|AChU7n;;#UZoc9)&A2boqgm(pwNGmnR@DB)J#iAh*6%qXEua;z)-v^&w7zE zjrzlb_x|zMpt}!=z_9ef>+8fpe8l-4;?KIKNv}2y_b+d!lgwfw!~Q6X)}|`Bo&%hG z=h9me^lM6PX;Q)i(HT0|3|Eu6nr++VDo%BSQ3e4+>@{q}2NJOSx{Ao$*|DO8R-4+VTzefk{d~UxUX6IvWn~{NQd}R@nKeZB(26$OVlfu zbTxYHyPCEQD;bv`-;ns)(Q;f#*dF6x_!-a&&1$jqv|szjxWjBlt>`aumt>oMT*zr< zTC9~X9@$R&dYnsImj>@9U@!6O*)VIhYZiv2z!jI04E{Oykwt4iq|)s`-AKe}%PpWm z&Dt*`z)e+u^G*Q@p5GCpJ>+Vf*YtH2&+(G@u+jXQhoO6k3S1oCX#&j`rzdb?6~Xpn8?C9a#|ux%|!2hH2jWFcVN>0n<||=-Dez< zR7>k9);Ockr$&}JHumfbiM|L?HwtyI(~jEBZVL4k>}Y|wweESxMC8VOVE}(;Jo@ql zWo5U|8%Zw#lwhw9HW`8eReKY1te5*Z!)`2UalR<8?;IcG_`j%TV}CugZ_kgS`tnc# zjxDYF1JPpx@CBlU*;`z<5S3dFEqA$Tl^l185aS2ornlC9W6MKJ`6b(fMfga!_{m@U z5kk>ycys$RjNe^MtRwLTejFlW(vEtZO$v#JgQpJaS zDzJ!#b*+MQHdD}t{WF>1h^p7G=igH8zDtu{2aKJUU&Rl*9y`$~1Q){!_jBT1Z(gf; zQjVFOO#jaVM~R)fUhb5Yx!kXja9ix#oi=9=!A0J)vqRZ~=_FpfdF?DV4C10jh6{a! z5F_G#XCZv^kuAaOW^vXNt)$4C?i%BTUz>^l8XMZ!UZ)=4rut0qA1gB@U4ehtYLY++ z&wAcL@e2J+n5A3uD#lW0XN&+qd9p-Jje9`NbFuCd4mykbX^{P}GW7oL@4}hSJq)9f zu{cuwI~GG}NpPgVNL7h=&lhDPU|OGX9L*NXP2ncae~cqEL^@Gmpx10lSU8uNN9UB;OkkGj@HdrI~E10ZaAE zm5%g2We@S3!|m=C(asn4F$P7_^;K79jh<2G26%<=P5`1MY~bh zJ;`#c+k|8o&8O`6thuq1dRvQ>RPHeN=fh zft4Y|Fp_vXW^UU~0zW>@+{>}0Z|wtJ11|7b0`NV$sz&a?d}7t&y;|gzlgUGBioKY` zq@&C3M~mp6<_MY#`{CY` zjut}`YJJ7KaL^V9`4dvFqATc}*j=sf%Qh|hFY@_{eb8DuxTD(ftAtYNab&cDOUf!#DqpAKxE7%QkP0OK* z678r_JbKNd7M`oA3z=TbTGIEE+ar*Vvgu8A7)QM%G!d`$nL4vzE8;$I`qTYg{GOP` zMh3t-9nbnrd@Hv4T9u>L{D<=zs=^uEO77Nbl>A^-Ymr%5!g0}#gl1o=-vwcjQMrIC zz(%>xd%wSuaiCo%7g!0@5g)%rnmg!TxLy?2PJ$2uj_llNy45p1{GJQl=OK5Mou zD39$4E_$WSdkARrggf#^A-Aa_y#XZ5*W7On`U47(Sn`xkbPbSQF8kqychhp%tLMqW z+HN(901OEkY6!Qv*(F6rXvB3WYC-4E4^K0u3@rlN_91B{ghd~37X16&#n=de<9umM z`XnNeTc&IXv3Z2I$N$tzl!Sc{27giX?@(OJ8D|J=Kvt~hsnPS=hFJm;kph97L6Wv> zOc{cGe~`-Jnm(JBSWpHj)?YqUfBbNB?=dnN>C}7n_xr+&;>Z1%RM5EpC%(2M5(0<^ zSAguO;{Lr)_sk8KSyS0gP3~YLvxlr%W69WhI{ZYh$2RK4JysDx{FpcIYILYYIXZxs zXP_N=(wbFk%+qg{iD49PUhP0BXaWD^AN*I-U5r3~giefhNy)BYsokwAQ~bdoQ8`ED zZEfP)L-^EE6iQ|K(9ep6F(9_)n0NGDNQ;N?^Z@ZmKfGoPRddKsGeNb!!N|l>ya6|H zqiA(OjOD+?zFTTzN+iNA#{jh^WZ*t^N8!I+vHNIA8{=RY$82Km2=+eAGbbD9?pfjdVGrxbt~0u_W3TV?HU0HJ4}_94J2;eRw&cY zF8SvDm0I&VOaskZDohN)Zs!#)gAdQYX2VlouIaFgm^aW){VOhc$?*(rTk!0VTTYhm z7HL%HD2ksJ@BJJKm+Xhn2@4MpCDLi7hz0$1Q3gje-^8p?XLs`InlI8G#!F7u#HxDy z!t*M0y2`9DxvkAu#TkD?)l+sq=lyUY;wpbJ@A=1Xy$T=RZ`3lDZ3|r4;m2p9qx#pU zX=D@R!|A-Z&iUOWD!JLROuhUwXH!FthBQNm>?Lp?P+?}3sMIh~)mhZdK>r11u%#uO z9*}9ffXIMsXCz~hIe)^i1gkM*d-N-MR1y^%qDx2HuXjm}R}Y~`9rbT=eQc3PmBH2B z*ykwjXspG@m99icU;#0m;?@Goa4AJqh3R~6w|0yyPF<;b6I8NgRV%h+=x{*MtIQJ~ zSx!~8hf9&rPRyHbWVnQX0|)HF+Undw&k&Bi3~S#2goIs9RF#%3Of4C(m)$C18jkzd zD+8U|QRW;6%VKt{!@U4CZ``QQgp9|kX_MwfMm05UC~&a1 zj@t?{6hu%h&f0r3B30NYHb8K|)iN_~9$_=>++^r#9$^U>BJyvyu`O9JtetODK?y#6 zxc%Qx_HTD^0jd6`u2LUQ{XK}|uU!a^9OU1T^wD}`ax9(5%7Fi-5AjYfEH!li@37}n zrFn|x3nn3J3WOZ4EsU4L^aY`moj!)-mX7DmqIMbVnV3XM=BO4vc;8}2UzkID-Sb5K zIh0wJ2ssTpfN#iK#&Cq)kwR+hGLbBr9HXY3t}8>y5uPN5c1xY;`E{IlS|qh%Nj@Ejy4sfPQl$V>j5g? zVmK5;TPn$zehn3#QZ5(!XlogSL?POVUy5Ejv_btoZghXA9hQ%cW$ioQwoTCus<(CM zv2(}Wcxqwhe34(&FZ+b@FmWcb5SW?XFIJJgWrd5LLNvD_ses|Y@2puSGrUlSR$J|- zaIj-yTqfAPsmIwQ6SWN1w5rT++Y14R8g#2bkB{aNJ=UBI2x6VEhABezwCNh@>fRKf z%M%}WCAauA>DaN9$MkcrB=PDTqj`+Ki9ozzpFE2$z2-rYwKhX}O2LvX8kK9ge|;%9 za`AK&vsVwG38?x`_G!=_R%tj1AGeNK%#2l;Z;ceo&@zHdmUVPINFncB(2cy5#vEg8 z(^Ekmai9P%Hl(zNO89duXSq?2)WtrT6`ogXj2qAn%=*GSD-%?t~q>gF|Z2KBQcRp>?cW zR#;-FoKb3unR3u-m_#Lw#H8P+XVqF6gGqAcm%nn9{@=@Zd?Jt0@M`M3D}I#35?Yqe zLml5~{5Z#07j*i%7%3tAEC#gcCRx^REUB_3p;#Iv^)Lfo4&CeYLF0r6R}W8?8VI3YULvjl@5ZEmL^-8%CKLB3-y7|PzSan;Sbc)y(2Asvw{kPA8a{ybttRjoTI_jjXX z)BO|13;Z2xR~B4#{zz>tE899~UZcoegr%8xT;ZEK8I8+}BWH&fBVtsZz`YaSg0ABD zLqqDwEq}N?clm2&s$s*fd60=fJ4|_-=21%+(!FTk3Z8S9D06s%OEYQ<6~#?_RUDF@ z-4!k)x%U}sY>Q%ppiRZs))ti#l{JYHvc>JC)A5F5(^36@+uhsUbM^-ImJ2I&+)Mob zd|q=VPc@#EWh)p-?gv9G-BH*3d%_Bb`W|LxdjfB2kvn%;Z*ahu&t)gaiG3!?T2|Z@ zbz5-__)14s9LzM%9tv&;n);zE<&uhf@;-uEDRHfidX+|%8cIvd_XGjrnAIBVx$7#w zqn;E;u0WdDqs=-)h1q9s2%@!*xPvF`K4a&-*go&sl)9r2%qwU`q_d|r3E8xPR!|fH zYEn-A-_ANdvta@<5F?YGkK2Xe^3J%Cfm?`P7~xd6ZCKqY=Q~OQY@d!;Xo4_5WBYgH zZiNPRmSyWI=JXJw-W~7&{xGE2?^n2h4{G)ID@*`l2s}*S>PZis*RfMwLB6~Q zW1#c|KmiHp%iaT|ykkUO9>5$ReSQhXGxBp5EPk?%B0H;C*!dWSKi-i|0w@Dx`JS(_ zeH;VDcH6$B0)2xF1i7uAOC5+r@;o$h%PG~vV=X-9tnRVdw7unE8VCBlhSz6 zV_l5dAq0Z|%ed7X0GR~uP>?+YO6;~_eiJTuv_o28o7iChOPbLI9|GlfZ@#rN0)!O- z)5F4I%SMNwuzn|;q*r0O+YMOrkKF0ufS-{;)qS*+j!1$J!Zf%1h}W1;`4Mgr_fB|8 z7lFiotW-}sJ~^QJGRWSPQkMk&l+o^dN&1*w6MEDKsnK2OK_#|xt9=SVHUu>HxI%#- zM3bS*x)R6Kn0#gbCWkiPz=BHX`+K7(+Bcc{oVKzgtNi9yw* zPN9!!sexVKHzUY1V!Fr8zB+6`57@?bo>-sOT%Y)@C(QueJjnA)o|)3uRnrd;t=UXn z+zAkq=H%tu=hkE9F$}E5`>lJ@w%luu`TOltQZHLL_@s6Jx%KJ6KdF}-cJ^^6Nwse{ zkkY@DS5CB+^cAEjvm^Dof!O#-mVZ;`@&gQt4VVBt(L9i{@TZZ?j^e8?I2f;Zi?iPZ zS3p;pJZP3+_y9k0q57mnmT%Az-VH-M5kGs87haSj2*d?ULn%|Y}^Ch ztNSO14<4Llgq6JcVGi$Jd1d;t#Ae??-TgP%Nk=cpNkYjXRzM{LpaGN%C{P7#BK6=w zE=ede8clv@fe(*3!_mOh&P$d@Ni7@{DAw76AKJ2=}#U|88cgubUcu()Y`I0nt z(njtdbMFbxn+N`+PoD6+S5RvI4JOz|SCUZfe>RYT8?c9;puGMaaNCKxob~1q}Ay%KCQgh%Vy}sp#?I^roDc-)kKm}$8P3@e(2G~(AvmK** z9vI3Ga0Ifb@-9>${CbaUC&?rIs`MDeaeVUa(-Jtv)y2&m+D;yY7vZINJPd??0AW13 zpAON1L5>^T7fV(X7(u2<93q7N!7xFJm2~N3o!u`ilC;b%e(P`D)qjJRH0op>*{_6= zgZu_{Snd>mQTx1tGX1D;C29oD zxqiHqHM1O}52F$n_$hV^`;@g7A_TSt{duttok*m)(+9N|fRRbE6a8NjLPO=_JjO{A z^UJ@!U^Yhy6Lo=MfR!|c=Vo)tnJO%mv;yMkt(PKhqyv#ddi%vDO@QH;&vg5SD8x9r z5`3r#LjqPOw9?@n%7t8_&l@#YOI7py3v1D9p zf%FT&KP2z24)Z|k3z()TMgpilOoGv2%|XuCU#`jT)^W&hNQ(mGt!{0rXLjhxE@o4{ zcC&5B&0ftSymN(61Q79|owd+;l?K)nfcfDCBq94oAtaV_P73@~fTH-=fSAI{I~oL! z%Z@v#A(8x?0BtM0B9wsqvHpNvWeKShFpK1;DKgskRW z&O!WK^)mqMSN^2!U*H^chYiRnU>fx#UXYyd1m~td7#AynY$g8EyYpk+x%Yx0mqyNh zJM_0y;5qYBeH36y4Gb?}0jJno+J9D^@BnP0vjESiJ7m~H;7+D#->=Uq*iRM6HBbQR zAEycH+{rJgpD%XgJ=ynd;4=5XEOVn!c=Ei%bi08d3wl!`@1X%x;0-1iuHeLMN__@< zjpa!xZ|eu>Ek4YugEjB0@5g^${ror#yeakz>jp*N4}5gp6uj%fn7{h-|GAL^*M&au zT_>t`g>m$k7oj9@rlTVOufRmz%TV%h74rF39*8^7YX16xvFB~j?^Oz~-yG1+0E^CJ z_K`Xq_*=+zDk|v-{^MVwC=hXCM5@VPsE8f zu%F(31o!2Wjn%dH0_lRmRVy8&g#6wfWS>t8Eime7_FV|U8-)4>`j7L#Odj)bR+tYB z>;e}@kTA<_7mWGG2dqq=%xf1o0Ry2P8_HNZcJ{35rGb1 zcJC7fL`L;RvX6f5%$30eM1nPLklvZWG_NEM?7zFDMq zF7UW8Hh>ts07&?cuq>Fb0>7PtdzAMNw!N`v;jlKO&rkW?HZT$UEAF--j=%4HZsnKl zf@k49{#Q*t`S)m|`l4Y11DS5YF#@yV|7wB$TMQ&O2ksVL-*&=;Mc#>JZeL*krzv?U z>Ia_aH!fpjQ5Ic*`<|Hvb5ULRF<){*di3lGHp^+lr_s7Ns*-)5}}MYN)i+D+Gn(o zC6R6HW6c(0FlL|gd&ch{*IehE>v7KW+~4p0xu02h%1>EpPe}=X@`a|<6lv``-)+7r zO#wJ+LwwUAOQ(aS!q@;wHL)LElMGl~+BOC5xmWX>4qn<&(2ZZ#Q&qKoF=DR}w_)N0 z^urQ`j9=Mn z&WZYo&(gzOg{AIUqi(hzidIBPqpU-pg>P$LjXfFh87Tsv+(Wb(m+A7p6tF)u76%{6 z+uj1afJ%6C2bIyBSZe@8ikD-Ci{IzQ5~q%hQp{O0;-J(nJ*3#<7@7_YJFy4{(=}~P zk&)foB-lvfi~)U2ozfIOtE3d6lj_*4Wk!fq5S9FMwF1ZeOtL=m3tS2U`CKjsK5 zl9DA(m4iH^NvGC@O4z^T1{C6hOLJWMNl6i0vQl9iif69qpXMGPl3iC_p%8$LKLqiD z7r;_4{a>Lw#x{=2{2PagC5vDLj=Ca55}?MOJR5=D!9e>|S3)wD)CCvepO64t3*`A$ zASqV3VIDysJhb|ZfR^hXEqRR-UX4b2p-)0i1rL< z*WAdRV+eJ=YzuS6Py~|>1{W0ueEE|SV#Zz+RN?SRY(R^~KM(kutt;l%X)W0;S3a$FGkTvJmO7{6cvn!V9mXgoshUNYAjU^Gckx;_W7U$~za}fr~z+T=c z`6(U2z+%)vHS*$fP+%z-R1v=7mwaJ}q_JcV=A2hukFr=g8;u%5vy%jmyfd1A5m*XF z{}cpl3vXS5cnTlspPq<0V>z8*O@a`e-dZS4)p-W=WECw46~KuP3HOu!<@5f@L*;Q_ z15`8_j0N5KTk2W^Uw8ZooVgsLpvY$Gi{ro)0l+}no&1=3`0Whkogx?A#Uwph67y36eXWb)h)GkWB67HV zWC!KBc{?30st1J7N~Laq^rtAN{nKy4)ZLK}0d>Z$R$|z-nCr!{*@JL&J#}|z#e29x zF-#W|w{r|LY=dh@ED8_d>M7c){fs7^yktcqGHTBZL)>x#rI)z}VMVeU^;@^_rr}f& zjj~BNxJ`(JHr;`j+qRNyxDR@DSvQR$Hd8NF;S2uTr=vUn$x`mm05(3ne@WU5L-ab2 zA{myxSVnA-orj!J3(N#!43~ahnhTFP$!q}?Z$fu!7%`gauU?S9hcl;6*I^c zu-Tv>bTK5j7+~gx4b0)pXp&)ri8^1U6xI~Op6i3-L~bFbyu~Frh$;laOVH>T;@l#k zQcRLr*vn-OoWxS}+K zF?f)xL{-ne0~}*zhSE@-|M5kUWy!7u^E0$!e)6 zvks-@jafTsS`<0l0_E}QO^`&V#cOSFYZa+Hi;u@0-*7F2N_nCJNjq9|Vue-w)$Fg^ z4n~?cP3t~H?NP;QcN;8+3i76VhrIxgG6h(wDU5ON(C5c(Be^*(opJi>G>%Q?A$A*? z--0UP2Pw-8*+~xiCJbI51>F|NK18)A4RKmrc6i2aohloTN3C#T3^sO9f| zZS_T;myya=$C$i@rIB_2I2sD#;AvnqT8nVvt4FIE zEkIpl?9ob)kmx`V6=1X)wVT%@Q%P>7#c&(WSx^kXG1XK%Cx&QVxDTn2q1ipsojjpi z$;%AMrnLaiZtdAB-YEw-TQS1mG0YPp?U~d+tXhG#PbS2*%WX}f#6;OrY_~&lcs09O zvcsDzakF=#!n2oVoDd2EHVVY}`vHh;WLe3>&TJznO?sA%6QjgM7ZG}OfB}K)Z>5n*+2R;8C zE}A74&nrw0u?Uly^4nc`VHL)kvn?fz8L+Vs$-{J!f_k*7tp$mG89c>-cXg7T0OLO3 z^G;Y&Qq2JZ)zwqXmiV?tFmY=q@*LPpc)*WDBS1*r={S(iSMXN@;aNs3ZsFQrPMZ*n zi!y>Sl6ih-zXb4ruD}FMzr9|lx^r4Z|HN@7?`lN`s zWzblDFk*%Iz-Eb+5cdlxbdaD?>R=d>Z89MbR0TqDH4*-_jkK_jiLfg1dkTxOu-m_? zrP97&E!bKD5SjmBZ^{NWLg=AOLrWw( zrc(@RODnpWx}KZ;)S75p*$hha_0tu39r_Sn9{A!i6F!_r$>lD9SI~|bKfvy62)h(; z#3UfYmDE9m(4*=ARFmMcdOuK*RDB#@6@V)IFbC*(8S6`Km**d$#rl_xQb}SFcEGaj zg+0yxP)X+OE`bl^NAuTtQ;j#i`*3OHJ0+s=7zLp$Fn4(== zXV=*>o_~<`9mgG<`4htf3B&JfpPzhtH9xEcifpVz>q`*t2OLwfA%64|V95Iz9sX6@#deW#R(^KXc(|dSwTO$`=99oCRqUnW2Eiv ztL?wJQ7oi*6VF$HZ}fY2>Pq`5fkqhAlzQr3tzUkpwJpLVr@*FX1#Av>N#R2==)HU; zeHu@6lDZ%S35d`|oO?&O!fTvGQ^pZ9mjN`cOc?IwnF{K^Cx#%8zTkjbW?)6u6BSvE zdpKIb90Ln?)@Z>|Q3pu`08ZSxge8MWGN?9JJ>0{eCn!h^1$JU|735kw$R z;T&&O9!AiWbiWcDYYL-LyI`y&C~Uxs3QRZgQp{x-Qt8yxg1`nUp41VW&Emw=Q`JLV zn_!9@(TaFD)UZU6mzKEMc-Quz4bXu8+v(yU=dF? zOPS=x$Lp5DR}SAVuV&AA(e{6c52G#1pb;2(DYOWge7$IA{7kiPrSN6+9NA_lYFTvkr=X0(N zRqJ#2ule7EvLETyx*^Ic3pRH_W;n)9m4nnX>+WtDx?Hzh|KfOfhO?8yS@fmCf-g~A zi(Rs1{y&F$q`pe6o+4~+ZsKz|NlOOJoJ$q%=nHaG8|j;9Im6@4kI~CVJBGGT+oaI6 z|JXuF`?HQUhCi0+LvQ)Af@<<1&K!t#;Sa6a_(W{#i6}U~;@Ezb;PaZJ=h}MUE4i>D zW#0Juz~%egt2IBbV4f1Ck~BLk?DQs7k!0|@boc-4E~Xq_(R#Yai8=oKN*N2#~4rue{|{)|JGaolx_e9Oil_3qvjAL*ykZ)$6{)E`MNnmA+V z>~*gMd}b7RHQx2pTw6K4xaZ<+GTF8GqOml6YgmlV=RW(@C;#nDHs;Ob>s!d6FKSOH zJvJd7Jhb-qiu0mM*&j1&$0b#J-}2s~a>MgCBsZVTce@q5RmnH$^z@wiK$p=r#bKj}h=}{nBwNgZXr~ADwt#`s|+ufqx^t4{Xj-NIqJ$L9g%6dJ| zC7vIkU&3BLnm6)UW)nm)5bIh*K5@}qbm%|JF>OZEoRQS|&pRaPbLB_=t`&aV+SPRz zefnSNxjgHTLvGjq|aoQCs=KIGV(QuoSqVt=e*V^EZQIG^rG=X#yB!`9kno8i(Xfu&cWKOfldyxxSI~W>a;_PawgTFG_*~lyU zZ|koaV*k|0IEFM{hdKC$#{QkR%@ulsxe=%ydGx&2t?{hP7d^*nXD^nm&bdIocGXXO zD>4inZ*+z!jQ8-bxT9N}gWWXle9Wkv{4(Qz*?D8x9HHIHV74b}aKEzqEf)!!J34PL z1O5T!7dCh;*6@eAP0)U8-y?aImCMysC9jhr--@dkvjJ{-Zbb&4^Q_w&TQ0O!xJ)|~ zSEEWQ{az)+Zw}O*yB6<(Wn7~_oNP+Dp>dp)U4QW1EVbUVzp3iX&7R&Q*?}8J+qBcV zJ^$n!6|bP3?EU9ZV?hyN^IPZBQu&r7w2e#L#?u4k(;RkA!U1~23ryd?970F>0f~wI zy)>72ssHBtaxS|4>Bu3davFtI|NRw~3)hmi`~H}jXbH=HKgT{BFjL~|!l@wAcl1|n zdQFnX6g$X0cU(lh^{~>s@#Ew1%+inGlN#^JM;zQ*TKH~Kby+CSn)4Sue7jfEuIKj7 zx_(kByV>$-f7^%h%8}XsoC`XKuPN!dpg(Aj_j;8R-KoJZ8{B$59k2IkD-No;5+vK38gXyC?ldZr_L4_mrl)BG){@2L={ewP#PNd(DxS*Ctd?z>d#;P3TmRL&p?6JEIz&DuYjHqaRP3yQNQZ$Z}Bb zy;M^!rm!7p*f8=rzdsOFbLD`X)teVm@1UTd4N$?ek&kC?7p=cN?C6B0YT zSC7`Aq0Zo2?5kqlts8k0>zF`v@Lmi}#-%jelrdY0y@@eIz>#L7IU{x2jd-doo|@Ng zh`y>HLoG=di`f(_^kkf)OfK?e>PKIUu9;am%*IUH^ptyOztd83yR4x=3S7=SB-?dY zSN7h%OQN?n!&(dH2p^Muufu&aXx~?six1V zbiNIj4-3NDX+QETkx2^+YBM=j;1WWcK5)~~{@bm9b60&^?!QgmKWFdRPJ4J4iF^VymcCyZHnC z6Z5CFzsIB0xl*0ay;Z-v4{q(>_k8f`)%USWqM@lWp;z@wGPHC;9t$C74OX{UCykdL z-(Q9-|G9L2Dqj5!V{Gk~cWChO{Ury=V_v!imE0+4$RpUM&ij7(3XdFGeloLCr+vd# zrlbAqLLzG-j4O2CZF+cvZK)U}|&56^$Ev0NKhH$3}Q z^Xn3TR7w%$WXB)XmQ)e@nNIW`VxZEUJ7pi&zsV|Ga zE-vd;+rGlEMamhy=?d!29mzcs6MeG-!9To!L4(=oUzN4z8yk};xsU3=d&-Tpvsz`w>>`r z7P%QwF_iigIvDWBE91Xs3GILDi+xi&-xc-gyvZF+ecR+l$+R-4ag{=}uJ)m4)x)w< z*Ho-@SGsR_MrV{=OnZFl--uXOVy=VMooqAp_LMN}