From b148a5cb175d973dfcc1b3ee7fbd05729cebf33e Mon Sep 17 00:00:00 2001 From: lensing1 <m.lensing@fz-juelich.de> Date: Fri, 4 Apr 2025 14:42:58 +0200 Subject: [PATCH] fixed a bug in trends that occurred randomly when using an incomplete time series, updated toarstats to 0.6.9 --- CHANGELOG.md | 6 ++++++ dist/toarstats-0.6.9-py3-none-any.whl | Bin 0 -> 32974 bytes setup.cfg | 2 +- toarstats/trends/utils.py | 11 +++++++---- 4 files changed, 14 insertions(+), 5 deletions(-) create mode 100644 dist/toarstats-0.6.9-py3-none-any.whl diff --git a/CHANGELOG.md b/CHANGELOG.md index 3683b4d..900d0a6 100644 --- a/CHANGELOG.md +++ b/CHANGELOG.md @@ -1,6 +1,12 @@ # Changelog All notable changes to this project will be documented in this file. +## v0.6.9 - 2025-04-04 - fixed bug in trends/moving_block_bootstrap + +### technical: +* when calculating a trend on a timeseries which is not complete, the method +`moving_block_bootstrap` may throw a `ValueError` based on randomly chosen blocks + ## v0.6.8 - 2025-03-07 - bugfix in time harmonizing ## v0.6.7 - 2025-02-21 - fixed output format diff --git a/dist/toarstats-0.6.9-py3-none-any.whl b/dist/toarstats-0.6.9-py3-none-any.whl new file mode 100644 index 0000000000000000000000000000000000000000..3364d3cf699f596787d59f996ea5449533a59f82 GIT binary patch literal 32974 zcmWIWW@Zs#U|`^2NXa}L!Jn$EHkpBe;THn~gD3+7LrH#OQE^FPNwI!>d}dx|NqoFs zL1pYYU#>$AJgg7?9=M~J5^LM=KzEn0<f9qenm-4tY3|)8wtP`u{629xrEkK`VyhM{ zVKV;7-dnTf*sMCP`3yBDy%f1$oaR?OwyDjgX&r0IPE%F&Y@ux4xCadL7R4)R9-cRq zM@p!<fLG^Q#?zYrOLc!u@7UohmE{?q?L5o->fhSyr=&~zl1li_6&x%(#`l5uL)KYg zRA030SP*`bfq~%&BLjmhiZ60gONug+QG8R{@9THOfT#8Q#3FX>*9+Ml9*J#s)+h>2 z&zfiv`sh=6m4_*h>y)>X?|-+oetyo_QEjb^)Rxj11DnUKY}sop&+j>}AQG68IU(>~ z;v<QblQc|RpDvr*V01UbYO#~U#<V}LlDM*+mU5orz8V{CDQvRJ=j{creD9UF3RyNT z35ko;6A0%#<@#>z&V_&TPRQ!Ko?*ez9n>pPwKH$R{`*mn=e&Lw;mxKSVZD;~{`5J& zSv2HbPsAU-{H^TwdgDJ8>)$S(c{j87`;1B7KQQ#advR*+fxkuH6E^Qu&$eSkgo#;L zMx<ZwF}+L{1_ovU1_n8-VUnDmS6q^qS5mB3P&qX`I)AZ+NL~Mw{|O&`!W1(vExTmn z;+d`%s<(JywQT4K*69mOB^pnrhX20*U7{mbZ{bX@K+F02KmWY(Tt=GbaD;(Z=V@j! z#czEvVoPP*ZT>F3|N48}H^(lgXQ@*<1pAsbTx<L<SX*pN7WjGU;l}&3*I(a#fBv*S z+xs=q>$k<vx&Q2S{QMjfhZ#X4n@xD{IlOGhbPKxRQM6X1-B2{h;^XECE_?28j(?`F z{X2E*hPkGH*cqjk#2is?`M}67!mT=E=3$2i%Q)QSxH$bE?7SbnCR{oo!lM4|Hvb22 zS${GW7VXnfp1*$4(@P30rs97OD=|rZyL{(f?%bzYiqH59Ww_-z6jd!+&vq4XbU)zv zrvB)6$lAJ(wcpR~%YB)5=ut@Y_wtAz;=8S6zx{3bVc{=ef851qrohfjncsisB}yp8 zUjO!A_}9hvia(5PBHrFK@Ye2++1x!fQLB;pBezfBCIgjD6Ms`__rPU8cSoLz<k+~k z>%*(a(8R=J@7uhdJg8QzQ!upCRGsD;6D8;OuWXecYsm-GhUDW%4bKZLYD|~<QGKsq z-vQA^V}7TgqfI)EyhV35dz)TWy0z==l0`C=vUa8oVM$@e(o0lY8YgcJf4OJ<5veTW zs~s%vvmUbUl3elJ<BexB|Dxve?UNI)&b#*J_pHAGo6|#YOI>>tc(qM!l9~6>HV<Pt zzbl)U-gz{!?SA(8TUS2%95+3C)5mu5k=AEC8lG#-3YJtgxB74|FnFCQ{GuYrSX$*> z<r9Uoi%n;R@uZh5-&EPcI#o8z_PgD}dhSyPs$AT}DxcVPbd(p$JX_#xILG&W<000c z9g_tLg$`L6suvtNs97y3Z+d*Mi=5J9_ls+K;y&-o?#sR>EYT}@)bnblR$^~8{{v}n zj%A^rnh%+N;|bt9p8Mcg%?=i&n+L81$lC~MrbdS9=Za-7;qUDHwfFz0V87kwDQrDY zKc_zy{JX(TZIh2l)n1i+<NF*nIg6}?eO3Pau@V1p#pSzN)dT|$?g};~=O2+p3BOAd z#MKRXxK)Lx?o*qkSw1&Z<G^~JlP<}z@A!}EH1AehX_%VucGGd&gEaw9j-4zo-thO< zwsTy9-?mLwTXk~c)sTz2hu@V2{CZmw_&WT{>&)!ZfM0LFcneylhE99(?d;@~6_;*z zRPB)u2wh^HdCPKZoS^P3*1!!ewRbk{wSWBaAx|#rqW~khX?aOYC(Ptu{bYG<*UTEd zy(!5xHsQ;~W5kxMes?<V?UU6h+c#OS+IO!l_s?tJ*_)$RZ;taUU3O3QZp>AgyEdY! zi%V`D(zzUHe|UEBeV35Z==a-RbH32&SY3N@)%_^D*v9Qw&Wrvs-s8^|!!-Y7^gWYT zw+<gYtJGK-6B8@-c+T}@LT8IEu4KDlbJ$Dy^XHINr#p}QT&vP3yC?Z+`21ZjZ!EkK z*s0SQt@CQ@o(s|YUtihYeM5TjrE}A_Z~T<)x0gjEKK{eK4LLPUj5{{YmQP~%zw+Sp z?l5-$wHZGZ!&pzBH)MENEB&FS@jo-7kOEadhV21|E;2AMMB^;DQ&Q6sOLM?wRBr#l ztiuL8ZSUFkH(#jK-aX~u0`|#EE^OGc)YCti<)KCQ!L8nF)~w$@<Dl5Cxh$8yd(Rgz z<US{P)~tfVwq)7zlZG{uX8Xi=D!Llq;Fx~F_gn+typz8;Y$RShPuUUqf8wo)CyQ$) zD(P=~_2<@>(uJ=M?`u0gf9^RqKl4QkE!gWWzA3Tz$NtJNMm%Ec?d?0WLT0Z_F18cx zK613`MF*!zY^BuM?fXSKX2k6%zjZ+R*;6rj2JO(LLDRIRR(8Grp47RIx6MLu@zu1> zKND{qxcsV8_oz4`EM${gqShT}<UTLJz>plnz@Uh=z|PDoC@qOk&PYwp29?=wWwP^c zpZK_&@jv^tzfP>1CW_=<G=9b8I7z$dS)s-5@C9v(`Aeprm2+LCl6su;(|zmB`seIC zZ|=*s582VZ^y#}ZZ{F-FE1On1uYUjP>9sen#9lMM`Ke96{@gF=d2g-mrKRq@`}OX} zhX)@X{ONl$b#~g;)fL-vUd@|o!DO{=_NuIRr+!&QXM3o(v!`y4Fv~qX{qLgv_dmVy zXw<u&bYa@ooWs=)g@;unu38_y*UrA^^P=S0xvdlWITDv$JRor6N80WC$1gH`TyrjE zxoW@})}@h`(_-&4Z8_Y$FLDCAyT!@qd*AqE9tdt`SR>kSPT=n5eQ$mh-MpE?f41dA z>Fpe`=&<z?-=59?Jxy2j)TX$P=N&$4?M}VgS1wt`dPn=3@ee)*mXCXMcRK1G4_LFv z)jq7@Me)rg)*Ls!Y|Pm3Gv-?Nw5!)UMOnk9Z_U`PmtQ%z^PO1vtg66+^Xs-hvg>bH z`D>y1ik0hy7k>#Zzp+Wn{6p@}AD6B?T)4}<&Gt*d4BpQcNm88QB4+2Tt|)Qu%ld4k ze>*L^e5!8Tao!KoA9bhY+<qT5yQV^8x`+Ljw(WoJ@80xub!(;5#!Ys~p;KS_NxnQ1 zwY8=>r%{um(UIekqGVmggy)t~yoy_zcQvHsY~$4k<&_Qpn-^X9W~p@f?6<XSvjXyc zXJ6v|k-yuB;dk5C9VK@S7i2pw?Q@J=xbPis)1wax&F4>-zqkAM`0VSWtL5$g6>R(Q z^6BpV`z>eBzP5S#`}&IS-wwy0Z4Q5bZ(n`ozi(F`&6l^|Gc!FsKW<*_pA8j%zg_(q zbI97d{doF?$2lMQ{>jOmnXuw;vuJ#J-$`vA>t`nxKjRByyuJJ3^8*q;pL}8Nxc|mw z5fk6OS$RgSy~h&sou~Cjty#omY{GuArTz%lnN6F{zw?T43oG^g*6Hn#RO+Z9QWf`+ z*PZSDhG;qWdL;q=7ydQBH5sD!>VBKUuv1nd<b&~{;)bcol9zh6g?(|jd}pdQpOW^$ z7W=#F-rQWyWgC(mEa>u1@0f+G=94GJ0k0>T2sS;e-OT5=<mrOHUp>q2z1yL$*Hl%a zk|HUp%>ClbtKT!7*|qk^#m1G|{FYN`mfkoeUL~kR>!v10q^3xLq<D2Ji@~Ic9$^)0 z^p`~$+B|7krr=N^dflLW*_Y<*Ll>v<?wNgYrO=Y*OSc=n7VdA@`Q?j+%$^;yoE7E> zs7rNKeR;KcQ}b4v_Mdw*Z?lE?ZIL+pwSAHAjVZET<{hO|R7)=}=3DWj<Iyk9!)>+` z6jdaJkNjLeEvI2YvatJ!c&)AjDPPT*)PL#U`&yE7x+*-b(ZR7Tc`4`q=NvaD8O^zM z;D(ZY|G5i40;aq;%eUr;dcZClySp*I%Hc+7FES2Ho)=_jajfrl(gpd|tCMn%EUla@ z8K!3??Q6^Ceb&mbHD_*mymH6>r`z5I%+j5%>Ga9*<MF-sBUJt~gqxh~>3b+T?NVac zoxHA6=?3O0N;&L-D(-h!y5_L@dGm1kt-0$v<8Jl4@^6nex5!>dO@2C&_0YNK=v4yG z^>|)*$1JtDz2yBSo@J?9&i1X5*yqMup78tG%sj56M)wZxy&v*fR^kr>bKjr22i`1c zy>8NO@7`RuW8dry&VY5_L^D0Uca`r5zM#jlrmQ%2di<rr`x6sxU%k4DB{T7AQsw?N zKdxRq;K7?Wg|Eab$yq`q`q1u!;uEhr?M%N~6z==oHKc(_ae~O}&O9w8vwe#nIaudU zUdL5Cldp}vuuFpN>C4Fsiq~#E5D=SivgSj{F?$x#k{Q-tv*e|$$~PW6cAPI@$@IND zrtn{LvYW{hq~k5=s<Uf;m;JWpuqD|AX6J-wYtHb{el5%2SMqJvxd7>RhPU%7_X<4` zvbw-v_1LCu?dSCtuV!E4a=VgqfJyPIuD14_Jxi3B4xWk0QP^korC#~<tms`cf7b|I ziY)B@f5GPTze#&N*8R)ov=DA^{O@4!`js5(jAWaW23{gwZ;z>q%{v%6fhAV)zs@ZG z|4Mp~&N3I5?fZS@z{=+w>Az!&Lxj)xJWk6!{VsNQ^W7<L_&NGRr(MvHEnKndyv^-( zedk%~R<R0A=NFvfky5-R=mS@H?3Z&-`A(Oxrryqd#n-+-S3|8~f4~I2;EG+#ov$1< z5PWf$_x7i5@#zW+&#N1(STpPC(UUgSa!om!YgrfGlj|<k)zi!UX{v3l&U~lCy+rHC z{(}ncHifLgCP6bpE*TltOq`UN!_8sx{&eq2pJ0iPu}c}1_8g2!IS}?h`;=<K>mbE| zu=lG@<*4SSa#$xedUyQJc6}XvYQw=d0`1&($LGy>TOP2IB|Ej`%9dB^KFNu`$r<Zs zTDWye?<$@-Z@Go}o47L^y>p*k)>ZHEihL(?Q8?|Xad=fqwc>x4r|xz|tRCBfIE%`^ zDJ!M$NL}9@+cx|1Vo{&H?+$n;JE<A1Exo!Su4vnGwZfQZD<n-9wNJRP@Q7sQxnq-> z!dq1Q%jPYf)I8T=kMWV^JE9JkY`#AyO+V+(jx2`N4UQWWZ`4K!1{{CX&Y<yw;bLTp zerLcXhKBk-Qf57ymni*FnYk>O<5t2C?ldtqNwrz7%Cl5nuzbDHRMGQt|HTObYnKIj zM{l`!R_xIX*2j5k!;HAXpP6jZ-*Ubw=3=;lZNJ<(m*73;1a8VnrfYmD&DsAqt*E|S zWXAfNXJ^W^@RUswTqJ(3XTh8YjzWvoCHOVFJ!gG4QAn14>S^a2F`;CQ@7dETYArc4 z4I?vS=GjWmm+5R-s;+fC``)9E%q`!p|MfN4somm~x8m+mF$Udkxy|dnK5n`4V5ihZ zfwBqNU-nI~n_t_{<i&h%=H_-TwrGL57tb~|S*033dvixt)c9mdwXwk_?$e^mp9=yk zQ`yq~wHg|Hp13yti}HuLLCogCEtL~j-gtg7J?xI-yjLr{-qc&)ShuKLrmAz<l;{o5 zC++{{`gU4{aq*j*ivO4XnPfU6^=AI&+|%pK->|swmAndz&g?$-&pSi(150m7<bh`m zp$>J5d@XAo?njGm?3hxlvU6RxX-1leHdjK%&ZhNmHvh5GeEm+m)aKol>jFiKx+PDY z^{X_ll}h2i@MIO^o8pWES5o%b&)vBBT;roR=a;ic{4rR<I{D0!drJ$XAKgp~;t)P# zBCHU%>894Ao9Y*R*!c9qXWv=zO*+`dbVuLAUq^mh&+|GZSNlC?#k2DB-#u0e`87@z z-|RT`<>ifUDm%^{l=@xQ!M1Xe$l`iCNrQ=&XOv`jsCfqLQPTOy9KF71)rY7{Oe?Bg zRRhX7x*h~f+`67;me=}tXO7sE$3Z_H`A>TOq1@l{$__cLj<##Zy&gGmoJ;EqzUeGs zYt&;RbnVE~i&`hOo7Q^n+`{!OX;Wq+_sQkK-?lIBiCWw{=TVvYzUc~Im;FAOArLuJ zAh41xO>|3atRstFpJ_z3{laxhCw0PKcsj^F**~S&nuA|FV^!`f*Izdj+f?V6t>4nH zDwa=<VanGQ<DT_mjL*I}EtvkW;p#&ri<zfwe>|S^s&!_uR9?TFqi^#2ja?h>_MZFC zHoy6EKx5>ou(p?zrK~cNT3i<RoOSV0aaopot0c8x;r<m}N{@fMz5k&;$6j<!z7FG4 zah2WACR}jtO;|B$X3r*0f7Thv3`cw>KRn36GM!)8r!wfp7xh`vzoqwo?=cKh`E%y| z*KG%%&oeq@cqi}eA?42niGTd`^iH-P6zzx+7E4KArBuyQu2E6GxnXi}mP*p>rb%sC zzt3FjpQ!fuOU(La^H`^(o3qrtWc&0Sq8R`59+uNs%xxgJT-?(C;H<bNleDE9s*gFj z$lqx=u_7)qS(5v5YJ+^uY2`=xtMAKhel_2Rg>~7|!%b7>+_+ma@9B}QOY^6Bo@v); zzVXX*r~m12Mdv#l3-0HL-{9Laai7V`&gGGY9h(Yn$I8cSlfJfbU;f*7KbuNBbiNy{ zOEl^W{r`wD@>&C<>dde#X?NqbE@4V9!<_al_|APKYsaI{ZyqT!HvR6jJau~JUj8=~ z@fksBGqxRKTXOjRk4*81HOY63wwfv#YdsM^tnw$R$-*eVrl4tBxJdLt{)v@hw^mxo zPPe`9z|(c5%f0ze-R5Nr*(<&?c5o$0a)hdOb@X(rye#pLmeo&Zys=>V<%HQhJ`5HN zoQ&JI+B6+we=fnEP?z;tHf)QA$)r!)iWVR4U8hvPehI72G{NnUZW@2Sdz$M+c2Vmu z>9q$-^DkFQ*BjsC68!635nTT3-HDg$o@<`YQ@=2&Z0UZN*mX<PLuS7Ic`0M#n!R?X zTD^43Ph6;d_oF~R@ML@7?8L7VwT{1S<#93aPj9f<_U~G>eo&6I#tPqSjy+<xcd>4L z*{z;@blW^BzRe}txmR`G8bzDh1thO-IktKeZ&Ju!|Etpr`VKp&M5u?}oYnYb+VbdR z>jo3nZHrzqs%}e)S#l!GrrM+Rt;5dp`)LYyw{Cm1tWfb^Z)DVSqfOcAUn8<o_x?UD zoV&(2Z+`7v;piJ8s{8lVx@OHgwChx4<hIalmMu$kR>!?}zY)s&Xkysa*u`6XWwrHh zpI@eQGC9mDD>$}PI_*V7h^6!Cd8an+Go7mR=+-3xwTqua{%}SG@}5kXmE`z!!3Ejs z@;z;<&uW;x$#*FXxzX}SILYCnji0A&yGOb1%7u5me^{+g)40=?amwu1EjInq*Unpl zc65HX_c7UDDDYV9^vQ3mo3bCfnRnXnYueU*;iT%+NG=J_`3~h1l9iZ4duAPq;9|e5 z|F6pSP#W8Ux0P9Q9_0P_VUc<NSl_hfFON=s`mys9i__Y#Ug4Ye7lc?i%vgP<QGDr) z)enn{tTXb06pLk(Sxah8b%!R2|39Ah_p)N!nvb`arOcRj&~CY7(W?IYVaW~2H-yeC zQM(csTPf>#eO0n`%4{81)yP{`yLW3&H|JwZuQaLctCQa3pSdZSVe0=CR!>y2-Gx2E zit^>o-I826%lr7^g7aPJ^%Z4`UY9)<am<K(al0x0vuJjG+=gZU{}!JSIUT?K`lsKu zp>g-Wd|8>-D0g3{?85E|%^ElL?;G@7i0^l7+$I_wHPvr<M78<t@|fVug-=S^i>BSp z+x$@X+{4Yke!ZQz<fq0Ft_@rh?|jPsS<G)z#CZ2+^1PVWJ)bxaeU>nuJa=ADjj+-F zAJ-VNw=Ap)x)o7XAhIZit1YgGWzPM*Kkmx=o_4ADB|LHA%<0$bS7`XkWour)7tq3G zwJE7e@RN_qBf0D|Z*{)>``UIfyxU{D?Kb(3hh=w18CWe}HqS(T!e#a;>umonnX_o~ zjV*7yX1rkeKIL2Eu7YzlM(ul&&lyY)+Lk!?w)(lnt~yN`f9sUre!r8w_}l9HQ_dIP zn{)WZYVI===f(0kPg<>}sL;roeEfu~&iYku7nA4DarHmS$G@jT_Y2={_dPQ@`)<9T ztB~@fuKc)M#1rGM2U>b&PWiKVifs*lcJbrsd!(PQyYZiU{y*(jf2NxZ5&K_!SeI)s z<9hzX)sL^||CJT53_aH}=g1G~yhmSpLq&P@mt~r-=KcNW(Xm?1vn}G6zMM4`zv(ye z+1_2Z8Y=g%d?z@`$mMdmd&buzTW4&G{N!hFpt@h($;evaxom1@`F+7|pZnX5K0Q*^ zmz#cAVMFNUGp)*Jc$t~=`6Tz7tY;U#rysdi|J>~DEeg$-Gud3HEDYNnocQ)d%ko(L zu)YB4BCZE7v%Kf}Z@ChpyP$NExb2$nm#-Y0ysTsSb~P*2<*#y&9}l+pG;NNq=I@ru zlF(;UbT2A#Ok(<R@IrDrThmKb$;4_;HBp0$&3Er4FZ~<x`S22fNoP*9y>wO2USaNS z%dhwR&I7h@f7lVDA$Ky;qg=PGSh$sgfgx8O*YI^_UP)?ET4HjlUP0y7$lH03Ere?M zC)6vj@R%{)ar>TBE%iZx-!Xo}UF*4$K1RQVSj3bsnhAPdnG_;(KIwn$>Mn85Z5iou zZW&8-y>!CnT9qnJO}%;7O)6~lrqwLkA;~*)cs5HPy!CG5B-Y8jqE$z^4~BVNW?sF+ zY+Cf;s5qOCTW8FlKX<<B(={QXE7#^+%eiXWI(1d#ocW@Aw03nxZ*=jz6cH)$t?g*y ztd&RWtrFL0ME<|@Gjdw$j(KaMm+G##Fy)THkL%Mc_FP)^&UDW8Ls>IS^^bOyi>`RH zy!O|hr@y|>e?LvW{`b@Cw<Gub{r5HP=gVL3s(ya_b$$2!`rmV=sL$TE`E-(U*K&=? z@20);JZ;)u-fHG4ShD@inl+1Bx<4DUg)d}PT71+i<Dy8z#0x<)B%`^x_wNkky_2(J zU8wSZ0dMVW!{qJd*UPt`p1=Rzy}9#al|-JWCz#H5y0t#J;%?Z38LOqQHZ=U(DX+2f z=|6Riv+UdBILpPSgyqhNIGxs+ug`t&mS1<=JcroB+2{AXI^w~xa=|Rs7uF$M?jf7a zgF-g7#UK5|FCoWuFx4dgTAsjxh#iFr=Z;VQqonuv_L1iYBGmkkxXg2BfBsWJFW^%6 zX`Ol@hfAv#O}GA_t~~ERL~xTF!=!#OU7in9R*LSGHI+H=IEnQ+s~YQ6@guulTw~qL zRC!eU*!67PYc02X7&SEP&wJ~wj-6M1gKwP;d)-amy-eokp6@9tKcDoLt?#UPY%KR% zwLLuf6_-o=gdR+P?wY=4)#pFUCr`Q&#-y`WXCB8|`TG;P+IBvlH@zX(N#NYgj;L8> z%nOfh`y$S6yYj;Msag#RsZrB2IJKu1#4(?c6}R|p5T>bkRC88bGK2m^hpb2af%9@j zazFIEE)DNKrpBV^&0#p>>INZUp(D$TXKvT!xG1zZajmAaeCV|e?s*Df465r78!$Ld zpP*cDwM8Xl_RU!BOVil(EcV7~shnct<DPClqe@Jqn$we^%jLP#zC^wHzGjuSOH6ap z4*c-?((K_a9@{e~>Y&CEpQM!`2~L-ee~CUD;Hw<LyKLj}DTgNPIKnk2n5lN+#}J(; z4m)vWuE4mq3CBN9wB`8I^sntzx2SoRpXbHKhnw;adESu-RZ#U3^NyIo*mplKW#h-R zubm}Me*7n<EPvW%dTP_w8nFrea#QaHn0>m+n8~Z)ByeHZi8U$(DohN`GAB<9Dm3<X z&zrufQDANCg3s<(ZzjI@&#NR~cVnH%$`&@$3}5{Tb&dacWDo1jJHof($)9K5H!fbe zaeZsX`;^qvpHGKu;7L>w+;h`cbnCpcoQ}N$OJ^R*{5Czph0m9xO2mE6E8~kv4mtBb zNk3S%P(<d{6;8fGPle1}=c)v`6tYFj2rjx{nkiDu=_z#jW$$L+r1%_#;B(7D3Nv-5 zvT%huhcoQsu5EoEoy^6?TRZ&%*9)I{!bLT+W-!_xm>9a1TeIVL)%GCSLMM&2T3e?H zoWDHj(4oqIFK>&7@7W(;Uv9szc6Z;reSe=mldxK-?pV^ZVTKc%>s+Bb3miStzB*;{ zma#9d*7<!XtS{@zmL6Lt_gnk@eimlB-F9cKjB;CXN#lU}UKdur3jfob-=qQ=xvH;y zJ8Bx+aM{&5(e~8c;NBfGpI@*RtU34hfZ4<BGA5=sv2Cm5AKdyro$(qAyV)s+Dh|_Y z;nq1?Q~w^nePG$!-TVLV2v<7k+MIR&vFVzRe+?q7mL8tLtbWO-BVp4=JLMajBu@)$ ze3R&Z^lz_l*6vGlz9~#_i+TLGtKflCJI@BKn;|>Ar`_ORTc+VOZJmF{LhYbbEz3R) z1A|9z{~WGf<gk`u{k<O{i5%g&7F$)y4(tB@@Q=NviaBz{t=4~!LaLvv6DyK$6`Haq zLQ~<7BS)WRY*KUuSO2Wn>NDP#N8hoC@pELn8fL*<b7A$mje6pXD*JY-rzHH+^!{j* z#P>#J#gt}a$1OfB^?Pir8&+LfB)rn6I#uGY|Ln8EK{9WWpSHhpTx+qi{XLWS9j`#K zOY4p?xQMcJq`tprmErFpWz1N6d-qpU7G6iDzPo~x_T2h3_aH-~PV~v^M>;MQ=^CDK zJLOQ%)vMGHHFaf;!Hy%`&rAd_&+YenE^=wh%UA9RFD9OC+3mF8UC;aAOSeBxF5@k` zWBc;kcK&rg<5~>9Y}`8a%ADh2G3#u;8@4Wab$7?3gJrjmoLw4!{-n2#*&WV3X?uga zd0+iaT6MQPu4(F@lqd03$tPc=9GQN$S@*}kNvkAvCCdcrpUzq=WFJ$zCAm{-P6B^y zvTE?Iid>WFY)wDnlfFl0FBR`$yA$;C=)DW?Sjz$vC+?7aX};(|e8<~lF6Z(U*D5rs zR37d*c*yWoW<bGp8`G@YJ-IuY%@Xc&T-tm+^YC543t`%8@)pj`DD!0fQ|V`Emg5t% z|KEN?pMN}EjK`<ciIg8-{=49EY2@eW<+Bq{b^Lu@crIeASOk}o&G9{=*8ApZmwefo zaDAO><icx*nR#Ai|JE|LOTT&fP<hM84cGtPu<uU~?`;0n9TyW&T%q=^@9^)*A0qdM zE;;t;!<Fg@)05q1A6;eS;k7vU-4Vv^lS^NoZ@hcT`^AG_!L<fCA8Hd4&hA_I`QDzx z?!kZW9E`SQYiF|2a8DDCk?>flF7ot^T2$}!6$?yn{mv>juy?tZy`|&*rA5jWTO$;` zm+e0kBje$u@@$jl(bFPwg0Fu>Dy?%mpTPPjrR=EQrhrSEJqsS?y?oV`TF5+UX`j8| z`XxN;57{M8&ueDfC9|vM{ng;m<Wo(xl9%Vr+jnB#_RrtfN0@7?J`su2?T!uX<avEk zaNho`hgZHIJU!vW9XF=Qr~fufA6(|UE5VJI<xA~@I&Hp%Z%aAjSFs5_e{{5K*&d12 zH|t;TK7Ny5ygkh4Qc2sLdA%#mFWSj*{y2H1SYKXkiO=TjKlf7Wts9^I5Pp7?(L(*q z$$vBdR@zy#-+!|_CVWLi^f~)R&-j8zN9z+G{^s$Q+uY8&So(zEP3E6*yK8%{f2mh^ z_8-yUVGQ$%ax(9Jq%6$9P}j=9poDe$C%-B`FEt*zqCu~q@@?$j{Kq>k{<C%X$9||i z`>Ad8Mp^ZZ&AVQFE4!_CVQ+2Y-Q4@)jtSagx!mgBtL&0?yZ--Pp620oCH#@^o0x*- zM;{7!KKYP6FTL%_guJa$^Z&iAp5w|m-<mu9+@c#@b2ZxbKbvOlA5-zkR_vGOn+h8v ze}6y!eQV})=*mprymLRFllmjU%`ww;h4TAWDvPaOGb2(lUGBNr_Q*HOFIm()Idj~} z-~XBo-`vV0e^`x^(s{2wSa^7`+2*x37Hz8glw5aq-<paen^#3vMP5ARA2a)IVH#`W z>yAq%i7KzZ9Ql-^HFMvc=)6N;*LXjlSAC}IMvmCeb)rF=Zub0C+*4Sq6S(K3Rl?ib zoTZK8(R^EHZf<?Q>))}D(<1+_TR8RFw*9kC8SQSJw?X&#T|<vYYbKjEbG3f($n^Jr zv(@xb<5!P)f#Ss*4L;iFRG!t<Kf9~8=-v~_rGCmX*M3-jty|aeM3QOE@6_Dyk>72$ zO#6LRIp(&>7gfb7t1GoN<@aKn)z&J79(w8dWc}Z|^;3R+Ufw_d!PbU8kJ7rpnOt$N zSN__&%WkE}&wwP~PiZmh9@!tCX>hu~{i#U4eXV3#W={5VtJ9%2YyH;PO)qM-zy2iS z;r1}azsu}*=!A>LztUXRoLG5aOURk_j|-<aD;}K6>R`OD_i;vIsK5W7->G+YYg~V| zy#CKu?tS4$19zRCRWiM`B;54nI<A*)rZ3%EOGFR5y*&1_<ZoQ5Rn2^F^A(=vCzd-~ zFI;{&z2^KB&5pmS?Pngo?37pJI?eHA56g=P$M-=GmbGYhKDazNY<26Lw`-W~duC2w z(=02rSA5BPrVYFw*#2?n`0V2NYkK7{*Tn}e%ht^F^bg+9oOdL}@9~$GgZ4pAf9@<y zTJHGbe5J_G#CFHzWp|!_PjOp5ZF@;&`(K&$%d{_K$T{{Yv}Co_bfv9Qe$raQS=Q(L zZG!tXtAd{;j~(PZn05-zs}hVT_-6CgTk#{8(fXR}*Dfz=Oa66vve~vTU!M4{(pl^p zbG7NaS-ipBG|7o8PAxio?W10Umc#UjkoXN7!&t)j;w<G%6Roy4S#VbVjtXJk_5Zjn zuiVYBGoQr@j~qDW<a^-t2FYHJ_0|^gqC8tT&%d|Z_xJw4`v)SN*-AVn&A9F4_voh3 z>V&1a?KYARXB|1y+I9R;M!VR|tQ;mwZXYwoD;_5<YI7#P=UgzknR&5zp}yHo>wO!9 z%=Ne%qhe!g^jRePPWoPxpXgFhJxL<8)F5Z}8)2=#Y&Du|C$DhWR>+`SeuFLAlqaP0 zt&8BR>V41GIZj<;T4i@V!1E=8r#Nfy;+%^;F;x<`v;WBR?=%nZ3)W-hwK)-^=a|^y z-;w_*b4ggt@_Nnx`HXcUE7m?c<l6dcl346%`x)tLyPt|*;hw#mjqgUnJMY~(w(2v@ zpG>&6o2|8IVaOs@_mq{dBszsxtS^vU`b*1a(F&m(A~^z_tP_@8i1<01Q!U#w<=Wq7 zv!uytsv1)lXr1(J47HOJUM?!fXMFGf`P*fys@}`_>R-7z>-^sf%j~}2ynddiWroqt zkUYaFf)BNy)-7~ND1Lw8W?Ldh@)SmocWY1CS)FS1bWSx5=wn^N&lFIivBqMOdzhGD z>y^9<lPn!<M56k4b~?AbEuXyTdA;g^+)e!4ZENOoZhEbeBxhx4B`CKc$YrxSzh(0m z(aV`Gi$8lz`w=wd=DQ~~A&;E(iyBxKdQLL9^zqk&&Y2M}bYxoudYT_ZKRVg>#A3q3 z>ZUH4*_PcsR=<@M-6zfX(vi~s{Mdv~vwG5eBZD6v{G3{S<qWTez^$*M9tCEn&#O+j zC!X_M-9__Er}U|6%~Pj$oA@iuWbf4P3p=MYuVnp!>O*lewm67JZ(Xuocf-AJM>Eeq zmi|`ridl04=P%O_%Bj~`Wfz@laF@E3raFgf%8mV}f)}}joLSekXT~MLV!;D_?030D zPfyg`5z!H{^67r5Um`*(?Q)Y;KPD$W(3PFOlK<q|pIXOP2rkh!c(v1SnZdf&i_?{V zb1a^AGUSQ3Wa;gISNR=5oW7H#e`FWbtu}tMtuJY^%t1B%_n~#QYMb}{|7Fv@DL%>X z<x8C>44<oy`xhiurq|7$uK6US@TvKwB?2Kwr`~s&9OScg^-rByf2wSz$=JRrvYV<Z z>Zqn^d+K7A#i#I-6T3L~U)wz^NSKRndCl+B=eFH<o9Q_JtkoMn1xFEu@VDJB-W$v^ zYWox4Da<ubxG5;0Q!?)2!;tRS-h2O}zU(X7wbW#<^4IiFsy?r-&)s(Qz>R-FCp3Id zl$7M~GxmRGliP1wab@>8%^0~!Uo0j~(vnf#Kcn_%6zhT|8})c=<~;~~n;Ezzg1eRJ z?bSrX4dri+9WHyautDHfuguxBESW$KAN}OKn_0&`6^j?Bl-IYu{PlYN{grRumo>M< za|?ZoZhoo3x@xnqR)@w`!B<BUH*WOz5ZZTs+wEZAziX!Xx!hsV3${?~5%hb!ROQlw z36G_Cb=f+uHg!DSX0vRvY^Pz>#PbD59j48l*;OKL6Kz>|YPFoxIc9J9mr=+2*7&cR z`Lbz8RBwbuv_{&U*Fi<otG!$FLO4UhBusRFDeM>QFYHcYRbQ6iC}=$;{^!QItPf4t z9xZUqmQQ9-%Kd9=t~Ft1vf{5*7jxLS7+PO{>|-lE#q`TNSG!~VI{C+Z3%Oc)n65AM zWr$zNZM?wUwy$Bsp;+lCDHfTy=<XS(muhbeYBdyVvb^-s_ps=qwu4y`XBIi!tbcI& zd*Tb{^%t}+J-*o^byRJQ>!}y-qIUc%srb~q>-^s%-3wJDsutf6`Nno-ukP_G*)NKb zdO{k?dwOc4On$B1`@{W7{{1g`sgeqer#NzilV>km_*8di97oV;7VV<TO-#|-mp#5& zlid`-95S=I?5g_*=?9OQ_FVLQJvV3jHMPJc^YougxqRT(iI)2td*p+gIBKt}z7ONu z;3#tbU0IHW$dsmpuG*{T-d}lD!n5MA=<y5uPVE}vf}XA0&TdI#j?`Vzs@ct?Q@w5f z%%18+`mSpdO#%<!Fq<gVv(3}gaM$6)-5PllOxNECeOuBRoOR!Y_gSeKTa=!r!Edpq zZ@xEp)h05Zx;w{Y>6uwJK~=nLZm*lwPfXdg_12Ai#xL)s4!r)!<mXj!<g=-;s0RDf zl(fyxq3>4o>AF6*NjY@Yp6A)TDb62nE$!QJ_MXDomXxFw)2{z4a&igFIQKUpeM|TJ zP51bo&NtU>+<z@s@$G`RGYh<U7M_XsS^Al~E8#f9UfF%O3$;xWcd^|~W-t#q;jE?7 z*gfscXZec1)ir;A-~PXiOY(rohg=>B)4W4=1xFujn_GOqY01^5-O7)qOkJ-$kMH%1 zCrNXS<=NDj99GwKYpKYzbXhVb=^pAds(#&4KFwm1gZ0-<shK$+Y&v)UOKP07V1vSr z$;M9{`&CbG{#scm&2>0oslQUJqGVVPW8b1vTk7>|7IYjxVWP#d;Sm3&74_?5Kh3=% zv-#b^o7d)^KDzJPexcLdb1T$YR$elfY2LW-QMt9&CGTtH7jhnQF?{c2mQpc(%Rf== zuABGig(<zqKO|l3-g)8S;q#ThY&YAq3LfVy@tGelBX{}b#Fps?Ohl$7{(7}?hTZ-Q zgS!d8wkH4jvxx2T%gRN2cG{i#RGPe1y2YHse%{r(l&N;*N(T8S``Nx9TVP(0|DaB$ z`5(s}^T#DzeD!i-w{G5iFiYY5O`b^%eR{jJztwd!L>Iqs7MK6S;rNO1`nn~5zRMT< zRgaWj@H4-z<Q7wk)9;7*KW^CXez{p*;`fiXeCEZ!#e(L3t==2=d)r|%iNESBwMUj5 zj$6;HY?-*@UK{t*DGh&LJa}Y%Ep8r52zSHgo6Cz=yt;LBnc+;41D&t!I-ego66B!s zX!6JMeILT*n3FcdOqOhp+kgF?Lvh%&J!ucS_m@7~_gcSP)=Wa;ujB7oy~iKt?Z5th zyYDygc-@wh1}FKB*{|X&wz*&<`c&(@<0H?5+y}cacm}H9UKIZN=da8Y3!~2OpTI71 zS;MPPk@?t6rYD(_J}%8}g;PGU9?j-Edd&NQb)3vY#YL@~8yE#&FSvLr$X(Jl;6RGp z>CJN|`K=5W4Qbnya&hOP_@nnc-j$^2>4mN>ydL~6uATkLMYXG&KRrm-TV`V0k+Wc4 zMBbBR{`G=dOd`{6CuZ$=Vw4fBtiIDTtklvm)4_IYf3tOZlWL0<-`2HHg5vgtTE|$t znBelY|JL(K3pa``_xZ}a%>QQU?c)`u-=6&V#v;D->*uGbTK=3d6YbJ`uCM#deZ5QR z?3<|%yxz9nF$vN>m+p4|&%?*k`Ua{$4@D=%Z@zmj^-ychtpy(Cc6y(z<qmJ_G@1Xl zWXox@ckF2~CGvt!x7Bvd;IhoWq3pypTW=SSwA=2w;#>nAlj>F5L^qw4`k1rj-qtN! zF60Lr&f6K2v{~DG*|%?>I#XtQU&?d&WFN4!FjoH9$>8No9h!1?B$Z@Omq{+U@TNL( zx0-f^q6y2t1!v9I9$c0F_5#1n>MhGBiY|*?m~f)W+N)~Sox^*}%1`Y(z4*3B$wOzU z^)3AErIWsNdvoOl)hnF&nd8<urSz0e-K%ro-LCJA+it&apTznpwx?tBrW~`Dknx#e zaPvdaJ2i)f++b~`-Cd4RSqh6%e5X6#J$-u4j8b;h&5CPh+tetfzSs2scdzx)!si>m z*c?fm@Xujhqn>=?%2M$y$vPS@=NlQYe)_CCNk3FMCBZ%+;ZWCV**!k{1Qxov=iKgk zq`cwdiR%W!mtXAH@Y&T7_ojJ$wOX0W?UkX&qTT1z7k{v^;|h2BZZzfo#9y43j4#gN zeKT40`+)~1EElF;&RnT7Ve>7Wy5?iPXA6q&H9xbcnR3oC?RsuB2V3^zsN&4|*Rq1N zOdB6XOj$hX@_pt^pD6iP`}=w(PE_U8<bU1hp`0gsdx71O?F$^c8SFP*UU)0=eOXcI zv<&g4P146)_NNFJ3HT*4ev_1n{_xr1>=Dk0mjN8%F<lcMnpLRr{?C}Re6h7ZOR`ax z(Ha%&{q2pp57q{B+56pCakfr<{yzr3OFQOrPFU8bzbd7pzs+m<<QDys#I4)=IlNt# zy=P0gH;H$b;HCWihBq1y$!yu(^WvV6!t~|$qMLS1C>QW;3$<?kS-t-f$MHM2o!bwF z$EZ$QwDQ`rTb?ttdUw3--r*Q=`Ae0$f~ey2euKtsmJ02MTfQ2*b1mA|kl^B0X0++x z$BTt$VjDMib?@N08Wa5K^%d8x`uo4`wUs$3@-#Pe14qtppXQ_hTdSy5COX<nXB^)U z*TCI)Z;7bZgg&qIw7$>EqQ|mU%yF@s@$>8D)?FX;7ju1mw!tK@u>Qe8Po@`p6Fb{f z`wj?cdU5eIN_{Ws64)7U?f!tZ{pCH@KeZ3ecE6FkaA-@7+2Q<`hqA3WuRnWGdO6?n zb9Vz{SfNbJ{JFVi(v#Vz?fmlo{G%Jzn;BSkpXfZd!S$I|-0B61tIpT`maopK=V{+t zyzwt%692CM8}IXq7EPOe?%k0Pne^gYOXA;gO<9++M7x;n)>`>rB7Z)LHk>-SKPKTu zG3$ecC-+bI@`3BX<H`Sz@tJJ4dsLCMliBarm8~1B9P1C4PI+}NWp2gCpBYbnzfe8H zTDEKrYudWgnz0*u7kqBLbtwDJ#_I()vsPvD-`QSOQ#)O}M!R%+fBV&Me^MW=Joi*M zHPCj;2Tw7v!;5$1F3ftV81z$n<CQ+MwKwBWnEsw$(-q;=-Tc_;(w64Tr2p>LKW@JG zzwyic$FnuRaDS;f>iL8H-0iEY>O9}(Rc$Xa*3~xL935I+)UObnxM;D~{XnDnOczd_ zKAo-3%Flf!jYmAxn5~k#p6|@`;|!mlb~pUo#8}MythwRmr_$!1o5aq9R~p32GM@?8 zXL`1p?F{b>ex_%uuk$`z9c9q}DM3z_`OIQ*#?LDtdMdx{I%u<s^|J;4Gw1obC;m5F zRXlKKJ@aP^=FNMaAOGm_-RW1PqP&1^y36KB_gopT?OggolljAU+-~Q7`>E1J{;PKL zq>j2n2Hgg)Pc6+ajdKd+{+C=H_QYi>&*D`VE}YT#G&bG8>Z{1AO2=K2`}1~5UcczO zfBItS3LBfdDf^qtLv|fZlG>N~<JtBD<{Y=LpE3U%r~8oSQ;V?$+e4nXwgNWZhbAqb zr4{EK)K;J4DElyC-D+Vz8@axMtooVT8LxQ1t-qfl)s+5!L*k0G6Gz%q?)Tl=@M^`^ zjeWDC%coqw9eB;d+RxqH#`TFvSIw+j9EICeuZrxyv&JOY#x~q={!i@-y~pIjlb>>a z*ZiCEeaR#7%*(4CzK7~6n#VoPjEb$g$g<?Iktm~G|6R4&M^i8AZb}nr=k7US5_SK7 z{P*9b->Y{oZ9n!%c5{@NP2TSavBJ-9j-Ow;I{oOuJ0?rE2Nb8SI<oSy=E@CcChX_= zW|V(>pZEKvOLcBeG@hj8Hz)A2t9jV8MfR(b&x&6ay(FJs^fN_VP{6QfZTV`qK!3B? zwmKn)tz1hlz1}@HP*Esf&DYdvPSd(c8BBkAq<<}a`exO|6IZS%eS7&ROltR+VD^vO zw>{$QNn0Jk)4$}vq3Fwne}A4W-CjRC{muUl@1XTFl^ieBgoHiew@K|!pV#UAk6ZG; zd;rR3C*?!c`<U4n82$?}FeqakVk^$iNi2#l$S=+;$;{8wE2umh;oX1ROkkh%hkAu= z?Wgo-PI}tj^G=bs$YWy6t+OjN#ClwA`)GJgUGe7D{(Z(vRJ=FM&505@ap=ZWOMbrn zKTVA9?B2cGKIml8q@}GFltNXKjI>lQ?@{=5z%x**=#ke8=0q*EH!71Z*}PXRZ*PA1 zAa|1Iyn>~c9*Z^3CNKQvviwQ-$&e}M?oHqMU{ZnajFRO?-hFfXSW~Z({Yu9@_)+jH z$#{PDzzHh7-X}e$Xq|1@6!<1c^Rjx$9_`SULb2Nq71}iZ`Wwz{ESfZRmUOJHQ1ai_ z=d5xZPCwU9Y0?zhoNOrb!J})<1KZVyj)b0UYo7D<$<7<{4>(+ph*q_p-Fc;Usl;1R zV+*O4Hs+ZkyOyg4AGp!B)4O<cNOND?%X2!imyFiW4dHyTD|1s)h|x4-&6~ykE92z& zUY#j@I=6<8UE=uT>*cQt8#oS3eP4Wz@z(C$);^0pHkn*!xt#2Gl1u3IG>564Wk*7% zKGaIusT^t~`zN|bVK+<NfrqcfmtPCuSi7y_)sYIFZ7p+|1+|`??+rUwI&-ydCF{vk z#bqk~6;6}+HnSQ@o$XQ+`(2>&{o3~GCpzV&zcToq^sKSSJi;09qd>7=+LEPK^XWv7 z88t_~i1(`1Wv`sXS-48Kg{eSq&z>LYLAe117R7fy8OH}+|FelN^5+{zlfKDk_As!_ zd*2rm_9FS1?7}ZAE`I8rCceg|{EdwI!M33CrfZ8Y+7uXmIGV+~DCOC1;}1sbQa;?^ zt!10Fpxy4!vEz4|kKYpB_k(->y6M-QHlExnc$TAddd7ugW$)PE&8)rJ{eFe?nb%X= zH3jv`-<s<xRV=(3_rmDev#%cZULVUfzid43z96c5?Ggd4M-jhY^i+LNoVS1Jl$wf) z+JhI^J03*uPQ1IKw5r8dcwhPi@%vwy?6*`noLt&%5YcupB64l{0>K@#3p~TQxiT92 zwtri>pKovIj+I-#E!bN0x8ag#Zf4~}EuL)IPMM1I)^67aYhP5YF%su{^Txew<|Wq+ z2lu`;`oN%KruV9{d%@Axtf{|B0$0yw?p@{cnP=PLJO^i~U(x~XANDkSk~|We&*mPx zuyBgl)WUYvBCfb|DJNzohvgM4yv{YtWWSyGve<cQ&wCtacI|l5Di-j4@s%wepRax8 zP}+X%f9f*pD4We$=0^3cR^oeCAMTrzJeRQ`LN0Sz&GJK>=i(phF7ymvU3A4E{7>hy z{S1fdCvtA$zH01v!v6{PWL0%f&b{xI+xzPt#Rz6eyiQ-lT-V3C?IqKV(EVNiT>Jli z{J7sp?2A!C!v;-<J=!;{B_-xM9-4V5FxXhh$MBSNK<)blALY5iPILXul|FB<f8+T* z=Qz6gQuOy7D8C$G>b)tCEhyF3vBmIT%$CUothYW3R?m1>`0h~9R<7%f@;ZMGY(JW? zHRf``<9xAYRtX(*RV*?rzD>K4AP{M!)^@cx-#Pm4p+09m1<uS0lNY|e^YWK0vRxB< zq3>n4^Wnp+A6M0G6MvlO&F*~T-YnPk53b89?OxNOz~K1j?$HNLn<B(leQI!9yz5dJ zqxpSaJ9D-voGWMKP06|?U9S7=`r|dlyBdqvUO4_^%U&bn_=Asa^vgF?eEwIiayjX( z<L<*u(+euJ6MiwLTXwrfK9Koix$)L|yM|5TEvF0ia<1(DWibEAN-yEY5?#d(zMcnl zd~T1WnuA}TxVq(U(W=<Ei**9)zP^*X@4PLpeO9Q&RhHR8Pvkq!y~#ZB<5OdW+pRUP zC;#~Wwt8*%^to@AA5oF<x?o$neA_wOq`AB6Pk8@Zt5y7a(zi$6m-Z{1y}r8KeWM-k zyFZ4$s{8y5!&e@&{&%^tdHSB6pF|7njnfUv%!RWSwcC|vH`QGf-FEm?kb6Jdf?rB$ z(~qvJZ7Q}Af90?;V)u&F=Tesi-ZM_NG?~%mKmXNlg}6Ft{;2_tXZ|!F=AU-YC}NjL z<T2IeXpcX!XK(g|Ja}7kapeh)WJMiw<^0!ux1^40$v+ca+t#mk<=@jQek#`%RaIW# z`On=nJMtjsv~}MvoHD9r>kZqvUT;Q5`C_Fq_p0Y6C+BUfPL`ZGTfC&}+1u8K50L7R zY6I=af?G4*rg1Vbcz7}}NMWx-K-*_Q(=#!-`Hy!9{>wY?zcJwAAw8L^*RG|XUcn!G zC$`!m;k(@%#=inAX37`M3ag@nydK?tSO1<z^VQLyM>z{}%vLnT-d1XNE_ybzeTnC- zO>62mcW_Dx%)B(^^pxr&iyNo2eh_hYITCqFS?%q^1dAVU52)V!*dVCSCwx$Katj;d z&D4qQmnJP?pCzOv(JAqHrvIfK7ml5XyKW*Q66z`R{)EO;MdeAH=bp^?VX}19rtgKP z!k%nrp1dYUdTp2*f1%iqC$j`+sW4u@8ahQY_Aq0w!s7!a=XiOI!gO~%H`h@Lne;a7 z7^A3H&TiF>4O0}3NtJqKHX3|0PHom$8dExJVyq0yN8SE!slUSxH?GQkd^B>(r0y=& z4K7n!ulp=u{A8y3K_!#JmQmNnYvvBe<V1BdIp0e>XMTmKeqqzupFTV0RfzA2j3sx3 zeun&zo`1LC?Z040WhbtQ4;51Emiq?I<+=H#uRZd{nd>hV4lFKXoNFiJFn`a>JrnK! z8-Kf&9LUviC1`V|bCy|dm!)s2U-vfG;GBy}&&;$-qqRzB-@0R#dwSW8O|EBhyLS2- z{*63$>zB3SXU#nc6<2nvW~a`Lik$zx%*P=t_TF;$o=LoWZ(cuIx{GN><BN>Shd%tt z1%DnDetLKGjrElNAO9w73k!a7h1=?Kx2hP+T1Kxeq05a{9Sw;2b5CHa!Bs24Co?-T zHhk=hyU=siX8MXd4-ah8NfD4!+sN=bQO5Jh+#|=I%<XVqpszLI)e_lFAq=*+US%?i z$aqYeF2l>t@^8aB->K4S2iCuAV}0Xx^l6*hZw8@@z8|zAn6|9p^65yZP6$?0ixLgn z<mvLpP&0YsOkcm?IS&(hdS&?jWp?%H9L*Aon*Q|7;?6vlC*7fM)Kq->IZYU?oJ)^R z3gkRvaqz}E)4DyEd>NJheqd<|zSHwG{DQZX(*N}wriUH;u0$||O1%$#vT&<*jblo6 z+K&~UE7mof`zxojV5ea4E2B#L@T2O+Jns@rCKg6Wo_OL}JvBD$C%e$a4_p2THz$5C zYg#pt%}^}6>H1vlHHs1^+_XPUV+qiXy7+vT$Lz$ScIWD^2WDEk)L&R76=7)MaDU#8 zy#?PC_r*>xuxt$68qDT<b&H(e|M-Gc4<9myWSV~e`%oq4?t&9)zb2gHef6q%NnT<5 z*DouoEJHNy?T-XrC^;JN`sS8Sp_#u7?weX<`c2+_)T_AfuJZGz6DI%u_vQy%(KX4S zDuLvxohb&(=I?r@aOGz5OVvn^pFeb__1uZHNZfvnm%(JC>wB(*m!ajyPh9e<&e*H4 z<-pPBg{@9=8@s0GEZx8Ir9?<nCjTmN7SVd1$d7z8wPv2<UGZgm-QA>B-%P{%*Q{^V z*;09V;`d$oPLsawVU^yTDJs1=bNh`I^?`SqA4{{j-~MA+n04l1$r%|xf7$h{vVPL# zI$19wN_zTS`}?jp^|?x0-MMw*QJKxA9e+=Jin{P{Tdz)Dul1e2ormA<e0+7m;;Git zOiOgE-HvZ^Kc4M)d=sy19<OZP@fQWvC40n5cIduXG3Ua{Lr0dB?rYik|Ac@4-1V(< zPk-8N#3;2Pne&Qgg^Sxc8Cw&@tVvl*0+>{PPA+KQRh$~m&v+q|alvI-hUoLXq6@U2 z3M@80lF@mMwXTRa@AYNXV!no7cGvrO`i|~cc_f=><=MS=PwDW?ypi|$((Sjg6Bl{q zZa3*NlFHdGy63Cnr>jm&1SN`HSep(VVOL`M_xQKazui+SUYewxk-z_P(u$9hmtIWr z?I;beyE0*;-F5C+aV~|cPQ7|omv=1bDMxB0-ztw+214)cCwXQlB!8+FxxY&I;`_8d zOOfCNdH>Yc!S@t@d=-%K{c3Vy;o`Lc2NziJ#m(vYvUOH=7l%X5V@J<Zidm(>Ig9!t zR=0gHV0`1bS$YYB^QU#0rt=;oSoci6&HSBDL$YzAU7F0?^|LBZUv6I4`ZhysNxJF7 z0~`|9c&)WK9;O$Y^|VMf{E1>^{9f%>QRZiHWTry4N$*+Z`FC&5Sg`qk{RHD6q5JDk z7Ww>}WdAqWH?2+BfYFC%(uZ5ioYVVyVtj?wLZ)tgKh^T<htTS@=@W{UI(-YeUYP3i z{o&b7Stg4YruXJ<KK^!_W&h2?u^}(r>}`!Nzc_qnxzFYooy*i`g*%(%&CPh`uP(LS z(d2Wn%jDy4B>Oa-z6Qp8oFfywdu!r~Wz%lRKDMdgU$T7Kyr_M)S%DMQPU<$6o&0~= zn@trWJ2iBDt5&`^^)t1wFYeXDRgrE>pV}m6{bEZuzM363RW)SdwU=j>73QzuSj_)Q zN>jY)*H7J3?-MV2{rkNA!b}0SAG{Nt8BCPUePa*0b3^IXF)^;^%xujC7ju$+JT-Zz zt;pY^9NA@_aQ_SA-#3hBnPhWx!{50r+x_a3oEmeA*NUkt3rpC3S^0fVke}VP@Y<2> z_8QB3Lafz9HKPr_CWaK8RTk4ry|SWWdp`ezstdglTaItjf26W#@yD6Zl9YWskJ?B4 z_%~Z@^VivW&hha#bq>ebmy3To`+Gh=i;ciraTgn}?~3=^|Ah1Ww2AUNB-{AoOU3o2 z`)e43HXIaZt@-+@RPgK7BYEcMv!?I-ol*TQqx$>5x4Y-tF?`Tp8-Hc0<m>o1Qupr# zeog&-+vUY_<$ZrAx&G1(sFYN?{p-M_X9r7Ohx6-eN54<rb2;_z37gkmr~h>RQIOa* zTgaaK@Zx`a1*(oud=aRSv|ak${g(+XYMXl-Ua-$u8nWc|+^;L=+-v^nmi+s~iXSzG za}U}Z&0V+ba`E#g!cR)$*UjDUf8f@e$G3jnl$w3*FPGYW*LUXcw5PXbNv_U~dl&4H z^z?M}X|o!IjUT7ZXIgQjBX!z=u$9Ss^-39dJ<IY=&98pT$eJFxQ1a%k-}Z`Ew91`d zH7oI~^^tpYV7vC}-y-1;W2AT|JvMzW`_J#p*NuNBN;@bg%)alq>8I+XW4eVL8jD*4 zBNo5D_b{AQQRjH&A~V^AoAx~T`tivoOQnNr4)rQstvk3U;s;wo)S6v>M?|c<e;n^& zT%ht_S-(y=h~<&Gx3AC|(V4zp-s?Lahso|z(o1`89jdUq#P3afPO{_0nZfT;*PfXD z{@I5YsYbode`mbU|Cvyg!PtN4qGsfT{j0vkYH7zFocu20o!3^=_H7LPYfq)UNz9wl zz@JcC&FL$CZR*-d>-rZfEsXgZR`y&Wd-24%t2Y$(v0b~AJlXTgeqY13NX>0~BIC5# zbt`8~XIfVER$|Ip-u173)+2iOCYk9`%-^pU$BQyBl%_E-DBy0L#g~?3=76_Q+h+IQ zcAIpc_lNzL!imzSCW>UfIc@1Ak|1z<W4Zag!+pj498Fi0x2e1fSTb$Nnn1SV{eM@j zx?kgA+*$lFDSbo7%Cz#*($c-VXGyah{2n=}!&Nd}^Yk>Gu9p*^Z`6?RSKcOeGVDpt zPVbkl2d71zJGC_-df$#`+8678J$X5q|MiZ9**dpWn{+pAd16>{#Lt^y-rL6|sf<1I z>`orl((`J#l{_&b#Zq4|?U+w>;3uJZk6b@~;XHpzoH?iA-XmrIvOOX!yKm$QB<DEA ztl!kXsaCXJMAPQ!qKI_QFJDeJsdq;n@wSWpkh}9VtHYsfv5OoZdM}IDEB2dl=s=7^ zf27_tCGnR|XSf_o8XtaXD@?m*`19S%$r*>^&uy#RCg3{pDR1!6la0AQeD^T^iT`%( zre69_J>Q8Q(gzc7P1Q(OYTT+MW;}1*B&%bt(>~0Z8P`*Fa^dPZ8hj@{G2QzjVRA`g zK?$Gr*(=>s@~?2*4A@dOb++D*q}=jnFDILxH;<J1ualO#$C1akdD+SZ;)PP3(swlH zsbm$KtX+5dYmV&tYl}3pwOh)4_Aj<|6DX@U<w!W!YAz=25_ODo@7bf%o;Q`rSFN1S z!geEh`V^hUs|xO7HJeqfW*`3*^XqDkJoBR+2Oj_J&Y3Km$a?g^BewgKKQO*NprNo) z)>OLcam>jTB{unQGg1yZALF>JbbIO3K+aqI=l9!bHhm9ak$2qwVYN$x+>Sp5A0HNd zdtN?WK0r9*Y;nTxQl2$w$MS!abt@;TJ-Nvg+8U_tKcB_1A^7<0+-cKw?&fY|KgK;R z&-b>Nh~$<7|F;R2#Xh^q*6=~*ddL^^MQ0lQYG(*YO%Tbt6VG^1UT^oCPhXEdjlXw3 zZtw4Z&aWSM=vaQ3EcyG@*V*0k|NTi+J07>KPuX!^?q(%5%do|sseg5X9vHjaP`onZ z?!81|4XKSLQpPTaN*b2%{A*BH@3$iK(ybfiET0aTyIdFidtKV0dJX$VS&pJSwr5OR zQmvBYuW3A5nXJqs_UYi!;$^>Eek>3SyQ{xwg^}_Dr&UKbv>m*CP*$;b&Q?jAJ4-7b z$6S#<R_Sq5zD!unV})i8Yu^T0m5AkmrdP!Bzn_=bE>PkozSeX4L@UO+_nS($z2D?J zXHrAo{q&U!=4UK8*8Fe(|LE)IbpJd2ZAx;Qx~xiCvh5#V^`3vMWmhK5OSyAF;@;ay z9kHC-v|U){I3z#+yh?2L`M`hAie)V3dwbZ#zcpL1x^Qt%+|2X(|2{oD%J9p4!|UBs zt~j0R4l-u%$~p3xgDoqRC)u<vHp9&2+Gi=<$On(U&Z;au`sDu+)vO8$#<@#RG6=?d zS#Z93TXZ(qOonG)Lilc>vW7F<CTE{-2%UGEF=L|ar77oHS(w9CGH8XQ3RxtvcxcXd z&UC)gyd!IE(cFIrY!<WCA8B1-!F^y|*nz2+N;UYSAC@w2`&|8_X<=37KADHYmc|RJ z);sV<#4%YhxRz}QKOmPibwh6Gy3qQXjRLoY!q<GCtKRkGz&A;Ujeale{3eMR9B7;L zO+(G0XoJGi2Z5sHHQ9mJ?c^g9)B9(IU9izNsGH{AB(g55^Kt5e?~m85mTXJCRnghQ zUKGZWR5B%Pr`94%9d1nq4^!r#MK?ASELNDX^+UW%Z`9Js>`mT!^8PlAUm7GBB`r58 z#hPe(UaHd;63eb*3BL91Zrc*WRX3+`pZ)uyATu@e&iU8w4sP|#GbH=}v1#pIuwn5f z|Go2G<lW!oDF5~Mg7<spOUCD~nf+^*=9NQNKejI@lfBOHWuaw=ZKnI|Ul+LVX74;< z=cIn~r}mW3Z7csOatQpgSRgoaJ;#*m$xa9MUrlP%_TT<@%JO-V!E+O%KD}_Bwc7u@ zZ1%;hgV$_zZXcLsY80KjCde?q&++%a5{vIE;^eQse!e=+Uu)kqk=94`KR>XSO=bRn zvcsD(AWq0<-$lkPsYzYuQ;+L;8}HTK$d`E0z42XEdcwtU_HXm{Z$FhhQS&X^tn^~h z<w`1!_dm?mi)P|XX9@IL9jDT#Z63XJTEZPq*JUO1<tpsoF5g=yksUIvDIh%RZ0EcQ zzjjt}C;ywXTstt>A(vB8Wu{;Ig%fto7VS;86YL9vyG)FBi!e_tR**D!7N+wkRmaos z_R4|)v6pv+zMp1rjZ@>kQTkm@YQ4p$tZ7;w1Ux&PJ~M6-DY)bo5~{3Q{bWXH){0mD zi{E|`K0T4e!T&+?q4+I_t8BxIU)f)b|GIRJ3j4N*#}l+zczy+HG+JI<7~m9|w)<#Q z(8-du3R$-_J<gn#`}x)3n_OxEmk(#*(K`ovrQa=NTg5cRWxJBat*i^{#IN-POz__v z;IsJIgq|{&1D1*>pE&K+-CXnF3xn&TZJYcJ*L>#M7V}AD-EqbzlW#k%-7qy!RjbU= zzujy<KX23~hh5t5fArRg>(?&x**s~>YT=jv8ymHotl}@7yv=j;jls=zBJrY3XN=Y) z>~yhR8^KW;Hno*w(y>hGucvdCUG85fx`!z&s?E{-g%8_hztb7NnYSD>d&n*5m2+&C zO6^bSD@+U1Z$HgvH_&;l=)t@EZ1UG+Ud65arY$AA9Zxj5mR)<qHOYjtV%EC~p_0f2 zy27^a&7>zjzT8_ir%y|#q$@*X{$s^eXPEyrSHAnWAXxBsa`zO5bE+-RrH=e8p6=1P zd)~(xtX}Ob{71gc^EjmEm8SE{{*lZ2%NlRKeAe+iu<u=^T<fmS%wkpZQ!lPO=q&y< zVbXQOjbedmA3IqDZ+bnsHz7^-_RifLC2JU)`6qPm?{6+$;Lx?C!1UnX1=&9*Z7?mq zo}gXLFQB*dqg=^qVGFVQ_Dc(srttoH{h?pO^49Tnc1$lfWmzq^^wI1t-Lqr@+seol zHZd)dhYJ~)LKALnm)Xy#e4nYYp5-co99yZS(BnB4?!VqXmQ}jinc#l8FlAnZ7E6?R zNqcvW@~urgUGJ2??^tD2pXyvtdCQpXK9`KW(vp}t^A9fYm9kV6%@w<qQ<co?%Ng|I zyqB%}m;R5}1E!u@!tU^azw*ywTXw%EPwq_33{I<LOl11B+%G(<N#OYdrS+S17z-X` zHCDPl_`S8^qlipmjoRf!H>8*Oui7{5SF_<1MSY8>cF`Y0-c)48T|N7F*{t)6m6hI4 z*?FkTOJnKr!`gYCYs$G|y1SwkZy(QLT2%h8W|h&FKVB*`R`5?)u-)ae-J~?`D(jQq zl=!X7xxSfxvHm~HX2!C-+DmS%zp@h*${gFbNLD8bC>JLDPugNyX*R1$%)eXc|BhWt z<E5A%X-fKe9_NyNvU;}Z9tGC3`=n2AT_k+?U#Mx_YUy6#g$;=-LfF~X^0%rRd%kPb z$~t@UY?arnVkd{?ZjDk~J+(s&lqQ5QeA@cAl+S*3(V2OkDg~eBxxGETkdMtVf58() zC%(w<mha;27|q${hzfsnWM06^dyC=6--r8!jgJ+UxL<YqQn0*iOHX{FiDZ4k-2#5u z&3}I^-69-)VOnHJPhsa=yHg@Y0;{G4ypm&|>=AVB$YzG!r<dI1h+D|bHko~`rD59T zj0L?;i!$5q*sNse@N7R<(iR;3muqjC16N`3U!n6g+E1oWdN=7znrOqW$cicHx2Lp= z`m+eE^lm@sC?e6ayTmbOmq}vy+%tX|ioq&^J|#VFb~`mLSJ)TFN#_Pn*pT*Hdq&vb zHCAi?lv-}(zjtZMpP1{ieSWXaG*`@1p5o1`b9Zjrq)7rkAu;ciW@w+hG_}fQ@%2~5 zCSr4M6gzI2t~PtQFnhx3NjZ)`^=AHelK8-^zG|0-)Z&S@Z09a*KD6ccT+>Ti@BM1c z|MZymron!lr5{apFWu%g(Z%l1om69=<F|`{2`ns&WDc!m|2o@o{?{Al|K0mxYW<Lj z#ck6;rA~7OK{dtJ1q-dCzc^Qy`cK^_wPuOLzg5;b_gYph&PuxR<pA&HCr=9*&wNxb zteE>;(qhpSfr)$xSs$f}4c_klVs-Xg=HquERm^GKTffPknY7E-W467?89V=*dyX@* zKa#i;J6%l8i%UZ5Y2<^OM^g({udKelduzah<2QeAPu~$>t1K1sfHCMy`t<zs@7Cle zo6VbIToJx!N}kIGhU?j3UAz)v2@e=pPy0R6KOZYQH>GsX1n!Nz(HG~kn$BgdJul?w z?XkVKE>z;|MMv|8J0nl8ve{j2_ab&5Z;lbi^TRR%4_=2{+`*x!EOE@?NWzJ^Qky-> zChBMvu!S@2n>8<Zo_ogTZDxluc-~A9a&7(<E}JHh{_atX%c|u+YlApuzj~IHIMZ{n z)z@UJFU$27U72;4g{jK&o<j7AOEzcPlFBA<o{CDlWH=-Gwp#c#n+A?2Wy{+hU*NTM zcW2`9$_oB2oOa8nwaUeL-Q$Y;p3)!le|}c`v!%<w)o*IR%ZSotNj}OQdDq@=Q`}$v zT$6Ls8kZ#tpRPK$`Q(v!JFT|wOT<m2zfR?x(b$%7e`8}<QF?iPwPWpv`uQ8Z7QOaa zSiH<_-rZQQFK-^ayt0OKs!MODpFVSKHn)rLM$5>%Y|g!(3Rf0f%hZ^v<jCsV6zIdq znA_caVeRC*-^JcNz3?dh^0(5C)h8Al^Xge%n&xEuHE@ke*g@N{H|ug=y)1k$>L(O% zS8S%$DqGePCLI@3-O~2mQKg!beJm`>M4w!d;@RIZlaC|t(w~%sxPWyVi=KS{*7)=s zFXzr%frZY~)=p|)Emv=&Ag(^$@RFct=u_p*Om~_!U)w*Fj}Hl$`aHa?YofWF|B6?$ zPPonI?hxJ7d)iFoPV(EVC4382TeTYbOINu45}y5Xfk54@ImK@CTV7n4;H326)6*OL zhIdcc`|aJr<=b|0k!oSUhok?$?7Q)0Uq@4)b;tc5%?9%?7@d%uCa-m4Y0~{UTJBqR zne==0oG#tb&Sf5A`f>h@*^{@uka!WEC4JxX)~<_h4}b5Nc5h31mEZlvmlf7d*|o7q zH}t-U=62T&-*<(4s7yJx-kr;NRg&AGKeL{c`}VB<s%m;Fgk}1ar==CaD{tC4-dk4n zR@(c_wN-D#a%_53jwe_YU)JBWXUbMJ))S&{{DtDmSKZxkX>GY}XrkwShpT^Oxz|@a zPd#1QeNU$DjP9akbIv{4l5^~Iu*6)u`+1X{cU;R@@yc>}&4%11&x0~tvK8mpKA#j( zzE|?soFnoVil44t#FZQByY<P2sHEV?S>fr617^(-a_ZNO&f?s=fH(imm7}}<-s#Pi zetLW7kI$7JAMz`kw#o*sVf@5sc76A!D_I-rKQ8PJJ3Y0luH)+V?YA#Ju-V1k@q5<y z_dCRQo%{AU^}}~r^)LDCe-D1RDfV#d(#*NWcV0!haQ{ES%yHWA|HT^4GxrQN-d<IH z7w2u1I4%6By5pCdvwd9KKU=0=eO|n$dciBfC;XDP9~eK&cvc)){V<d%XyL?t5(~Nn zBB!abw`pqL@|wu8%S2GCTkkcO`SnKze2XN?y3{uxO`O<~VE2AT=-xvb^3#2f^KZKK zGCS5)RDW(o(DT?2+m2~(6sU0N+HQWyWBvPnnGfgqJ1W$&Ei~I5j{0{jYg&9UzG>G5 zbLS~WTXu4#*1x>RV!^p%L-T{5_FXK(#k~>&@p1i^jvr{-&gINzlaZ@c_@CvI=ks}H zHk0MH@qFnQ`+6rdCO`C-rubX#cWoDBc6aXZH2nH|#;S9jwue;sULR+zILWdoSblNe zgd>+qANun+_ub9;owQkZwcVzpOY)S}i>xe9rWmZ9e~aP4L(8S%nFlpiIDg(;sm9q` z>nZ$5(?dn*g2&<M&tq*8twkJ{7tV2f!svWap5eUHlE4aS1#32CpQR^~pWMqWRM*~F zl(DSX#cG=L%A0*JlD=zXtB6-$58D!aykVKS@adulZx!Fn_*uAe)k;bB`nsOe@A{Th z=DvNuesA=Jx}^`+=YFq}*B7#hQ~ANZ&tpOA@{3bsx`VXD=iJHmJ{vx5MMS#3vF(j_ zYE^Gk87g}_#82FNXmqk7^GhGg1JQdK{)R71m_ARqSgU$Pz`-}JZ`s_mtzqZ4My$|g zWYJscb-_X1_A=*d&S~5)B;((DO-(si^<3BDqqN<juG*rPx1W7fwve<xq_S!EE*3$% zCu}nke=QYLH9y-c8M@c;fdGpUM_9S(6S3y6(x+$cmyQ-Vdm&@X5x=r!>HX{9sPa7* zi<GK#T(xiVxBWksN!eNHv`U0m95~nD|IXF?KUehK2{-huY7V{5+J9ungu@H>A95_+ zb<6RY`P3<=7XPd?v7f*HXW&G|u$}kLvq*;CIR4sBlfU|Y&1c_u&GyiFZ?`c`6?!G` ztF)BkVEMnsxF?EF_w|}xNp1I@Y<$?@2X7Mdxz%@U3d}`AW=^%AwShz9={()^RUd9i z%D9`G{5uzA@9pLvcx%R?eap8M@4t94Z%dPB?l;ZXb9-c_yi(=}`}lMIvY%7FzW<fB zHU8DTxZtqg9Dc!x*4Mt*szvTSzh%+gy_5OVD#E{Qdb(-h<NxWRKbQNL-guOm@*sLm zg{k)YJ9p=uRXLP?mFb3YaH{H6@7G`XLj`ZEFt~qlvUOdnSF`8e$+({r{#)EyxR)z@ zmA1oG|Bd|1ge40XN4qS^C`k?87;d_Zp^@!m@Qz#Q{BzC;+phTalKHN?T4s)0>W%kO zAAFbY?`0P#%iWr?$xBhgWLMQW_GOc2p3#_GTw&6#W4JP1?AC`YfBnU#a<f+~V6}Df z64TuNd*Yj$n_t_W4|iR>Dd&I1<-^4j&OJD9x>0NYt0#4*I4`hI=h(~kY3?c6KNdf> z=Sr=3U)5=(pu^vv_r|_mT6-yHnYpb&p8nnD1#dQ1Io8*`GPg3x-lbRkrQ=i-n`mrv zu7pL^T@&HPkgx?m&3cN>8BWT)x*uB{ytu5LDfv#>)2u1e>OO4I{#Tt-@Omky`}4ce z>hu4%SyjH!c)IZa(XF39F1XM3YGK!J-#Ow-7w<YVabdMlo1$aqyQ@XhC&w@Icu*ky zu5MXY=*!8;r3U@S9(^|!JM5pPYd3XgMaNzVyE6$}&a9YTb}MNMLxFhBos=yM8AdAZ zLGoU={%eZ=ua@#TbN1AY_?tTF8i%td)bU;Gh&ao~>#voWn>j7bPJg=JL+_`-PQKe> z`rbu3sK2#JYv^;S)!jD7;KRFLo;LMceXrMxGR>1WJjxp{^v@(-uXID@3v(a;+I>F^ z6aPtM*6iafW<T-%XRNQP%i8#BYtCHyHrKb-<7SLv^`j}DuFQOPG1xQPC8zhDocBus z@o#gF{*bVH7IH}`?Yw&G?OH$0+kVnt)34Zfp31S0$@#x)b6~5z=!~_#o42q}=H?4* z=m>v(>92iO!ozyTUC+{=2lm>nx_CakcV0-rzxHPBtk>=W55)D^&WE1+%~Z$q{<r-9 zsa8){7yNkIf8OE7(X(c6Vr>pDdzip5&-3Z-{dJ;eA1`m3>(7e5{HjEiw_-H|1A_r0 z1A`38;A}}zYF<jQKGKoXxqXgZhYfgKzq8$UtgyE{;E;bYR<h^S#Aw4VQJq%5{QE6B zQBK?DfB*kmbLr=~2duZoGtXWxJk5U3MFD4?=hpRU0b4AKHuJLQX!u79go`$>lAYtu zwD?>~JNFm+%<oG#TCDDR(R`ct)t#=*+uloXdCg><`{~(_rp$d@;t!r((H6U9<32yx zLv!u5Tct66iRSLx6Peq;GxfjbzOnOK8js~-k8=f0-nSTEs<$DX;t4v@blzniiHR%> z3@*IrM^QsPiF|Cu)G*)t+ZF<Qr9ae5ybzKukWttCw%}U(B8@IqvFHTz8>f>JFJC&L z<TP2!Q}rv`&-Zf|>YVx%9(B+_XO{8jXJ@7PCaY9xEN$*RncX@;@@?xCzM9T}tJ<o@ z6K*Zmb!iKbkpA~Utn}T1tE;_I`jdJuKTNQ&c)i@fYi9b&>j77+K3{rKVsWD7RF<lZ zUX%HbwjGy_Hr(1Mxj*_(*y7JK4ztx}etM|jDa02Yv6(Tyr8KU>d};mKGj~E>t@_rz z=E$3ijdwzhRdM8Z$5a%iTw7vds+G5F{`0ch)ve#SyI+^($5^BXR^P0A7bJcwr)AO$ zi-`{vjx+I9NSfaCYY6n8<q))R_Ww*1KZonHTZLzLh3BPjzBywCYgNZ=|NSoMD-@#J zw6>Kp{cIDOc3SP=(@pEyzJ1xvt}l17F+eM4oop0G=|WSd<crtKmJ4>T3v;Qu=zhH- z>j1azm1~~c*+rFOdA41g+5E&?w|F7vGP5ame>W?oA7Gt$&fBQ*6qlDl@}rx=%Rhhk z`Z#{K#0w#bS}oSX)_D*7ALf;n+h=~8b}w9S?@qyQyK81i?cTTY)1zna>b~Vw{eJiE z-*(NteGeC3`yp%Rz9v(3&5buF{1{c2W_Ip-@>;%L{m48{{e^qJ+}ix`)bfYBbG7vr zF}R*INUU73@s5HGXK3<+ZG{)ATZ+{tH+jqbQM$Tmt3ua2@0n9o?$k-o`ykwPky*lZ zrjS;cge=2bDfMfPz1Er+IIlnK{HV@zKrtuFC{D5VkMW|mNP+WjtDV+vd%h>*$Q+|u zy^{$$b28@7e^$P}J*X~HP`_ns<f+?-{j#iX&ztdcs?p8IJ9jW>C-3<j7P9e6n8zd$ zn;h*wIcqm+woH{tx%kfVXUzMk9&dl4m#f@jJx{5{HEVFism;@`U|Mi(^4eSaER(I& zZ@&9|e8HLx|CW`mwwJc**Axy>*>mn>YuS}WS1MbL4#l@*{g|-e#kH>7`3LK7{MQeZ zu*?xD=yL6z8JYdleqPqB;x7-DeY(v#<&a}sc!`-%VBf^oQ7gM|zV)cG4pW(&(V=L` za@vmX#_5d`XZ_i8COy7)t5xr)u#;Xl(@ItwjVsyl`yP6{(A_Bedk#nb)dv-N)lL6= z(&asq%OtB_eBAJTLw)9%V`6sf$vf}7&SDZ2jg?g2TXV?Z^tFuLYP^Rvc%7Hp_i(+O z9jv4NXU_@lvtr_3u6NhI+rYV9^q*q)U6G9HC#q)}IN$KZT|SpA@N?c9ch^ssB=*al zE*J7YrCn6vI%6`oifV!Aqgd6)Y2x!$xtomd?^f%b742Sp{lWW>S7y6f6u5Z?cxK;C zHJ!hG<{Q;ZTuFjz2W8x9Qv&t>rWR>s&b#nnS)yInY&C-<b78A1FHQ&le^_y6yNm0d zztXz@UjMs$+v8oc#uTZu51h?6*-iNH!%aZEV`?1J<1O-USNL$QULjqU8=J%})3fQN zT(_f0?4@d>c^lG>pH&nVmwcmokVU24d98EF9RWe_J+t*@-8^{6V&9IwvQ0lNB@3RP zdHA&C<K=f!wpKj%?$%V7-)mkKuM@TX>(-wqr7RD*ZT}~Hu;1D2IIA=-yFCM9HVssk zz51;3^8+IT10xFqgE+PVJ0CoW);h_ze-7weW{*GY+UhIPE@$<7%&A(trTtxPih@qy z=dFb^ZU$cc_syzt?zt7(rtB=wz9#=JZ4vM<yu|mIwa3ZJVZFAS>jj(I-}jm>`X#ok zh_i6JWi4^BXZ67&?>DV`^{VnB%QgSifd*C61FW02Xw2Aa_;1Z!?lY^UH}6}mozEe+ zjwiGw=EZ%hQq`6dHM!0GoOiMgvrahBn%(0jt2aqV?U(fF_&a?@D(UjN(<e%loKvw* ze;~8B^@*X_N}s@;olj<&Y~SX*#=h|W@q<4;9xwc07uY3HR21u*#WruzOoNAt+e=lp z9Ms6Fw0Us2w*K$`nxngQ-ZQyvnNS#8xo`K3xz2u1=Kiq!pEzq9mtRKMJ^o#p%y(o} zJ}pjr%hI}fQrCOQJ$pVcI^Np8|Na5?Igd`oHk(f_`<a|z@O^9RGTlbrlHHX7XMQ*A znQtBLJ#)n+saW>&IgtwwMlNS9QF`!uand)H*iUx^woP1Dm#@Fs@4|E4(|hCtRrYVr zziC#e@od|qP~O{{eykK<AGl>xkGu3CK4xXRxO(M^iPq2ZKK^~Hx?@|B%25sRCC=`d zSrWAjHJ7i{?-Xnm?v}H^dnaY4%GtBg{miCoe(YHowa)sx+s>l*x7v@--8hrY-nG!D zM6!yZF!_d9GE4j6jWag2A9@)c`xoh+0Z<}|os}r?i;;oBmIX(JS6G^uR}x>8nyyz+ z**e*`?>6Y}1-2S?C++1YR$PvhI<ZP!!=T!E+tbc~zphIrb9Q}uvoz(n?*&aQmd3LI z^78VMibu9bRvLvlU&{AXp7kIxplrrkrMhXQq2cS@n;#sqD840Nqxo3y#fn6(Yi->h z^9|0{o%op>JJ)DqvW!ixO6(>7`s=H-6<yXP#_)6=IIkDysQ6>Y(E@{yS`$Sc&X1@| z_Pn)dj#E*L%c_h|z9qqIewr2w76f@tow7wf@LY_1V(X5Aw70SOhd=Cn?X0k5Nz%TU z@Wt~IUN4apUHNxe{ruFj6){_+%I_F`cxbop-}>8+(){ukd!2GTw(aTfXW8F^8cjU} zmi^)0+;h1gtn3=csebMU4>gV$NSu>hSYZ0HWcR_xdt}ZQq-e`6|H!z!;)}wEYtE+^ zP1qV)o8B}#GIjk=XYsYRHtCfgo-!7shs=(z3Om-TyRdNOiO(Ni>nu3)h&l6hr9_bC zWRtTew|_X*CoUDkBYRMQY5fwvn<2ZzQVTz(NMHYCAlc*`ClzUR`m&hn`|vXR+4|p7 zU93;8Jo~c1I6>XJdq+adedA}3d$U}2ZZ7(_)qVCHKGVfDQ{pO@Y)t$jcV+j%!=Gb= zY>F?PUcTGDU8OVlb6tR~*_@A+j&FCoD%o&ry6QdKy=!`oTKrJ#%3gH#f^CQO<!KVi z7BagZos{Z&KEL#<G@@kj$f%E!d9?n^9A*XveJ%zDNo*wxY=Zu5*y+C8W+Hc`KiJ1) zsNUL?v^ztLZ4bw!lF%)zQYFa}W~zZ^f?Zc8gc!LQ{{0p+sp{uM&Rv>X62@Yc?{=;~ z?Cl=-<>eXPkX0Eh%`u+dJ537G7Ntp@a(gySJ4RzFTXkB3>Eu?%%PFR3?Jo;^FQ5GK ziDG;6bT5vtQ%?u`Z=LB^r1i-y;51)b;9?cg_CJ#rzhwHQn}0XA`M#WV&`K>W3*&OL z!+nWT|9#b`{%T^BF<SWPe7lck_L^m;Qb}zy!k(UG;;`Xh6A)Ejy3;oNvZsivUTEA^ zujmQBq45n;6H51o1xG&Ms(ZC4-eg)-^U}maOU_BV&AK$pddf<#;}&y8*0j7ARn=Ru zlw~#Bj|86g@hd01Z#|rwwJMYKT*>4AWwC{pK>~Lo0+!x(l!`O_;pEz^dYFaj;8o|| zKuh-=#(C2tr?_o+A)A#Zt!w__rQ_{eCk{_%DA~E(Yeq$Zc;&JBl(|xRkH7Nlzfu#H zc;~tKYSWpItJj_R>t?cOhW(bh-?ERlt9Y?J-sNGRu(Ym_^H_a*#-f0|vA?(ZD|}k^ zw(3;s!Fj>^zc?|8dV1$2uelg0di#Z^yUJQq?TWm|()_ZMmwh_&D3odS{1e=Y{VEr` zf0dfv$@*WvsQtHc<<gGXiW~kO4CHXzQ{ndHLW_1`toX06r0u@Gw^PK;C%(Dh{^D)Y z)ke3f&X}z+8z<<OUEp3^wJ}a_ddjw_!!lK$6FxfyH8VdH=y-In=E}$8i&upPKRWF9 zVCBL~MF+R-*E}g?-pSV%wlyJbchNoe(mx;d+3u_jYdX(saOcQ3x&Djq_y3cW?O*Vx zo1ZEEe{pT^ozm7D3%G5*s*Ct<pW<+J!TFwt$Lmsk4qdWr7X9$4(fzwj|KEkRQ*N&~ zBmI3+qR)(=2(4SyAN~Hl|8*z7G2Go@(KOb{$&Kkdo-7l4tTb)&O{3WVK0#S;mZlj_ z-{$Ih-{1JyzLRnfjv1^n-Fw7ZJVGXo-@oC!vC|*j^Npv=)K0pZ#ZRm~y>R37c?{;< zp-VL)ZgEr}_`0KY3#-$OJozx0$B9QQLN;Dr$$wsyq4Lt7H5KQQ!n*~GJa@^QWGLrh z*WaASqL3t>u-oU|WufKmu6d)DmBSr>i&blB4kcb#s3J-gRRi~WQ6wDyd$uG_m; z*&jM75H#bz|82j{C$j|V**1jN?(qKD{^jxORE6a8maBFDJnga1TKlC}xm;~$&W%#V zjiNuAm$-!<lW=`?^`7YN-D{Y0tbVc~vh<pimZ(cCOes^C7#KoW7#OsWv$U>(o|&Gd zUP@+hiEd_ITE2cxW^!s?acZ%?kEgS%U$ARz(Am7(20Z)NKUfP`Zf&@jYx{25Ew43a z>nBW2v$WbKlwGjqz$EdK^V@F6<t&-J?REo0qeor+661vmroD1=_{0t>ZF)ZaqKq|v z`B__=x|%;M+Sa?4&2C%#dnLQQpZbqazr*uq@m*V9>$D)KU;Xi%jm#44bLGDNW9>_} zt+?A(AhDILx?$Ps<BRzYMz3fWw|IR~@xZ&<TsC8wg>7s{n_2sQ@m0>T{o!q!H79Rj z**UXiZcZ7w#&a)tzV-Z&Q1$k~xpJv#S9v$?nI9Sz`~Bb{9*48n<zjl?zj=~)<y`+V z#YqBZE_@00){DEn*K+1;mV?jd?nr3dGWoDTifPN_+_R1UzRr1L^S1i&E`ccLJHJmR zUkog~m+#?ecSCT_yNeP*5?kz!*M2wWpZMEY_v-y?PgWGYeDm_Y4x6m(vO<SI+vce5 z?AfOKUmpu@ea64-{oB?{reztttwr-+7I4k}dhg$b;(4#{vcJCg?kelthkxqsM4v7# zy?;3F{<gMLJ0Aaf@1bm2AZO1I{Jmh?W1XM%SK4nn9a^XrmUUW%J<w6JJ6-Uz#I!Ys zz8kfy4e;C;*dRRVX%_oYu84ij>y`SQAKdpkxtfbhQpY81SyOW0lfYe_uKU{5S)Q@h zJyH_piVS63Qt?=5kJ7fD(;5fWtA$nu>^6v4({ie_G}I#@z-Xn*7SEY2MariaJ<3pU zyBG6ByVJ<#@y;rT$Tf`h$Jj$Sy@M}tiRCmiszkP|&?paDBh2bK@t91vzs%l^%csO| zH3*+#G_h4+bI;nJJz5jk9$&n3RUt>%&on6INB0zsA}=A)os)j>7_JD(zp-{j1&2>p zSJ|dXu5FT;%Q-aD#8{sF@pAs%R)4NZ`OCeKKfQuU=YmvLbJ#4enR{&YKcqGeC~Zr< ze{?#Xi-AE^3GHkqjI{0R8sg~U7~*&~GB^LWg%Ei0e?j36vkju^FSl*Km0h@VOZD|d zzT&eV?UUgX)V*kCF0@ir)B4r={gV4mOwv8-$<I9}V96#IzV`DrpATJ3Ps?!Qo%^?2 z^X$@&la*hJY<4YU?e3fylIh$l*FWp}=NF-09#@nEi?W86`m9TtxB1EeYo+Nd?!~zk zt4;=KC0*;O`CIf;+cPyVq-l1Q&9CJ;H`XlCp0M9+*_DqkUv|x&qqpdOiIndCvvsQ@ z&Rs5@vN3qsL;fs7?bWO{Ho2K4ku$f3818&lZh7@d+tczKFT3E>9SY_zcbz(G{-@U8 z&T!$*h*nMGuz-2TU0kG||5>~F$rO=^Ws~{ipZEoQKUz`3`Z{~x6P_)V$IW!6ryguk zG+&Uh)Kg~i^c1&$d>-Y;>mF%xq*y#SeU`l_a++6OMB9G*`X5qp<<Dz0XEJph?YeX2 z)`mmAYVt0}yx8hnUiSTWJvnueYooDU;Ko%96>~m&^9FG2*<JC8W%sfZuhso~{(btQ zX>V71>gFHoXLEYp_iwtdy>7<5o#KX-zmKi@)AusuqmZTD?n#FjPHbLZ+R1n5-C4~a z2V>Q?C)s=woBj7pamCM<hp+yAJ^Q=v&i4mQlY)XK2LAnaw$jADB7JA<iJqqh3378A z3-4WR60_;`e798J@mGb3oNe~+=M`T{&zc6VnKDayWmbV@qK5dafE{nY2F_|`nQ6Yx zN-ogxb5yg8L|Vncj<hagUjrwXC(Vi<yrx^WODC!rMthcQ)yS-3pYW=}X4b^Xlba-2 z86P}Rk&BWj5t97)*nypC)sj#@US**KmAkK>r08l-y{@j#ZZ+|F)T7ChmL|LwaXpzI zw1HWYb@8qsp1VQ&3$mnG|2^806ZU`aqC*UO8Xj}?EL!{MQ9@+*KjQ+c8-aYnmRY=F z6BxQ{Q<&$sv(9VeloB=9n-IlPkn-WSNb71YmWOv#8PvCXMD)zDC|jB0xb4mwN0~kk zqXwgcjgrCF;$nX{8%|P@?DsnTl-o$#<K*eNO^sHY_gD4aJ9p|=k{ZKaNzu|RvCre0 zjIGa_#W)u*b(k72(Mk-KUgyKAlf?9H$p**#gAKkpyjt#ZYWdIG7F!>)X<F>iZ)%#L z#5n2r;$tC_x0W8|IeuY1gHLu=!>OW2*Io-wHO=9g+YxZ}@Z5IS$TI;fg?TSTly&!B zn5ohen&K?Fyj$iC*U!w;pYA9#oze7K+4J0oRdc~%@!ykFde7XRE7!TyVZ+*|F3Be~ zET)KGFrROu93(rt;*tB$`RtGXysFa{;J-R)*_oTud0(DAY@rf-Y+4LMnC1TDX*&|# zYm4Vw=U>~>a6&mOf89Zk%hPiIFq_QYuUiu++}8Bj`{8lR#Fdj$XK++L6m!47Y>SD2 zO4%m~t2&8o3oDOWIj!Ley3@2Tcj<>tk!t@-bHod#-%~pJX7|Szt-{7Leb&w5XmJUc zAUwlO@2y$p%I3vuxplQ4#<XU;%SOv?HQ+lMKFQ>o_oM?GPIxlyH~YeO^Ch#GKwx9f zCL=uuZMFiAx22Z?gcoK%df|3hOn|*{Z=j`~=92rz8fKVHj1;hV!FA}(DJ2ut7u5>( zo<<tyU#ZNLb15uW^<+#;&TT)ypB_7Fcf9nrB7V79>lap}p7Pz1|9^o~2<UuNmCU8R zI+LC(Up(o6F5|RG%ge>qh8wR|I%ZY#UXtSmAIBlf?687Zr#lHL+5cw=G|rh@JTs|A z^4ijyoF`R^!$9YwYE<40&2;DrE~)4h+cS+NVi$w*mM=M*>sa5fhzf1hb*pAo^wvo- z@OCg^xYZ=SOE*}j*}G%T177Xa$q$-OWo=q;>nZbk=i65urVH>ETxIlJVP+U`iTAxt z!NHevwVtgKy;d+|Cx_#LmhV?X)~hhIE-oxrb-nH>72dTmtKsbPwKgqVO76L?2`sMQ zvNPUrWkT|=F75e?W=*K!x-skQI_4m!Ep0M0-UuZ;s(GTYU2X1Pg%2B7U#;16e9^Lt zcibY`+q%~?yEP@nExYq9Tw7dHIIZE~_H2vpQtp#=+;13OZZ?S*Qsj;N$lF}8r9mfe z@x=4Y9@z|na}KQQJ|o3)u4}_lvBi<>6Am1!IPR@?z?J*?ixNfc%VmG<o1Y)<U1?+{ zazWwtnl}uq;}_U(JbLD{br|#1Id^L|RT-GPGm6-J|DMb=75%n3OL-TsE_S^;{YqJm zoPdhN=@Nav%l0Ym*A_50EjS@<wyNKyu=|GHlfyTAcQS{a2zIyNn9hCU+Uts7f!(4T zW+Wcu&eyhBZ82YZ!V$x{al7|Aa%mpyjWBx~e6)4zLIs9)ffZaeh22kBo@!iUR%E&T zxqYpWm(C+|=AxPzGnniTOsI@&7jUWF7wu%*zj#K7(R6K_#oxTSxYbYY?zi9l<nPZ< zn|J^G@p03}FK4H}Zoj~3+G%^kXie#D&%10p!Ws=Et_H6Zd1EbcH+iRbPNcf~*V8wY zOPqfn*Szn0cE@)=r+3o@H<@sKaLiX=X5aJsY0KZNHuoKS@4lP#x`0W}?-uiJ_np<N zBJ<QZulR}CrPn1$f3~b<U|JJ<YnA+iU90t(*D%XOc=J9~d=fGL?S>Uoo;|hMkXc>! zFH7GwN_gw0zkKV<`tI9qvwnMd;ps-@M74?ycOUOp7u$AnijTm?ll9+1?$|3W+~)22 zR+8mVhjvYj5MRuSn9T_rRuqUI?^BJxrqPtX&QD{Z)})&O1#bctZ0)T4@{8TO(TOQ; z_gQ|`wgS<fTa%7#lZ^ZSLA~Wxsi;$J;ToHwaL$eAcL|kMJ1t}Q&aU{wM=R&D%u(j* zPJ{dNx2-7KA=AvXNw-+~gP*>3VvNp_x0cW3CN@4)KYe4_VHpe08zJwaL~0KFoqc*U z^NLfvQWuwHuIJhK|MuZ$f|`8a)_uzTst_)7b@9t%IUgGm-&ArH^EgatZk8xAE59H& z@k%Pwt+MoM#|2C*I3>1i^SD&wo!!gF6lnIm|D)i_A6<)QmriN?CBo&(V0ymok>l&6 z?-G|Igx~Ecw^|(G>l&()qx^A4>6*pYZe7<reEj3$-&u2alrOn$Z@1%-HDBI~TJGBk zZ-eIO+kQ7};d#0{?^yHWZFzjt>ulBY6E=SHd|>!mHhRs}_mB4#-H+r9`!u;$xJV-G zi^Oy@{=Nsly&paAy*B9~%SZQhH$6YTy=tvK&8=Da&BpNXlkFulebm-&Q!PCEUw!M6 z==m%OQ%zG=yw81p>!|#jCstqiq<4O+tmofuow}h_+U(y6zIEYWSDuZwe!H{zWNF&+ zv)iX$@NaD>&eqngC^DOWzJ1riyz^7bJa^rdE&a5Azj4;W0~anw{I7rAQ?{je(wt8Z zH%ED2y{pa~a;^5=#gtF|UFjhWv&@!y|Jin7U77;t6tl%k7Us6x-WY%5u9T1eiO;J| zS=UIY&F+|U{=9yT3g7(OD|bDg?Xl`==1g7x^ZL*97jK^aTztLKkDpC3TnFEs{QdXc zT4rmtp4<7-o<=`Ey-X82b*pphyC?s5|4i<^7`WkK%IZIP`xwhsy-$0>{Bu`i{rZh; z|IbWwl$#m-Y0h7ksK?ho*J<AP!qM8zo5ApVc3n-hp2&Nl9oG93>b@&7GEKXq@Togy zS^43@QpKeg-7Xf{r~E^jdjd66wSI|f#(-|qV_;yA$J+%9_i%Of;Xd!-d-fu)x31Q? zGv_x48C)@bP~?BsTh~kHr2eKL4L@B?U5%6H&ufNVFxY;9`BYXB>(l4cUsQ>`W}70m z^jV5pq~;WnsmqF5RhC6-iiAkLe)jbF(pPHFLW>F+8BqI%`OaRyKzESpGB7Zx<8@d` zenEUrYFTQIUP(oX_DO#oEj>>i&1YJT8^06>dv_U4p6Ke;74$i%YbTl&fm6<&-Nwwo zFqaK&=Lbd;J;>GBKggx_%E?KqJVe+Y{H+yRbJ}Oq3swc=MJpa`P!tGan`tcg@#XDz zFQOAm1-IJl{841_-}8=(?IA(8{f7<K7t6HCzjaxlv&iz_`sZ(Vo@zOB<ld?`ny=q< zPUSFrp36N=?$hJbLBe~DR(uVz7c_DIcW8bDYvKIkn=?}8Pt-s7dzW4)NBQd6)}=dt z_8jZD{eI3x)>oa$5^gnX1$5lcFDMM&FP7DF@5&kDW3P|R>zY<ya;)*o)IA?=&xzXn zNTMv|Q_6|A6I+vxD>bRQ7Pk35FqJx`|Ei9AF{dq0v~N<)mWuCAmuGL*St7iud+wdQ zJ=#zD&TgCL-jvuXX|nz{V|yXn<J;CTvf5dfY<5J-^=#t}3O=<+!D`o+#rI;0>IC(Q zYA)~o<F}gS><1ZB@8lctizW8$?7qc0#qFMTQRj<^KCTlx*94UObjs2@bZ#HN_Bzp* zwGTJ4C`9J9<<9fvJ^X)_`LBXJuEHgA&Z-?*Y2d<G|MBng2xs<tmnW~(F!VaccxTzY z>b<P08;s4ePOad-&GytH|CM=+|CbFFs!5Ba9zI-tx4A@r>&5%|lcooFZoCs0dNrzC zME3{_cSYK{f6J?s*Naa!{i~ob*|6yQjm_I*7r(in#Vz(-ShdjfDSu<FZ$t+N7ssEN zn*VDLF+P7feNCjprdOIfx%-1NRa~`B>{8B7j=p?U>5{q~bDGkH{-fT@DreR4OTSB; zd*+;L)r8j5Cmw7I)&3?Hel52^yU<&XYwgo_JF+7)_ugLSI(IdT^QQY}ciL<@UHqDJ zrjq2!qP8|g{sXI&avTnqFZ#EnYg7C>&B7O|^DEBE{EyAu7ZV%yCPr#*wbhC@2lbNJ zvTkv_+G>}ur18NnE5@@W2Nj>!G%01i+kKsP{Yl>H8T-ob@94GN%$faTi&YL=!$!Wf zS?X(wr<8sZ`=ML)X41j$vhSBRNq+sWE5#xmGxccGx2==kT{FMAxIfjFD`-ArsMqy% zmaJw;EEOj4x=Q!e${p)I$UkzQq`c?E>+H|F&m>1Echp%gTeQpi>*vhxHw}0H=uzto zoAJ|#AxA1}I^WC2YkVg<567`z7U@WG4p<-jf5TC)yWZZJ`hOP+u3Gl`qu@s00JT}k zza81N)EV5$BIZRfPyN2I{zX-^oEsB|&R6DIgOIB~A8TD*GV{sXgGW1e-Ln39>+3f0 zyVBg7Zw5`QoW#1qLoPdCXOcjDP~qh^)jfwlG24H5&Ys9SDLRCAW#K+PMwy;9wMq-t zP4(hb{Z_TdV=2d6xqX|y{rg;VdHbqO_w+a)UB4gT&B!FefcyGt1}NCp2x63AyUZHh z2=u$7K^j4LTjMe0l1d+X?y^48A_Zil%)&A-u3QEg1H#)HZ!t0;U9Al>3i$$NWV1li zDu(T-H=}~g0pV?p^=#P8Lb@jv*(humi-L>+;cboG?AVM#yJi&GIBYj_fQ$m+ZH+S} zup5WCrvuq6Y=>Th%mLwTja`b^%tAi`8`(TisG*<m12PMQw>7>s!Db%Dc|XVoVms~u zWEco<Yc#aRY9RFR2V`@wZA=820m9oF&%59>2fj@a*))(h(O1%dOakF;jX6=sra_mS z;7GC{lhBu|fJ_15ZH<A^SWH4*yn^m1^x0yNF(ACHQ8odKQINS~bVJZ*6F~-m@U}*$ zBrJv?&n+Swg>8-rWDE#zYb;N}Viat~3E3D>qv{gcxHiZL5Z>14m4RdoMl%C_lpEb# z^f5({nIOEaQ7s#vxv&vNbo0;~e;~6!cw6JyQhertn}X<up*Pq-27&Ol#<VhghM_dy m&<#azqJazq;Vq3P@fr$gumyOtvVl~CivQoN3=Dg#Ks*3!<Zhw> literal 0 HcmV?d00001 diff --git a/setup.cfg b/setup.cfg index bdd4fbc..b1e5ccd 100644 --- a/setup.cfg +++ b/setup.cfg @@ -1,6 +1,6 @@ [metadata] name = toarstats -version = 0.6.8 +version = 0.6.9 author = Niklas Selke, Martin Schultz, Max Lensing author_email = n.selke@fz-juelich.de, m.schultz@fz-juelich.de, m.lensing@fz-juelich.de description = Collection of statistics for the TOAR community diff --git a/toarstats/trends/utils.py b/toarstats/trends/utils.py index bd6167a..650cf99 100644 --- a/toarstats/trends/utils.py +++ b/toarstats/trends/utils.py @@ -63,9 +63,12 @@ def moving_block_bootstrap(method, data, quantile=None, num_samples=1000): for _ in range(num_samples): bn = rng.choice(len(blocks), nblocks) samp_data = data.iloc[blocks[bn].flatten()] - if method == "quant": - mod = smf.quantreg("value~datetime", samp_data).fit(q=quantile) - else: - mod = smf.ols("value~datetime", samp_data).fit(method="qr") + try: + if method == "quant": + mod = smf.quantreg("value~datetime", samp_data).fit(q=quantile) + else: + mod = smf.ols("value~datetime", samp_data).fit(method="qr") + except ValueError: + continue samples.append(mod.params) return samples -- GitLab