From bb58173b8c47a067ce67b9282c7e00388ef0f699 Mon Sep 17 00:00:00 2001 From: Anna Petrasova Date: Thu, 23 Jan 2025 22:19:07 -0500 Subject: [PATCH] r.sim.sediment: add simple test (#4978) --- .../testsuite/data/depth_default.pack | Bin 0 -> 25906 bytes .../testsuite/data/reference_erdep.pack | Bin 0 -> 37150 bytes .../testsuite/data/reference_sedflux.pack | Bin 0 -> 38550 bytes .../testsuite/test_r_sim_sediment.py | 106 ++++++++++++++++++ 4 files changed, 106 insertions(+) create mode 100644 raster/r.sim/r.sim.sediment/testsuite/data/depth_default.pack create mode 100644 raster/r.sim/r.sim.sediment/testsuite/data/reference_erdep.pack create mode 100644 raster/r.sim/r.sim.sediment/testsuite/data/reference_sedflux.pack create mode 100644 raster/r.sim/r.sim.sediment/testsuite/test_r_sim_sediment.py diff --git a/raster/r.sim/r.sim.sediment/testsuite/data/depth_default.pack b/raster/r.sim/r.sim.sediment/testsuite/data/depth_default.pack new file mode 100644 index 0000000000000000000000000000000000000000..579863ecca98c84a9fa757a0de79d72b5b76aafa GIT binary patch literal 25906 zcmV)$K#sp3iwFpnv<_zi|72xwbZB2>WoBV@Y;-PgVPk6m?0pARRav&}RS=b6!hj-* zT56e)jB?Mu<{YcckqjzH1TiZL<{V1JoFgKNVpc@N1Zo)&!GJk|q0ITOx$Zs`-8EkK zAAh{AH~tgG=-!1W_v{tsnk(#mY#Q1$Y}(ahWLqB(Z=YeB|HzlE{_DS5ZEYPK|NL+K zxt*=OoxR3-q~>q<8Xg+vF^oQ5^J9M5H?j^48yx5(>|E>{+uAudZe(ZU;^5%q;^1hk zDfs$d|K;s7G;C1A|LkYLsa;&0u-dse+2vpJeBRE{#n!>e!QS59nf%}0*4a^G?eus0 zzlZ-YkFa2`|MKtsyPyAm`ejpq{{{G8fd4-@{&)SgOMACgUAwjUkNLl|v*Uk_|8@@l zg#Qi>E_NDg+rQKQe^C7Q85-*6YE2H`sKDKS#IJ(*Ul9Kb;(tN>r}*ElbL%euBk|wC z&hfAD-@YLJ{~cdJ9)Uit)*U?pJ$;6SS+@ud@(msq=n;Vb8ti2}bXf2ZSL*;TFJo_y zu;GEO)C@@|2M76h7vNPv{4a?A1@XV&>;EnO_vqZNd$<3Q{NM2}@!!d@Apiec zzJ>=44s*5c$o>+JH+@!#3jsUZIU9bd)ti$NW^I>@!DTsz6NuU!4)Dl0?=6XiNbuFK`RMXr10 znk3f?a=jtfCvttKUkwf9T1u`~a;+;@7rC~S>o0O0C|7^E%GhBzS*{D@x>~L=ay=l| z6uDlOYqnf-<@!m#7AY#%vU05^R~xx%<=R%RJ>=>s*I>DhmFsl5E|%+hxyH%$h+NOg z^_pDo%k`C9zv)*a6S-Q*wWeI{<=RB99p&0buD)`W#{i5X*tF*hK&T?%b*RFCMAlD&s9U<2!xz3mCD!Fc#>wdYOlQVo7y9H}BO|6?WUub&g zVd1En1x`^A{|n;(j~xFC=Kp@!_zznkj|d$4cgBCazr=qR2bY5Q|95=Z*;)q<4jLZj z6I$Tz1@XTi{{Qgtzu^9lA2j|C^8QbK#=pJ)+x9Q_f7%zs|G(!;f5W%_t8ruF(ZNAJ z>bJOIdsx^Y`K`UZBdN+bG<%7pXJ_m3*Y|%B z*cZhAzvJtl=l_L%HC(A*)wyzY4!f9~M>n?AuO=b-)l|+|m#Cm$OHR?RrR1bw8Ts$c zQ}t`Pdiu4zJTG8TLcdlVtzRoW)UTE0{Jhm6{aUq_ezlg*S6$BI*9g+DKgs#?T5^uO z_FVn?^9TJ}SDtLGcS^rDkn{RBo204AxqiD@`qe?6M{tyL=T7(ZtBX7r(MX;baFyo} zv~o^gxawEkQvKTGg??@3u3wv<)2}TX>DN{}^lKY={-UjXp7v$*YX^Cbp`)Cq?;NOK zyU6nszsT=@y{=!o_tdXFN^=tn^`gOnz{pwauzq(8LdP+Tc$#W7u z?e(kg@A`GnuljXxhJN)QtX~7%$Mf~=BL{Iytm<3_jwvpd?}5lkY;Ax#KP*I?J(AKFU{(&b=5|DrgXGod$5+G znHa3S_$*5sZ~Rmnp*^78(X_Vq&_omMjJQJDsV&-SM~M{e=57nL`x6FekG720ei;+2 z-17cOS$Wq|n`3iS(FP0chAs!R+d4Pb4(aRbdUs$GZTm4B?oKgP3Lo`XD&7@Jt*F6D zgU@3XmnYuJ!0dI(*kZBDw23v8InKGtM%y0B>Vw~vHN}T3mF_gsrubdJN4tK~UUWSj zw@7Pi*FbA!`N?&*Nv5_~kEyO**7U}&ryrGhKUT3$KBuJQ4pCC7ZBw>?Yo$ys{akS| zvQ~^er)!_j@z)+NG+%rF-9hb%eK%Lt#gqMd?147&)I{xkx0TwsI??pd?8ln(i!+uO zUi9@U`KA8dV==K?7pJwYZ1J+GLAcQuQ%#|MMKpFDiu5-dxin&Vp~U&+?zQN1w0T&L zN3UWnzHdCQ%6D&`d>d)Mx3D=0;43Z+8!ABx$+ z3QE<~%F54v;Yx#&tCY5_Mk-TBzffX}_EG%$n2@-%m&ad#NVRiHOwpe2HjE_M_q?rk zdGXQ^2+V?y=hd^h--RQHfr9}&+^^$2ywZo$oi@UXzX5$|yotjluT9zymdz!Pk9(IbVQG1SAL(-Cw8S*K781t7cc3b z{h#E)+k#B-{InxDAB+<6!(L3I+9(boE0uCt7SKEF3z8&R!DSf@NEUc}ayNhsZ2=Nf zig%f#^lW&Ut*V|=U&a-Zvf^K@m#%!OzwU6UE+T0_55}mvt)J7jmFn>}kPCY^{dE=$ ztc+ZYL_!bB4hLt0z)9|eTlC3*IK}YI%TM;>*iA%(Ka{fRL-6?YsEUzUCiX{4gK}X_ z^goW~nxf_j4(_!Qr~8$xWI8Z%>$&=~O0`(OYTDbBwA&fWwjF31P&@p!$L;WY+GP!! zX<}#i`Ni$<^eyaL`Qoig+2U2Ot!_)*mAc|HDiT~kUrOby=SttYfiO)iiJI+=EPxyaMnKPs4~{p6sgiv0 zv7~9TdU`fSIPx^elVObn3Tr0EW9)^slF`t{5D@73=s}q@X}X;5eRsj$jjxjD@3hl- z-GAr&JZ6ls$&J$=zDBR#8)9G0@{ql4U`+BC$44W#t=xKOaa@W`Xq4vtvw9YXO_ivJ zK{JlKy)LdwcYv`bnUC!;Y?3zlndnI&D;ZT-(AWt*JR}T>!jK_ph!7%$LBWJbH&79Y z6p#Twyt%}Xe<$r32tg(UwIFWNzZX5btrs!z&6Jy)i!0q1T~oq8tP^$iloSS6=PEU} zJVB)LU0g*m$Q;R74_t>ULH&gLBqTT&Sv0gk;)9l<1bRf0DRhh=4H1)3(}U6PfK6-3 zsOg_$-_@*}7G7_?-Lh!u@1&4;2!w!Xz|!;v0XZXLURvEu>@9Ig>1=+0 zjhm1e(835_zxor!!^uyH-0@Xu6O_ifAn}ntNjoGD88*pf5E1ctP!o{~wrUa`sg=A8 z4^FZMwvkxzkMyn&yJY|&qi51?U8+n<*>`?r`qfEHhz~YuV|vRhXK+?;_n#|T-n!uw zm+|C{kIpXmd`e>Ih2JLonN+Usb<2Hz_$=-D2e*7C9P5)-#p(5;14RQ0ClvZt!M5+V zqN+@Vv1U=pyf6l#L53t#JT}XQo*PENvVy!}!Z2GH3yBlNh&CeG1RI4uFaQH!$M8x5 zNcy$lS|X{|1|=tbIh-1h4cY*B2-g`ah;b8^KM+MLUQ*WY>*Juh;t%q>%tr69WGI0 zL7$55iL-akwtMig{CNNSo>fMdw|QCekD?|+i;quil07`pqiV+W4CnZ&e|$^&vTB!h zfHp3~q)UC<X*-$1AS%IpM8VrFV2kRBlhPVocDVQ6!hVTLJ zhaAb)V6mhOG7VS=1Pd$$hb5n7jzw&vcFqOYu0mfl^itg>q$I{2R!$41N-nvCacoa85fOy~$U+U}Jhh^bb3 zgRnJ97x(A?qAT_0sj&3upg8u*R=Q5JCm(?7Nj+*9CCj0w6Wa5qmG@toLnWjJk}0eM zB!Y}x`SpDY6b)hq1B0Ew-AD=KH>5zelG4YkbFDyAuim{%70$2C-F|ai;IwkjTR+^c z=@)$7)_ri_ac6F1<<=WrCBb%2k&%@uw#XirI)6-jhk0MGeEjXp7(AXCqKBp6*|9ZxeuxTNf(3$WkRI4) zO2+IV{R#7yo5WtL+R6-vCaMnEGazv=4RSF&4a|Uo5%aj6Nh!J_b8qM>22@qjcQjY# zPIpqqCJw1NKN4O*8k4BYmJSdkqYxfFGL^->J1Hh%mDse&p7iV-)Cjm`^Kc4<00;pw zm&^(eOt48akcxYHMSoftjmWYSIVgQ8&Fa!i8h-yB~pre&`YJ^#iI;_a4<@05iQuf zf!8o~$N@G*iXvIdBnZHbhlc#HHxLd`3Wyt8+m#$DmcDaPCVe(%o(!###vo>LBgmC> z2W|x_3o{?$9wlZS@)c`mG*PxqouMpwu~Zqc`=&C`;*8?)c`%@fGz&6Gu|n}o{2)TE zTu~aG%LN&PA(4GR{V*$126-XSa(=`~WwB8NC+$eB5a&ss^cQGXz(=rf1XQwf2&$~6 z!|9CcmQ{Dmo;qE7J#Dj-pV!;xsqZeni6}C@$5VHYQ$xMBzuVr^y-vu5S-W1pi#2mt zyEK0`cSWBrla7@-z1Fi4f=r?`;auHXnYggK(#I{HgA`DO?1#f7873#52y?I`SP-NJ zo5qGfrD5Z=eJ~D^BIFJQ0VyCtq7!0~+au*d+Eune=mwrlQl>2utbjI?af3N`y%QyV zdp}*-IjW2@Iq9RaW>Z<^_@WgsYS;+LjiLsrS~xBGCVC!l&NE=L9O#{7OXf}5WJFBQ zr5`gPK}>{!!46<&7JQyqO74oIMkhI0m7VlAMdus zkkcW;4XzHwM7xRUMy*7jLF+jXkpAdtfp3IR@K}WIx6jThr<195D-MPbL7gBi)c+P$1SrNU&x05SQ3NCoXIV z@|u$)ipHn^a;C$yJ;sxcwHj6ZMyW@>=?^WuV{hD8{<(C^oj2_FJvUWFZl=;ct~w!A zn00psoQDGunJ2J{6R=}%iz^cweO9`qEP$EkRS#eXGWmzd3AgE4VTS-^dU}>P%5=~Z z%9li`L>sFJ(YN?Kt*pdA7H}1rPy&o(-4xsgcm%8K`DGT88j_N-_EvR~w&S$8y((6D zs{N+yS-w&E&2y&m?EXm>sx)x;owD)FF0rrqRZsQ4_jelHp%-V))xC$KVUe`MLDA3zEu4!M$u zi8n%blz%9{{u~PEQgc*@gG3=FJq#Cdft!{67Y8D75n^??GBV&Gm@AArxOjhIzTgQd z5OzRk$CA_%kTfyCD-T2|pHy)^3bM&JF%hRgWNc}k< z5sAu!ZgwufSaMb zP%cYyl$i+a@E7<5k{YCnT#j~)zYZ$|nkk)M(r?IP^N=`zD&#h)1^y%=$A##s57?nA z{Ki~YYT*}M^+H~v^x8pCuc&n9gIGVp9FR!J0oQ=+$xTTA%66+b5!A6ZT!V3tdlOaY)Z_lv z_J^T#ht zSuAWiX_OA)gfLrM~8Yqg+Wp}KS%~Xv)DIjsePiCyM-A}~ij5@#Rv2Nmv{Rf+8MRzz;HqmJ> zr*v?r_v+nH|7|t*^*Npg)CvG=?2QyNiw_cw!yYLge`_s_r)pVHq(cHY zarg5yUAg&R#N6&TS;KHt7%b8=w7*&%+pP3Doy@f@)S^K|$Z`c*x7lib;iC?m8SDnsWD5>hAX$6&J6n zvpl?99)9>xB4eJ1)ew?utY}?WIPDkfrd?-30ha?-&@%v{fenxsoK#H|LC(=ox-c_W z`D5of<ph3PS!^4p$8tUUsE7BHH2(6%&O^zlgk zoLWKXAeOS^N`I5&j|4>o)oLMHH2tb97=45B36l;g4uf&mh^aJNY2lBG5N^)6Mp zlyZ9PuAms#2~kc7-3QHxYq7siW}KtRm;T{b7OeG8SYvz3(O`RzG5N|O-OJ?lPw+_) zt?Z$QEm!=+;Xh7svj@zNLKUJlOjhk!KumZ@k_*{B^KrEBVC|33ci_tMi9=PnWjb@_ zH&J#;VI}inb)i22#zK`Sh(U!Qbd(gz6bYs+B9^=oi?@Cjn?;h?v9z2x zEXIl}+Xsp*&htglfW6A(AyXOYMhW`8qy;D8!O2g!&RE(ksx6|X)mBl>@H#w+?1kh_ zmIR0b*8=5*`3G6g_B7OIF7}b1#aFvIHeygHN`hi-2Vq^uErF@N(YeHDDc57WRY4B+- z&?8wznGU*!CZH8H?}h0>Q&L1A(=zYY${2(=01xgM_(qAiH43#K1V7j~S*Tu#fHN~FwEya;fxnkA%iIHIoMa#@AI1{w6WYqv& z=d5|csMlj~HTs<}Vb~^wMkz6ZI`bcTI6XNf`NA2bmvWsZ@uU@(y}a8WL-#D7kS*PV zgL;-mOVF>zso(8rY;H&iZH&l-sk+ znq$6Y^Hrac|17Jb2u1;pyQlPsVvW@$Q7)o9XXB`xP zZQ*#kn&{oRl_*(jq4MhMAkk)81JS^7s3_I&qw;W5jxzm50vj(WR`q$AhA;$CBCl&= zQ%AHaatop)O)^%hO*sS%QU$h6jzsQ6{{_yA+yPyBxWbjk$C;0fx4EW_y6vuK*j|Q> zEDAR~eyCWRlOqq<^;mw>pxen?wXz>XtPK!O_Kv<&o`$&sfHv@ZzBA@`6tc--< zuSmp6GGv3hdrnjym%l3tzi|OsR1_3Q?ymZI4XxzZ92iIy^fY)dxLG^P8KeY2 zaYDymm7`{jc!FU}*Jyp0yPHAfkocC7(F5%5I&^99q0haO_FoRKv>UoE?VVFh{iQ|C4m7UCtHG4|r$r(?NDw}ok zIWZ@tx7boUN+p@>e+X5OF9ZhpLY6QN+B*r4#SGws07w?#ckouq+4b#1h#dQtJx53o zR)|*h_y^OGLXp_$J5nEMOH`>cLV4b_4wt`oCZ&j4CCdx@@57bUw^xwNlo;9;BRVHf7*barqRvS9hE9=V~rxMjx&mCp3 zxR80_75jH>@mp9VVd;+B6C$IH+8B8m7zDmE9DA~|cf^4^1Fs!U-W#p&fcB6J0KW++MfK-dqGhnNzBV1&?^fDCk(yrHf23Ko@Vw*Mc|>OB*&7`n zSvGX2krL6X$$tNbeVg7-o3**p;oT{}+|_zq3A^Kv4~bj;NuB{_An7C`k=RM6WHhYJ z)EKir=$`zvawlg@Y9SJv0TM$77`1@V0TSv!5OxWPsZcz4bsKme#0`Cb$;o;I8$5Ik zKOi9>^deE3wAMqJlii)^6f{C+Ko3U>B{c!-@mnu0)GAN^Xv<`h6e@Kin#4X6nV~;( zV-Zmr{0?>MeY>w?NDTf&&xUA8P$iV4IkX!%0$B}x5~nVc+Ds9@INIsyvx`^E`B^)h zSIu~lI=SNVLT0rFMegiSIdykGmv_F`yKL-LKHJ;n-NlsIv;4-!FCTY&(&lA5>R;$w zc}&;5*|{Ue8gD%k!=hnmu;?+JTZ&`jM>6!+iW&@3fkp}X`A7kw584hmj$|Sc5n}xt4&Hjm07ug!On*r>W}jLcrx8R$bL$(_dckc=1?|o%d@%!PBjk()yZeIMgThgzO1`hr>vsg}D@4Tr36OEz1Rw}27kOztfen&AD z#U_gy9AX`EB!qFb=mH#;qjF3@$yBv6tzbA%`{77kob<2>B;w+;DW4G)MguAJc_46YY>g}!BfwROy5Z0 zF#{Z&1&ZQMg7yg~fN0XMD*U`6-Ivb?tf%JQqrYmH$2R=NxG`OUgzEFCk%Nx-0WDff5^XO zl|v4tAV`4omaDcg4jS_sL2D5A#d0qqMDosf?TZMg~|EYGW`w_&AC;=qG(wiYQ4{R?)}NZs-FbY>1bfTy`MYV@PCZS*nd4 zs0QT{v@}VAC{ZJZLb&iT_!__jJs~0%Pz|}IsNTQ1&Sc|wvHa;5b|5kXH0MR+nn%j@ zVEvFV{W&a`P)WJ5O@S996(%zvNTDC-nVvpLaj^G2{OaqfB@MppO03}dC3jo&TOa@O zFY2#8(%aws+%La0+?8Q9FnUa@?fvF8eqGsdUy+ddM`m^EIW2FjInYECr3VwF*^eR? zi&cFpv2;P=P=tY$;JVldDjv`RNC=u?l6bJBX~XKkBFhMxfKdXBkc?MNtxV~N(7X+2 ze;_%-2BA-NfJ_~b1hC7z0{TPemhyd8GpR++mypQ8*U(nzy=V-OmmY+E3k6j67*ckZ zea*!3tlGMw4?pt+4?-~dUZvL#5)1lT@YFQ=nG3DoD@nPvkuGS}K`Yn@K`1>2eI)v+ z45Vu9mK4k<9W;khPZ%Yp-O?_tY3UFim6lbha^T0B-m`M1X#0+9aU<fJ_mykS6E8Av_67VKD1OvuK z4g(&C2_6c@PztIoFc1JTqz~Ig^nnp!DsoxPG9qi*EwR1nB>k+ID3eyyB2K0e)D>S8`D!J?b?d!5WNOHIp%KQ~z% zn|U}bV&>qtTQZu|y>X%Nc;oj2GcD$wYS_ytU+gTa7dy#IEZMV0XISbje>O~6^;t*? zaSvtQ06UN9aN2TV-$NWXOgf zP<0%Uj0i*!_Nx7&Wc*w;xdZjUM)e6vcv=9j*9 zE&8N269;;%!uT5}1u!d>aA7Nu5&#)u0K@=@*e*xduKEK^aX%U^eE(Da^4d;7lY1QV0DMjt(VYYKEAPIu;0PBx8V_ zi85m%Ugo7}P*U0`Y`5HJ{UWJ^1~6q3;}k4bjOv8)4q-T{6E*~oB8SB6i)=0F$+usz z-j%Yw+Qg>kxEmSuw6s2Pxw3ZA4lSM-BkPMIo_COpC{pipxuf&wM-9$UPL2Lquia;&c_U7Ajb(&R&#TuirtSLMy zMIW@aQ3WI)g+Za8ixD3DOL8ndhP&LY`oi8jyWD_&F*h?^2j02+Ag7b1@90h=XRo%O zb8F=Kk7MQ~40Gsi@XqRS+?2dF!A%p*6+TF%jH;#OcF=*1$CaclZc2PodoJKa6?Nh= zF7OzM4)iUGZxWWQ*M0!uM7g~L8LMbUVNUg|h3(D%(&F>&Fq=oxT^ zgCq_?8k6}H3c1dTNC8uWRsos-6ty=&;(}JymJb1;9O7X=RC_K6tl+8;8{!N3A&3w_ zhya}d3(cj|phnjB1o0un^Typc?U2X4i{wKaqn)lms%YKKOeZH7?DuBo$ALrwC)7c8e zcmfM{j=T>x3dun#=pP~a;W1$EYEvEITuxzZQ+g_kw|-?8!2lQR2m>M*8X0xlokyqV zUuljWz3_T*SlA!j!lfCwlGl{7%#wkn*dmk(XL<8o(L)&k8!Ori{w=S*T2&R1sNEh>C8VZv;2NwvO!R3DX^yA>LVMO6kv z(uhm#ENAerBY8d}u1fIGkEKuLW1N~Amf08Axs-potB`5PTZ7A|dqsB%4Z8X|xo^W+ z**6a?dDwHo?Mip8U;O%uD*fRmnl~Ixm2)+Uv6LYWAUW{>r1~hVLjf=iU@i)>WLP*i zfguv00|W$TRfkC+C^dgqM~LZ3{j2p7y~DSQ;ZNG|tO90LkkX^Lg?1pI46;+lKh#4a zq+0kHOqKiqPy_7%7Jvf~y6o&j=p=5~_tsaY+&Tixw|w+7H=ICKFy_^Cy{jlSp|_aq z{FD<}KqGi6X;t*DnJJEC`~*$msaI?`%T^~3C2@g>{LyWTA{Jq+{8rt(^rFzi+B z8>pw=jfIx5v)X-*TuiO|fu@pL5n1F(IdVSe9PGa9G<$>pu3XRW z;-Os#+*F>SLTN^1>>I|@F7SwJu?JA5T28BO7z31hkoaNUpn0~Xc^N6VUGJDGMLZ^Gzt^BR-n0$wi z^oJ99@Em3fnc~Ah4*@Blk8%!?9v#J4cs&FGkdea)vL?ti;C5u{@I1`9ps9~M^ws4Z z;#)W$N})#(Dea4hHUp+`(v>S%0#Y5!1fM#g2%a zT=9cX)_(hS;eBf#`e>VGRnDzXhE?@&qot!~iz}Iv`;=aK;`+J7*g+49+fQhG_)_Ka zEAj}UkFn-HAL&4Eg!3#&lO&62PDD=ie3UvZ2T@8+s76vq8jlD!mURv}8PD(l(GnjK zaVN`&zMtI5{s>>dcOhRAl?u^7AQFkb`9Rwx>5|OVq5^^eYGD-okfFi3;d>OuM8BsK z#Q5hm;SQu-6y)|^EW`uG{&$Y6Cx%hHZM|d~3W*e^M2Gq-VG(@hT~l=Cfj%08Lw)-s zln5MN{O(v07w6n%)4d=+pB+_S8BF~;%AmZ_l*H~X;ek1eM%9hj7#o|p!Xl*Kg=6mT zeoqf*mT}&q_krCGet`$_K*C9FKq?21^D#D#?`j7fvw7gzl69~_t1%Jgif4v@%KomL zZ2~Gw!v=YWokP;{u$}0;%|@K4(VP=0k^u}{tw(~l!$C+ndZ2-}2ok~Z5Go=fQ!&8J zM`Mb2;r<-R4S*CT3?c+Fh0s3FZUGx=TB8oi#|!gBWSOU8b6PdVH~j$=SdCbitcWci zbD=p&UU|G|*uv4Wx5opHtxubQY3Wg-s>!{0UOhLE^E<0SRKq#<`WxC$Gnu3{SU!K_fgCfJI4S_E_6810}eL!Q4Q zg^?ixiNRJN2TT_Nbx~&kHNY7hwCVx!h}#4nm~r(?{J?bv!f0hd+DBz`n=w2uzx#1h z5nFk&aGA1;q)We8SX6Exl09?fxegxSMNt%tNS!EMPLOm^q$P1AJBp`NpI;1IW8{)@P11w0WFgt`(cmpX%ZOSDZZx9P_ zWJq|zuSAUyoD~p7W-hA-w0GH0mILYP)Hk6QJ`8##y~A}#z0fuY8Im8FN-=E(O9;B% zU(8fCA&F#|BM1byplnE2XxJU`+tCeh9uVCns|Fxt0Bj=rUXAe&Q|~Mv)!f&sknUmXj=mW_pHEzEb7)mudUyNi zh82pa!tPvLQsw2PxoPpYR zZmw31F=j}ZjH$l70EwO^hI6u)epF8FUnNSs>a8oT4aAfvGfGe;Lh41C6Lbx{&g9KN z0aHN|gRJK=r38V2yOZpd-+R~KOHnl8wK6)*D}G|ilhxOIhq}}l^JMnIy≦u6EvM z7gF_5fOk*7jK1p*j45^KYNt2*>KfFURCVEgJhDMVzDJzxt=7_fsn;KbZX91Czy!rAO9N>*S! zWogp{9+tYj;jy@Rb}3YWau6w5ELeC`jGVMl9=l|N!jJ_N3K@aqBl4mZvMV`Mk#qB= zn%edDPU(GWe28D#{^n1A8}Oi9r#W#q46gTodb;EK)Vd@8dM<%>BBBn(@o)glA3Oj? z3DZ}5FMtYptcpbpji_V<;+1W(yz~R44+$?Yi(G=VB*)N&)}gU@Sg{Bz03LxUDHuGI zCm{!lAm2fm7TQLx1UI4tg_9CFu1Q`%a=`=0tQ6%h)P>OMhTI{4R2z`sqIye`RPMFz zC3ZHq65iqngB-anjBvo|J(v{YylcYp2V6nJAO=*0c|(dSq-6cM1QW}lahhY%W5z7( z<8^#T@#y9m^%AcJ9qwq^u-(pDNvXLP{O!XwE*TkWxT}fRz`_?t)N#mHE|i0SC>0;< zSV7o1=*K-_#o*hLT|h2mz-n0qqJil{@UkkzK^};Tjp8U6H28o~vT&cfNiGC)gFZ}+ituvkx(;J8W_OVlZ z1(Tf7heLg=mmfN7p3l4`MI<^X?pNL=B7Z5za+W0-mN!sBhAVtD!X$*P0y`{3ZN8Hi zz`22`7%^6-)qph!%5Vlai7W{M0K_+8%}$Sa0?m?hAe4Z$V^{=m5AQ{*7syvDYA_13 z1PAhfh}uQvQzt|^NI;|*JRD>Wts_bRmci9g`ULbqtxGx;(bcN#q+6MBPRX=+1k56} z;3P)n>o0Xf3)ROw9Ql$ADv}~D^ziJ}B_k*`=TyS+i!%|b_%eJ9!T=u`Xt+V}ZVcZn$Ul%j3XZzAZ(0mgB)g$YL`!^QSw zrC7aQ1=5G;AZ>Nnj*Jdt2I!!pNrsscoGBqgLId~GrcANXv>l6q6a&MC5M_gdK^-|B zzzl$ZL0(w|0U5zrNN{lyrdo1G|6`ElOsI$c6S^hK03to^UW#XJF=LC7ILxA>^YN)= zF;T+uzOWoUS~Rvx5CQiVh^lQPc}Bl9Wf3wN@+I;ZY>RP);X8Yh8n(_yk)iP7(eqBT~Ft0 zaF*o?AxCWFfHT&e7RBek;GhSCQ%|vDdlWwrCkb|8@G#W~W{I$Kgf>JhoF;;4;aJq? zOs-#p=KyqJhOjDw-3?Hgkq1>di9zlq4}bC)17ay$1vw##S8DPKV?$)e@l65>k_8YO z3h_a_?bp<)+^wW!t4K7#fag`%bvL|3!4p%oNo= z^#G7UHzF`|sOVxiN<;;8MFW)MMVB?b5v@ToW0#M?w5V2Ik5NjqFTU`E+({ds96I>S z-}!w1&I-_*9+YyeHy%>}$pg zb?gsH1|gO`F(!71dEhbX9U~x(qzQr+U@v0ajw3&d{+^w5_S;|Sy3d%1T9mA=sLB#? z(eC2OVkc&$z%VEuVnvG#@Qo0N{2f6RDr(+`9Zt3l7#=9i1p3gU(o+Jj34(-K$%ovS z$5l=k-rt#&!t>ei9^6%MVBl_L)dd3zK;XpyF;YI_2#UsAKSYAh(DNa*)jwt{#t*O& z@2s9HkE$I5ssV);z75e8nb1n@F_+S~Esomp$oyQ)+8wbSAG&+($hpxa{f*V#TH()6 zINBXwW+atJm##)OGML5G=}9=jNx(P#JZQn#>Y=!#_HT3q7^g^*c{MR-=2_D z4-7%H1$E#8;T0^39AESsHdzev^UxKYU?L_}9KqKDq^By-# zkqLpX5UXQ&MDRs<46cYRBIN|sAxH!CK?N37Y{a36d=0h&+6F_zlz-~mtGe3XmkR5U z3NS*1YTz>cHNN&xQ`<6fZR34cQ*3_EzU95DeB$}ofjeGXq`NPfSk-)tn`hM(&bfyT z^3fQ5F@QZE1#nCtUznN8LI!dJZ^*0wD{$0>q3H}q>Qy@WtIqOfl`-vSq|_8oFPbQ) z?3~bR!wpaf)?!lYN}_Jvy{@vlsun^-t6e+A*Jj=vtx?Vc4Z?IR@CT_lD)vAT^m;}H z97KRrLoj>56w8*%x_Xb;2GkqvU>fQjaO89wy5ytzs7PR<3-Zr5ZHmSeppPsKK2KVq z{gG7Z*U>-L7X#{uUACcAxPvVP4)a#c!0L>ahlvgOdwz)tP=>aeR&q9ymi!iw8uj&ZopjmNRVx znoz4+&tF%8p@o#mcNz!ppA&+Z*(gVvzWVJA|N`dBnAn;0pJ$(em z@hadSnTu97_SN}P^};rs0zh9MDuaqh{GbHrv=KaUL<3AB&p_~(7qno2OkUhmQj^v1 zqeJ0U53}nUOz^q4e)fp?@Ct|GY^U^moifobHQ0I6nHKG`cQ!Z5C;aaE798#$0S01{ z?@aDaF?a4z7Armw_zXn=?lc)4VmavAjPy~cOYqlk$8azVI zA?Nk@G2)B{*B+Od;)zp$=$0~9luKJ6nhtm%nt3fnM_U-QO&8Pqj}{Ig@7ey5-9y2^ z_xPmt;%;cBGUN6`KJkM(2#zrz9YEkk;@J1iRvx1U?jw^x2tr(Vll+tqM5CUdj+}uj zf!3hv0SZjd3U?)}2I_+$ef+I8hOJ>ka9o+BaS{1Y(J!1+0JPLA6rc@&7xKciMjar3 z1Oe4~;!J%oFsVoTqUQE{t-DwbE;;jdzlgb3A#pK11Jd3+A74c4=sBi$o8%K+8M8v~unjp$q0ufmM{+ zV=v5OJ}Q?g&*x&zs)(AvL6q5WBMPqQJoeK-<%(b7{!@2wO&TJSUE?(az z*6$f41}2;r({en8Yv~}6c;&ZJ-pbTVW6(Wemh;SiIUlULG~gx2a?CnZTiTohiHgSZ z{#N!jus*U}(KPgkIA{2svnd?=herU#K|!y-uw+j`2B4NV;7MRAFc-|OAYp?(P(FaC zpgoGvpdO1+IkdbuPFc~zURX_g0l_3b+sGqch+rIT&iS$y&TWa~bDygNf^j{Rl3w4H z4~gLxh1)kP3sb9kqFB2Mh*+?Tf0+J}#wgO<(H1E{`X@HPP(R!xTOfqYbE& z;!H#3(*1Bl8tWF-O_cTm)+X^*&}WB8G2;`w@&wNr0UMT`+$ap|q%n&DSMIQ^5rspT zgT-Sf+@-|42y-dwoV$9rB%lL0KXc|cShz^Mxj{6$*q^H>hNw;PAS5c??S2&`)f`|*_z(O!*V&CvpxG9J;hH+$d$V6jQD`%GN=I|qNMz_|@ zuCVaanG!o=4h(lMvu|k1k)D4Zqa&XvtM>^A7jr{Mix1WfguA1>x7SI*K4(^1-`EuA6kvh@**BH?vX)C8VXk^HP`z#~< zED}RM`hh4EL|17o*5%yb(265O@2(XUPNo~!bL9mtTr5Rm!ox^YiXr;IYjI^KHtwfW^tLKScI|A)reJ&-iifw zf}zWsP$-vS3MEgIaU=tfI0+BEOH`d1vOC_zG$1hdX0>}a`zPMIpLJq}`Du@x(Ruet z=9;LX8>-kit|VW-f{;s0hY!QT`cNHDbPZ>Ae$Ziy^4WI@Uz-gJnpWp7SJ2NlTOexH zwoyJjv}Y#qz3@C?6uL!OZI#00lES`d`zi&;G4keb^rRk8PFn)N;#AtCUVn%Sr=of0 z2>d~7w2xM56BR?3i|Nik?`nkQ%M3o2Fh)>_1y14q?TnazR3gE1CZY26=>Y{AizW>L)wt{ zm_m$tMYK3+akuVJXRCdEUsNn}>-wE>9cwhLy#GP)&5C))K1?;<>=IlJ03PGoM`TWS z%6;Y2M?i;zqbDk-z9=}&rX0V%k*DB5Iv!o$syu1ipQn^A6n!m<42a~;93~3B2K$L+ zCp|?a@7p}8v7*duQO((iPvO89)HzH@LZ0`aP0}tz^%iHvjaG47mMK*Fqi|_aOc~lQ zg^?5GQ+b72-pn<5ov6GyQ#racR|!d9%bXY2f*{>M14F%@63l>Pt@2;a?NRYxb1_h@ z3NVYHCj`)Y@P9-OMBD8Y9r;H(P{ zP2grA9C3bC(uiBAm{x_hfZcbx2MN1p*M+8ian4P^LdlSYPwn2qdioKdn^^L;9AwH@ z)2nAB!E7KaG9iQu^;SAuri1Qm)~PHduEz_puxMkHMEH;^>fd1dC> z2yod8#V_h=RGlpEu>!V{eTa(j)5Ny=Y4AHlCHm!hb#$0)QR8Xm`RKm4N1ezVQ2tQV zw4UqAj4$WNft_`X8?H8!t-}P)7+2yl~4#y<)7nt6Z2vKIQ#PmP=I)_cs=o@ z7-v{Sn7DV#o2-DEa6|%I#J~x-1F9g)ea?tRi{f~+4Rv(L5P=6-0&XzEv@t5?INyze z7!Fh-h{#tE04 z2j^7fbYweLAZ!(BFi7}Urx*?@P#Oll5%^K;M}k4N3)0O&8*Tx?BT106#GGV6A@Y6~ zD2lJp_)z8;1O{2*P&MKVOd4pT-s%fnM7EAgc}%DKvHX!c&6_-zF%^|X$Vd*#Gi?VR zqY#3=iR4#s2@*A;5xfE1j1MqU?njt`G#+y-L>K{cq%T;$y#1TY(v>!!7gK7O`RMs$zcWp`c4J?B|@sqfJ4Z>m&JDr%Eg{n%@wS)?yLft+9} z6&wg+gqV@rpf!y*5Uico3b(x|FKh>~D&7u%%_#&zY@kPf#BLN|S^S$kWBJ%5BH`A) zj>`0gPFnROH{u0mP(jLJXMiU)eOHec%3_o_A7LRDHd@MZ#1$Tr1*6skssSN^YT%#_ zZhm|7{T&2`*LiF_vx{rmKp-d_9tZ{sS_agY$Bbpcj1Q{Ft4*N^RN_&A0@Q&B099cF zWb33#bRCcOFC@Y&LfB~NdwN)X`N?^4%u0(Uw*uTwTrJY3eu(|KDShpS`^_$5(=6%z zI%7`7nr!5yJ2+_cfVm6JEuD(gKBmAKWu z%I7uF%4?5wo?HSQ$F+$#G=Z`x6oBbDfEAQQgey-0xS$7ug#80thyZb60>W3j&4>ZW zoVJ`W&}sa4q5z;0V6D6+12}=)29smBEw7N@5FY+PrvSM(#{7wrqZA>PzzqwG7)wZk zXCU*Jx36$8K`IEog^ZEXAo`n>UF6-E16)NO39blD0u!MxN_)keg>O)o0KTELMUJG; zCk`(=yuNthJGa9lYi7<4h^|`8Yt}ip$gf+AC4SquK_4D4?g*)aWB}SB?YlGILh{0B zu9q+$nGO@irSELsZ{zNX(o-uk5sd7rb`sD@wPF|!sO=RZ`fLVddj3k=$VToFbQIzbrX9-o5CiXt`l zC@=2YF#QMVz+ng6WP*GGegpV{l>umglK3Ck5&j&GqV~z*EC?sS8szE-Gnjmq?+al6 z17c4eG8M;n^4=!^RmtHlE=2>>kRSq5aEU&y%^GiWja7;8gD{A)N1L02^M{9Wd0bXr zHo;*6#C+cM0CswG1Dy{=iIyjODTlJ+aFT;hX@N}vo3dhdF!rE6L9T<^8U~Nn^!i;4 zuAPjjxV(_9rgMl2hy>VQW^6xQ*K9}KwljgcpZmn>N?+`QlTjjda#^iRf*Hoi8?aSG zdT>&h$ciEnj6Y!7l*tIgR+;RKh>mapQ~;3b!BbEjfl5(&MO6Tf4Aex#M}HfYJyH#_ zMUW4aa#6b>ONCH?QN)>uzJRbn)QCXv2wAq(ncXtgSs&pttvXx-6`>m-RxLTECr`oE zJHN=8F?(m)mx@nbZE;xhB>t+gZ`3FweL)abfRRWvp@EU$h3E>by;WWL^7#qB{|Gfn z9M(W_4{uh=Xjoeqlp3mJ?DXPESOmqheXHOkyeR+pjHvLWgqSm|C*BH@AD7g&8)FS< zGi^oZSi(iOqwNA+rJEr-aW`7$^S-F4H#|j6rMTV;G0L0V^F{_x9sxih5rT*8G5(#~ zthmaIvS{Y4YQJTMLlYiI(?VYEc*MEwugMMfC5*C^tXh-s3)h zKD-P+!z2jQ2TqK;M-WV}J%0yN;6Fzd04+s10w-2{_Z9ov#pJy{07!w$jKh0R5U#aP z@O#%-Avl!`izky7bw0Nik-Jxm3F}L6!EXO_3oZh}hLArYvqPohil07mU%c{+a}4|r zC}n@Ma{L$_>f~U?P(OTR-D>&}MFz_5WzlqonKY)!Y|@oI?|+)s$p0BVwv}zjU{*VUOOT zqxi(|1zG_pdsy$el3V7zf(s#_n2Z;JN@mbzjz6$4U@p-v%0z^+xB~_R)e3oDmVqEm zscPQ+BYlUht6=Y4lhi@j0K_`Z2ex82S+L#xFpu3oR2&%3cVf2sWyNX&2|ju z5XqHSkEkyjVE5n9E}A*+ky9;ryx8rvPW1G9r!(8WK-cT$Bfy0IS`*GfD3@isfXoqU z;g|?Z$Gc3z?F?c__HN;}L5J3%nEN^b6(52^G4)e*Q7JA#z1P-6)50)gqScD)&h=_! zzl->s+-!nlzonhrFaLHmW0HeGzQCa`M&-S|4I)aⅆOz+Ls8Z=_1PLz6z_KzK9Ad zSggv1#mjggV$Z!-ICPF-65KC|P9Oq4C^)Qa)NOZhFjy@MJM>f1X|u#^QFl009qwqkk( zt~tLteoMI)6r^0Pc1Tpc@liZ7x+(lK-mw;ZzPu16KYiz_&4#K=MELizoX5y?+eFjN zaOu7xN@hqChj+gR7^Sx!o3?6+%~;#$2`2{M&8%1?Uy>JAF7JYD^`&X3&a0Oks+Wu4 z=5C|ymBfm=ZIt&NgZNINl&y0y&!%Kwo5;-xRHiT=2bstYm9n;ZXE9~d6*0+OpEby4 zF_T#YEut$JQ@Y(2$h`O$9>x9 z$Sn?OA%vxdFD5oVvRHXFr2`uMV&9q3Tn9t-7I(B_{zKlV$uxBRWxn)TmY_^Erz)0l znHjeyex}pV-Og2F_8eQ;py=h)&JQn~ESFa3eEIc-@)c+$I?1sTE{G6eVm61ECvP$1 z={qu8(V^K6@qI(E$nrAd5*14Q2oQ&lY9{rCYONl$m9qmbkiP=({1_IQ&;NIWf3vPq8ZR^$haj66RX^D?&NJyOcedU+e<(#W8A97NQ(( zPpfgZGWT#mddX4|1j&g--F0P!cECw~r6suqSB(_KIKs^>Etou3z5Cs{F-HGQ!!vz`aRX8d&>t=hj zq#;Rt&BMRg%Ncw2wS7MAVDCjJ>B8klPKl+QyY z=o{98{;{kxpyp@uKVH`L;;Mt60w2Tsdp+-Jc3<-W1>&|sbw)fuu956 zw+eSCDd(OA@g)fG5!$#oSEjW%^(qcfN+u(UHCW1#4a^G1w&S{v6>U%birbw7Gl!xr z4nT*+>0h2~Qg6uIw}}If*4Xp4)zj2YFF#ln&#JT`Z_e)@LbfDNoLu2C)l==*47rJH zCG)LTCHp0b?RD*$w1KUIIG*^vNi=;&xIL&b%E-#MxifiMg2b*~M~oN&&;w!s!eOpJ zAxvp}8t*KAGB-jT7UdoHb0=!){C(U4K)Q|2Ad&%8UcqC2eV@WjcDMZ+TocQE2(n<^ zpw*LgBU?s_gPZ+D%544ZHmI0@!;*zy3vdE(t>~zdqLJAp&WhC{Ehq!I4)jj`CpzzP zVKxsgf`S1h`z#xIidQV0G+#{JlaM$5GHbS_aOtxY9TM8CK4b^ivo-Csbp4X+Nhy}~ zcLrp;4|uuz>V(~U-<$4!_2A84>VaDIVpNsMtFJE954SVItEb}8anl_36jfaMb6}=? zR5SqWNE8(chilal0|H{E3b+o3P=j` zrZK50hxR6D+YT_ZtTCzKjX$oYuX4Z8{>83-%MaAISyWq9z&Q2pbe4{Klm_Of#%lUF za$SxiJIMDAzzs%D+Q@3kU+R!=2hN+)M%@8#(?-N725QP{9Nk_CT8jgUf=uIBNsImwl0hkIrhqieMnO;qX{iF zq0*h-``34R=HGYa{OfCa*A9EtG*mmTMz@6{tIjkl8mhgVUMQYF>&@gRB3on#8pqvPMb{vQk<2m)vORzWH&Zk*^T7JT~#utb+wwjc6d zPYAGzeGn)Lx1Y(cM8Mb$+8M}x$%1f;3g*cn0g?@-PUR~qLHIc@Z93oq^J4d)m;6Xj z5^_Enwcu$m4{|A8iFLN(`kghR@aj-uz0OnD*7uT_HLyJAcoc+%6=kJp&`9Znipa~% zxFnm7^AY<+WO7f2S|n*uCdfa0AB9g8H&H0XDMEFUUVXJ3NOLfQ*y%%{F6t5cMW7}qNnH)6F94ADrZd7}LJGZF+>=a>#m5vJLU3B0oAC#4&)1dh?W?o$TcWF7I+x#OgA32C zrak1>h%ITcS&Rt$$Ya!a>l&IA6sdK=5Z{3r_JqL_@9og?>uXUS^Clo(7e5-vv`IjxH&JAegt{8sc z4C=d4$i^rN!~EsT6!~4hYg?}6$zQU4;a2f?zRw#JAJu8~)+Q3R`f_Le*lON=wum3F ze9F=gGJt#t$2VpT*;YogD_A|kH70ly^=2b&QEaK|(T zsFmkuF&x2LM`EEGkJ`RO%T5}sst7a9n0d!<9$6C7IqmGpfN#HFFXbQIJZz4ph4s)w zk!#F4|J)_twUYXy`W(#E;R$&(E$_Zy#3%q0bR_3bSU0fnFQ1=??wdDarbA3wRZ*0= zQB4h+fN}Z876cge%sFjDzJ?VR50k=K8**oP@)PU|RR}Kc!9^g2BF$Fu)sL8k#wZMq z+sb>#wekf%&?r;{?4eN25d`xMa2fIsxB^C;F(C~+M>P;Pc#WS;RRKILi6o7IyhHd0%VMiUqh$I!s4+L-*a9! z;m?H=z}8>Wc=Vqn-}VD2f$~O`H+f4X#)?Q_k~(ZH{!{Bv>{N< zLvgv(`DJ|3IBiNf89hyzc_*DOjG_&bKBSz>Y5Xj0TSqj!5jQ0D(ot8JGfkXFdl*#i zU14$fEzicW8J?!kI*e>+k!|?Y;K#ef_?#-}ig5cfPVzcgF5+>q%Vee|P-k3IE<7 zn3a?QW7h5+Hcy-N*C1pmj@B}zt)aPWbWXnVW=A0sI;Ovej<9y_@>>S%yWb9{eJD`i_dHu z_x3>E`cc7WKOXM)4Xo?jHfS%XE9~vFY#f)>BFSVPCtyMgwQD3>h>s5DcFrD(KVsJJAaz479uzmybcnUVLx#V}I>eYe_vkp?1ci zTQ6-7KRKiCo#Ey#se0;Q->YK|muttWiie+n`ut-R6R&h=zbg)p=yAgsoekGEH3?8} z?T&@KV$US&qo*vL_3N?~b6sW;5?XzLG>Boras>aSM=i51?@T?x?xxih+8Yl5{fL|* zHa5X#-(X~88-ar%N!;ybS!)ut`7Vi72@(|4>ja>@f9FTa-tGbv zY+IY;h-XIFB;y?uvTGZjobujaAC#$r0Go zcCSM`Gi(z!(ubHDynUnBsA=FqA|VHm6;f9a2O<&ts1Ob@(?m2 zCUAsJDp(y@gE`jQa74FCGeg|EC0lS^bG^fLpYq_O%440K{TYZgL=PR@F0K>TE<}60y z-KGlN9D8mSG$jp0%A0UYsAakp$uTe{wnz}_z}yJ95K-U@5mz%@rRO!6)z^ravKHhB zQUm!IP>5BRP{CSj65&`5#4I=<@4&QOn9XYDDnUJB6jFsGN-S=Bk055WwuuF=w6tcP zi=Fm8n-yPm{*}hR>2GdsFZg!!w(5I6$Qqhh`@xp2e}8Cn)1);oo@zgMqU7=5a~t@!c0TQuPn2bJ znpm)_Hen*N1|GAlF;Ik7q7va3E2Gu_0IC!khwZ@w%}6@f*b|m>AzaZpVF$vng?o=m ze4-a{k7N@jdXe;kBaz1<_D2|j^b0Hs)bg-k55wMlz}~v9NuS@^@{g%o#c_$%%UcSX zz42EU?d#HC=m~tjV@g%Bqs5nwy%O!^V@b4cu(N7Trw6FTYRBmS%7e79qTMj0>o> zrpqwrBt2nw_KW(||NM>&E$nrh?D@OZEl5IK9ZCexD)V{2E zd+f76oF3cw@%~3%syeQ!m3x7WkR(}jNFKlF_PS+R_L2O{ey`{d#G5Jq>Y{{d#6w^DyZ`k4|0q13 z(u9hTAxK_9^30S=XVQTz${wf57|v8uOGr-k?9POmG_-^>qBo`tp)L~Dk~0M@sx_%b zB(4kLW+A-fi(gLHQ_*PA7hjz9ny%6Lw6G9dA|zbNNX)fZO`kVn&|!3uA_oGJJVS`7 z&FQyblNmxni>dS*yD*EYATVwgXu{v1)YL!;Z(qK4Y2h^|t<*X#Au)ihNXnzX1(kzHP&>~ja4=10kRDiPA?9THsWpDNXT;=_}MUnQ!SjgcjX?CkZB@=#O}zb3~G%1~pNy(I+^hv(TGDLAQ~R zR3uS|rwk%f4Uw}^k?Bu$b*UtMR}+eBb#-*wfY3n5rvtapphg=PI2%;OLziv{YKBT5 zd+8?ia6DF%)U>c733_899-=ezHG=~O4jede;J|?c2M!!KaNxj!0|yQqIB?*=fddB) Z95`^`z<~n?4jgV0{0l^7P67Zx0RS0g#c2Ql literal 0 HcmV?d00001 diff --git a/raster/r.sim/r.sim.sediment/testsuite/data/reference_erdep.pack b/raster/r.sim/r.sim.sediment/testsuite/data/reference_erdep.pack new file mode 100644 index 0000000000000000000000000000000000000000..2e2be5e2ba05836a066ddf62defd84656266d992 GIT binary patch literal 37150 zcmV)cK&ZbTiwFoPXp(0F|8iw!WpZV1V`X1ua%5$2E^uLEYXI!M30zHY+xOi(X-1Nv zfedLL_FCuKJ5h!rWXznRGBl}32#rD!l_XK3M5w5B9Bb_qGDJjWo(h>|$UM}OfBSmx z>$xuXeZ8OOeXjrgdGCKmpZczIjeGB7AII-J*E-h9+``;qn49McXK&91-oaA;7PIWI z|BRNFHn#sU{cFCJrL~o{l)(zAKf#uTgn9OB)rbPWvBk-oXpJgDn2N$M_Z06lM2|R+OFP-}C>8c`IAW*49>FYiajK{0nSt zqzvr-%=q^V2=M%`*85k>|9@)co%z2r|99s9zj*$4n>1mJ$4IxyqyKI3FA&)N=lS2t z<{$FEwN*O|EdR{-|BL2-@1PK0#=yqLj_wTizoK>Of2aO;>VK#Hx9h*l_)!!7TlC+? z%J!f2->Osp{|pQ8^!H{ATs{50yn{mxoB{)U0)zcM{eDd?ny(TR9JrV<@S8thWr1ht zGJnP(z;gj@qvFXJ2y7_Y%7!wxRPhdB3@pumd6?MR(3Y0g)}(#7Pg5#>o}nIA?R^V@ zxs~AW?T-c4?F($pZT{YD|5#x8t2bg!*{S#i2K;K6QAB0MuL*%A5m;K8TU)Dy|7u#Q z1O$34@b+mxdi&o+0p1Hb^Hrz*cj|wq{&(8{nf^~5?=ofbf6Msa_Mi0M#=3L-|EJip zfJLEZh>HIoU<|A}^FXKmck2INO8*6&`u~^G|M{MwA^*N({A>NU z`v?7RAO2IF`u}HGKLbmHcJF)p1_m$k4pC82`8AX`_}~9@CW21=@6`WJ{qMB@vHI`r z=l9>J|Nm6~)!wq{)c?->-Vx#`QMJT!%EPZ|Nnyd-^b&R0+7%@wf`@$`-l8*Yh%-S z{qIk)(C{E{#=yrf&@)uUD?HTOBe4Ao5zH8@SP-n@w=&FQexQF)uy;tvqQC$TuRksi znsr_f?9~5G{qG$Ab=v<({a?DwGa&T;di?*F`fq3bkNR)bzIvzr{|RPhWh1b&wY9br zF#j?2yP2KYpi}=l^}kd9JMDkC{s(&o_NT zo%-LY|DF18Zb$4o&;S2L^*>;lpWl4{py2=8@&8}wzpbT!`ls{%lwGI({~4y#0olP4 z+DJG|!f_JLl+af~$r&QqNC`Jfm?7Z-2}>kAC1H((cO`r-;inGD$x5gxVIK)iB%~x9 zA>l*`=SsLlLWy?B#Y(tM!aWieNQe@im9S321_@tF_^pHTiV}8{&`?5i3560mOE^_R zF9`!BTrJ@S3Aam_BVnP0WfESH@P>qsC447gYX=p&NT@5Jv4qwVI!Nd$;S33VB$V6+ zps-HDLUae*>Ub1|96i6|A&tMJJ0|Bb;tkh(eQu$R^Weq{@>;w#{V_~LFf7ZKgArC zUa@U$Ig^}DcURL0l9G~%G42eCPW|uH|G%jIcijI=5j#ukPW}H&>Hoq7|E^>F%jbVu z{?qfnEUg5c`u}HG#}hs~j1ZNTfdSrsYyNt|W^m|2Nz>Zew%t`KA%V;OZkEljS>B!@ zp?}+^y#x1aR+#tS&l3EaH9ycVY zDIM#ZS*iT(t_1G|jDcNep6b;9PW|uH|4#cq*Z)Aj;D6sS{+0eyme&9D{4eXy`u{)0 zj1B&y*x%qkiv2Ckt^W4ES!_wWKLl3RjAi?ewVmC6ti@RV@w70~VbrhmVQoiSTeYue zMXY}P8vGagQ6c7%4Q<;~4^7!w3B1iL?HDWSkHvo+>R;?$rQ-$t>m8|>{U_^4mNS+H z9k0pk_ecD-ij_6hzS_Us(XU2Y}<#H|ET{0yUyo-{wem4^?&vfO6vUNzI0GQQm3Y**Fj}T zomCf62URb1P+d}or}3tPTI)L4RjY&D_I6N5QfJtsvV%P(b#i)7I@nuM=hav8`t;K~ zXeg-zG~zqhUsA_re6NE813EaUt%D{>9W<5H>6sOG&_YsYW?3&ugOWR5Y$SDvcI!Jx zN$Tk6{0=fq2Za|qXg|AyLnU>B!`5`rQBvpTl-0qJRvjEAse^QuJZRrVQpY!TSqH~S z>d?lgb9o|6?NgbT${0@43 z=->iLourRW2Yn@Vev9Z1F0SgJpQH}d|78aQC3TET)jAj~slyAEJVdLa8E zo?~xV{eb8fQV>ddg0YSTWi6cq?`4(g*nP9{^hPbOn-+(h?? z!}LD)Ab2%U3kuHGv)jJZpozQ+o(k|H!TZK>(iM3Gz29KXs?FT&{ZARqb=`2E(PW0% zWJo?$I+OT58_Chln`p4QL3cTD59fx?g9s}Rm?26a0bgGWFZKxHQ$k$$G4GAp{Ot!= zV|iEBMS;*?`|ZWtt_z{scPzS9ub?j;o(ZEgx>6pkmw+r;ibjb%U7jKXBYqA5muS6G zTg?sNDz%G!V_pt}cP^#kf=u9CZ7XU{t%r42SHXDM5#oLYvv`qFP?^FPO;MlvJ{iPW>A^!v&LwAyut z6h?TFH9q4@wj?)`9@C2L_Jn++uWjo=YQ@(XIn`&3)W;#*anKP_+k=H$pCNbmkq-NH z1_fd9CNSXeahh{(Mdu6aNc4zc+!*_gcARYi>wY-n(sxJM`@_#;YI+S-bTSAop4Wkt zpO4^4VjGpJL}ALIQJ}QM4(17mQd7cZFkNdGyHnV}o36hgI;CPG^6d^H>tkodUF7GG zZPBf8@TxlbsB45f-c%CBX%*NEj`X$qSPZB=&EDRt0bUC&(0*nHQ5~0n3i@d!Z==4z ziam<3`tD*F;A_B+-4_P>0p6g##sqR_eL(xjfxmSKs`lIm!!v83tHvUM@Ln-2RGY)L z?u}w!fAwPzH55ShiL2}t?K!A(E!!?DaK-VM%l^1&{3A#WjE9?icMxy<;reAhNK8)w}MYQd%uu1fT>p_@O@dEZJp2j054NyDw2sS<0V;6Z`BHCcW7}_wcn4d+c@*)p>srzYa|lt*UV1SIwWyXSck< z(1fpSFx7~w)`ybYFW+F^&&$LNgS^m4f77_pxl#7#A@ne z+!UTe&vZKoZMP@VcU|Vg8`~#nJ|P0fyDf&O9!jWY)P)|r&kYmKZUc+-6sSCMt#sN? zFF4jP2O>p7z$B!o^zgfNa82Am_jBrt3LAZX>kx>pG%J`&>#g925d*Js0(ejL5h8F z*p1A2{Ft^`kyt z{q)6bp?V0dSvLpQPS}YveOc^gvjf?s$q*VIgR$3kz;clVyWH~yxEOz?kB@s@8am#E z9jIGP_3>3e(-{xIVOJLPo4Q{RAM_C&KPL;yzR$$1yLJ;p3l|)JJB4OW7(#_Ok!|RI zhJ9MGoc;1`Gp*m$ANxC=q@!l$;+h``zjX+^s;?PDn67%=u+MNR^$y6mm%;0< zQ%YYa6~at>!k+#X3USZe;drJhl-6~F`8{_*-@9R;-XjJQe0^ch>!H-{_`a}v@=mJy zfdTGcT1(#@sEv)zZ|QeVukmih#8S#+C9awN9X*G1a&XvnsI?Q( zyT*(r7G;4gFom=5CK#I0s_`YJLW=?7Whx}nXWOWbgO zDH1wsCn^o_1xGFs*Ix?5%q_8)xwRMEv7U_!O1ERQNk5!q*oBjsu?%8*a;O*W!p3eo z1LDky@Nrorju=5xwY@H2vXVD-YH=wJsOe7%oZr!gb6tpRN;LwSd-#Ub{ke>s|3u0mK>i*znbDq}I&4zLMyJ*XE^5FHM z!S;feI`sHn3KIo6FxUAz7;kA{s|tdl$ukgU`Hf&-UHS^Pi}#~~@pQK2=R)imS;ESV z{eZ;xD!sVb9v25>)AC=QV^zZmNGN=ZJN=UB`!9L8r^HbOqvP!Sy*{|?GwgUFMJDkcZzKHj9>On3p0}Tf^f$9CR zI5FocuA^3=Kr@Q#**1fk7W9Zwt`X4RE-8?}9gk7?IFJpm?uTB(d%%I?8(4dfA+X2g z4y=o2;ZtoPdvAXW9HL|3@n#--&d7vXx}0R*HpkiD))D1fgE2XUW;13qLH8y#7-HoJ zLsl~Is&E$@YxD>nULn6vS*}|qfOHE_Z8)X@k48noyWR|pU!8+aF*0;yv?e^RbB3FH z7#3$`v2XKEf!sA2FgH2C9&oS`*oJ(8K!HDNIsE`UY+_;iJ3kz`~h8uja{law>^iCeoF?f*WyE+` zYibCW!(YOYZ5wb4^oJVXr#SfR5IR|K1|A&@q(!<@AVY5uEa-9_j75`5?)v^HTC*;K zooBrkE}d)oy+igjazIb6LA_*_dibYw6pSfgiNi%SEn5x>MK@@@E)OBIb`E7z z2XM{Kv~*dq3tR2FoV_bIngy+EpgeRcD(xw;D^(hU%aXn5Ax0W-t7jnw?EgUD-I$NF z3rk6Y=q65IzX83f&f%P^ub@(Y1YFaRX6(LC5UMtg<6nH}CVW;?%U`Xv7Dir7B`;6< zfNR=TveWf2YAkgsdG~fUdrxaMENfQ9tEUDqUk++=C4KJlvlD5NX24l4^TTG~brebT zaU?Ho+(_i8M`Y@oW;AoSMdZs@mvXMsn18=7J7(hp6bs5}zikOPcI*<&>0-}y`RFe+ z;%d1TM_=L08=ts?w}R=sa7&c#8HPIUXQ8}d96GmtWM7|KjViD1L6BDi+`Ev^a^AOr z<@Is1cXw3X^_sqF;7sd9CoxLgBo|CKGc=!Dn0;%ld7oY9!?v6mSNEaduPgmT9YE|cF}=PCm!I6LoeY_kGq&)9DqwM&xcLJ#rRbDE}7V}0Y~ir z$Vutv6IIh{E;1{K(c?#R-S$`umpH%XQ*RmxE02xlH?LU5l;k&&Pf0$sZgU{qXc~(7 z=Z2#G@d6Tf*8vx zL|prvWAyix5&JcZNbZcIbm@^YsQns((j*r(S4Cs-eKj`wCV*^k2#nf#$Cu> zB<%m$ol7@;L5*Cd1*L`$;lUy?4xAOn*1tWE`XQ|lp)v$DT!LWj{gLQ8^A2w4(Z(4m zO~u4BLqN++iZ(h>gK_Bs%Eq|_%1_><9{6iR#gj-dlzR^g-;AXlk7%?HBa&$guVue? zNLwLIFL3;X1Ijh2YHKxYxnhS(R&j9Vo&u)JuLpxY9@NM97P$3FjgA}0pt}5g>e^y8 zSd>=7HWqU*{AnVT*o1>hzvnpTXJ1-o{Wka<{Q?(s6OhRz4Oq|JB4>=|Vp4)H+2y5( zZOOaHf(aKmL-Q8lfN(3mPtH~0;`Qfw=NZmStj=aQN&A3TOCI`GYD1SJL2xyjMtN0v zNO)I9tC!6qgNMZ7pfp+f)PnK+gz$y@DqdCO*mEnNoBde0I~}>0S!2nF(~IbBXR--B z(2V5in~~iwtH?_7i9R3Hg66ST;Mf~Q`p&u_IJKrb^qz2?t6_DSJ(&fhR9A}JbSTF~ z<|#NQ?Knnz6tnpYP0&S{3nNX3l|25oAChj^3)~H-Ver;k*2!Iu%`w>k8+7%^g4;@r zfmR;c87-%yHgGTw&%#5&Rd&J8^I-AH5hjQZLywnNaPUM^SlL4nhC7t~Ueqgs;&E_>v% z$?fd@+%^zQ``$is?*dtysv!0G6S%m$N~C7dpHgZjpm zaCjVn(v}58UrrB~`HiH9y3OM9P3?H`{3?sWo&^NQ~FfyPV zls66nli>4sdxbkaKpaV@MZTg+#wekAMk>^wh^KEVE`>$=XTpHFZtRzfv*5s(5MZmy z!Pjvy&?fUCMOGf39n)nmZnVZhJ>Gy?;%Ho=X$`YStj6+^WQwg*$7lD(!G(i{T;@s? zrsRF*Yg-Qr{p%0G(t{!N9bq7W!*A^u8Xd zm};61MZ&E>y=PH*ULwfkX$lS%wQI^NMQ)G(RpwUFQu1(EPiDjNeD=ZinRslT2T)7f zfWNmM<=1FJ?6H|_owYPbofd-svxBI1ycc`wM>)G|td*3+XaOW9qJ{SQ8b%qFUm?o=7N{_MHGQClD z@^|q`ZWaqCWHEw*piD ztQ$w1zN7N%Buv;@&CXb?kE@oXp)t1|SL@!P`weaY&%>FhRJ9Dkj9u`aa{@7yEyGQI zZ|JNxV@mVv9OCurrr@FA40q^+P-wj-o!X#H?u(||qvg*)GbuezS$qm{!-_J2Ez8M%W&#>5ycx{U;@ z_O#=9!4&20RRFuJJ>nA%(M-rmvg=y}9_^#VEPA+ve${UylGo z)3gaVLUSn`ZI}j6YbjJ^N@3m5=g^k9oH}iz`TGd**vew?yh|{iVNfpZHBmXa0UzYd z!5cEun0{F`w8IY(sylut6;C(?6>GktW&B5KQVSm(Mb0)%^~`?3Q%6F2((1kQOQ?OnrFTQ_btn(Q+>nn=C?!C#QLQI zvsE!`FE2zPn~a= zedZv!G3W$6(WXD~o2Le)8*JH}(76zmum`@1GT9vu-C$R?GHkUv4NW!$Fi_kRbxv8( zqk5mhnEEQzB1bTOOgbKQ+(5j#E{0vxe=<`InmCIC*+SPIf&7l<>q2H+66&_eqmx%} zIK1d0oE7(DKU+A$uJIEwTW%?AvD}Fct1r{t=I5ZENhF=>9f|v2_n_BLdj!vByrEwQ z>2n4jUXrB>JJHSNJE}!i<448yd|aA_@cGpyKK5IlaLtn`e0rS{zp%Xn5Lao+r)3}F zREGDY56n~J?u^hCPCXgS^&4x*Tp72TTYT^Ww>mWrRSR0dxT^-}> z@;UqNW-+VRx{B)|XAVakS@Ojvo!LqyVs7=+_5tr2nB2Gl5;v@%Wp_V>`rz)=0ZRjL zReKG+rrd&!3<53-Y$_fa?Lh>5|R%5aq85uLdW<%?o4U*@zLSpfL$%A6Q}RvR`|RUs<_YRBFBn^1E-0&FK`so%!+YbNe)JSYfo9)e(=yE50=jl>d8|*x_t&`+mE*I zSA~IZaRjT-?2pDEk2;1nn0~dCt)I5(_ZsplwF`-BvBYe>GJF%!hs#fjB}I7}Ot-Wt zjFQVDe810-EL7FSoH7qVQK3IHo{UCi#U1qKI9U>zd56AxWC#`vlfn7RuA^OR3Vd9g z1jjbX17|aVS{`}`2009cZfEzB347jvThTC~!LnJr+Udp2w(YxNzsnJr&6}~u7Ji0J z!@5A2*ADFdZ4!zsns6>&11hx!$5sYFr6NyT`9Fi(ieKoUz?Ud0{e(J6VrqoJUQC}X zMa_al>>2;j?&+>#I5~SO+|7QAUdbo9Jl_w()A>8Zt`&9svyw`YV#QL?;;0H?IX_o; z!0QLoyk{|K?KMp(KcbnJ-Y@|hVG$&1DZ}L_KD2gjCCXP1Ej_o^1@#@S!`r8<*l*NB z@%o9a!YX}F-hRe+;WFNsd-chZ>%0CkefP&em~M0#R0;zi;+{9{Ro6sayf}keB-Fr> z1sgFnXCn;Dh{fTjQn=7Pi*R$*q2Fr=bKwy^?9~>sPb|yz-Biw;(yk*3edjPY=B09T z)Vt%IQ_UDud;|3q#zUp!eh|6t0|Tjah)s}#nwGQF;?EBmihVRka@k9Tb!k$u-Hk$?G<9Q3S(ylKhQA@9SqUi1NWYK0go z0vg%7KgU6(krcdJ(~Yud3C9G*SzJvr#ZS1!iUukaholZL!Vm zr+nEm1BPGqgUfL+Vw~i|8D(<^^!hd)Zdx{?TiZsQ>{$=TD}teTxdz>>{wlxj^)P$u zL7zCTMkw?o=a?DJzF2)nj0Yda;6?U5?o-S`{WGu8<&ZtZT?}bX>w)q zL;C#)J32UY4>zs9CpRr(G~*QKz%_@RWfH{)OqE(vMEz_bfR_1=KXt->3!jN2_&ckY5hlhrP$hjY-h|{CZT$-Uy#`mXZFwb3k$H zb@)^eN7NEGVyU_xG>kg(dxzl0VKAXBoI4XP5Uz>6!MnMf5Eg%Y#(lr!4pv`^QLa3I zN^P`&^}_~0s*WGz96A9R&tJmetG(g6#sQ>gWz>v(N`}=caB*{Uxc)=?;+Ck5kYzcT z3z>S7Zjc&5BJ1K9iyP~>=)~tx*dj$g`k{qJQ;VT|O;1eGd=ZvitCZJWe?&ZN&?50RcQGHK@`axZ-Wa`o zDvncr&oGK1OdU3GigSbM+NB@4ndL0e>QWD(FItFNe-~<_?4UABykfuC-KcCx*k1bs z8C01evEaT|ZNb5{4G=eXJI)=g!LIMW1#QNyf>cflMNt9R+vFzZc4>j#XKm1SlsXyu z#EFdDd7K$r-JNVyPyM|^bkvKWY-0ic+;O-_%P3E5QP?1I){WtFG^e6`PB!+mUkk$) zwb4pti{Wd;a2Ox71IC{S1iglK$27;oX1QI!8yv?lwJgw*eZ{2r*5dm2OJoMNjDznQ zTS!Q18sQyMsg%J3nbEs_g%6cuK0K5XTK3r|l)V=$yf~p&q}e4x)H`#Pu+{#9P-%-$tXlU~oQo~u z8m5alD>Ow+&kG<0%{zpnihq&=Fqn%SFU<__J}A8V_&Mb34N^4ey<-suy4eQcxR0Lq+NP?ida-QFHd~s_-ZCf z{wOb-FOREVu7QKmV_3g{XK<;B)e!ZpJs7~N_iHJF8pv> z|Ku{5JbWH{ejLU+#&VE9{QykTse<|6BS8dtIHM#Dyb=I5m7hTI!Wu}r+zUsq(8lo^ zzNoq;4ehla;piFBC8r$@vflNs%=X1|h01aVgsy`&3HM(dDZHMaDtvj}Q6yK=BvSrz zT-0s$SCMq9jgT6Dp6}gDhL`Pqke4++!RyJN6k9nsi}mUz^4~0W@EYOwgylmy(e@`g zqGE@cBAm5YwD=e;Od9uz@h$S^q)bM(cU$c^9gkD|*{@v+0;y!P;m`PR; z=ZS-E*`w2=Yi!o(X}~V;0y;h4LfPspbgBp;8R-w%$=Wa` zs~@b=ybqA5>Hue?9{&DZ)|4AJ;Kl7cKE-OVFzMrWex>;+VRG|AeqsMBSa~9tUaXvj z-Ad}|Sp6+P%eKOIm;0q7(jLNCr)+RA3VRq*tjgx8a?f=F zaZ(cdbc7NPaxf!dwOW{VYb6mY7<0Rn?{jbS2J>DrHhj*qBEISI5dPD++hVyczG9t> z?_v!zMyy_YkZXMOka_kdM!2bOhHz-tL?H~(70JgWiw1A{L$YC?e;GBW^%ihLb>NxjPE>#7Ow`xT#Pq5FRG)H(L?~{@u{L_}amx{^6+G$Z zg)u~M_6<95RXuy65C8j=m7Tv8PYPu?cg1JSv*_nSZMigVgr^KQ#`rRPJjdgjSpq`A zO5Eal2@1q5g4ET$P`kY!eONmI`UX@0HQ_oejP!&1heIIo^cJZ1wq)N7>B_#`GL!n) z%N2Fa&(QXDMd)kw1ct|4zykST(mZAiQx&^HXf}MZFjd1rcvI6|cy-Jdk<`d9BCX0a zkzQk!NO!>m;qvA{_WI{=?(rQf9Ojw8RSyj0xn7(2#trxR>X~c#@tfXp2J8!5sN;fB z1E#R?qt9~kORLa2{5X#JJQAIU8PjJ{({WpKGVkP2&DqI%(K}gREbcAE$jZ%Sp1f;h z3eHVrBKwWu^)`Qk!enhWwo4Sm4eUX5p5Mbor!Nwfo<=xNumyJKnq%L(wb&`z}4n2MpH`9*sH9H}&f+((k{JhqP``@7NR9cQyL-cMSJRs`v8ZxUk9ZR6zdwYkM&wi!IGQ^81v;6c3I=}qQ?}5^uO`O5PJDkt9MO=ycLGI|@SYGhOlb^Ql96v+n6u)u* z3V!nkYd*Iil$VWi;B=nf#ijGUFoP~n;s=jD#aN~+!WkoHU=C}F*H_%(%;dY-rQE)X zFT->>Q)79ke)o=3^X~#}C6(wFJdb{~Xc~iMJbBG*K;>=moSdIAqqT&_s*;=dd`f@1 z`A!|FjL~C0?RLRSUkWk4YaUp8TEq2mS|m<40BTkx!(`v_P-u9GU7O)et4}*k4Ry~Z zT`wT)$q1*fmLEpjNr@0yzX^VLYmnoGoLxb5Ex9^23aP%$rTYZ8a7lOoJ@fHkc$z=I zRAc>8a{T;BlpSF4dxyBL_Trmeo{9`=pYZRm7K-$L?&PnepCU^y`C)U!db*ri3QuGo z(2;NCactcvy4%Jp@F~X%`dPeY^|DmiaobDSZ#(V4{QCr0|FsHsyv%^Md*A6TZKiOs z2}#$YL`)~ibiv2=v4w1VFY4n?ICe-$eLxmmES15E{RwQ+`!DGB zHW2&yWMK4_Iao7h1?GMn&sC1F8BEwu_xPdP&Wyj6pVEv|UR zx|V2hlW_A7D}mF)8L%;T8*Cmuhv=PHL&``7t>-EaLvC9_iLER&=WD|#!EQ(%~=^a+YIG5_3jXW*7Imw3ZW{wxD;!kHNgzEQz36(Qo($ zxappU(r(}2MdV;m?Y4=XtjHZtmotTRkp|eWfug%*J7DCS)s$dO)^XR@H&7{Bj%3xR zaEb%okRGjXPzLnDy?GQ|wl#yifeC1SVHDcltH<6Oc{H990zdVa;>HOZ@XCzo@ad2M zHKe!TBK9mvFOcR`EE>3UA8AsaK4q2*(z%h%-HDGu5}q^9JLaL*?; z3R9L3#w{)VI4#v|_@J9dVjG>v2;)9w(R8&tJXz&cN?k63H5o72>!Nk+)68hPXrMf7TGoS3==qX9{YDIp#Q+6X z@$}0zvEZ5H0d6((QE_`N*bWax7573|+_(dt8p=bC!%PfS@JHE#pJ;OR5@s|!L)|_P z@W2GRiLY0toen(<`iV4>#uGq}d#AsG{82dXV@WWbFZ+$Q;2v16s7h-=#>-t&B= zC?H%;+?$>b8wLfU(2kIOy`K?lA6L3C$c$9V+i+4ww$$uw1;NNlOEistgp-7!oXos5 z=u?tG7`l~yysrc^LvMs zJw7S?6!DYan@}L?HtPd_*T zp26r%Ga&kdE*d>J0_RqSfyaXF0&|(Uu<*bI$kf;iS5AB@&3SHx5l^e=T`RNUT;O;z zwfQS2*F%p}#Pu*!E{DE!YdiH})El&ze1=ZW9fiWOM!33u7oB>VVAlse_&ST@?oRS& zK4&=bnyYp2RLESWq_<+}bL9=#ckOgC%%YKV&Cn$8ruL%Vdc^QQ_v{s!*m{ccr}ZOBcbp0N ztwR`_bDM>%N-s2W3Sx3!l<)(idHP6_a@HX9QlQ!0z8R2m zW(iwt`o~#JQ+ZzbSih_a)HIZtt}`y$AN>8b0^@U zwkxQyNEU1EJ5Y0W+OjRh@40H<8j<7j zbNn6qYa+SwbNsc(fw0%qjZJmaCpO)*NV9e>PN)`>37Gc!0X0sRmKrhpAb;{h*@1h;4bI0H4~HBi&@eUVGXT4rly6LWHM! zm-34e_6zxgdi?nBzPP+sI6XYU8LPKOvy0*eveLTNDCa#Ljoxm6n$xQA{%`}#$d!VO z59`=(oDxjFo6RPUnT3wq4r8pHI45{8LM>MBks zZBoWP3VlTKzQl3UEI;G0)TxlYRfr>}2NUNJpD|^ct>DbcAV^I737nlhC|=8hf-RdM z?p!&0IeayIbh=Y|a&|7{&^t;^=ZkI09xZ<2haXoj0&h#%E@llyC=i1){f zRPBr)m>y<^8ctGd;lUrMbv2f1yC%la@2qinW8h`DxO^Hq^wT3d*RJ4JYZZv+R&Em&rR^5S^8-YI zWmEV|UQc0`hb@NmwIsJH<>~vmrot|(oY5!xEHAZj3v<%Am1`*o!|w$H;8tIEU?URI z)v`mqS-jXeeYuS z)pHRfwB@7p4jal`b`K^L3?mxp`|uq%o#^!o!vu>03>Qqp14Hs*Z`^fEzrP$ewE7}- zs(|#HzLU|IWlisv(ZXe28ga*^5Ik7rP9DtEr|us4!9{Ak;e+)2ghTI{3e~^7=O(Y6 zhI`Km*i|!caQ!0OnNxNyT(9MFLem3x_&(eQPCCn+lh2jLwXNN8Vf+G^-%8_r!$%l$ zE{k3L%?nbjGU?E5(^2Nw7_POz&=~^6i)cWvfDOkW1sIk;BvQ~aN`oZ0GA=HNx(@F4$9cbj<1jfPTP*gh|RSze#+kJ;Z(NtYjimQcni>I^4M_pnK zt|+jPHs2xdP1o0(as^JHwK{U15?cF*M?`WBygIupa0e&;vSgI$=De&GQArxaN=V6mGKIc zZmLJKq963IYk=XOT8%TOsG3zuMBeg!X?@K@ zl)CjD=f#YohqIZew&S_rytxS+1TmRjrNk|C-OSwjw1C&wsAqE5p2kPT&3HO|G;tcl z(oKtcLjSKGY{aBI2pO>iO2XfxrKLVR?|KCsro^*vBTrBz6|L;LA@9Iz=Q5CSbAgG1 zuCQj~9ky(p0`_f;qO2XCVoJmzy0?BN>KtB(>jnhjdV~9@eyR#DoVkSc6A)8H2;q`dmM(BdRt zdf>qA#Bt1Q3~{MvTRf-ZiL1RCg9Cb$VQ%nm9m0{g5YnV7jh^zzz9`g0hxontyy*vC zU%nXB?oP(Yn{M=tBYSYkm=%HzA6CLWmD6m9dk>gC`$b9PgDRA>K7?tHhGXm622M5p z9hvWEkMm|Zk^GRSq-)4cVjVZ1u5y2b6*Fs?rktDncnuxqa>gT$OxK5tJtD})q&z(Q z_yIlsYzb~I4~3l(D`8KaJNr1Z0S?>9gHOAL?E5$feEg@e{zd&jq1qU-y9=n;3CEy# z*(0hoI~^{ieukC|e;mTzg5>xRIJfx`eL*B*+1lJPZlMGo36a9j@*Q*+Y>i27yfn6pu6ImtBwdW7i zm2HKA+&!+zdWX^D_R_~6BTxShhw)Qh*GOI8JpisCe(IgBg4kf zGu{g60F&e7;ee6Ar4InCwq!@PHvWDbqG&3^=|=Vh923mR^r|D>b))F%KYXdnDlQl; z^R4vN*hQGl=%UukBWSZPg^p?rf_F*%+0PXoC21EPp4#!2I1n_pHWe2vapZQdcK=!JTt}m0jFsZOsbFiwD-bjGRq$2-Zy*jlV7f6sE1l*Y{tL%vIJ{cwsve z({&8$zF&YgFXwUwGkaowSt3Sh4}uB84O)m64u43Qt7jbH{s!fzz~m z*!9^Rev{u%QNJ=Lv90VHk@UcBe0Zh@CaZbi=xf)BRyS4RHf1c`d*@Ypv|=i#OxS}5 z8gtpx(IenVXa$K}{)i#*tAFnh-R!4S%Vd9;bqvtE?F&})E5zP16wDNt!I|tOjCMp0 zH>{$Gl&`pkSH8O7>?_Yf)?^WmF8M+Yzn^Rud&U9o%1^|YQ-vrk&_rgzW>~zvh_zfu z!6&OA>TbeisE~C++mKeFoo@_xUvq>|F0eQDHoHn5*mcL0oFizj)r>Xo2D3%Ca^T6) zee7y22KHa?%bu;f3{gG;a7tSYN|7(vD!LMUG`oQ7x)p5TwF)S`vXfR1Zh^#NLXG>G z0N2L|U`Om^IMK%!LQPu1%5FKe((x<2T{8uxZ+=G87tLH3qviN|&TKkr@>zQBr;nuj z^*k=lHk95K=7EO>USkApamoDW^ClL`gi+qN?V;}YHL&6jLrtnTu=Dz(>di!u zQ-2P-d-kHU=FA8E86oWbB0m`3xE8KYt%3r925h{s7q(8SDeYBQf&EV2r2=B?G3=8$ zd*IRqs?4Lk8(J8}#QNOi#%)^6OS$%8o-dn%8@ir?u#vsF?)4V%vDg^5H@+!lrp#a)W-}a`)*VMS+rlKMgyx&8P<`1MsFrbul=)GxB+n3z=k%tmlCMGB#G9oHeLe_& z#9o7ni^t)e?&$9p_K;FOh zT6>kQZUym>atIw=0a1=0;Nadi;BmfSnKmz6_p_56kVX`qV{N@x_xcrJQ9+EH;+L5*;8_U>@{f+e?q$~Du8#%|308XKLx!oPkjTZ zOP`?LCTxLJkq3f1-)B>uD{aBFs~_E7G!DEUDkRlE;UN_^!i(H=Fho1%^+k#vi0de>2GSoaIcK_BXaW z*n;8k6IOAp55)D|gVNO5u&~w+*SRR*)9ZnWn0kz8MdU+Pcqa2w*o%vblYx7RL*;Ny z>RgV8+MOxbQkS|w->&lM>N;6-FH-X-j$4sV5&A?hFYV|JlC-2Ac- zrVmyyZo90(`ItL%XNL}0TCPMhWOE=XHidff%>)9bI&nv~6_XQlo-+?mM3HqOYlvs7 zEU7x<*b@HOG{SaR5R^#p6 z8p+zE?1Q3DUr4q|IV@ch3iHgA(cI6gs5ce>J|>>%3fT{LCp;t%ICbPDGY0eSY=#MT z%gFM}t)Q@`gA^B9kU5U;h>EI!95OQ{{C;D`Foz*CB7dU0$5-M-h3}C4N}9Y6%OGN( zQpu?EQ?6G-5oNYda3J-6@MQvQ_U3>zlimTv!g|{eqAc|gcNCAh-r3hVb zjPq8>b_)Hni+Ia->#&#VZDCTl3Uy*y6}v|618>G*FLu=4n5SsKWB1+M0n_j*W^T|6 zqGvFGnuorDVS6v5CHhwkg*_pHlu}|dWiFMtPJoR$oPYQB#=m}t7vAq7f=vyqghVRg zuFinB&$Ch4yUVb$UMN|t?_>w~4hy;~yO{9ByJ2AQ(ArxP?oAyl_dyE_P*;kucA5HAH{E4{@AHaPX}VwF!5!+;A6K zVOljbj(9`Ul4r1VmNz3Dser{V3gE~QE_|jCtm+>H)6FUPNc?%M8!jQVQ=G{oo9_x$ z>_d3GvaPIVm=96=v5TpS_hn}o=F^{USF(v7c{HjhV`ZHRnLa#5e53XdMgB`P#o7Z2 zTiZcN-32wTOQGVoY#~!d?OCx?pQ!sW|1Q5HkG(ncbWau5wkVA7ktRFU!wuC|h(Wsg zee#F17S*@R0lUR{=(PSq(EGlI%t`x7lGJR-^1@OOO_iih7Oo=YhupEHVGXe^vmnFW zYq`rCkCO7&1*9W1AEqdBLG8+E*k3u7wQSU3ozKX5X%0n|Ou;REeEz%zO^T6U$7=-8~F>a8!M>QNhX>FR$Z!rh2^r?(EQ zRyPviUwLjj@f+E%Gl3NCSqwSN$8_cxL##Qqv;-W)|Yei?AB1CmHeK$x#<^4)(I!K4TGrHGmRlQI1$C& z-bX^8q!X^QJp6?NJcpav9R>d$N5H9JpK(1SVDvX1WF^*Tg37r>P;IS(J@)ZP(e4jS zk_0?u=s46Q+(J+89foZa!^rE_>10!T4C&PE;SUAAgJAASa-by^d@Pem+_?fWWy2fr zmN@zM~K+KwyDQ9EaI~-<(&LFGBTu9@BZgPA^7l`To1SK(NGSqI4wx8Yu>5C=F z@uDwGeApyu1kppGC#*HW0j{4hl>A^Sn*_pMhylNNTuQ$%Bqzl2fZiqa7 z_7vzfYA`$NAZefK!0lFZCz_8BgKSw8o>y;5nk*WyW#}N)vOym-mFJ)=-T`=a`zy*g zs>KK|ZiSgC7nvGjPIC4ZQ`6>@!L6J!(tqX=`4Snzw{j06{q9q#K*JU`NUMxI7Yso5 z?HkziOEo%Sc@n*{Rl&Q$)>20KIkeSJd-(d^Qt0Thfv0h1KYP`=fp_{{92@uTxKPUt zffCP@O_a{(nMsSF+F(N5qxxE~Pw=N1(J{Jw&&mk{T6v zgT-c@5Rmi#Oi&Fy{z1!_{KSuv)&leKD|en#K10;}*(ppEn(?12TbF#bsF-bC;PLj;KkgCm7IQBRmuD?A8$sAGe$Y@~|CbYrk z2wx`n#B4}9kVzhjsSCZQzqV8TIJpj)uHcn#uVcfccd>I_m!bUc<8aXWE(8!e&8|75SVTw;jYzy?s!>_A_}BlK@%tTcqnODSfWqOn4itnPI9P`26#{1CCu9`X%rGC4L#NxW**C*_XkB8N{6pPfQ{ay)Y z>KoZxm6G(p5RHSECA)&|bFC4zeW z2ap;itmwEF$gWsLCS3M_cq<{E6y{4Fd@zC4ee2Ng>NZwcBmp%CYCx>~5|W3_AVR9fTVdrnZ&;xgNum12XvKZ4Uf4R!^6ixs&&9al?93T$s3 zAV+VWfaov#uwc@DH2-yxuw>J5x{Tto3RkXTpU)aN{J91@lX{3|`XZ?BcmpdxjX}xr zAiVsmCaSJo2AmCYAUS>je1W|?uh?#4_Qe{}UYNrQW99&6X3z)!ay$>ShPV?P#*KaW18U??+Eg|+`Ir(yghv?9K5H!u4$fef9fy*{9 z(<>QXKJ5FOodBabQhmhSI82pU}Lnd{(WrkIc@DB)88JD6_o+N=Mb% z$mxB&2eSG?H;o>;!fHBh(|Zi+u9y}7ejF8fR6s^_BupN=fmQePz*_yeFnsnC(z<(= zov*!^^sG6SJAwFyW#B7>~AcOMQINXx2$= zd(f3Q9ehSEjq4Ly<`!Kk^s}5@%DGv*hz|v*{kss0Zw#l~-9Pf2EPvzr zOHdxB$R=6C}RdK_0yS1xt3O5lg98B)a82aZ6^%A(=F=T{1-Z3`@c?4Nst^kAqiE z2IPGRB+;?C?Ad8%^p?24e^`knu+}f>vh~_@(m^?Tc|Zrs5PX6qS6_q49|@te#4ldG z(F*qIQ(58CX_INym;a$tHqAI$BZ`%NdKt8qABM37&)8Y^p3r0vht+PcAjfXrBWO!0 zm~DTIV(4k)mG@%UV%h-WCv_p^(Fr{9V?7l=*^#>DvYRwIH1c@g_pCt$%3KcQ@JN#>;p* z>gKe)X9v4K*_~Cfb*60#3V{Q=?P5-y>NtjHX#b|T^GL_XY<=P8Qb_%|R8(*CT&(hmAgbu$~AKSOYG z&un;ScM{B6zJhdRFDZ@+hpe5kWZ$f$>UGc8g?~qVLYyieHBN7rlA0 zF<;pK%9H85wJi8PQ6SBySD>tCqY$nM_Zy+7p|G17eE`E%q z-nddRfycn&`EoYV$Cho{sD$UtEMvo0KA~H-Ek%RsTWIHWcjiyteVlRJ81mld(OMU> zNNwj;{6fzK6bcAF4i7mT*c81SU5=f~{(5Tg?-N5;e@!9BUYzG<6v4iWPX;;51R^{l z2J_M-$g=J?qE5unV3uFqvk~iWHNOis+IrlyXT@#;8^7g%f#FRl8 zseBEY`e)cwy>zypH^j;WszX7G25-vj`E<(uF1l@8oHY$m09Bb!%${TGk@B%JKnGpG z(0M+&8|_YsWd@?Y;$&>_@*?(s*aOMaIzTSWn+=~YWY6Ef$HpGB=NZY>(W~c{)KNFm z=oR&q>|m}PR4hHimV5iquVy~M8ehcz%@7~GCM+{pi0^)%2f}ezQ1i=#l-qMj)$UkO ziq2srChx@5qHnl-kt;LovXzXK?Bb5-tbq$Rrc>{HM2N198Kgz@kwe_RzowC;({^vf z=Y3r1>{eqszM4XUUm;}uqItM-ZV)z>)5B-S8en~@7I<&J#@gigSGKQQKzq z(+Yue@VryUaC~hkbNBONSe2v#${Tu^v8(3fjM+yLq&kmyKTana(Htc7X(4s@cabmK z9<#H?+F9|i4NzR7kDGP%vD)-2f5Y)BJhN{XoJzD~ihdLFU*KMp*7=ZXxu^x!54@?J zSL%U!*oe~~>M`jniW$k+NTy!(9%tl25pzsDOu)^&L{57fu@CN9fJgHSJNlcOpsUn} zestLi3gyfJk-AZ9M->7ta51pB#0sLW9W?;O4XuvII11Y^1}k z)|tr0Eepg4H{S&F{qkI$qFawSnyE_z3z~+kQ52VlQY&YLL2FHbVU!lXw>zvTzT+Cp3^B z#j#JnfK3vgwl&knV>i9={3T29BNtyVdfW-gaY|VD!W^aONy1$HCX(~m42ygpM+dn- z@Q-bsIIrEA;qiWu(k-f}IaHs!E7AdYOM{7Q@Oe_ay%?Dm#*wDuNr3wDfDUm4uhC$# zZ)OTK%sokK4QoL6ty4JTukvUsN`}p0-Xz&RoaE@lA@h^dP|VG4vT@mxTCw`$+)=Ss z;;+|&*7xZl@9m|`o%(K~$w=0;RxDu7-Q{rOlqH$j6P=(x(oh(<}@T(dh?Fk0}&K2NSD5e?@~z z2knx-7)qGwwCIu*;Af^nSVKjUUf0bGPk2L2L+w%Kxue{dhij2%)=SWl%Oo`kvc&F_ z1BvavN~E`+5~%*Sh_F(FRG)Y**=ufvbQDfg@-1qd*9o>vWmzh7YL){-xpngG{Pl=l zMLtvQ@%i6?*f+P)9*%=>$8!s<|5yuRpE4BKnXoF;4CWB@kAl3bn7*BZCu@ zU|OUo7$)C<&Eg0}h`zuRn=u`!ypfbEw=?TT?&6GjQb>78Am7Y&H79?FH%t1#Wk?;| zPp4gxW*3PG=)?1tvJ0JOl0Oz}f$l-f_zyRssqaR+i9vpKB9&cq<8ME9qQhEzGRTHC zxZ+EzKP`Z$3qJ^N>rHrJeH1j$#zNf45!n2qfy}ts0gv}q2t^jofw$fVv1CLEG!(@{ zy1`9sbbATWdUqa0DBnS9#fJR!w}!}GHVcaH$AYL!1k8)#0qMR9j`yFC7~M=-MzEX_ zR(BF%@FUis{ug9c*;7o2FUcLcgiWpdnYegkl)mB}_wJ5$Bq>h=rfXCZ1@A=0BjSdD zj-1Yv<6^eAZ9XK64=$FJ>u(Eh^zSl;XFkaQ~; zl7;e=KqCvQhCYPc-cGz%x(?S5HREquz2tnv0)hSOSetojdA3mlUuztrFT)ArT>9pi zF?)MVo^CEoU<(GO6Zjz@!E;TS8rLVR{2?Ry(0CZHr$a)RQ*(kI-+LaS1}Ko*GM^PO z(I=f*XW{w!UfjIhkv5y(Mg~jFKxv!x-*kD{e6cP-=2*Q{h-mgtOszx zu@4+ezQR2zJ1Q==7Y80FL_O-K(asPNtgJKKO>~ z_(+wU7Zy-Ou}6r`2Sf?Yhta1`-Bi?9H*j-14>dxk+N#m}T>n!>#O(0^PFkn>Pe3?P z6YK?|?%$ZW9TiwEK^i@@UBv47h_Y@4N?0;Yi=J)7qi1wIrsY3n&{0-#y_F=P;D)_2k(8QZDhnJihirjuVBvt`!90MD5X=T^^B@BLC++O!0;N8@CZ{cPLv*KE3LJQn{dPrZF|o8nlj z5UDwLNa@tm)CJ#hazj-OEpXN2eu+6qy9MUZF3KOc%aA*B&^w8IbCSm7g2%rZQn>s! zr*}cAVAgvRre4Y)+gM+uwKI>Q~xven6VFIX6?Gre8jhSe#Zy16lG^!*eGfOXs85 zzX|qjxiu`Tenbq$KcK7XCzxYXt=I&qR9@90N+{WK4W5Z<(N6O@_=koEieE24I}Yrm zwO)16241^yoKi@4xgv{5D-2(`!A}3m zC-;{xW)8Ij=rX->HYm`XdXll8)>Zp=QsA)fPI7&)QvQuGh)oltwQ)+IT&!`g|=t>ApUCGs0t5 zO(_@Xeb(YHp4D9Yh7s@=Ot{7OE#@s3^dgg{R(Nm|?AY}4}sI+v#(^PX>;C>xuidTT_syt>^-YYfhf#{KYC)EbiJUW9~~ubHI_KMM2| zqB#c~&1YUe{aYZJ)+LzTT*yC~+jB{7Y?Hue=n_+KX(BCkWJFN4>@oS^VnB9X3n4Aj z&XD`Bw-Mod8XO`vpezL$5EXpIQ~lJ1s<}42r{`*fH^raw>}W@j+9VB1?d5ntIt0XT zEW?*lG~n{30yL}nBH37Y7A4>Oh>e}2;q2QX{IJFt$>9}(*0HrX>yHwo>261NU$Q9u z_Z0HzrYULeJ&X<>;iFKSn!l`CjpJE6=&y~0gEtS8;}k_A-@N0Wl*uA()+Oj?<_2U% z>%m5*lt|>7k?*1_z$e)k7yY#X(xC~$9}CahExaOMXVRrxhs|7h*RhW~9uGM1v2FZ@{{t1XJmixos*qp{u zZ~VsB_Q~N#`Yx?iudd+F&`jkYdXhZLWpfh8O?C<`3W93K5bif)jc+8u&k9Z{rw zt3EST8)8Vbyf#*Re;*zB>tj^+uV8e{3~-f0B=%o1O!NaQNtTBZ8LsU`;yZQ;e!Kc% z9gWp+%OBzWhle4*ZVs(>+k;Ju@uVHL*Wf6V5X{dJp9}RQm?zDh5nd&N_}vb2FK`{yp>ZZR@`o@?)JLCv7B< zpWZdZ37xvH)=+)xtazRb|MR52nzy2Tw$qTez`VzttHTy^#I_x8d&Vm56D8Zh%$q;gENifPc{Rqc8f|=Pk3x8M`j=clbdym$hZX#tp zTyPA$&vjlr1lficnRz7+_-sZ6$?bj1%5zt+dOvyWqviWxdFpJa;#WXX-*zlz`jCW= z+mMr1{)nnFBeAh9$nn%F#`>o(T*>6*ya#cym_GJIA@{d6tuI z6;@+P?-ERyJi>LrpShFdYPd_pZo#5QAL)rx{=kR%_pqF9E(IhbSKPHN9B$f{ z6IZz+aGaV?b(S`RWV;={vJ)bo%e}^x08ervVmmxN|5VIZji0r?Gs@E?Q zln@bA@@@%pA}yV0EH))R=1(Cs`vom?D3TZ+^+GSsxq@4*7gAX{l{~vxf>KJ2Nxj)| zQh2@%T|QlbqAPsy=?_2#mOUhgw~Lb1&yq-4$RPFMBM(vY#^{Zr?zGDB16U$-Guf$9 z!CbX8C99vCQNiJlVBH);oWsAy2Ib7e86$-C&oqZ)5?a*pnf?C;WaU(UuHAo`f{gk4 zTqBV`9I=`}{=zq-oUJjdYgnYade&iIR>^BR{_S><)j66^ef+$ZikAJtYESrr zih3-`mz--v>hujJC|-tq+PanW4jyK2e*VL57iDpa#!;T9nhL~3@4=VqEm>%*#PQf7)BTR{n-|mk91~Sve4*x7Ekn=HGJz# zVcPb)F%R~Jlc_k!I_q6+nIK;l9)xiL3Q1bNdZ zaGs0H*H*2XXZuoCtJdMxCGIkZD+05aCXUx;%rXBsgI_MATRV~;!}S_^DUcsI!xjC$ zyVgadk}qE!aB1HjJN{A9PWF8IipyKBLO@jtvnK2wq~sR@JNYb2##y>U~#N6c~JL-oV`ASeJ5TdrM@HIfGp^%W2vtOTtQ&QPzm9@;_+p+SDT;F8lbO0}vSoy?8Idd}ID_SiDi`al!r zJv~bTefOfMs%oGeTF4yx4iFXpN^0lxNRPu_GSIw;C2uep_Of0Ueqd9$32~XMp zsqfBVE|#4VXohQ%ko#ZAQu#3^>A^N~sBk;hd=QRBH;18V1y7)C&>DujFHnwaCAjI5 z)&g)98)-rK?+>wWB%L|T4Ve_Z8JHM`bJw-bYUxbwEk9|w2zR-cAlbVESbjJeaJ^>CU=9(U#&(itC^CW ztrN(J+L`d*>@YA5cz{xfDV~#c3ntY*U`qUqNmQu=xuf|T8GFbf*TycgSy7KH+Lp(i zy5Kem(y5#EEg@YxQ$b8(1K1mx6O~OX$OUV6X370b!4v(xL`+5rt%!X_ zVxLu!wh(=)@U{HkXH4DlEFm|%($V7WE=>66R;JNsE$RLk#R%Ot2u5-WiB6p3V#;+I0+=v(M>tDb#jA<2$d8slU31j~bS;k&9mo9%D z6ioY`&K1+E7Kk1EkMpEb*Y>T|Jictf2Ep18F@aL5rC?*ra)Hv0T&~VO4NhV6FaC|7 zWZQdc8w843%7VjHmjzdXTgb$cSS+(;gbZrhVtI`W(qXw1R%`!8{swCxLCp*|FW5zE z&a0*OD~r};E^@0+JE~ZB)O~#&Cz1&F$PV#3PdpZ$8+Er+8NDielTrbvBVQ9$T}u*U zu!(d>Pa};-AEQ0nqexG>B0EPo3*N5Rf%G{VY>dedV82+fYWuv1>fBzq=z19v&XmBN zN@? z7q_mXbD9k(cVG{SIAKbXq8>AYJNJ_bNtxVOxfN*h4JV>^c#_}-N1Q2YeaKwQ{={GA zcZ7*DR${KomNAE|9Y|1y1 z^7sk~$($4WyttCS>0A@Hqk?)rP3B&H6YkPk2#sm^IKPX5OPdyhz1&Zjcqx=H%QRr< zpf25hA%;yUwX&b|sKh~kp^Lrt{vo@{JG*%A7q!vMgcrPP$;tG|rzyhxN$VgteXGD= zr7DCw9wt}b3qWdr6J!_Uz*Uo(*z98mG5k*(4!T6(@mPP>`D+>^-zvn>aROH9ivn0g zPG@I{R#Dg9zem?5>=mR<^hA@4BuL`xV&+ZfP10)XLSo+RAldIb$-x0te$){|>V)oG zZoi5Z;l8NlUwd*J=2rhl;#mi3^wdM<#jj$fA@MtjU-}4Xv-K!?F@+rdDz~>~;$*sF zE?MLxLQZe?W%T{tGnd3}f~HLbEIy*kTUk7V7O$BH&!Ym^?LK3mJgPx9Z7xL21)YDd zJPPc-LpYsBNQqA|XdDU>JZVhgZ1yqZt9bSZ&=)s;e1(eO;D-+G;fl^0!*g~Uw z5z4+=!-lw`I@9r%y5JKVE}L{_)t!!dTF0(FM2i{Bth0>e+lf6=6uu6sAgu-t#G}@R z^gS#^qiJhMRahq~@{GVf=AY0;2Og`F_5qd5QG}CuMcBIFCEoAw5e_YX1ju^@1f)Mk z>r7sOX^kJ%p!g8k3rEQD&JofyV2qUmTS01mGtBVsg^BZ7EH>`U?D8(BB$OwC(#e(V zlzba@@>x$*aoHYM#?HrErvF4!R_YShd#0!^n~>jg-jYVsG)g}_nG_m4KqXcwOl^Z2 zF~8oz9FEH8rj1(&u0IUrXJ{R$Qvwe=huQvhVNIv(HY}VftS!?b3Z)Ckw=PFq*0Sf{ zeOaz^9r45a6s|qhfy4I3K;oc3iMFd_=Hz4B_@?O`;=7%rJZi?*f5GLfJ{K%d3c0~C zamnE5{I=m!=l-juuifQNG#ch79J&wc;vZ>^Rb{M{?+l!=rWDGB=}@9wK)UX20R7Fb zSmHn)1pYLvquS!_vr3BWf2eZpv)|VX#SRA|>-}B4yDH_Z(xYeWtlD%GIjNW0aw8oC z?Vng{{}QtAnLb&#beQz&D^lmj^=aL-GOqcamk^OM6@qQMVS%bN@td>~#mbeyVV^Ws z>yi;vEPRMJUwulVKA(e`?)ON1=~iU8d>f{Uny{DB1d89YkW7(wA$G;anE7#;nt83rxP+esW4ukx5iaBx4_M zqWA;bh~Av}%+q8UlpB-}MveBuh4)_A=`1@cY~HBDyS>Un=$G$AvebLYk1J#NeAvH7 z+L(-E)%2~d-C!ua0-x26;WBAy5xKJw($E|Azo10-|aH|Hf_!dk$T)u+tk3LdNe@D>{-*KvnCqC~K54=Jz*5Lg$ zHhT6JHp2KTk!^cNr)Mo@MGcP9cIrv&guFk@alSeuI`RY`t`BB)E?EHcCJ$EM=a8do z51_V^8$`Op6w9ySppx=7sz-kYTQA$qlN`wxuG(b4(|>fC^;vzJdnorTx~Lw4x1a0# zH$!ZSj??A8yfC|Z20o_0Kxnxe1pF(psBE?rwhTNDvTJP_m+&XHz9N;jdh3_-trltt z!j8M}Ph9NeG?h-Y^=2;epR{o~>V4^4t#3TwiN2vFU-ghH<~8($j8Nh&Ooq3 zIh2+x!EJpOtoXP*>wd3B=&2CF&VMit>Gyiss702 zEydH>aoEmJc#*eb`~f|1K8wzjRu*RF{}6tOxXwxq3eocEiga2s5630fqZffM=!Nd0 za5o5O8!=Om36+O55i1f?u!-DkX$JYxOVl}!WTZJ5gYp%2kPFqWq@tu$utVh%{P?WM zqkr6k)KU@rMR7LN1UJE1!B8(kWhOHXtV!zsro*TY& zkQ01CkK?dlfUjbd%dgyE$hR*Js%@E_$jWY8N%v1F5sDZO(ec@z*z6x!WLU+Vd}}j+ z^R1teLr55Fkt>N$S&Km6Bppb;9!VPYG;xS@FwFJ;gd93bnc1pCSpM@X{Hh+&aRHXt zp9C&tdp}XyB`8-lRRgjbz2Pup-~z!!iFkAgH)NyA;e}uOB@?hLqm3DtnI*%|rfF zr^_6iEaghRmM*1^{Rk(OE{616v1br}eHtcr55t$e5lro;w|LPBIkazdKRPd(@$Z41 zH-8Lij}U*7b0{AIOu~dl;pZVOOdcog%BBsHi{Q#(E&B7F^SDfCNfP`_1wl=_`9H)m zYu`AjbNo(aa}26y*q-lJ<4>Ehi+|{-1Q~C*h;QO%jgfdk!fYv$ z!r37xG=4~y+b)KgPOc<_yA35RJw;DwU4pY_AH>>~3J_;}nLgU=j{DD-*-$cs5j znQt1!vAymf-XXwftc%;cBoemyy29Il)#Uv7v!tJVLZg2^k>z*gK;9*b5#0axNE@An z4+zytR?uRLdCcIPU^rbajoMa+Lco$ZSju}%Z*e@y3gmv^@aS=z8&LrCr#GBSB7t0y zE0wjAiw|>fORs<~)#9s8*~1r^@|9CD_c9qjv6r?p8e@N6Or#}DwnJ~XbPsEwrV52J;^^Mb18ljo zMcrg$O`+EBsgN5oP58y&r`@%-L7oJEi=7o`5?nl`1-ddc@1%)~&|%iAx>%`Gb&I#H z!jn|JNZr#y?#QWKg0j*t0uM3~ln1U7by-_7_0R};|0tVOC(a~q=lK#zH)Rl0(ue*N zGsr>7pD41b8Ehha}j zEut+#m%#t+{0q9i1GX-@1`|KMWCrdz{eQwhU8Ii14)pGLX4xc=Yu_l*} zc`^^zR+Ezcg;>e(3hasZAWA<(iKJHs9^1SF1f`+ukq-^5^D|4(SDOnrQfA=R{yFrL z5<3)D|AU?S#tDi=f7Kbv2Ek9aKcs7h37d16(7z*_g_k|^>|C8&>=lbk>`r7qfl-NL zY?t6Yb6T+&M4mon7D_!uA<~xC9Ow;HqvRS&4HvD=T-;Q65)xJ+;e$_t1CL5NJ;y@Z?>jqR9#}r z(ne{OyVGHoZ#~4mae!+EkIAN*o!}~ej}#b=vd8}01QYyE!n%rXC|9V(A6<08NJEd` z^>ErO>bI5cy_hb6LKSr>UrGvPoyo)UH^oqb_74_4Z-uF1N4VU_YSd}r%(nDNB>UJxBKFyfl9-hbX{8rg)kW=)YWEc< zZQhIX$F}3+O#+Oss#D+RS`*Z@1hm(0#f@$vaLrqQ_uF=mlmRiOH*Pnoz8gnw^WWCK zKZe+Yc8MA$0ara}`R9YUt>H~My)Qb?27gOJ5vWGP@r4?c2gyyB=aS)y#!Wk$3F0=kKvs-BfLNeBU^4z9mn1*wKN0 zbn5}TO0Wo5UTGCh$yz1M55ELrN2^Kd<6EGbB+Z0;-^xsRT!hy!yJ>^KGx(^*Wm0-g zlibyjAx)1ik)Dc+OhT|H_t24zM6vfaDOTvmYWMzQysS49M@}2sz91E3RmIV!x?wWu zxj8#NvxyfsNDFiOs(7hxSA_;QW+VQZ5^C)A3UsT>j=CmvWQwOv`8OcruO~oURtkUa z5gD*_Qy{a}JOIx{n;><4D6;mr3A=0_!*dY%V8A2;{$887_)wvS*+xY7VNG*Mt$u1L*8BAi;Qtbg?+rQ6lreuVjLaRIGyj=!1>9KXx_ei07N+Hp8y-1#HmtZRb&Vz$ncJyWOuKs-v%X;k6a4cArIVk=j9ls@c2*tCt}D&R_RtKbCVw6GwwVcYa#f~T`!g_W$RXYPzvGeX|4~6JH{(Aq{)0Pyw%F-I5y{_XMA}UQQ1*)R;AJn4T5Zna z;GizrOlyplQ>mp#_G^LPO))xA*bnbId5@?LU|%6g_{&?oQPYUKpvz#yu7|F`|;bR#F1Ruk_g7#g~MJ z4{r(1E*Arh52^yb?kG3Y!-`2MnlAX#`J8cCy`1ylmmjmU+J>|TIx+ui@64lW{N4qg zCWK0&H0Y=Zp+TK@@Anj;kc@@C35iachf^8)ijb6wG8aP06oqu&=Xp;u6b(vWAtZ#% zQ!<6<*7rBu`&-NJ-oJjoweG#=bN<--dG`41=Y96t`>ge@H*KyM)|9M}6~^qKN|i6r z{eC_ce*S8E*rk|@j~9^vW2(59(beQ*r61sv1aO=f{cDF5<{CrNqvtTVegfo9_zEU< zg%Fp%6~h*1gH(_Y=LQXcOKz2T@BLmrLnWFGj%WG(mF8^q@1<~|`xa?t{3|>xO%e6t zPRi8tk*&6k;WLvn&?WXBymFC{3N>SLJyjR)b*ZE7pK3uWOJWDr=hAIA8>wpg0N&`^ zAzHTJC_Pu+2-VvFcGj)rHO>Xmgt%BpdH0q3fx1EmeKB%ku`qU70KIj6D!;nRKK6xS z8UJ`(1v@|DINK|(0ltME!px==Cb7y0qpOahrkOK}((J)wyALs59*8sF)l^}VTVh89uEB2iI5teVWX4z6DsZUP-!@W zD)TMSe#aJESFsB+dYj;jTnS%(V?2A+-~oR{b)!iARSOOnW&=v0H<_GMIih{aePxb2 z&WcVHx69T=IH2JvRjgUJ56gaw!!gWt>BJ|m>Bg;qrGs4^1 zEuIJI;sP&UVELM0j3aH?~Kh|EmM4Vlr6To zE>aECpwIM{Qg^|5kVlY?L!9-_Az+>LyaC&==XIx+q!lkfdZr(*g}GI%{qosCu6y zsvfK9Lz6jls9`g3)n_qHc`SG?3gwOppTmaqEYvxWN1h+&Ko)qFggLCI`oTV|>4*%T zPcOurZ!_UtwkdTUbOTlIUtlgx?TwM$bs#1laQ1XRc;9D7>GxJ6k>?H_cC?zmXsp*E*|{Ol zL@_V-$y5u*$#fO_%G8Hy$pWWa^ZWZo(G`ukRKS?R%5V96asFGofTurXsh_h&yG2uE zW3?|zXPi;6ebjuLlsG6s-O+5$GV%)}>ubv0)6DtxO+B#nEHplN3~1JIrrw)?t_Q}D zHEQ0tJ1LK(e{IJivjLp1Fb9f<#mfv+=ZJi6h-9OcVniNuKJfaIJM`g&B${k_fF2Im ziN(F&K%sXl2CrI$jDsc`YpIvGJ&o;^aAJdFUbYfrLfRI2@>)cF6qo<)QIg0d&g^&#EG-0;|woq&V~|N z`(bv0b13`oFB~5gt%bS;TF)c?CY04@J%R320pF zghSgu!Ht;7=+^BKZ*sDhvsvZDb4_R1r?ngS>58M-jgmujbnbN8dN+XFovQ@{j0>5K z4WH1>{4MN@oQCUvT*q}KPq3)19hDl&Fuy$sK7DA%hW7Hb$~Ie zvZcCPPtii#Dn{ga3=S!@!}g7tT+y~BuqxjIv8(-1p?EhGaGR-tZYgUo^`V$%$b50l zBhI-j$JsofzYX2QHH|eXZVYK)M?UuCG5TH zHa-_j<+Uz1^QD{hc*peqFlBWO)O`Plren6^$s?th;RW`e8RuHXs^dd=*?z3R7*yYF`xTGHJaKko`;JhLGEUZ7MuWH=T#dqPgR9#z|<1>2@1 zfbqSh(r*K9L#PTRXPQ4jmdZQGJYRwnzRjWPiDQ@*%WTR06erkRk-_;7>;c3iMc6k! zmhjRB@IWEuyXcYH8jzbNJDI7i#Zy!wuS3vG;gi4DwZiY1uk(;P@Kg zuTA6>13l0vRLmVS9|x0np2c1P%bCu18HWu#!YNlZg4DJjG|h@&;&Pifg{&)NYd|PE z&27ga23`3{w|2p>53lI_>U4hkX%>e$77)XOA4&3sF!*5Yj3a_Bl7gMPNanVIxSRV4 z?pgD>{sY9qce`3(!BG*8PUwb1oOwK`AHa^NUQD%>!fDEj$<*_Yo}{70ma5$bU9^b<$e%AD|Ug zLnh6u27~ga++IlwV}2}`oBQ=M<~`NsRcAC(vcG`y*vyjOyR~xb%sxT)J!MkQsAo82 zwj%h=oAhgkWO3Kw?Zyh6pSBw#%PulQ&%8#a)|HrF%fXMD{V*?R5=n3LWj&WsHn?7o z23Gve?l4Ye&Ev#uEP0PbYu(^X>Il4|v<5%=mtvsbG}w6VAnaYlLfdvNPUK8Dm%iiL^!Xnk)fq}zp1 zg*%F@%2IVQ_*y2_I$1_9zAq;U{w^qZbB7#Qk%@Op4tE?LM$iQdqe)rZUTUEH8oR|R zqov1H?v7M~8GVW!wmf^ZV#R}D7h zTbP!gRa}iuJXz*j$bMQMN;k_InbA7~i5l!=1TT{5-gyHVsiZ3js7mLmw}x~37Z&`J zFMWo}xIXzKQFq~b(!br)cKd=yFka|HOg0<5ZRK z+H-PXaue$IoQ^u{zk-HpF}bj34f9iJGg-gbjx_J>L-sB9#Es*Z;tHK;;d}F0@N)k^ zRLn4f)8?lk`%nr#IOGV47j~oG$YD6~n<*V2K)TB6G99GblZJhHYoq&pHBS3&IvG-X zfxIg0N>w+>z<>N0q8%e5Z^fsW)gzsxRd=g1hO>qV9cUWWpYGR&r+ot;%_d zGos&chqOib!Tc+Ho}7*&99uB15V3ojHQ=uAu%FJM)vtu8?p> z<*f~wnW;_|#@9;kE{#Q_dOhx4Odam{e4ER>W+wfX5Jrv-&qA%+=G-5?L10l+3$JTz zp(u1Mn#`GoiA~|?U9F5LITM3CqxTi8!e6k+Dd9}pg-y_J{1L*1F=h0Y61w0{R{eTeYL3GZJ!q`k7Aa7d0L z8c%!n+p8#0O@eMaxPQtI0kH=LCvf*Fx)YQ zl-PB{J|k{`@0kdIp9j(Qnm%lQ&ok@t-D$rgF;J18iUz(3xFtaYC+s*4>u-mEM)qMc zmm9_jo}40W0gupdPj7PBBZpWU^nr?FnhaHFf||BO=C0>pZpOHK!lJ-CkYhCt&IS(S zzBjHyLyb+S@$>|`3{5~cH3e9fKZYt5w-W8g5!A}<0S!`|fF~wTr@F;)sPsHq_$F@^ zDDO5TOT=$Mf5c#+&%vjKC53&U#~`tE^EoZX-y#iyg$EeX*Y8Zr<|VKvQ5&Yb)#7V= zexj|-kE!LF!wmCy8SR&L8rPi@Ft@MIgUnkQjH#zHqwzKD*A8iV-j&glUPCYUFfux@ z41+y`F(a-Lr&V6SCn-DG;WgVtsSmR447Ma=&F26%_0l|k#@Bc#KWl_B(}$pyX#pfU zKY_q959}rIE$lXX8hHG0fyU5_g+Hw*7iBUXD#v&L*RP4$KiL>+hqQsuCn@}Kp^8&h zeTzC)4~5|m`-3Do20q3gL35oJnDFR6xyhyizfYfq?vU`~RxaR0Rc11%)mONyE-7D#-b_G16kBg~JT*<&(tTK4Qtg zW(>=}@Xu09%TYG}neWtFSz23JD+reVI~E`+BGOyZG5r5i|63D5aHM~**xrg6#RzSP zm5_`aWy?@oOKUZGbo}2(Eb$SCS@`)!ME>i>=nQJ1@Xr{xVr(t{Yx*Zb8%Dv(hOx1+ z5!zUG>aDD7ZEO?-mjBN9_YMj1{_pDjPv!qhndKKHnCu-QupT9_vSKZ*SZkrc+S1xq z?av28Ekf+mu{nf9s&$&LCh?Al6id|9)Ru@t#1ijFv5#QU3V~!~gn#hJzt^1QCl>ro z5cmc~Ef@Iu2Z;s#5uMsjSxxO9TL{EK;-%h^{-GiE#(&Qm3zl{CFY)S_M5Ld+@t;YD z;OiX~6lpB*S#Cd1%{VeR%*$tahqS`}FZmT6d2yuoV!zIAy%vXxeSJHdM*4Sj7b)=$ zi3kgo{9WIv@SmsFQ54}P_Lg`>L`uXF5gmH(sLuMq-XUH~y@R5}_Ld_pZU1hOSEmK` z#(y=%SRj^!Ms)NN8qzt_j_hB}=v+KwL5P23=Q4GuBJ8a!1)))qVNsDC1Bt`ztpz^* zzP?cr9c6!=kmWqJkwnljDWA|_uL!Z&$KFa15+(^<+|gFP70Zz$M~)mha^%R7BS(%L qIdbI4kt0Wr96562$dMyQjvP612O#4y* literal 0 HcmV?d00001 diff --git a/raster/r.sim/r.sim.sediment/testsuite/data/reference_sedflux.pack b/raster/r.sim/r.sim.sediment/testsuite/data/reference_sedflux.pack new file mode 100644 index 0000000000000000000000000000000000000000..c684026035710c6f6c500927a60ed61a33d193d0 GIT binary patch literal 38550 zcmV)HK)t^oiwFo9Xp(0F|8iw!WpZV1V`X1+Wn^Y-b$BjtVPk6m?7azCjcwQepEPLJ zJRlj$RB3Ks>w8{Nkq{z8rqU!DD3OR}m5L}0DwRwnB}x%8gd{5SoOuqndH!|Z?dSFU zp6B;_-v9AC-sgUg|9iKNL!W(K*RJi@>ssr3p4VRM+=fF9hYoji3$gcf^Y9E5`L~!+ zhyCAZWMpdg>+$c$jg3raeM92KkjDI(Of46`2y#MU{|4+@ZGyiwy|IYmXm(Tx>lO~LrGt$v%)W0eI zEiBCb@A==@^dIuSiMdI8{2Tq5@&6ai|DFqjyh+E@v@_iQn%1fRo%-LY|DF2(hyITq zZ#Ut;LH|vS&Hh>cO*{4f&oF;CKTpzeaPyn%8Ms)-Ho)I2Akfdv_xDFW?$Qeb1Ll*C zue-anhuh*{Khp7c^AMOyyOEBCsinZ!)Y8yM+B1lBj0}JGFlS~eFfuYR;o6t`0!wLM zx5aae+vgS*hQ=0uZ-1OHX+L3RX!`eN_QwgM-)AF+mgdsF0sg-a3@kb6kl!Cz7;zRx z#)c*)(#w7y8cF*H%<=H_YJc_ipFaMc9-aBBQ~x{lzf=D^?f+H(r;Hyx+3CMw{crY9 z`fp*`x&HrCEZE;?G3ku+T-_IOIx|71{&(tsr~Y@^|6u+1^!5F3)Bk^}|1$Y~O3dH?snE&u<+ z{a+)2g`tI!sgc0gvQq&5V)=iLpHGnc(0@<-|9`CiHMKOcY>$6aV+)hc`rpp{-z<Lt9&Go-|yZ(260Qgtc|3Ei?Z_htd|4sh!{;!#tWvBlCsW$Ap`cs|GaOl+kPW|uH ze}SM=|NpZ3?;q^z>+ZKO@V_?x|AqdW8Ch8V)BXSU)1CVNCs?-*bR8t5sgT2k94DlU zklsQH?+sXLY^1$rjSpC{Ln$Mu0nPfvX_v9g|rlMgpd=3oGs)6 zA%)r@79-?(AvX)TOUP;=PYQWK$W|d=3i+jj;*vtD2-!zSLm?rgy^vFcoGWC2kSm0Y z7jmPJTZAkV@{o{cguEi;10mlC`J;mpGD2zysVk(3ki&#@5OTVZUP20+03^bNOcF9n z$n8Sz6Y{8#=Y+f=oC|j|5YkflXBQF2pY8TK{O`i2d*`5Y&^hQF{3Qmz za$x7gxpV#Bx&Hs3y8iFH|Nj?W|F_4(e{n1De}4b}pFaOL>HPfv&#~cCWx@&u&PX4w zzf=Bz(WPag(Yl>M(W(EP`u`W!|IYfqzq0=Sel@hmgx|tHPygG4{eRbgb7P~<{XhN` zGt%+%@ef|?8Ppl=o%-Lo|HofX|2y~p=+u9q{?GIH_r?Fe-2dC?pZ5PWwe0--|4*@w z9lkq^LN+sAO4K+m86U;2G4?QanD_rs2s)Q;yH7)$@HE78+~ zbUOQ#JJMFN@D2QTy~cn3{U1i1`+xp7nXZnJ&OdhmN*j@h zx%nTPe;WU-O4HJevotd|HvRp~AODP`|C96WzyG(h#%7$6v9X1*$scEoO)UTT?SFAr z$C!*c<3Ao-m`eZs;Q#R7rOqu%JN3U)|2y@+)Ba!e-|OGF0qDQf|Mt9Z`p@70X=2*> z{U3jd{bT(<-$Ajz)pO73pro+kPfA$tFQe2!Sz*1pyjceog!R^nQ#+{ixPvOfT&^b7 zL3Lq0x#o}#YSnbGr?B6k*X<7W7S_uP%j`w^3hTZ5_Z3cG_I7Z9upeOHr49}j*0T@! z*g=E14i426UXn$^^e3#>H$B@yb78%?MOz02t2#)+egx!qkQdg|TOIG<@W~yte%L|V zB^?|o>=&@h=%BsueI9+FgJZ{b&_UP_G2X9(6NL5r6W4aoNvDI8_jGWoO$VKY{R-2C zx9(ho_4YGkJ2+dtgL8!a7;eJ;iMhgl3=iSg{dfxd3A}?l=p(F$pP$}A-+>+U6ZT64 zj1k^G64v_%3Uf!$cVScs`$3lW?qGtObNAsdVbRt){_x=y710&Tw*x`?b9Tr2uVf}(XRzf;50SNIa~5@|0~>E@D4wEM^% z+MRrsa;v)2+$b|nAxwfhk+++)YCclt^s&_7KN%(??_lSqG}bib0m~fm6$9QK$IlH0 zkodU|mZ$b|cQrq8OEho6O7Am@lH!nZaV%5Poy&%deaJTW&t>;qjal=Gau)tdiRsBN z$J|k+YaPb2OCJ4w2vniR6%?rZ!ckH`08yNKSlk9=Fxy8X4XF%60cL^5>C0%sWvkl=|2Kzr7M|}(UR;_?Q=r6LeR&rO1ip%}$%O%tde>%uw5FB3?q>=mqA zzCkc@=Qn{=GG}S8QP2Ihau|8+{!W*Fl*8@LXoRTz5f~i~*?ZiHt}RZ0pkX?q-%q7e zj(g~N?pccI^@YYP^5wRcFXPHyo~Mv6OSnY;o#<*UOAmvBnc|zl%xx#J*i>m49N2~) zmu_R%eTD(Yy}8IGKe+eyR+Q<{LZ^pU!Q13MGf-KGD|DK%J#wtai%e!z>c=imzQD9i&K^&f?^|cMnzNf1q4%?7WPu)f1sN;{`PxK47)>!&9pP7i!tTsd~&-KZ6saRz$hC|*34BFijc68D`7xD~$~q8;x*?{wF}CsLg5`4&+51#7PS z{uFMTk|UQeWIgBU{){solt!}-&!(NM6?)6_5I^lGK1?ZPW1J2nqx%~=JSi7(KhGn% zrG;`9jOMNnSVM-P!Q{Mb3M^VIu&MYkVx?4Z{i z+FS+Mo$2V7e}$UH>v9LT?}6^n;dCHmHI6PDhdj-f{N^)3R!16~@Os!#ycl)}>r~yC z}6^t*AIb2J!3jZ4nc;p;q!UHVOs zSR=tnB^{^06X!{N>S-F&)tm}P1yc6v?}*WD`VWT)Hto!!^d>#p@qG$+e$IOe8a|4f zUDqg>F;88fq@6A(RWuT8R!HUqY!;nNYk`R5yd0Y;nGKprat}zje~B;)3#-O=k&E&A3e#+Ymx*1?wSMo+U=ww z&0*-Z$AI$oHZj?onoLFIJF|Yh6^D(j5wyPw{h~k*e?}wxxeaZ4I*#s!e56bA=g7pM zmI`_nP?Gi{dLq>xBwP!mtgm9^P(yUC#(bmfUapWolJFl6;goc*aBFpUaX|@wob7>ig4v@+3Pj$y z3C36)vs@$g!7|flpT(1&7r9B7cGH9nIrLO*J0_gS!{qSI5FNf8<9h6+iVGL%{IUad z=2Qds@_;jU$ljWIOi_VS@9U5ra{_}*{BeKL6>RI((w>v^(A@tO6IoisB!5=0p0NQ; z#qBJcsC%9by*QZZm?R+ld?wXQm;)CLcU~sIgx|cSiGMsH$x2!6tyRM0%T}@z9O#q$ z3QFt$gtiTtir_j&rX=sgdRo1~tu1dbxp^ck6pqoB_H>eR(;540RzNmxCgdu%P+d_j z4NDn787ng&e)A*tcKZbTy^q?Rq63)W+clsK}4oTtUngVB%4Wm``asEC3lO5yT$JLK;FjNi&V(5AkPi7wY?L%WY;HlCwd z|2K!(fbKl1Z&l!MS8ptA>W85!vHW=_b1RjvH(A?)>A1e-8`B(ilI7gr#`G(naw+pe zsjfs6;W`#Q)*kH_%RPz`&vP+S7H5@evs4A#_YMTF|g_?X6;SF`~*wH zR34?rYvmxee;&-{o}~}p^yuLnP}z$yRQE`WYOki!YxT`=)IA5C@@(4Vw}Csadyuw$ zxlbc@W^><1ouf@MHgwD29(3mQgVd2rQ2YLO4v7j9aFdOKNP3I8rO8}%N<6TUl=p6- zCF?9`*yiJ;Epv{u@2@S8YrA9Vbh}6pvf%{hJI0*`N=1{i#(M7OODDnh;y_XluHg1h zcrCb6qK~d-PtezY06emK!uj|Y>~&axc#*r0zU5AP5HY4V;J2Ls_$odwCg(aroW=o*(+&rc+ zE1t3Je>8nRQ1naas+P zEIdXR?gSvhkg<{5eb|`b;Y@707Bl|*6Q64aFomu5ID29?&feb%J%hIp3k`r8%Ya61 zCjES1Oke#PE1+diuTy@t?P1)ov*qJ$0B#^y`M% z@kX$JQOwd-U$N@(RE6*HX*vul2f%Dc02bcOU}7@kk!UasjR8HeU*;?`@1Mf5zQ1NG zi+Ur);{hB5PcXg92;AIvmnl8Jz`Pz=Gm(8~*y1&-QPG&qpmrDscySao7D4cN2<lp$?@HjdUjfP79N)IS`iwY$Dh ze)uaAUB6Cn%zYyixM-^C8PiY>`rm|L?MwEQ|(TJ)Lv+g6f$*(;Kg zvZ7IsQ)r{VHRL5S(6d?sbG!BX<&Y!IU%1MML0ryuKd9%ng2J-si1rk!L^D0(T4;Cf zeA=5)Nqye<(158`WOedfZ9A!^Dy=GC5^Y+)k}=*YXIK+TCC`J)U*l`GRRm{=ml_KNx8wvzaFXnZclCtXtDIW@Nq!JvA>ty!j5@ zzS$i{lY61{PBIJnnZkE{`i|Fls*jgz+*!A~GjLr;2InK>uxs#FjIa3u({P^3lRr`A zXaSi$e##j<5D@2<)vik^h*La_l10^cyuB|>Gtyx52yji4M(;PDDD%oX+FARWx^4qy zoE`+VxLH&_Vk6S}?x6dsWxuAqKH)dGN3UyXb;2+@x^xP9T_-qgRfT2S3s`gwfWur9 z8fg+kP0cqk#jKb#vIf$?Ys#E_f-Kj4#zsM2WP-rt>pSlFy=~+uvW}iC+$=5a-XZQeDLWy~vK1&xo``*R=IwIOA^%E$ZcS|2<1G+UgXJXVT- ze<7If(Pu3xf1W{Tr8NXjE3x(KVbq;0VoJ+P*z6x)@#M}nq?cP@$AYeif33#Hsjc9H zw6ypgxppkr?lU#}njqFr7V8`xS(p4s>{gRv#t#lN^}bE4dkfFh3^ z6|{3d2d7ij*zP)yR~}f-UdBgb-K;=dTht|;*=UwzGI9k`MHP>kEZOl*XXiZks?or)8pBXNP0#* z`O6-l%lB_lv|$iRdk3IjTk9`}%vhcTt#(m&Q@( z!H9`UE!zzg;U_2K5JM^lfA42Wc;KrH?pHk>er zdtw{L95iI|{c~}o_i?5$>50zL*W+zRlzvYb6IHReBf ze`N=bJ0kz-R~k|nNx_*5aj-oV?tDG~&x;#b7qv@h9(4`-OtNtCyao)6Z0YTiixgD7 z92!HEpts;M+k9gXTM`|~EWF$GV;Zpa{Vr@?*Mq9()X}lvHafDzgp`Z=k^ZfCa4CRAGbSrlnEadE|Zt+cO+PZ=UOcN)`Q~hcC zg|o<+uYvyA&Tu@)emP`Xi6d6es6p4nxghywTG4Mo`+_GGwn{5;>2O!nF810-nzL^Kxt{+{JNDm1k1QV;k3S8& z)+%U6HX;6m1qz@1#1*B*2sb!~1Go3!(5p3UXz_J+bCogwR__*Tm?H(hO+BI0sD|2I`Jio6rIP@hITy@Iaff}d`sY3 zyBZH=Ls@m{CRQV6%apF_py>N)q&(~f!!f0FqID4!ex5?F*DfKaS$2em3bKiM2&rMF zh@WXik8+Oua>&Yq?kIP3flu&z8tHKlelBlBvC~<*T-{9B` zCZ+}lBfnsA#75K{slk4$MQpTW2@CoBg@2cSksq{TC{FHui8zPpsEcsIu$Mt{(L#kt6ES4FB12NOFw#o>1+YE5g{umxN7vZGyRhY%fv+#>HENhM- zwLb4h@%{X{jiXxW-opq^&qY8VcKR}ncXgS{2J0$(&(4IhM?SQrG$8-_DJf~3;&e2wQZCy{B_XERlXV0xwpF;Y?>T$( z!jsLCbHSD5YdEcX4bvBS)BVqvsaEbM6`XUSty_887+Og3=X5D&aw`;W*kiV40rpHh z_)FBIS2!T8OD&8|_F&P6NK(Fg8hzJ_vl+(cS+L9(9LF-&U$X?iP40=11XqN|n8GF< zbo5aH1=!x;{FQPD`)631x_Z-&Vrj}6eg&qbYv4G01$>-sU_R0U@`?=BqpFd*dmM61 zY?0w3ho=h+nb;TNVA2c;h&OIV!w~`Pn6i>w=K6DrGR1UFVg+A#^E1EmnkTj%4!}*{zK9Go z!fNZ45N}PTnp0P3NJ9j7uJi#HbukGB`g`EIP8I7c1|rFQIIq4#jlEG}%;)VBZ1@_6 z@&|@s>-DHL>ofTsiKGt|(lq6YJv|dQpw~aUQnP#v)%#q7MpOZ60#tuFnIxab#7sG?%=4kzE1`5JQIu50(z((} ztI6v0INFmio*WNl(+QVrSb97e&pw*r=pYU^)uNf!tNu(aLynnz6=7=kb&$C_2CMH% zvECOh^0i;3*|;v2FkhIA6NN3%eiw_ur&S^EeHkO90tU$cE}Iov}$1F?9zU* z9jC!L@2JoDW=e4>guYQSeVh`C^FuAz!vo>$j+Y`+=r#*!y2nvm^b>xcts$~Qmy7P( zn-1E%gPP|@O79a$kF)Ba7<3Z)tHQbBCw}y4%2dSnZ}{a9tGhPzLS`0|{Oxzg z4_mtLmWCdlX>856BW%XtANZUoh8vEfuzuE2tWwOP`k(2v!ncyl>tv|+*m|0w9ZO!@ zmT;*aYbf!-7pfn<2jNwrur*gi#l|4kYef@tbk|~0f&iuxy@=U0XyHtv10P#1;71!& zA)@XbRMt$##1HlKX~j#(Mt#F*^*k)ms|OeV2qt?@L(3%!i!x5rGusgu#v0Kt_!th0 zMbgY#F>vQrg8)yv*{xHC?=0JkL~DY)l+ER(ZUWp zEs8ql$-SQbkwkW<(imAM);dLnpWbT)+co_j3`=A&VCWqvev_b0U+0olgAd&d+s!4J zc$3(nL(udKL>I3~NGYF*tGbcAUb-n$`xb%lK!VoGB1(T8Ognb#)5QIQ=*zip^kIV; zBJWgT)OG<3M>05j?1SN&Jje?QxoaO3NUv!2FNYNG%dmJ9*NrLn-^F^5vqi7igAjX_ zfcAwus|tS3GSd^7OjQ|fJD6a5q$;M==FW9ix@;GxCSH;lrcn26R^T4KIQLN|mCN^4e9-C+~mqpd9vkz)keE*FTv7u=@yd|H* za*7vJV_caT#f0G+;DF1_{hiO9c;dcu8nh&>2^WjyR z2F+tfFnMwwWoZTQ;a%SF#$(T-waYV{e0DX)s{z@Y@tZX%WZa@hKPw>bT!5&x*Kqda zb9i<4pzE)u(#r#z!P7pftS+OraCcZPvPWK%Hs7*md{*rkFOiig>N`>j^z%G zAj@yWHhL)VlUmD|%8M$L+Q}g8fG-9Q=!5RDdnr?VJ3UT%PJQLN(pvu^82T!Svi8Jd zlEHkoaoco0uip_C9#(BVRozex@p{|55FsAT@0;>ZXAiPkn7iCcbg7zEZO$gErNdkJE4OFy zLu2!KZs$InN-Ux36)TYbaRE%8_rb+;(%7jeie3%^%4}Xw0lnfW&UqD8sprDrnit!k zl){hPC&Bh5`mpp*rYvN{Q#P=CHWMGmqe`+5wN005zp!el;Cz1yww?t^`O^^HSB+j5 zEzpv9k9D!0!L%LYu={=~8@&F-FNbIh>PAn;ZH2kj7JBW!hOIL;VAnF1@bk;sc<(bJ zDAV_#d!Y^J^CSV15)pJHdOS%;WC&UUfW!k(4gSx*kwt{ECA(pPYPkEA{zl2(-n8BC7!{sf1eKA)t(q4-x7yby zkH6|a45t!~V8fBcsI;DeTMa&_NO!}IL8kQW>|{DRejfSn9#4br4#(U_c}zt4C6k`E z18vh)nAXB|?A&WJ)-57|_1Cccb+`A&kXmY8bO?5Fg>)_V5S2|Qe$?s=KJM!YT(EIK zO3yZWrK1jm>nBN+1 z^1@t!<739vxm_k-$!~OevppRcGzI*aC+K-c6q?DCsVrZEYo7j@#%gW149{2NUeAc3 zp0iE49s5MMqf*`^F?%f+wyrlVvYd$u_cDI#gAgRF%ZAZ39gKO_pAu`U(dA(Qt@vm?EZ~U! zuw@S!Z@laSeZRXF;?n(~;OxWg`KU-WCabwh&+*)>cCCYSB1JFnOZ&ZxFzcu$i#5t- zG5S~VqW>kPqoT;vLJlKL+6=QU2mgwYf#W=&lq!iHJw)jAzB-7XF1PBR6T;u?KA!cc z*2mS#S?Jf-5jHVj;Zv?b$Dd55Q$hXOV!5KB z^GMDu4q{(tz-RL!8f?0RioPC4+t*gM^xb7PvDKAzb!bD%;;$Ip?F`J_ThV1f1FvRm z#jgs83xwe6J8^)6VgjgKp{oyROPdx3>ib|;)iH#I{)VL=CSm8XDYSmLEA;Y5k@eD}*!&?L2Ai*- z;6^JPHa^5C?^bNgZ($;P)v!sanAsMO;_03{-yM7S)y_@)v61@x^S1-}>SN`6;`BoP z+=HQ3)SzlLS3X#s{vPeun%dyJHV;BccI>|IPUoKB<@Gx3CK=0N7gGR4E}r=#^2sk*{M{TqamQ) zlZH@^{sz)7kmn>fSyN1wC5fyZ#`V~$N?PUBG_`jSqANCGa%vh9cu{D}`_Psf?UN^v1R!6QvboD(nef9|YRF9+H?_WT>WjvI1-=Wu! zRLpgXL2vCz^l`Ziwfdie&i9jWc&3VSHX4#D%1m6|jM+>uWR?5-@oEv@`P{See17mS zzD#@+f7)XMU;9#uzq%}z&t9{XUw_$z_cn^BF9QU)Sxue*ybZhvj)YDJLx0LZ5tf9=2FY?Jj(k?yHGCNCo z5s`vNmaW_s^Xr`NxINIHvV+2-TDd(-cXF$4ui|PaG;39K%!rwJRbl7&LP;K=n-!L|eGoWVI$a@TPnlj@7K zcG+Osl@LPL&sNioyVZ2}Rh{6}g1hLc@sXljjcCJ@csPfovyfO3_U6elUN-kLuYT_u zA16}9r}sV2?-JAHPmHnR8@j9Uca~@K)(^h2opGs5#^wU=vObBg+}abz3$HT?&ksx| z%!f^JyUUi!+A`yw1*}g9W0`8sC}9_H#?lQ5vIL1nROlVj1{^4 zWb3U@F`2ftl>B-tX6S7|adtY%cb`rf{5J077jIg%(S=fD%#f%4m2!LcL)!kau-$4* z_Uqgb^5QZbJy?R^*%C;umZnWlH`D6Pr)gu@PWW`$i4tW8Jn4Q0i4P3~;i6JBYqkw_ zb$>zoYv)q^^)6_QpURBBUPjiA30SCk>6fT0{^(2f=Y8mU{~F|`68|*tyH&S7dzsPa z`}lrH3zB0#!X-8ZiwY!Rsr3P7$9{wS@QhnKW~RVbgIR8sUoLQ-ch9o*c>}SK zI9e69o!&TVU}O3Wc#df!{|t4C?w=|MKl_2Z-D4{0>#m^-`QDhByb%{o#v{&2A4fiR z$Ic<5NNnYp_-qSYGJDQ!ZVbllqm4|`zmXYPn=@k%M?R{wk&g)dh$Bs{xVy{~7t`|b zF1rC`kKZ$?2gRsS`;ODq>#?eB35DG)KwTeQCNHs;Nww_3OXEg%$Lj@4s@l#rZCK39 z&3wal%geJO*+iD(@|hhUl+Jp#r?!!skJCGQ28m0fP#kg`Q3>mDTx1!xolC@~E{O# z+|gQ?23$qK{NJ#4Vh_xJJCZKd9l|6*Bt`4Aa9z*3a{95o1o+P#{V7Y-VvM4$k4y^~*VldG8>mw1u9%JB3AB;fPh( zjX7(kBCO9FOjmYiTIbF)DN7^lay5m=?zf233u2mAO4*OXa=y3m9qe8w0@3z2s+c(; z_+tzri!VW0B^cM1m}2Z19h9lcVBLn@8vFu5=pmEf2^jHiurO-NBMmHCS%yh!wit)ROm#YyPgyrA_}p zYu}!yx4sjcJim%_7;#51r+y4wDBOa` zx_k8Oelx}`yMlGA4?%qVH;ix4#`5H=gl}sJKk~TD=p1g^JqPZCoj-})x1o(9kI32c zE1ekg9)rK%NA~qNdUBvQVkIARopCrX9Pyn6~} zuM)?k?&`3abcQ}X=mWJ+&nT_&E2(SVLgIp}n4T0w@e&pkQ_#wV=IB%GVlNU|yqwhB z4Jgd^AbD(0q2;~qP}0{;$a(Mi%OP7OLg-=p96a6i7=LOh9 z5I0^1nQC3JWb9jbMIS*({ZTj#x2G1b`P}Zd;oSIX2L<7mqPbCbZ&FLPHvRVOCK9Df z5&h*mqE>~_7ujIwYd(j^nqa)~`%G=iyHm2kR4A^prZZdQXkw)d_bk?ij72|lBeVu_ z6IAjkW5sKVwtPm15AViE69r6HiNgcCEe!U>%+1SyOR~!^hsc>JP{K4xa4Y?-AptH`J@Zm z#71meu@}jvc8D&LLAs9_-IK4UtuG&NP4c5S(`(L_c4-*`hbt~L=Mlq*O@k@@%zfNY z9)xU-W+Z=`g&>U{NRnE_78q^CNS_eQXSwiyv6W`^Pv;Kulev3J$|RPREVy^hP4G<*6`eo&meRj#n|y7i60xxwuOrQ0#w{IVMCg>~bRQ%LOz+<34wAY!kQ9C7WAzB#xW6_B+WJ z9H*&TZ^30uA<9oV-XY(&L{u$<0mbcg*W0SnEE}2Kx2K( zF`0+!52~T9xDTsM4X9r07Cnv5Bhx+Ixi@~-x$K@5f*Z+41XJd{<%;t2sdASCd>qdr zEn+be8XmPTEE_RlT_#Gmeq#~tZA?t540ER}r#+8mQ>*7DD3mL3#c6ePbw*!`n%14m z&OKo9zNwplg*8&Qi&hXSj7H5=|$juqo4rqj5+UoxN*L-qljj30zC6 z*8@o;%o<@!2zOIIv&Cbx*hN(N{Z{U4_DOx^YDwAUI_R?j9-jYJG<`t$nvXC7# z@?}eMT@i9*Fda;G;8l7wS;%cOI9LxcGz1= zJb#QkVKbad(`e%!yKE82=Juj~Vcx$SvdPqh(w^4y`xOoO2bb-beBv!;-v23+E;7SO z_s5W#q5)5LNBBmV!PiLzGG0XQ^;&4>*~!$$Rg)_acynT^^SG9w8gzD}6Lhom=vKvT z3~)(7s_71_^J;^~h6?O04`un=SF=#nCQR&o3AyK5!R;>uuWyY$JL8}~ho_lM%2doh zBDFz5BVE75dot8($x6^1E9)TyX7Q)D5AO>7-qKFxebZyll zy5ln!v*qpB5`!Pi{kKExwqpk0qiZZ%-)}UNe;veBH+^LOxhL`FRu9O&QD;M=4cQE_ ztxRd=7WmbD!07@vG%hOQlZH;@tJR-d8GI3^(3=+p5sho8i@{W?@39%pO-+ci)J1XF zZhCm^4zlf2V0>ykEAvQYIjO4HR<@t^TYFGdst8we){(R=wo$ypV3PXm$gL5zp;^vP zVX<`r$u3?d_@K(Ty*^90Pp_B$nq_S`zl>}q=JWeoy7BqD@3Zco&N9m#+nGseGnPI( zg~CaCs52af>T+GAZnMVHHy1FX>=wFS9L9ZG@s%ro-GiGeQX+V*Hk9^6t)v&jbP&B| z8|z|a2k#pT5I18XS}x`yyP_UTScO|g^-lHRbqtCv`DbqMPfGYq#M+A*ouia4y@$W&K|x34*FVxR02 zL{5E*!k|QE9axKmAKjTup(LAG_M6q@;xsI!n8*oVsC7!*!%ohXiPba{+GAYb<;DnFJu%%?#V&I-;pvxMv}stIW%guXhO3UUuU;eygHaSrWk9mBKFe?>^9!Bxn-YvQBM9Oh-ah%m8vYB?Mp^O)_wDsDA1Lf1;UdkD#-58`+~JMJub2=w!p zV87~cmbrNia~u)HZngg8qeuI)u*C21DYwC>qpr~2!_k(n;*gjsOPNo$()SC?In|nr z^t5pZl`b}fl30H@4eCt=uU^v0y@A}j;>W)nvfg0{ziW*bzg6QR%A`+YXJt7qeX_&C z-(I3F*$dNFk3#v*P#ofWpt`OP65FSEeSU1G16Sv9=Z~J`!hfFUVlTLIGfZ5mE@lFf zf+NxBq{bvi1fe44J?pLziJxCqFy~{#S-i+&UO&^3m5!Rg#;X~zzQYGGSr!PlnoC$P z`X%_|7U;LU3zP?|a)*X*C2`Xz8t(O(^9zjR9(z?#L+~?tuN6jt{?RmakTTcs(Uy|+ z^=M=42iUruLi_+Dv>pn9tH=Zq2yHIXEaF@P8aVgebVCY~ z1oL5C^bJa8+b~Lx!}D8p$QpNt_B@TI8oUZg_hcwT*@(B+gzlN&ccetz1ZNQyDTX3CMz~}Vat!S z;Y3njig+@aqK0l{s@~n0jBN^nw04uMMh)edh#{%Ok(X0?#2erLOfTv)Xjb?vPHEf* z?uGJl3eK9#584*Or)e!_8pGc*v$hhZ8#5C|>(n8|-vmGVG4e-LVfVo`IQaG!He9pE zg@bdjApZ*MbMX`s)LuLih81K?L+*S;>|3&Nw9*xquH_< zlUc8uElkn04SQW=(Xe?q28=vL%_0*Zm2rl zwHF#{TVQeKAC%-lC$EN%XFn|9O5yM|#XwXeG#a(B5qb@W1@p)|WNiAkKYdag; zo-#9RpVRBfBIJ5-3;2Yc@NaH`fwL3(FLZ>pSs`1zR-8AG$m3^i%Hz_UcUe@`c?fg{ z-Q`>z{h9oIMSjaNTdQ6JeVNhG0mz<~3D=qLFnPpd^k|-h#+!R_qi_V)o!o=mfCpG& zriLc&BUAW%2QF=4c;53Q>z=+Ae$VgyheKFxeljaN=)z0xdBx5&HnFo;MA_x$uPl4T z88&Eg4qi>y#3`)ATYXK~=q#gqw^@4lr>kPYxr>XIqFO)XTLE{-4 zHqQSx+pfHXjr?sDlW}igIt^;aH?y(+C zKiQy>=UAWXCfGSAi4M;7fXAm{5bK(-Z*3_SsFId9Y^UV?i@Aa*xT zgWFbTI@@(0PLJboa6HchFTs4LHejk^@h^vzO^V{Bl;5x>MJ1M$(jRZ`cVosMPBF_L zbxfyeJd-s}V?7&XacG$bmT7)M{~SfA&fm(}+)N|Kf)(7U7&WT?)E!;M2cfi}m}xvd ziAPmoxK-ea1s_Kv(<_CAGDlW$PJ_87`?8c%O1z!Y4c>V4BW9^Jk%dc)WuqpDGsm*G zsEv2Q_VSU`JWU!@9*UIs;gED$PwV7waxF8ZxTqjE&dEKW|L^nw0EgArHIl-itL1h5S5WtMzU9A;B&7 zLhGBBDCu9!E$(r_p`QEDWKSnPemsq)Xzrr4c|rL1;2L~?D+>zSm9Q;h1D16! zMTZsfXq;M&`)>MT(+o?DueZkReoriWG#RBoPK5eFXGkZ&oSaAR!e1wy3B%su;P3xF zX`>n9akR;!kaQ(>kaFk`a#^yQ7G6rCq3{MbhZc}rfga7h*?>)Q{p4<~fZLZ9Voup$ zJaGL!L{+=;3*|PjBfF1+%$w&x+pJ(bdxDaSq>*3$1S4@Arul!w@a_X7cIgbsX#K(8 zsU^apyp%9z@;O>IQjY4P^eBEI4lHaP?f@@9(7%44i`grWp=Zzp9QN`d_>v7IzD@;? z=xyNZdY8H1Y~epFLO!Y%Is6SBN)$vkG4!U3RLDRNZSh%qVayq`j^|(j4cE1g&-oJs*Ep$lM?=fm^ znf*T?`20{fxx7D5!w)6l!z1zN-tU1ErG}tdNHe;LIpO5Dh^~W6aA@NTD25Xd>hTyQ zQd1%PV>ZZ-xXSN#{=@Y*iZMS(fWW!^aQ`SlQ;PzeX^w#U8(z51+ZW9?PsdEtR=j#x z5%Wy1(D~+N!lI?}^e1sP#rAxrl%1wzre%S*U+6*cAV+i)d`961NpuVUgtlTg;ZgV& z__8SjHhSM+kK=wYW2VGc&+}k*Su5DeLnh4FelkqdS-`vuYQR}~46M*t$v0TXvaH5f zHreeZ*chLKq_S#ozcCUXpLE7mVLu0=x1EN6I!tDTlSu6HH{2c4f*W@KB>5$!R4{)& zS-FhG!WcI+XpuzDT9n2&m)ql$Ec-}N9ZiJIK1Tm`WNIu zOY}5!S8k_a>8enBKoN{BmZQu|b)5E73f5hikG$ky99i)QSC5Y%H@9PWeWot|aR1VO z8M08#fW%t2;w3>g)}+k92M4XN#QQd$yzw4A6*AD`?{dse3PGLu+t62|7p=6%!NKOq ztncbaR%dC>yDyO7b631$*8p zJ%F#olX38!J($>|h_9{&*>7L6!rrPj1y-Ke%6gX0VwuhUa1>nF5JzijlX4Ym=)R(j zb9%@uBo|qi2I#B0L&VieJT>ec#*18p>i2V4rKAel?>+@Dttzo7|2Tfp&0{Cq#Qp(! zDWXp*;Gg7IBYtat$R6NWSH z{A$Ku{LGA_H!yzREVwvb2+7O~s-9S(g6>$jkWzty{c{GKO9;liFvOQjeDG`5G_0{7 zf$wVW3Ew0h5IUb2gL?wsqjQ@yhR(IeIO`mo-s}QB^B>{R&Q#3(UWRv4?x6LKe6-U( z2-mi_z}c`)v>WS)+)+amGf#vSJAOjO{9b(GAb$n9$1S)qiG}}pF^*ID~-kFJ`-v7p5(gL=`)aV}&kHmvy zEx1jmNS$^cjUbz6;zYLJNMu_DKHv-CHZOu(?E^3@eizEOPK0a8Y9MugDU|3GvLdrs z*!9VVm3W+j@274+t>s@-KJf)Ou_`pUJsQ5-ZA1z4d$1?O3TJfg#X6MB`3*j0XJFif0#?&l1j#jTpiWl;RR$f#UB@)>kozf& zH9CzSuLjUmw{_V0bv(swy+IT0Z;|oteNTN@d0-#iRGUdbRsH0VynxoP zdxqzhAENlp7lgCh&kLu2JH@8GI?ooUIeZ+-z5vz?f%y9nJS6L5|D6x6w4iD_n)kflEgp4cp54}Du<_tTNk zc4R+}KW>E@hvTrSsS0PSaPY8pKbmF_fvcVp=(M!~U9b^qk)eiOBKq3^mhRtRKWNnz z;f;e?v~2$=eChg<3?vuRkim9%E>#zLSG-15Rs;>lOmSD-4IFfGJ|2xI!sPH$?D}{d zPaJuGxf?uimGn%E8*-DC1X#e&{M~50NgK^ly!n8D!7So=7`QDd1Izr`u-rNm&dr?$ z+Wcy_+8JNUi)53KaRXY;j z9-2jL=0tl<_Zs2DcU$;IU1LE;gB!@CR)bsf73jXNB~)zoqqD0w!f2Cz)~b}lX>MN* zJHG`W-*X3#Sm!zsZwXe+R!@upJ$_OgbUx4^M0V`z$&Elrx| zO6!bsu`{THoZSv!`>Z{GTi6L1QjivyZY^+OugS&cm;?CzPn0xydu3RGxS0`$~4clzop&5wV_D#ns!=>nG z`UJN=+KgK3wb5!OkMJc2tp^ul`uA_>x#=lefBK=VW9U&nR$Yu2duRlWPhwCsy9u*w zT+w@{3??jnjZLPOBym6sn`W0$as6GoM0v_z%IK=cfGo*U7Hf2&Ox19>ur!QX;X}QY{f@zjnro{;gP-c@d-Zj^aTMvCByr6YrYi`Wf*!+PoHRn;&DQRR(rw{uMT! zJ}t}|0J5NRE-Al#M&fF7P_!r(v#rwc`0OM2E(p(;v|wUd2c};* zi>pG1;Et~{@L)$eS}2R*z4Bfxk==_M4|uS=nw!ja9uE=)uTb(~JSJ{&#TAmWXmkAo z?x}i+arfl$*C!v^5*JNJtK+F#5KZT6pHRitVzTX%qsQH!h0bqkDI-ap94-`+*h7Ds zruUb2EFEv}Gt1N7?#ovkv_zUYbcXY97FY8<9;z@od>tLI@1bjFH^J(tBu_OS=szGs_!+3F?MIs(Zc>@YPwGml zpr~EjNX1tJFAs=n_JkbNY&Z>Pdqr?~$7iS-@a6n1Z+NEb0mBZ7vKM3si4TKObnQM2 zyj6~cZ-I>3gtXW~3cr9H@^&2FQt7~J!*kGAs~<%_?t+4a4w&us2H{S>&@S(3XdLbjca|+?3!dG9S1~KmVub~E2k*fhVt?U7gbxP$ zC1L5;EhMe1N=nba(YEOAv{^ij6qAqRtf1#)JU^dy|LrH2j~Q6wv>UQqZ?XM0!XX`-S(}GH%YHi#e1q=^RTi|-XCF)4oNkKQma1@f-92o} z-A4+N{kZt`V3zDJ#TGf4foEkrzvSFR<{~cDgKSa z*QZBfYr{xfk#iEqUe&-60pW0Emntf2Y=ZXZg*f!^9?;G`28C&l*^5&}@Qyoz=`&Su z+|w#D)c->JB-hdSp4nJ*(x1fVCsFvCog~xy0EbNbgF9dE2iH$6a6P!0CN0=#Z?=0O zzEmp0r*B+o`YJt=D}IaFj`wI-z7Y4`dy3DW-NBdTH%U8fE546^gvCSJNbc`>d{HqD zkKVRJU8!@Zp0|k9*ENufek{(q{{h-lXX1!Qjj%3j0;EJH;rd(-=jd0WLXr_w(kxV# zxQ)psw{X_F^H|h!6+e8hA=&jUSX}Om*|lkqul<*lY<^>nSR*Pc>B5?U@Y=U@1*}Wx zf?dAdwDCtCbsk=5FD5T-U$X|n%dPfMzTBg&69wVBw>k`cm&ZPsJs|OsI zjWDLZl$=DHU`E$R{%y22>;HU-zn8ubG$I=Q0WsLzf#w?`$ZT3VUY^!Yqt-{z_|?O) zAEPm;;y$`vd5H$f2;EOlp`_PY95tgAHy->6$u26;Cp-n`A8dlTVq;++e;sabD?PKV& zC>f8Q=oDV+%@yu?xD!`tnnRXYFg!e3yu^HzHtF?L(M*dm-A4nkKAv?-yQ5UPLZG zCX#EeBmT%J$26;NB)?XchVoG~XT&!ew$+0ILk5wF{5QPhFqV}5Sdh+eaop{u46RdV z;3+wAyllLR;^tXl>ts<}du}8udF^3O=ADIA!Y%M3`#j8t%eXJ{JzR1*%`Pvv2c0E3 zXef6A^;56ocCQVztuBl-VrJu}pT?xFds4{V7%a@s?xKgcrqkqTMG9N;mS)5i(w($m zq1ntyl&!Lo6a=erT9PP==6)u5i?7&eA3@vAM$_wQ68PX@9jw)pWMk}7*vQSuq?T`D zgVjq|X?G21$e6RZ)?ci?`w%E6{Qe)1|Grm`Zz#t5egV|H+M5Kb=V-~oHoWv!3oW*v zfg>01z^SQw@!;c4G=57M4cR>wO(vD2;E^#}Z}l1Q_X2Es)`6Ypeq#5GdW;_ANg^po=EB+CyqvAs^NGEs)ZWi3d`}gB8 zvU4*TN@`JU{1r0o8jk#nBsl8t4L^2$hI^-Y(4PmFNFr1fipiq7d^+*=6mMp9P0kI88UCl@Q&7eJ+oQl+s;c)$;^Ht8xgLq_N9W;Q?iDo@ zKBw&i9dY=&uaxh898YnHxKlwH2M^tj`HCBH!?L;Ta)~-CH{sD#;tin=O8Rhob3)R=ALkuVE?xE`gGhfK(@q8{y&lSx|c3-*jk zpmD6@C6&cgG58B_ShNQ= zTFpbj~- z{(2G{Cm$fKM|<%^*gD*PT@00WzQV|~R8)CzA8&@u!3zuAN!olIwwxSFh7Wep{JS~W z*Et21h7PhC~hYPut$gRByN7D_!S+N%a>=YnIK^HeO>O_mVvo#WH@$^0HF6pG+`T2t@QA^CvuCgp^JPGm z($J-N8mHP#KyB~We4_oo?RSBkb}$ze8qaJG&*vroGw-Bj6`nevjNVrcL)43TD0lfT ze2~{9rLJyryyt_fr6X~iP70Ju$>7Sh61e~3bX+s{HU6DyPpWc}q$ga2+wZl~wx^p= z^Ya*L4Bd+J153edW(04PVajq>{$!b7l7;+`7U6rh zi>N;TE6x)&$G@LnVwR*k*5X8p`PM^%H655ZBM=|zJtxbwZhW{s86%XAVC?quxGnc2 z1-!gTp+m3Y-K$Ed_p}Oaw0Gj^L_Hke-H%>p#39~&GD~}7%A1S}=S4QPuvxSJ!r0vx z;8#`(8r*P0gR=vg^&}q-${9n(B@uqZmkzKLxrFv)O=0x3?du>JF796s0y z&6ZoCdh0}RiaQKbBbwL(-Tln*RV9cOUSc1nzhkc#_cPD(zcBQtBy0FOhYfq1$%kE2 z7p!z1#IB~bu>`3aCZ+e2Td?@b|A5$)6!dWBlq9I0R8L~1mX!Pa6MAGu+$ayf3{LQv-CUAX#Y5VZB(#PddHu`9<9uWvq1QGR`RXZLJe{^c^e z;?ux2U8&t^n%4wQ;u$lCy*^r<*3-_JQ!~+v{kd~G?4a(Yt97(ChMj}zzB>NI zwp*;Ox07%8&gZu**vKn?QezPTmzZ|`9k?X+1j-+Tz%9DMMtN@IYU@|9!5@p5t%e&D zowb%#WPai=yjN$>6(_S*J68THS!%k?Ttt5>SNg37<8!LW?Z;h=-+m93HmTu~v!S^3 zO&s~D-4|NOhtk!|K)f*11b)Puqk4KF-1udOAKY3=_0Mc9kvV~H#DHA#qoL@x4lf%g z#=TXov+>~mH{ssSwDubwjk`h#p@@9Q`+=xZU{$GuKz-3AgzOAdcBbUFdL79^lj~8-Z!x)f2a(wc zNw^++0>#RYqHm)LoK;Olze-JX*u5J?;-pdK-V`=*n=Vvzu0XFlG8hnLh;332NVViE zsX3=&c&@UVgdqtcN@33^b z_<9Rol|`&vF%EAZ5u(Z8R5YrbfnT-H;OXTP$-2N6mCCO26;H%DrDx@Wf-D!i!nlwB zJBLDIEaKT&)zl=Efx-_P!1anb-y>ZQ&9aKnrz8)ZP)ix@?}Tq zB^~M@_@z(E{f2ZiY8++z{3MZ2SMhb$Sv<7R2IXwmqTi$|l%BAkq;vkFcK2b_{p5zD zTAk4P@kv~7xe@)tF2YlzZk(L31Ls6k!}y>}c)({SDJ&?$it;YZpS2K@MvetII2>KK zKg1=*S5SLH8puj-K#lcVP+L6?Dt;V;F27&6eDNh{bAEz$!#yGP_*&Lm6AaqZWmx{* zXs&6qHPbjI0qtgaaCq8Da6WJpS`6QSndxWN_|k*rokgyrk2dWq z{_}kP-}rGOS&`~tA|jk}cNurbVzOYB_bNfaNN+~RPoe#)12kSMmx}TPn0DEiHj2%m z5rG4()6h3|h??NRpfEPiLSrtN}r0=sMCw?#O;l^=1Q zlO}bxKEg2@UV_W<{hd#*62M`Z`b*%YI$#AW=tC5ZNI zH>K#1cJw{(it0;a;OC$Ke0h5ihNUba`3pDj>7LnRur&b0#>(Qtc|q86z4t#H z#6&c@14U;4KOB;>BEbVV1RWbjardkr@rD-f1kZc2`BC%4Aa(vCNZe}xYR#T76!<)ix7FADBm&~@jj+pn zK5lul90xu7#^+ylV6Q*>qT%nyxO7A?)VkK8-772Muw&`MS2O2fLq`+^{Olv-8)F~`cyIwkt4I6&&a)@h13qn;H+Q?h;IW{o3o!gIrKQw z(wPdn)5gJO&o1ygtqn^?_Jhx)v1II1MSFZK_)q1-`L_voKz_|dc2i7`jJ3!7Ulghn z22lLpQl(S!LLd-|^NRPMab_nQ>?-Ud1e#sQC%TFWlE&I{ZIb%DRMczES?WGZZf!#U zOTnZj-GHU;xA1WBGZt5s!H?NqE>Kw3ZkM-bvt3@Xxvik0&^EkN+WP)do&Rda!PeV$ ztsx)S*TE;STYfdX5(#H#jyChNM~;D2S{A4=N)ki*51`~{7-$w{@l?nxd!3!#H0Jn3 z%$OZWqA?>dJ6#gLMOoA8b4xL|rwzw%T8{fOWl>Y43j?PWlX^q2P&O$TzduXFyNef* z%*}ifJzPi{Cc$J}dJ*ds?ctW~H(b_bi@UnE;@4tnJo#-59=e~0@d33ce`X__X+DlQ z+h0Vn$QYbk-a^V-zK~JV5;E)CLIFFvX``z<$yNJe3GOEC`Oitlex%UUA)J2u>*CV> zL?(SbfNz%?##>}{Fka^f_}y)XfMas3N^u6vP(Dog$2_Q{T80mQX2vHp0b8+R4Oo|? zvtOsQ|EC=`2kNo~B8tno-1!z<#l5kda%_j)lC>H(`X{Wps`0v9XU`$Q^S0kybcGG4 z?%m5d@z1$X(J9;nJjks$J(la~9Vxh;ukzn(b+z2LGkwz|kop>FrxiNg=7;HRJJHxV zf>~)vFy!qX94a{z)!Nda*6%&Y)^B8!uABsiIV0io#08k=?|@d~r!cBC8RLwz>?Jct z2?yWaL)u-w82ot-`gi;#`}xl(w51DGZ(8E0bVc<0{1=sn%3|=+SJ>3*D%9WkjkJ{d zuxoV-Nyf_Hix3@}R6i5@8hX*GupV8PD#5705BMfI0WI|&LvmIc#B6-W&eqoO+3UKw z7Q>;CW4Z!WzluR)(j!cEe249`K4bB!Fl?SU4X@ZtMs;&T46po-QSp;Va6-#|c=QTk zwAK=myXME5Tx0qAeZRTXfrjweG9y@Iy93JWcR^*&9GJX)G^njU!he^K=e*{mGIL2W zc6YKYXs%iWs^upC?C#U8Lv7JU(N00)0gM~8<8F+FW8e(wBE zaObJ;e5$H2XJtHQ&K1G#rf_n}2%vcJ9^4^gf--xzqR7hyaPOKwBpl7bfNdNW`W>d{ z`h!S=JtwDK5}3U9JO0^gMuR>YVr}{-lv;2UbvNV>h*%j)m-<0w#Q~Ug!3vg~@q#(= z(R^dwApVi2KWtc*4Ucay)SFlaZLRt^S8NvwYmv}Q6;F8T;mI}A!EU@A)t&!Lm!?Zn z;`(vIAx%qB^V1u)GT|_*KUKrE4h?1_7R}_u2CW01$>p%#&I$}BW`arnE&jtCA-CDK zlJQxM{1@2>)~uTdS|zjo0r~d9L?9d&#_f4y!TT)!z=bW%w=-^T;#B9iaWCvWxWlRQ zxX;q7_?U0~yq?H>-tGG~?#-zec2NPxmCB&hY%}mJkAck9ny6RO3v(+*k@B(U z)TR=G3lo3RFp)L5!&n~^S_?^KuQ(>AU56XDb>QWT0+jnQ0@M!3vlU}n+0owlFxWF0 z+s~z7#I0~PXzCMCTVf0MEWX38-wWW*l?M>~dnX1JxshVhM3S- z<%dMv;@jU9u=L2cux(Ne?7Vp3e?T~C*K>j{7-9ELDw*Fhvr=&FD#z(hT*=#h)8jSH ztmVAjy7-FPmE6No5nSc<4nAKwkT-U)7#K__+Aqh{&Z%~ zL7v6NEy9Gz=WKS`D#nl2ha+Rh;o zF=_l-mSLB`r}*|WXYz^i0mF6UQ2y{#JozaQRhF8fuG0$m+yL9_2lSqNPt@-=708eUT}g-PtdcDC#7w^<(+$0((AAbp?A_qYD1NtRT_*By`!P zqP{~lDQg_$wp<8iIy1^(jfpecoYTmh+pPZq>3AUEMgF*R-ZQrGg%-+OA}RAWPc`@@ z9zEQTL(4g_b<=p4cu8*l6)*1N)mrZRkT|Y(?oQv>jk(kiB&B|UX+*Y{d=U?K>e!|KlvY?oPd_>_VhMGDoEo~zNJ!{5! z#XoV@fZHuyRZdAZel$+j6s=7ppdjuKhK_2-iqCQ=UwICuzeA{YvSX*-q(WQ4YB+P? z9Ar3-fSq%mLyd+tY|A!-n>A0+$M+x(p&-bR@WVwO7I6N_K>&jisF>ylPe$fL!N^Gi zy7K}hXZ1otjTlO-e##ZcsVAAhGikm@0I-a<@~*B#=;a;l!OHM2cI{#G5hcy46- z7v19HpRQzHBj3==8x=HuuO1T>wcw-P@B7!GyPR6g54W!1vRCU0#{72XS_(yZ7v+(> z{n#!1oGf+j_S#ka`sp8em4GLF97yx3KFVC`q%FLL)SYWWee`Iu3?w+^%s5_m|jfw(<211&-%PD@)2FOP&`z zV6`a=E{9a0gs>ViLj7>?;x`0`6#b%kiv!j zyvJo4hS**1@v^zNqLNj{D-5`+Bg}NzSvI(|1`2 z(uD^uufsU`7npvz1l=?@qIKvYdw8ZG3|}M;_fnU`-b+6~RbG_s+kYRns%4^~#78u_ zY=uXs&&B4VBQ)N-kNnSN(5E-rlxH>%)2m+NXn6$`tl5am%0{4#-35>xk&B8|KSAb6 z8FYu%?+$jiyIS2b-lj%mj z;juj&7nuS39X;U17DtdSKh3U8SitRkc!FQ$zm+@jAWkrsb7sA+)4=7zYxdVChs~3d zVzs*OS)lwv=J_#OxO&GMn)!MaE5;A(bZM7e`_vVDbK}1mCjyOS1p@#26TI%Ad3;tG z$6Ktq&o95WnUC_^#%Bdy=gvpz3+|`%@fxcp@MWsr%q)5pziZAuyAYLC{IK&^xrAUT z-fYleUL$%1Z_qM>Tls4qKPJ=&;^!#iaF;T8RI3QFde_AoxV)!$dgTmuvf7m=@I+D>u?16?~?x2RnP?UBJf;-xkP|;n3qFemY@|v5yoO6jGAx2kHlbJ)G-ADJ;%aL=YJzL1iE#rxgFJkT(NQPL6ird=x@o><1$|2?m6zq z+(rD}ql5T*txxGejU1+0>Dc51rF-8 zN9HvdXH_ZVjE+3g=ysxl>6Wxe;~nK%?w~u%j?+z}7E(I(6-!g4@l2RIVT%?iY7fJ; zdl$oAkxL+VEFPNTvN0ku2}_>-z@pDC=xgtd&Yv%%{(mM29g@c(kDQ<(rw(p?p9a+v zyx?x-a(ECH42B>4U{Hz|yfhySCJ)}hxC#rjTyPOG4d$_5FC^GUKQ$Qh!H%^hd%@t~ zL7*q1#xB&~V5M2R?d{K%3%hpw0h4PFn2PEpZq_SfCgybm6)mFSzE{M*fdP87koVFH z5Io(X}I{H{iR^Q~U~^auk!QEVL#_XYfZH+#Og*p&0s zIKr*kP|ckitk1`!Y~!ySILIemN@uISsG;XI8yY;W3ynqIGMYUfG!DPP)IG!SbH&iHFR4a7e`h~_Vwa2+b)oJ>tm$Yil zJ~~kS9t(`E1}8UO&G2-U6^FXSn&>HrmoP*J#LGK4Gls(toIbl9WDL=85Yvc zd-T3!(LZiLe3=GZKfVTC5}KI7!*0P_pD&zy>QKH>d=no(`2~NwJA?a~(Zh>m{pDTb zZu85P`Z>!#TY2jh26l%&HS^2--thHbKd{Dq=Ir(wJ2*Lc1FBd#Ls(}bYpU(Yc zkpUB@F>sN5&hNK(0BM&CLb=*_+AF@6?Fww+_~Ef!NuLP+_4@`qC!)Ye`0p?7GL}iih!;!7voed5<+8YT&?^l`Lr0MD9?U67QqAfu9#Xf%j|Fta6QwJSD|pj*b;Jm<3Ck>YSp#2**5uZO#- zk;wTZph&|MOgW!QTJ|wM3mA(j+xNobxfen7OYOhU5JX2! z$C%a6aQ9y;q0z>97=B_0S+C6zt}s!+C+6j_Kk+WJRz1jMk zaG(3Iyo+nR8N`oDILwQ3lLbBcqq&#+s(2NhAnv@{bFyqUq^1oa@V(Xt46eNe-GX{J zUwa!vFN)C!sRuN4raLD5(PJsgJ23vuM4GEyLH?&|F#DbkPRgtU+r1;$zWd)_{RMGx99TeDDi~H9);DOmU$?>WPUWgYVIy4_W_H07GaT$20Qn#9Wc$& zg1^@LkXx^ClKH!q@~&4O@W(a9d2gfpAU$t4G$w7wm?f?hrly9^cKpV&%?k$BIuQ)6 z=VNU0TiANr0-By~#8kU5d?dYt)~%SxuDqPa&s}Z765r2)_L;k}dUQJkdHM0K=dS+C zkPU**`1{UEe5LFp46M975P7yVD7X(V35Dps?kTQ~XvNE8_M-6OE>EOTTrQ;))yJu=|EEE3XE+90I%8ff>}GWQqdu-rVBWES#RGH5iKZAyT$%b#$9 z-XBbUUI)#-8(2ri4X9YSA605v*zx_3cvTVP9z~_X?+;IL^i4OGuYQL$ZmwoEhx+~j zxu!A!E~b1W@1cixuQh}xaTvr@oE8+xtEImyu_gnyREy%pls7^+PVD+x=#1t zUa1zsNu_FdeR?Nsc&!hWBXgi_5`)y{8JMp>8$%6LFw}5BTm|z{LB$a@`~Fbu_02++ zXLE#BRx7EmG#DEsBS>PwXNn50A{w7zzsp5bxR~2dMqNSB6RJy@hbCiQ;0g3^dX3wn zJ<%`E6n2)*h49pP7_!Y5Man)y_gpFPNO=SnbBqNwX6>L5l*PLpKhGM}UHPw$%VEsc zOg8P!CFoB#gO8VGVey&y$m!VQs69K#r06*oz!Quay@OxUvzdvxuH%LHf~6k}VA)TM$k89_ zOdR;2mZ#j2>PdDhc8=lhjz->Eu}AQyZZVhPbd0M?Gv}8(0Xx}n1-_`&(CuY0IIhGJ zDqY`0Xxo0)ox<2&ulrE(CiaVZ3X^tb;D9{PdbNz!DUc%OAR8Ira^U5Jb2WnK-i}QSby{lJ2`SGToDKU ztez&ETD=X{1Mh%4uhmdgk{6J(Nt+0jNK6kp(=k_iq=EBt7fTu!Lae1 zen)UBitpK4jYT-Mv_WVnQNxsHF6L50TiM=i-e{1VfEJOX*>lT9th3_| z|1iA$|2rl6@F1LOx}2^bTS7B$=i!~dhjFmQDY%(ki@UT2llG%Wq_yq|>D+ZBg$pk* z}u z3{LNc^fT@3%gz{37+Qx?6F=crXhhSFW=L5p4IkQ)VNm}W@ZOk$KH`_@^otImlhh`m zway~y-*XQ8J}IEtRZl4D-GYOgGvR5Y5G-x>;k5e$dbdA-Zp3BKe-05iYmp>8@S|Lq zHwW#jR(ab;6yVHouQK&V(9|t|yj8?8!apb*+D8Go<@AJq^FV?Wd&0)BqcQTsvP5(dboO?J` zUHk75g-|Lgq?M#3=fj%gzQ;o8sHBsmgPfL}&!^7H`4}E066Khvbkc#uBasl1N~9>I z5Q+{|ThH_ExA*;H?|tp-`u%>__1oL$pE2f`_Z;`VzV{qsuCdmflX-D;7h}kN$P9dF zhR3oKNq&GUkyIXl_;aqb#yg2a$48K$pYBXV##eIpLN*Hdm4e}&b1;5}h|Y`ZMP&TK z4Zjn@l&`@o_GB6>8bGCf)L@i>8R|7pLXrLE_(**pTAV1t-V#-6{&+vNuINU~kvNW~ ztQg$M>16ZSh_WqAoY>-}*V$dmJ=g(_`Rqv^%1mnOWg@F}5_wLh;L-{8DBk`MG7L+} z$y{?7vY!Xdt5-qD^5w|eco0KE7oz{Za>!q)g&~Kp(Uj;yj)KDuj=%;kXVSH78YCWv z$GXnpdhtP=^L8$N5+A@VRS}TbsSCl-0_|cBw7E-~Hdh{_J`X8ea-M?Mbq-F#SEcwY%ARdrn{Tb4-^x$V4FHV;Af_k$ENcWD$-3Nx@tDPm(ZcK-ey(tiH z)sD5I4mf_lGss_gNLD%!lJdxdY>l=d6PEs9s`3eQudfI;erbgCdmo7nD-FikR>2(u zNpdLMh*@@`mDy2q3}%-papoB3a+c@SGqS{s`MS#=-mL6L73t}q@4Sb6Jsu5FD!OFO zf4~lZ{w)Z!I2*(V^q=M;Ctcx1#%2 zb5wukk5)1I7z%Uf+^=%zXE+3v+9Pn+A_`iR>!_ZpHTB9)hxm`vAvA$Oe%>K;EEt{Z z#iQ|-{9Mk3=MOo2CcccKz$2?+-De2IraC8_p4FPSeRh2xTX9j7QUQ~0WDwA*&q+2{oUy_ zl^qx|g7kLDV`SlRcEGbl?aMFaw4Q`b)iM}8Tk1u9nY1)XC~^@&w9x&=c!?l?MZm(vB=qs(dhh5k$69FX_+UZI^)W=Zqwj)^C@(y{bKa#`f z#TY|Bq&X&ISg0+v>fA+VDwon3yPi|cUI^c2_~5CuB#bHff!vj6sKVFz7@VdH z*Dlt=x0x=8F4EXiS_}8gV{qJ8UleRUiiH_v9JAQ{9GNOgYnKkwiTN>*IU)gTj0j%V z-H%oy4>5W2RVZ`hgO~5Rp=5~$%%5}`F3oDw)r-pEr0t4^CyKIkQ(v{NOPM)mtBE9v z*iFO01TmV`qXjK-JvgyckmS_m66e#=Fflup$@DH}nqJIj9zKp|671!fM`_87+DvbD zM|B~IZ{J0}7v{hO!Hb}6(#+hQ^$4|0;-Jd?0ox+Qli2v(M81NhxUA*}Bn-(y{gq8H zFL6FO(xLw^2^sx1WnU`>eQTYeAY2+$GnHWXG7H#WU~vGDT7|uTjW&Iw>2%(+Q=EqJx%3d^k(U8AQV=g>YrhEF3Cm4lyogas=Q1&S7csY!7qUMv z-^iw6OEl$CK)k|2t?PFKed9{%a`82YSw-xxsi*cn=4djUx!kP2a zL3uRs6Xbb*C&aJw37mSYK-xE617@Et$R&1@Z#Qe;R?$zKQss+l#amGLVgVjLoQp}q z5p;@51SL_qoH)qFa(Mx^k%QP@{}7(plM7ZW%4yKLVwzecK{bO%g$$jGl1p0QbF?r< zbKgF>N+ zD};U)k?0>Nj+^UNK(7TCY0yMkA$^FmXm$aeJYr4*hN7rb(JM@MkHd)VPvK_zd1zSj zC!p;!NPR8}bdCi46kUjxcf+Bo^d&s2Oa7Y`V(~4!ISiby#ltZzRO{0_#8K4%E2rOM4q1)8T0W9L@YnsG0=+BI`a_6)s~l3h6j-ccWYqW%g5k1bvYcZFo2AVYY^bL7j$ndgS`59 z(Bu1;bqJ65HV~@!gpLzeLFc(5+`fGcO<8&Dz9&X$?Pq)=Y;`@MDx-;-FWz1`ED0K*9w-2Ij9xsl# zUcu!0$9PUVmkO==N##BE;w{xuyzpHK^S!p=$Q~gKy|bQ5x2EG46>qGl4@dS>Yg~09 z4W$fQFi{xb{>TwLkuHr52Fq!EjVk4TJOz)!Vmk9?AT?C*$LZVDFeG6F!MGhxA1a5_ zyvY7ghrBLmIW;itgUccOsT8CP&xb^pVTj*r0+SP^ zeqYqX-UTtTdOOH~L<^{HQi766)-cOd3D!?OkGh+uW1JJg1A7zjqF)D{_R57O@BKul z>FCnz%3G*9PG9GVnFVMvDVXSg5bxM8qcf9>@IkRX)!j7%-%Kq<)39U|3PF@^kVAoD z6KGvqhvllD@%TGce6;2zlIWAra#))#yuX|pa^}!U>z?B0q7XbZi645s>*3zK%ZL}o z<4)I0s5f;79*;jni!Zp)wv-omMr;GFH+qQQZz<8R!SC34*A5TQ>_uqthBuLC!68ip zM0r<$_ih>F#=ge{X@YTAZ(*M7Csem^L_zU=ByVs&==hkSLy!cX6?skWeUbo)#~iZZ z)adG-y%~82m(yK>ui)}XJ?XCGP@R_nI5%!4?A&t(Wi+P(pR^+QD=!0s*w?^r+YEWb zA{ZjzP9>W~@UAlpjcT7l!Sgh-sJ4LIc9w?|t8-vK?`F6htgyK4RBYZV&Dn1#$h3A%*VtCXe^jONp>E8@d@kj_tR*R9VZBJ1; zb{fvke*?qq^>Bl6MTan7RMP*5%?r};>+=&>Zybu=fzf!(aUZVNI18VD4r0>!8eC?q zi`i2|=-fjKKmqx_b;p9ixMd`Fc&}z_0J#{Xa;)zwYWNFAG658lBBMyB#weP4=zI7jLJ7GumJy?%x`}N51nZ3+2qa)yBdX8*sx=$+m9}}Im8Yr4{1w?p) zz+ixZBKUx@avdl$?q5Tad%`fcUw$z%SkaW!kbK+O?Q{q9XQpLPAb$ zC#nOUp!MwswCY?&jrT9nn&pEIrsb$h#Bu8-1zd4%A)Q(rOqJud(>0y3^x}{qrnBdx z;f8GZ-k(MJj`vXI&6nxi1Myf~X@@87iKBbAEfnY-z)82wp{r0CK3*?@>}oN%J#jvY z%kRYr_8L$q>yKtJbr^r38S`Dr(DHRWs)o)-Jt-L^FU-(lp)H#J)WxQUNTst1=|$^> zG&84=s-4QlweMHJ`wxSNU%N3zbv)hfehfR4o?x7^8d}a`!AH8_D zwXns=Nymx*GZoyp9X71;J0YRMhRgt} zfkW*LoLV;?gBMT6MpFkmeMb>q79GKhR_SqlJBE0>GoK{%ptC8G zmRmH@(7<+@(l1RnPMnA*Mn*q;rYn4$^a?rGk-fR`JKW(;gZ%hI5YQBXBDPl`vY`-p zG!l_}J`NS~{=iFS1GvmE9jD4fVR-5TtS}41su~;0J2!`FC1p}S_B(p0+JWvddq?GZ z7vh(-O?c;y9L~=*fi8DZOyqu|6TR7(F-H>JI9^m#K^NnrBq3b)J6tIF6L}iTAfA6N zUcKi=7aNVkqXBQ|)OIdD-))6cdIyl#cNWYqJITC#@Px_QTuBZFTEmv4)p*gR3R4u# z;dJO>ayhmfW*7GW-@qJbdTfr;qEpG?jLXc`9v>3yoD8WAi%C6Uk>~5al9&c#sD9uC zo4z5$#Z|+ZWG9H7{Sqg?&4lZF>QL9?E*|dgMe|RW$=l8!jC8^8pWhk9mru#oo+QxP z?T`M8H{piZ*?3gY7hlGdP@!E-bb5O%cCC9%efFHAiY3Q5*W>MR^}Hk)ep`TUi-+ho zMwSNkb<*I2MpWQ^I^IrQjlmK((4sjIPU$~`jmCPg>P9$tZ7T;dI0)fb1Z_!taAM>f z%JJ9Y@ltOT>VJS&edprvo{RWYuZFG=k)z@R{ZvD06<+rZ!H$jR>9}?t++kjY1rZzZ zzV2lBL{`9KXva+^p?GV?AuP?Aibe7rxOF)P6bwMU>-Gn;mHiUP(-9RzPyrK{!@wM47Ee zIf|d$AS(AQD)q;pymbQ{m+i(G5$QN9x(;}kMF9WP77|U=?GI;X#p@)i48}}Zg#)pyBuQAKwAPiqf z0JGaKpe?iyV|m1&ksFD#-qvAo+ew^phl__T-r!VES8}TJ50DtI58l1WC^%CXB}^V; zq@^0h1rK7rv=OGVL@}&57=_#xkkVONzzM&Lp|YL$NzZK5}BnN|U3@?&xCyRR7aZW~6Y*`SBI3GUjY0LAVv;h@qf zbW47PEQ8HxGJYoVmPexL)YZ`SehcK-m_o_MCy;kVALNY}fzUF2IQ*m$`WvpJX1o<{ zm^*~)9sWdDt!0?O^Apcpu7JS}AK;m1F1(apjLWVaBO0p?5U6+n!hVwEBJ%l~eCI`>10OpYZ`If3k0<}A#MQ8zTkFcpJ<% zxQ%Y`5K~9iL380L2shEisca!Qw%ZO)&i{R*95cRNiR=(+fUdVcFv|TtiYyVqgU3_x z!pf&u9CaA$UB=Ui7J)QLteDQsFU2si{WN$zgPK#nVE*kW)Z~nZ7ejTZbNno-TGT;X z&R!%u%Fug95mX%}fd1EQFl$a6Olqo!c>#xEQ^HMXyf+h?A4kJ;fp^Fs_#9X6=Ae}C zC-hvDgL}|AXOy@Vpn;4gO5hVcEnu zcBc`6*oscFwA$YfXEynToW3-V3-s(SzW&NwM%JD8b z{P7X5f1gO_D#=sVqawOeRU@1YH@vYn+7a&?X<*EnG~CuZ2)FKwLrUH(XfrT@+GLMEH=~V3Z zAy*?_!kptD$m&CyXW49+@XXDDbujt+viW)18;t+ZA1ZCq&>|G8jx;4rl+^iW%yK*si65M{|K{+&x|8{o zzodb(-;`15>jj)OM%B6#!xd4@R~SoPC}VIqFM9Mv;##k3aKddPIF>yG|2ju-@3{iz z+aH5~(p_?io5}P_J|r^T2S{FpDlxbo1P+8v3g-2arbS{fkyjaF8;YR#ObaB+*}xzB z%AxvdIUHJL3h7N5q-#J76!OAhNxLrd<9#bJUB8$l|BM7bu_rJwNR%9cB&PV2988=n zL5g~Ikd4M6BRv z5{X_}5voWTqkP&_IO2DPSytD@9zw*?sPro zjp{qtx2hSskN=6nvogR{-kTiTCjGCELmTn%Oy=R?<%w8W!pAel|2-TEWTx)xhYOG`&fO;dez{A;MIYxA(w{wv47rHhN@|IY6JwEw>;o9w(+lysu@5J>4mJQd<%WkI) z*VfY8$&2OX?6%X!@9$#wzZ))a@pf`j9JTkdcVT&XSi1ZwEWPbrgm-aWxE_{XTpO0v zZkC6Nr@gbvKmFRf*n6q?S~@v$G_=(jW)aKT((muhUnkEZmIv3_a@59U7wcCm*TZwv zPK~WCY~wxJqrIyOYk?=%T1Ang^UL_xRHy98^|0r%?7g^u8*tQweLd{Gytpo-p*eY3 zx^dK0)YzID8XBX%xTD(maGg{}+cY#Zb<{M3|7!Ad_4cslj&^OY_{)0KTvJ0>_}?RE zahUO8K zE7#Wc*P){_%KfFIr<<$CKimHjcDwNCE>AnIrN>TBFAuJ#=V<#crR7G|bJ^)*>EzA* z6}L9)e+IkrSF$;B|H!%=i|gU)IogY>%dbdBtN-{I9-~`wESG<3=gCo5W4U^Jjb?e& z5!a2Q!LqTpwe|KK-S>}x)V2$&Fs#v_Y+RjpdUClo9Cen9n}@6Q=&@s$m@&o}V~jDz m7-Nhv#u#IaF~%5Uj4{R-V~jDz7-Nk8X#NLlI~lbAKmh=!8T26l literal 0 HcmV?d00001 diff --git a/raster/r.sim/r.sim.sediment/testsuite/test_r_sim_sediment.py b/raster/r.sim/r.sim.sediment/testsuite/test_r_sim_sediment.py new file mode 100644 index 00000000000..d26d2d06efa --- /dev/null +++ b/raster/r.sim/r.sim.sediment/testsuite/test_r_sim_sediment.py @@ -0,0 +1,106 @@ +from grass.gunittest.case import TestCase + + +class TestRSimSediment(TestCase): + """Test r.sim.sediment""" + + # Set up the necessary raster maps for testing + elevation = "elevation" + dx = "tmp_dx" + dy = "tmp_dy" + depth_default = "depth_default" + tranin = "tranin" + detin = "detin" + tauin = "tauin" + sedflux = "sedflux" + erdep = "erdep" + reference_sedflux = "reference_sedflux" + reference_erdep = "reference_erdep" + + @classmethod + def setUpClass(cls): + """Set up region, create necessary data""" + cls.runModule("g.region", n=224000, s=223000, e=637000, w=636000, res=10) + cls.runModule("r.slope.aspect", elevation=cls.elevation, dx=cls.dx, dy=cls.dy) + cls.runModule( + "r.unpack", + input="data/depth_default.pack", + output=cls.depth_default, + ) + cls.runModule("r.mapcalc", expression=f"{cls.tranin} = 0.001") + cls.runModule("r.mapcalc", expression=f"{cls.detin} = 0.001") + cls.runModule("r.mapcalc", expression=f"{cls.tauin} = 0.01") + cls.runModule( + "r.unpack", + input="data/reference_sedflux.pack", + output=cls.reference_sedflux, + ) + cls.runModule( + "r.unpack", input="data/reference_erdep.pack", output=cls.reference_erdep + ) + + @classmethod + def tearDownClass(cls): + """Clean up test environment""" + cls.runModule( + "g.remove", + flags="f", + type="raster", + name=[ + cls.elevation, + cls.dx, + cls.dy, + cls.depth_default, + cls.tranin, + cls.detin, + cls.tauin, + cls.reference_sedflux, + cls.reference_erdep, + ], + ) + + def tearDown(self): + """Clean up test environment""" + self.runModule( + "g.remove", + flags="f", + type="raster", + name=f"{self.sedflux},{self.erdep}", + ) + + def test_default(self): + """Test r.sim.sediment execution with defaults""" + # Run the r.sim.sediment simulation + self.assertModule( + "r.sim.sediment", + elevation=self.elevation, + dx=self.dx, + dy=self.dy, + water_depth=self.depth_default, + man_value=0.05, + detachment_coeff=self.detin, + transport_coeff=self.tranin, + shear_stress=self.tauin, + sediment_flux=self.sedflux, + erosion_deposition=self.erdep, + random_seed=1, + ) + + # Assert that the output rasters exist + self.assertRasterExists(self.sedflux) + self.assertRasterExists(self.erdep) + # Assert that the output rasters are the same + self.assertRastersEqual( + self.sedflux, reference=self.reference_sedflux, precision="0.000001" + ) + self.assertRastersEqual( + self.erdep, + reference=self.reference_erdep, + precision="0.000001", + ) + + +if __name__ == "__main__": + from grass.gunittest.main import test + + test()