From 02228634c08e440c99f0b7d0a15ad03ec621454f Mon Sep 17 00:00:00 2001 From: Claudia Vitolo Date: Sun, 18 Dec 2016 10:56:37 +0000 Subject: [PATCH] Added description of package (as in paper) to README and vignette --- README.Rmd | 4 +- README.md | 4 +- .../figure-markdown_github/example-1.png | Bin 20719 -> 19548 bytes docs/articles/fuse_vignette.html | 6 +- docs/articles/index.html | 3 + docs/authors.html | 3 + docs/index.html | 6 +- docs/index_files/figure-html/example-1.png | Bin 68878 -> 81574 bytes docs/news/index.html | 118 ++++++++++++++++++ docs/reference/fuse.html | 3 + docs/reference/fuseInfo.html | 3 + .../fuse_hydrological_timeseries.html | 3 + docs/reference/fuse_parameters.html | 3 + docs/reference/fuserouting.ranges.html | 3 + docs/reference/fuserouting.sim.html | 5 +- docs/reference/fusesma.ranges.html | 3 + docs/reference/fusesma.sim.html | 7 +- docs/reference/generateParameters.html | 3 + docs/reference/index.html | 3 + docs/reference/model_param_table.html | 3 + docs/reference/modlist.html | 3 + docs/reference/modliststring.html | 3 + docs/reference/readmd.html | 3 + docs/reference/readmd2var.html | 3 + vignettes/fuse_vignette.Rmd | 4 +- 25 files changed, 191 insertions(+), 8 deletions(-) create mode 100644 docs/news/index.html diff --git a/README.Rmd b/README.Rmd index fcc16e7..37e4b42 100644 --- a/README.Rmd +++ b/README.Rmd @@ -24,7 +24,9 @@ knitr::opts_chunk$set( [![AppVeyor Build Status](https://ci.appveyor.com/api/projects/status/github/cvitolo/fuse?branch=master&svg=true)](https://ci.appveyor.com/project/cvitolo/fuse) [![codecov.io](https://codecov.io/github/cvitolo/fuse/coverage.svg?branch=master)](https://codecov.io/github/cvitolo/fuse?branch=master) -Implementation of the framework for hydrological modelling FUSE described in Clark et al. (2008) and based on the Fortran code provided by M. Clark in 2011. The package consists of two modules: Soil Moisture Accounting module (fusesma.sim) and Gamma routing module (fuserouting.sim). It also contains default parameter ranges (fusesma.ranges and fuserouting.ranges) and three data objects: fuse_hydrological_timeseries (sample input dataset), parameters (sample parameters) and modlist (list of FUSE model structures). +The fuse modelling framework takes as input rainfall and potential evapotranspiration time series (areal averages over the river catchment area) and returns a simulated time series of river discharges. It can be used to understand the variability of expected hydrological responses based on model structures. + +This package is an implementation of the framework for hydrological modelling FUSE described in Clark et al. (2008) and based on the Fortran code provided by M. Clark in 2011. The package consists of two modules: Soil Moisture Accounting module (fusesma.sim) and Gamma routing module (fuserouting.sim). It also contains default parameter ranges (fusesma.ranges and fuserouting.ranges) and three data objects: fuse_hydrological_timeseries (sample input dataset), parameters (sample parameters) and modlist (list of FUSE model structures). ### Dependencies The fuse package, as well as the examples in the vignette, depend on a number of CRAN packages. Check for missing dependencies and install them: diff --git a/README.md b/README.md index 540d9ee..c7ada02 100644 --- a/README.md +++ b/README.md @@ -7,7 +7,9 @@ fuse: an R package implementing the Framework for Understanding Structural Error [![Build Status](https://travis-ci.org/cvitolo/fuse.svg)](https://travis-ci.org/cvitolo/fuse.svg?branch=master) [![AppVeyor Build Status](https://ci.appveyor.com/api/projects/status/github/cvitolo/fuse?branch=master&svg=true)](https://ci.appveyor.com/project/cvitolo/fuse) [![codecov.io](https://codecov.io/github/cvitolo/fuse/coverage.svg?branch=master)](https://codecov.io/github/cvitolo/fuse?branch=master) -Implementation of the framework for hydrological modelling FUSE described in Clark et al. (2008) and based on the Fortran code provided by M. Clark in 2011. The package consists of two modules: Soil Moisture Accounting module (fusesma.sim) and Gamma routing module (fuserouting.sim). It also contains default parameter ranges (fusesma.ranges and fuserouting.ranges) and three data objects: fuse\_hydrological\_timeseries (sample input dataset), parameters (sample parameters) and modlist (list of FUSE model structures). +The fuse modelling framework takes as input rainfall and potential evapotranspiration time series (areal averages over the river catchment area) and returns a simulated time series of river discharges. It can be used to understand the variability of expected hydrological responses based on model structures. + +This package is an implementation of the framework for hydrological modelling FUSE described in Clark et al. (2008) and based on the Fortran code provided by M. Clark in 2011. The package consists of two modules: Soil Moisture Accounting module (fusesma.sim) and Gamma routing module (fuserouting.sim). It also contains default parameter ranges (fusesma.ranges and fuserouting.ranges) and three data objects: fuse\_hydrological\_timeseries (sample input dataset), parameters (sample parameters) and modlist (list of FUSE model structures). ### Dependencies diff --git a/README_files/figure-markdown_github/example-1.png b/README_files/figure-markdown_github/example-1.png index d9ae682df0a708ac8fbc3025bd2c9ef9ecf0973d..4d2a22b8d7f5b4caf5c2ac6fc8915dd79da19c51 100644 GIT binary patch literal 19548 zcmeFZ)mvOa^Dc_Jy9W~7-DM!STX0EmcXxMpXMo`D9^BpC-66Q|$@kme^X#*4&iMn* z#az_%>Z-S^-l|orSI>ki%1a<45FmhofgwvtiYkMFK>)$Pz_;NbLEp$JGr@p?L4qmD zs)~U`aBy%42na|>NGK>MXlQ5{7#LVsSU5O1czAdO1O!AxL?k36WMpI%6ckicR5Ua+ zbaZqK3=B+6Oe`!cY;0^C92{I+Ts%BHe0+QY0s=xpLLwp}Vq#(v5|S@pzL1iVl97>p z{rdIWw{PU+>L~%oSd9oTwL7T+&nxyyu7@8e0==;`~m_3f`WoVLP7ulKv-B)Y8(@*4Eb1(b3h_)zj0{*Vq5~^XIQ$zYGiv3=Iv9jEszpjZI8U zOifMA%*@Qq%`Ge}EG;dqtgNi9t!->LPJBt!otGC z!y_UhA|oTCqN1Xsqhn%XVq;_D;^N}t;}a4R5)%`Xl9G~>lT%VsQd3jY($dn?(=#$M zGBY!?va+(XvvYECa&vQmKww^8UVeUlK|#Uq-@gkB3yX@1ii?X&N=iygOUug2%FD|u zDk>@~E32xis;jGOYHDh0YwPOj>g($p8X6iK8=IP%nwy(jT3T9LTie>&+S}VZIyyQ# zJOBLo)7910-QC^O)6?7A+t=6E-`_tlFfcecI5adgJUl!yGBP?kIyN>oK0ZD%F)=wg zIW;vkJv}`$Gc!9oJ2y8sKR>^)u&}texU{siyu7@!va-6my0*5qzP`S(v9YE7hsVdqr>Cdq=jWG~m)FjwWqn)Xg9@7Cq8rsJU zu!ZC%;%Nc!$o4^x{vu#hM47YmcSs9COzLCT0TLDMKpZe5<5Gv3|Noc#fAq`X00cu) zTwiTq|I!`0+3tfS57AoxRz)Hpwe8{2ewTU^y9oIq5JP1y0)#3z-7jk^csGAp~;S zfpZlc<)*HPe_KT_dq9);diXl!YG9qAv;?0D+|ZPoAzEQ^OE#@wKToU!1WLY5{~n8y zX!jZxXaU@tAiDi-hJOIY03t&OjwHMz<2~*CPbb6ynO~TLc_B9T)$C7BGMgtWN)F_o z-*tc<_26%X|0c}R`?$h&2b~xF2M)(}#6LI|XTa(zVs5*8?8bCNRmhyjICCBNNpq7h z%#zN;yBmgQe%ulmoneZ(79^sjaWpDLZ6G$orfY2i`eS^|yI>52yiuyT4ibvRzZneE z<`)y6_a{K~$6Vr9ZzAZo9L*wo@a6J!GI?NXy<`Vj)w?$l1Rv`#`y_NbJ6k`zD)wO3W*A~Ku|-5fXK`Sd>p%wi@5d*+{ zL)h`%4zgDr~Y3VQm2S(5%Hujnu{#{-{prrl6V*nJWTu--vQ>Bkp|Mqss;#zqk z0Yy!gHoeQ{B-ip^*h;9TNjnl3+t^zXK@3VlS(mh>=gDsaxirQCA{weCWXU7PM@dK; z#u~_#f~PC?4@)PX*`oWgoiyqH%o!$6glJd}iuf01T6lJ#!O+i7PiBu|eJ%%FP&!cO zSR;4V4$uGnH(mzE?sYa_^`m|eHSiVSEM3U+?}`*liq6{UbJl;{q^$}R7weW#H%H6^ z@lug|%fZ9YKbWEc6vVr-g%J=PgAp0?xDK`TME_w~)5;!tdq?uS>dk)TD^^rWCu($G z#F?JN@YHxiJ*jAFED_(1MT+tNUdarY3}M@G2Wm;S@8@6VXlFHpA&?URC;V%bONW9taYLNAB{BO|l!>A|bMT{-xq83J z8F{F&sO@4+L-;P5G7u+b)B~(_fNcroMY0$NNJq+`s0fgMojguLxC?l0MP&D<>~3BR zN$+h1+>@Jbgx?JUWjS6}`-x4u@=WY<3m3w?QQirj%%h(AXjnj@Us zh@;~u7s&0#nl(R2{0;AF?V5AAXg-7s#})9P3)1#YJCWsXu$#}J(4_H52=K;9gIj`1 z%`QdBL`&(A&jeU7;wpsF0lvgH2Fs>gv7-urh>3C;sPubo6#B(kKNJgiZM6GGm7e7@ zWN$K*y1fS!ueibkNe{!a;>&*sxxiTV>c5`>l_Wfipv8*Ti(Q(KK_oM>;ELKz+YRQc zoULZmfrOMQCe|+TZL5gM%tSRKRT1rl`KPG1V7&PJbM&^{rD*vO;aHBhQU}PADvk@j z*_bL2yPmjV?~C^@LepZD$eml6L9EIAjYq7cF#6QEW(0Bh7brJfoOJ1ZhyM~n5jCCt zPVmZW)7&(|nlT{e-wJq-lS<*9&m>SX`#_vJdATd#qI}-|l;%^hawKV>%ywDU+&@(< z9tZi;uK_Zwk2hX4S6N_b5^F6Bl0mXK8Ij#%qx%CQ7&41B6hs3E-%BZ5*=vex{R`RJ z3C-y+C#_b1L_mjLX$0lpL@5Dt*nf*M;ifny+#x0*SD{w&l!r8^qCS8?K1Ens3Sg?c zbQ~oDEl5&8x=23YXo=p%O3QLc3LraBpaZ<&>Mo`++@_bFrQrJqNY#*chEU+E$YQH=Y94F}Rt*xj zRL}(Yun^|AkpAORu-pN;QtO?~ntSPg{`&Ci_&tXJ-DbDbyCQhL(bsIXkFYB>gy&I-q|4L)S*#Iu=j`}UzumYUvUTmqC4Ijt z|0py;3r1JPzUR$&qM5FM#abx2UCy|8$;0`$lx}m_(!)3Oq=kKnM_P4m;Z@OO)F*NMf8`B#M8`>!DR!_U%LzntpK4USgATH z@MOUix9{VxOY3zW&TKt|lVqWwm5*r~~5cRo_9H=S^2w z(;J9xwmDKQC;?DbxLy8E74e9N=*2MQU?0TM9xWozR->MdhaR2q`q_T_{2M&mq*7C$ zT{~`y!T}XTl(4SUW6~h^9E9G{$>+Pa6gjf~Waoj{C(pT2Ed1Nz@jgPO|FZ>yaU+7u zYZ8!y`~|61sxR==ScZDi@J%~=Dphq#;5tfY-gnk`ut2xJ+`-=>;5GiI#afp?T9amJ zI!zONxmZQU$+q6xtA3_Q;)xFXyj$djcGP+xi1+I;bC%8m7>DBmhh}C`42w%L0Nx8f zz_-hF5uN0dj*-gU**PxTY)z%H?7y#|mH{k_vQ5K)1V)(KDa)m4*~r?^fR*h4Vj~SM zZ1{eY`0E4wirh&zbeSNF%~2R*O4=(vSr&Doo2a8&C`MvB4X)Wsg2ZS4@T0HuMYn(X z-8h&I@D%gH=MQz`jSUXmQGLTwYgsa(`RN6~XZl^5aZHBnRO0v0K=1tE+6uH>!{Q`f zW`)@qx)e7kzKzDS(%Oca2rPrD<}rw>-hT$VM5ZoX*ov9$P0ZPRG(D@~dd&|YymS6M z4qm{uZ<=T#+{AsUUMMBj7pTVc_G9rq^ta>NSdsEfaUnHXQ(>Hj?v9a4kb5Fl_5Rvn zvtPuiH-WmdX+fZ67&(xg;6<9H#yA~?oGWHv)PAcR_Hs)B?h&KXqJ!4kcNa~0p}9P8 zg}l=M`2mEray=|W<1h_#Cdr(F|17#H(xHafp?P5~eV{+ic%#^9aPWUPHEd^N3NGO> zt~i3q#mH+UM=%S&Ga%z&P#*5M?01SAg>GZDZk)kHXU%q}5uKprK zk|-#;ki#PYDS1K)s$})(>60B;`G;6JKb=p#l^b%Uc8u`tYjf~T9;H8pdd;0NxUQAM zG7_u17D`ymo4CsBz7vUM>s|WpnxBQDw?nJ-I}ZLkLIMQ+f(uSwo+J=GX_!x*e&kj& z+au>xz&kFSSWKgggUiQk<~C?w<1$Y*3h(3P!(BWT_2l_H#x+y}$*v?q=8%xSn+o#m?Co z#iL54wYpXRN2LesxYMbD5E%O9P37b}lK?r4w1HX^x)1xky$$7#bVzi$)JvNe?hnBa z2o)XA`<9_%(MzRSN7O)47DbGDBU{TgCU@C8gm)BtF*)*I$?C` zc+gVa&(7Y(JT8S{j-QhsLIz5rOoadr21}v2RGj1otSYqAGnNy0mGJ;3e|@4q)2O{; z`-_p6k8HhAYGoXCXxJho(91OVLS_zg`tZfr2=g)leTSg+eeo;xz@q=lE;su5`{-i? z3RCI!=hSy$RK#wfGOfZOa24%}Z&@FVgTP{N8L6POeWf<1*xyH&CjL@ebDPZVxG(96 z{^Z*;`GjVl=9OxgT3v|#B(U%QS)dSpy~YK?#t4EpO?oJ4)(q}D{cyMRd*jB& z7<}bcFaR|$r>fIqeecR@aWg3Z@LqbdR)hi&(i~RyWyJJ35h$PgDstvZHofKVXnf%V zX0u=S2Y+)CMWqAS`Rgy~9b1~ez&P7_Qm(2ciF$HAi=#F;ICrMeCqqhp?3>D;U%mt(IVRRdSl0OB`i zyprKgeJ8wQ&aN+|EtL-b0a$tExxCiNU;S;M8M6NTBs$WuKP$f1Z~eJC1e!t;?!&qH zyGZGw7IS+eh!GK1p0a@!_xMtnVlU7O5=Al9??NOVz}3cK|Boi7(rWzV0!jlWOB5wn^8-F$^2t?A3LckJ z{BdcB;A}{^ipM{*nI#;1`!6bDWcsh7pd_GPt_E|UHw2`Yy=3;Yk5o+w#2|E&$9a!4 z4ST}t8eOEObf*w2aV4VPqBbbn*|vvT2)2`7P-F&_LZJ>ItkcoS;qcuht5VL#9#YAo z2oll!kQ%g1BB0vnYsvyPg~*2j&TJEN+EY+XUNav;gx$I&vQo{)S?!=Z8?*&|sPh*5 zr7%ja32Ej^@c)m;p(JM{5R9zC{C_GW#`c~kwQpGt7T%bTn^6ut|EY7`UAWkGKv zd6?HL>trMYb$~-ui+^ET`6|Cs%(D1F69m;yH#bdE`GKPL>^@3!{QljcMI`dRcZ}yY zC_)_pU{Hr5`{f)`Fan&bb_Gn~N~0Z z-tu*T_X+;w&WG;|sS8Koa3#wu%bz4snDN#r{l>}R1%oe7IG&dx z=o(L7FOYR{@^j}Z>LV;BnZPV-owH)?KZ5@cAdz;Yx55uQ_MFdp>EtHbE!sHg{}b86ye?iNeFog$ zX6yklC zOTYI&MPOSma82NBfRR|UuZz$Vm+1h-L}m>6ip*L5wF%|7s=vP32%L8?<_|HS^?xTq zY+4wS$&dlmFRx0NN{bk@eUqh*QS5sbb1Dj?(o-zyvUdC~0_wadF#^B$Q6px#lW>idoyeR*&AyBsiV1sC({quz@VXitXIeU5 zp$miV*jl9+mgC?w-tk)J>7n-2%e)YbpD&rhmwBSVob5E{Q_crBbcLY|{4?!siU}oB zV&($_ZatL6!YY`E8Oc}tYeTEFFR%V;TQXKZsLCAuzDQE$7ApAz4A;AjJ#u7Ro^ZY9 zuXZhcwPuyv@FCv(%wOL%n5_V}L&5pJw`4GHcxSi+S~W3n-PBlSCT#m2Tj5+kK9*?E zj=8{Wh8%4GdGM%yx3xhi^+gnVvqam6@Urch5@UUXgyyzrwi+b{#hddqY-pXT&)Lgm zia^YXxC^~eXNY&GZo_K0_V%;EEkKo5>1c&T=L6MFp}Gi|Y$04sP`3%8nJwBJe(s{< zx)sm`I!pBrYwLr*i2nx@kzmaMZ2%9WK8;d1V5=oi{f^%H+zwJvfNkDXOj`##PX+-@ z)-yYlu3OP!-h`q<AlYW3eC{Q*@=g)A1tDU;u3}iq2z|_mYhV5pz|LjCKm!#oD&ePyH zGm%HOsuwO)tDKhMsnd`AYg|+|?<6FxNII>}1C)&+bKyq9Id;7X*1e8U>O(`DQdP0= z{Smmu56_Z3w?Y5UI0A9acfbl}zZ&eKFz9)sm)kfQ!e4I;)>yfvrcY2H@?iYRCBNc& z7RYv0-su9Ye~}fZXZd#n@$XT1NPn3qij-tV?*~x8DW=f0L0RWLQs8TlRVVEa zm`Vy_Ub8N)DuU1K;eKZnY(Em)@kBafhh32~|M@uY$BhysNv4za?oRt}<>|SW`MAW45A~kA+ie@P{>z5gJV8&hhw%-AN>CO;x zGgK_dUEru6U~7N00xa8qm!v@WlT<(%W1!Xa(E`b%1_s8DBM%d5#BmiD*{7F(Afn`m zUS58Fy>h;LnUughGooSn1L2Ve1qAg1w|AXfHWL4F9fAUFA7r?73IuWBc=2Ofuh}h@gy{#0|;F(_8fAuVy)5lT$;l(%}AqD zzt2Kky|SvUk1^`K)~(ZSbi?6|jn83!^$6G?GG`Pemkljm)1dMu|j5fZUcsasVOs z{ef3048s2_7x5@3q>8NkrhHeV!W;gSD=i5t7C;0pmPFd2=5c#T)&N};l=O%d)M0wm zb4Uje7Yjq~nON{xJiHNA$X0BO$5fzb`5b$yPEhQJ+Za0U5lwX;MI3hqfxISn>X_66Bkz% zg9~V%q)GV>LlHSE?nAhIV4#Z^$c!XIxdawTZS57DzqE%Nw?q&_P>bDwZdo5NPwd7> zpaf?!q(Xp-u*h%JPs#Sjg8b`rm88%WiJB9HIE2#{WZ{XYu@#W==@{wh-S>jgDIaa} zyDOMZ2oJ}k5VXyBs)Sb52BS`uaAHLW;sQl6FVpBS0IHc^*YjxlYB!K_uPzjtDKi+q z5;3@38{M<|wlp3vzAMiJnU8P&*@FU}G*jc?@A6k0!Lq?qm9IRStG~mSA+cQ29}utZ z_UH}{ZT&c(Dj%4oP2hsCqn&l60C%e&{3v%ayYvASuGvVYlT(#IRplQ&2<3r!y-eKs z+f&pERL&kR&U(t>h50Qe5&AcJG&{hq$JxVh%JKlfjl+UAvyu9>cbM#BX5}H#gZCoC z%7adn+#-Mden2b-1KH%^&x!o31AbA{FCNLz6O!*pmD?hHlm1w+OTt^kz(sfQECqsT z^TmvX`gZ7NMIs(-IIQ4Fl&plHXZHIG%40mC8BNUV?~TkoC9{im#JudVsCWlz1G zU7Ov#XbhQB3YyQp@MvnfIu!W3wLrUP*}#K9WbF=gXSz5~P)N%pSd0SH9on~mnASy9 zb0Cg$!%=XJx?_fT_T3tM9sW^CCeLLsabEX~?BjyuwnOe4j!9TZ$>hXGI=zs-(Nn3m zcR-%ulzwA|OHl*>C0W}iy5V|OW|BEz?0K^A!Zm`#(LQ4Ly!ls@=lCx_*Jbr4Va7ai z6V;JRGLn;{V9Z~4@6r#euoRrER@PEVjML7*gQ}Y|>56pCg;f%;Y0y8$u7%t_c*M3=S|bvKZeK7ZwW zmXW+i02PsS>!7J^93tvecH>0tfx=?jP2T$a6uyH;IuNTvbEnyAQ$KPd|FY{}Rc#Hl zGj=D(PAjX_jad%xdycKnu~I8o#rMXkc3z@P$R^q0df}iZ{25Yk^>uh*098h+!N!gM zo&~5uM|3K%T7B9eGgw!KYp#Ml;kqvwYB_?;%)X)KBTe{puoq}yTCOwc~DYy zJAmXW!#B-J9eWYC%uHl1Dc4iMxfJ#6gXt@(jWiwm_6VXFgg`AoVtJX<{eQOMvYbF! z51M!>7v>)4` zR~Nj{o*a3EX_XC5G5}B4#vFd zS<9ZoWD9(*;Lh4Y&>LOxl8v%$+*h3!9b96K%%aSfryBcWEv|?*_>&B4wQ0ETa`tq5a5Y{z-kk zjUK*tK&`o%j5x0@g|9)2KypKE-dDySJWA|1;QZtF;U%f9ekUWN?2ZxzQwQ;$IvW_@ z6yq^M^&!>6C6|Ezj^5Uf=B%PIe=_7;PxeVd!)v|xUB5mDadj+qMN-8i2LtcHgd|?x zb471kHK2A6%P;!ytVDat({7Ab@dC%&_8!94@J`REza0A?MMtKpnXNx#+_otXqw}Bu zbhyVM5yJ13a@cmG{GbY_$4jUJy|lyv*`cetk&Eb8&8yon;L%`DUgcssp@!hwcH7ps zwl8e|=w6CEx$vI-%EsGc^qqIX%$3D9it4$*vAbY1Oy!UUjT+%`Eu|&mn(AlfC@ccl z;}S^GGYrvWxO)5T>(X)P%MwaDFgaLRylrQVKk5z*k++{Jl}nj% zy4j#E$<`tGmyT9BCUp6^Q}u#R`bAZn@{XNST-5n|iJ+YCKn$5%W4loY4!4?(yuuT+ zK&@+N;mTN~#5V~` zTcL=5$U-F|?V^}bawY1{s8;sy!V9SxOWb|Uq{GnIP%(L{YHo=3brYrDMHOU-frQ5| ztxtVvgR2zFhy|Qj>Q;Mpph8pC)4ihz3j-_#TbPKy4Kgycvj}?SKGgza#X@+Bd`tl6 zH{aS<;upn3qvU?$yT6~-HP|`yhLVl*NV4W(lN0BV`iZ)#+7(P|iu@Tn395E6)Ac~S z=|?hn7dRkAfa*sWlg$$ISa2h6e8D!&w0r)$r^1Wksnf*10kkY}Ni%T{Q^Z4K@Se+0 zDF!DSas4Mjnp)Sy=wu{#hePM3Q4yMr(XburD=$zVjx~_|mv49?2m>PvveRG=fmQm( z@gCChV*X^eSLLskUY~)1LQN|MFIS}y!+>iaL*r_gL?oXD$ZeB9?CLczp;%$@&yg#3~dWG?pBXhy9hGjVp59{ve}Zd z&I#xR4NrZ(BJEcMXkO|Su_Quk;mORwk0&vDiaea2HN)3$AvK7q(I7XrQ zM9EX|OAhj>hp74kW(?0L(i&okI=fM*M|fXKob&#Iv!`L=V)m_t@nt2#U~-|AW8-?b zvUz>1otc{#@(A7I5^LrZ_1s?zDg}JSM!*{I?5;CguS;@}X~G0$F_(Y!&%tn3d0j>C zxT$Wi&QWJIZX_Wp=E{XI4ukkl-QFoR-M~4-&ev#+&)W77H5_SkLFGHlCoZ&u3*mA( zQ7i=ZpDaPu0@DD?D93QDbKH`rJkN=)lFa6|Y~ZX$-_WVk&0mQbk!br$@pXVm_z!rW z(*evrhoK?-jx?Voh7+5)5#gtMCEg^ygeKA=RbXSJSoesH0%_Fb*}D3akH7dX+~vNE z?P*%Z7-O8u7=mU3;c|HeI zk+2Svq1|;)npgsLXcIeGs}{msl~o5K(GW8iAGNIH+i%4U`xtEq0Cg+~SX9`sHTM{F&)k}A%r)8)|_c?mol``ja6tu$Twa9b~683AYhw(N2MJpFi-WtFTIdivX=rJuF*;vLc8F7HI zkRJ&M@k}z(Jq`Sdc&pDL7kpzH3~^BDIPH6H)y1yXWag{-g_cCxpUSS66w-LsHz8y_ z+S~tw(2RSSG?Pu;R}q0;nvwiocuiyA11dS@BUK&T{2-aj=w1*IaN3eyH)95mfIPv6=dzAv(CIw+^=cR@t|tnQqTEY?AjL=qSJVNx1vD1>c=Y} zX5o=_<+2#TY<2GEn980k@4z(_i{dbBU!Zp06T2o%SZ8lmx4dTo5BR=w~RH zG*0Oi07yA@FVk=GWw`z#Sua=+G%*u9HhSgW?wUjz&mvp+S?f>skMV>PDzf4O)m#Ym z{F*f{?!KX<_Oh6D-FY__2dm(mdo) z@X~(oi+0?`g43=}c<4~vH}I&>PyyeEusuyYt5ziHZY+jbmCZj0Q2POlcg<0e*nU8q zSJlh7?{?GCiAyMubDrj01dbF@@bK$>VNPae=X zm-g`e0rgV%ZVXp2NDf5`Ss51e4k2KUbe~qP9RvQ_eFI36Ju3@rsuhoFwM4?QY=?+2 z&TpSeeBjqTMoyA9*j@Bgny!PeO(@;i&q1@`HzJPDp5PK31Dx zw-L2rUQ|TDA7+P8&vJCS*T+AjsoD7!s}7Zgo#OqwcU8AHyd-j5eiSe#mwaj_=GvXS_;Z^9m#&1XKTtVbcwHx>aa{J#VN_DV(goxXzZ`B?$K}=qUb+w=KmFxBy z#fWXdCFeY}VC9*d*#G+|*ei6~T-*qTUZc{Qys4<^uPhr~*y+ocmqSnxJ9cLrW70Z$ zFm@32Ta?d*oXh3nPoK-wj-8$9@b*x41}XP#{;#8B#dH(t4*GA2O8)Sx-*^$#d4W?4 zO38WF_X>om`-(}Z@{X*XxeZ5ES~bxtZVJs7bsyE793s`j+AFto92(*S&eJ7uPsg;@ zWb;1$Ha14wPq7$k18__?r_J}9&743eHP7-~X2JHGIZOtYiN-RE#sbEyEz*bFXf6vw zAhwUK%hJ;wQEG>n4N|#we0ynU;=WhqByxyqdE8Bn{LF$}^=`%R(VnFHTmn%W!%kz- zK=yBe;GRvw`a+7NOcmUOKSIiXjJ$Pf(B*DYn!mZKXaiFkdB33GTXj!dBEgnx5uzCE zit(BRa)!moTBn{lIyXc+H}wneZDaVBD|y?Tz2zqK7rSR3Zdmn(yRa(QIrw8STUA7sI85CcQQRN-vQ5n7gQRCVyoTgd%&+_dSpjrMM1DK*hK4!O^I9=HY@M4VjTxOcC?TeOJc%ZEnwzM%oMatHx zveRmJmEyw3^w>+1F%#=qMRf}|WHb}&0Z&A#CqGS@f zC`_OPHFs%`+Fo}bkyMD}Dg;l_o`pKT8T6=QO106-9`}Ie_iFadBsX}_*=zQBZlP1q z37tM4|F6~J=?3}mezyM7@DCdy50U^Yx{qm;km#yv)YMb={tBii9Rl|7wZ%%CCdK>SXCPQVS29h z!xGor*G*hM(l^7)UXuJ?jp2Yoy4meX?Dn4fQ_6H2N-?EfHDllk1w1WoBMIdJv0?7( z-V9N>w{Z%ZVX>36DkJ6eToPQKQ<~BC`~d{=*$m6wy0}<&Uakk`HCSS9vSK zA9Mx@<9D*}Q&t)Hr_E)%dqI8PMl? zE`F8!KIe#73@KrxsV`GM{d7!;S%N4imnbHx(hi!{!pyzGX9Yru&`~gaj#t#-LB1_YfN;QS& z1aYJ(LoBE?j|x-Ht!Mu_e6_gjDAWqEvtM!hL3dtRM|#$rOZ~o#Yj!U%_$`UzMbt}; z?kj;y=e>O?u5Kx>5$lEaVV6Z%D53}H2`zj<;@D5^H`x}~q;k-fQY-08YbWD_mVhyK zr)_|@|86`MjJ2TE)TBiZPasviAWHYHcZ~w;oCr}wG_5-f#$V+j8?T>c#x>~aSKSSZ z##TCk8o(1t_hn@tHa~Cwlt>2iqgRxw@?Hx8rhoGxWAq7 zDk@*8nW!DT5~K(m(qd9%4YIVQBCH}L{L%-Sm# zK=~fitJ$e?oTM{|V-WRKXG)^915aqBL}_}P|yiz-m<6IRoR z^7-^CI%IKu>yrE(L*i-G7f)t2{S2%AVEn}oLOvAp(Oh@5h#PH zvAFx>bx5?V6GfahlB({GPY3WS40#%cPuOlW4(-0@>PQ)N1wj^}q^twV&ZlUJPY5Uz zu|w9D-G1o-u_s(ia*F(sn`BAg&d%ZH>6n;_4bOwj)d>E~8F#%(riUk`u|mRBPns8@ zBDy~m^dhv31EJ-}TWyZx-}t7V*ixTGe8VBb?Fdfh)hH|&szG)8ApPr+-CGATzo5iGX=?7iMSRUzf&dxp~65P|WFrcGM# zKYDS0AfPI}TjNW@d?=C?=Ya-@9mYq_cj_lI^OE)~Jx>+kuzYx_|&_7D0a0RBCp zJ}(Mmm(=h-3%YxAt~@N`i0pkrCZCd@HRm@13%@d`3MyvxPGfE~EmqI0z~ESx6m}VE z1kDpaRIJG*bm{IJTT*b1w}VWAU!o+3vJ40`@i-?=Q+I_UnfGHU?&mOy$vVkv&rNTYn!B(QQf;y`Ztsp&ZqKI9+b^=#`!!H`U`Nk29)qPruSU_S0I zNz}>ikhH>c#l*1KTN3YVl35+&b7Al-7Iha9fURf_G57OL`z`;$R{>Lzr`5HVt02e%-hLCBGv$i#o7cD5*TnVDNOS)fsNTK(!U5WT zzmagEAWvOtDE;Ff4skIoVxcs$tP7GN!_WCI4k@_#QS5UcsH3(m8lGTozJmv;Wb{70 zb;22tTCyuo3NO$pQoDqvM4S*}5T7=jE?y}JjJA-LPVR17Z$5u2xghM`I+aW5mU0{NIwUf|umjr!~?z9=k(J{LmD+MOm z#%q3h&MGdEQ+_b^RvjIq3szT75xAXaCUsb3Z`Ze~q;}U9q!3eTP}IP?S(U8Jt6)?x z;Pk{CP`Z=~>0Qcb-G*NH`kYzbiH?2vf_}c>PRGZug<{eQT`(c^G4qX4dtPOI9neya z@OiJ-GuWXhp`^D#Ts!&fqkt)jgpYg?py}_Uk>KD|a_vlS;p5;+w(~2#%aaf8R?i-U zZfq%AYn#N%;D{cobwfIBY>`GXLF`nh)_O<5&lcjdjt$2cK^g-q>KUSc?;Yk zi^3Y8M$63yCe*J&iow(SOvTX>!lIScIaTJXg;7pQ3Gp=SW1=Z5Tm(SNRQ{hi>uM{s zMH*JD-j|2fO8mk3Msz(l>~@zok@(q*M0MAolbCUY6~CLj*|+AinQ8#r?_W3=Uo8@YqUwAgGZ*kZ{irUyk$xePZ}kCec)VoG?E)>dYE^}1RquY3QO0GL+5otZ!l)@ zGB9kFIq55p1#mdBbJyEYxT~OqUjQqT~REFVnG{F3ayBeWJwmjpx>bGe~ zXpIcPqxp5lbAakGtw#SD5pP)=S3z-4pS^_;oTGy9tBkG6*Y?SS+%1CMiso&aIY&m1 zM?O_MsKd6+nlDJ4NiFEO&}?{Auk{{1QiTi;PsTqv-gj-y$W!k*vrKseF1nyw4-i~q z2WFkIH%cvZW>Ca?l5#Pk*1*oRd|G@9>0H>yvYqEC9=Qa!E{p2!`tVeVCX`SRgeKFD z&^)l*XmhG}9-y{fuwE-Ole-R%y-;nK(Mx!wUM%(`;F}!MZ5eGZj+>V6?#- zQc+(YAFHjOeD5$052MF+oa^{TKEMBQYi8@U-PGSH4jSYLaRgGj9F=uZj+iblY`4#1 zqmq;nd^sh?>6`1`47L)VW71;?`ICpt*~_A=4-1$33Lo2<=uplfnkYSeDu0Df;aH42 zXg90nf5eL6D3K=J)F^gF^B5TPWsm-m3bz@5i)-W+ad4k|OpUlk)~&6x?b?gIFs)-w za@wL$P6}$F!UlfFQlE`(3OjI$53k(LrmD3~Yb;sWF{K$jGE__R*8$zM!GdEKYC%xo zSCiYRO+b-k(hLo#;Ki z#3p$=?uv*K&Uc(P&`G@ThNr89XZ==#$oe?55A11o!q02`mI7Nj!!a9KSrIyHMzU{Y6i?Nxt*9moB&macW|dN<*CPDKHtcgRR=8^~dx-Dl z(Bt&mQV+6GHa|n`G>^U-IKUMCbH0nR_~Xtb}TLUorU=UYgb! z-8#a0x(vft=;2W>i{RN#v}vwfWGPo@S>sM)YO7R#MaBoE&a>UheGP0M&R*^}WnwTi zmDyXFL`Q6X8Z7?;E#k(a&S13eY)GIiWSA4SbnL6QTXKzOe?@;O=Oy>`W%99TUzxH` z>t%p}=7V;#z|N-6uCs?oNQ*K(cv<(%=8E+|&KM%)bKa{L^V1Xa7V3A(c~1A#LTGp5 zij?6*4NPy+54R@oYtKdoHGZB*>cWED0$O4H?j1x2{(A!rW?bJV`pGqWNJ=ao%S84Al_4zS#tBNFEZX zxn05Zyo_OV1#Sh|Du_1_T(hs3&Rz9BmOg%b_NEkU(?w#xX=KQ89I_Dpis#zeu&{F4 zaBD_lZC%_Lx0x|y-%QOiM$c{F^}r9~zy=YrzdqaF6!+A+FRnEZ+KUe&?C6d1WDOz!M~wJ_u&pv9keQ)!jC)g! zi5yBpMmN#S8ku7xhmhM)NFj2_DP`Qwv>c{}Rmk- z;&h~eNfa=q4<)BfyeFOH@QuK5j^y0)Sjc2vaqo?4gWNv)*Jl21TFrF+vZtbz+d8E_ zr_Wc>q#wm~Q8$y90Lv~1QZ;Cwe_AC84*6+BSYY<~Fq+Q4Yz-XdMj7A-skQBaYrw}x z*oDwdNK2&`+Ri@QW^_HhbUl8rsJ*Drl$t%6?314yb2hOO)d@+ZEBezA6?4JzdQIpr3;WKz{_?#kQJ4kn<`@O0JK1;N>0;iaM zR%XGVzkXhf-+Oj&fnI!@_;zg?x9rSdB6jLh+1y$B;qUB@9`68wP5q@Q*eBnc0qcD9 z&8PelTJoWIQt{LXl6Wpltc3|?P3Ltz{76eyRoSZeG-1bE0lP6DRpE$LNoi5Dh3es5 zBf=QAd@TKelI3DL{rv5uIXIsU<@!deb_?3hmG-5Hy9Ym09*1n^brJ|?sExW|9QB-4 zkJ!RZX2=Z5hmW`;?Q*4U25K_3rYdA;90?VFw5jG++pHmGa|!GhQ}*^)Z=%81wwA-4 zGvV7~t+gbmu2SR*G$*@$D~=ICx|aIxyG<&7=CvzwFv3SCI%~NgL50wf5n$XHx_6SH zNDNFMcC-}L%{ggkoic1Rr%!B!3YfySQ4ycL1RfWEe_?$pR&V3P(VXRT4lwW64>839 zu(o%In`19h1N3`WZb`wBXDzF1Y$x`2l?-A3h!?+5RYYCcv)w8Rnj{3%SwsN7J8Y zC#}*@1;<%U(ZAW96AUdM?6h)x0==gYf+*aH7fSI|b1il#^zx@3=}bK-@ljuLp)V8Z z;VKFl(ENki;9yYU+VX}I(eRYK<(#-xw$VXzf*I=5-NOTi_N4^6r$2&7pf9xiZn+%Z z=VS1Vc-Qa-_KV@-_bRv`mKGe1(1aOxwQJf=rnX-VK$bLh4bNA-UWbcK;$Ou|$_%{Jvky+J&HVtm1AUfGMbPjs2( zsk248Dj+YjEdVVfIPZQX7^>2aUJsN5xQlmD7yAwY^hE;)NQ6OBJ)k!Y;*a#@)SnB1lnG2 zLMRAM$|b)Rnb)bs4&8zbdQlipb}tBVon3$U)W*D88(nGvdlY+J&NysRT~?~#t;A6w z^8&w>r#bzpfT6zyfQ^88vmtEm{=`rJqX_aJOvS-8*J0o~vPT!n^ByOCC=x=C`_NvC z1{UH@UOB(6v8+F!(ZhPSP zgN+Xolo?$aWU@OjGNU1*_4ZKxeoD!dLA?`PPfuQXvnN^W5l&*~?q#hb`^y>G#qHf8 zyO?`W9%Nq`L`TMVIrV+l%858XLE;4x@PgVbH92#c@>3SO!7%b-41*GsJWxB==LzXc zxjKbeS?krD+X=;V&_ly8S$P!S5156oF>Up}eQize);#vjYFq_|vSp!Q)KGl8YA?;!2D~JsdgjR%$-$B_EVCu&tqzSg~o$&`hbYjO^&+ZS|a%09oL&{GZ?X dFH5wes4*l3pUjfF2oOoTEG_KJ%W!THe*+|Jg!BLa literal 20719 zcmeFZWm6p6_cn|M3-0b3+}+(hxCHm0fuI9}OK^90cMIBnqlbpTNK%!Ib6I zB|xX(;NTDt5Rj0NP*70N(9oYgeS(33frW*IgM)*Ihetp_`26`ZA|fIZ5)v{pG71U` zDk>@(8X7t}ItB&?CMG5p78W)(HVzIBE-o$}9v(hEJ^=v%At50V5fL#lF$oC?DJdx# z85ub_IRym;B_$;l6%{o#H4P07EiEk_9o?5NU+C%S85kHC85xy4qobp%tE;D{r?0PXU|?WqXlP_)WNd8w?b|mK6BAQYQ!_I&b8~YG z3kyq2ODiiYYiny88yj0&TRS^DdwY8a2Z!(9zdJfQIypHxJ3G6$xVXByy1BXi`0>Nt z-QB~(V1IyxpMCN?%UE-o%UK0YBKAu%yADJdyAIXNXIB{elQ zEiElQJv}2MBQrBID=RBIJ3A*QCpR}YFE1}YKfj=$ps=v8sHmv8xVWUGq_niOtgNiO zyu6~KqO!8Gs;a8Gy1J&Orna`WuCA`WzP_QMp|P>Csi~>Cxw)mKrM0!St*x!Sy}hHO zqqDQ~_wV0bU0vPX-90@$y}iACeSQ7?{R0C7gM))ZLqo&E!y_Xjqobo^V`Jmv;}a7T zlarHEQ&ZE^(=#(Ov$L~vb93|a^9u_Li;Ig(OH0ek%PT7@tE;PPYisN4>l+&zo12?k zTU*=P+dDfuySux4dwcu)`v(UHhlhtpM@PrU$0sKzr>Cc9XJ_Z<=NA_jmzS4US6A29 z*Ecsex3{-{{`|SSySu-?e|UI!e0+R*dU}3-etCI$eSLj!KXp^45m!e)Dx}*_%8e;BO8pzYmeWibm$zv$2o!yCkYBhWoyz{CZ3!I_-EM8D-uSkDVmnq?a&t@b z?!a;=k&Pus`3)El`}GUG*(#g0naPw&LYCB%FM<%n-FfH~X)?J7M>9!Fu!@6*57GWR z!z-*yyh_<$MPtS`!;0|=;Slv~jDjtxY$2HDgX7maUNH4%`3n(fSrvF6C^}NPo6D`- z2zuzY5@A}|?lj?dNIK;(Rb_yxI3gyE#GEqVglW?p!480x4KZ&>)B^wC5f9vdM=Sp~ zx`(v_!5|D>`g^VE^-?f57+&4C#eu3GTF39>$b#a(BZM%1(!aqUfEQiJISNw*bjol0 zoL;CaP~$lup?$B^e?qKw%>RTva{mcG@&6~R>?QwC$h(IAznPWY4F3t~^!`uSMfsnQ zHT^#!V=wN1LbP<0|1D?iLH$n%LGXXVE~5X0{Ad3OC%Uoz6VluMCu9{~M*iCOK^tzh zs3ki42MIioa)F#8&c(uaNwd}B{|TYP{}aYlOaCW?L-8e5b2iZe&cadnoktO zEkS?>f#>O)gM6$iNQsm&={w|Y3qe?okc;XT%1ofe7_@s}bRWv_BMz4E zm&7cBXe7LAZ|OcvM{GXWfowe$ESo4>WDoN-OAzdS=!3}@H53O^g_@lpsVMMs%?bkeRTQRhxK=m{t zlm@Wjs1x@5FronlHSG!2va=DkC4cv8q}Z@PvG9`x+vnY-({;@clg54e%zdD&xscyj z*P{#>hACHl`nnb2L(X2xe%V4)?%%&w{G*9vlJ1IpDqr@uvef|)vDb8!ZPVODx+=T9 z4mIfe*_~w$n^=4QeX*1({_!5cOja@2lvJ3LxY`*p-$ z>{lKPAn?q<|A)9(SzAqwY+rQ}p}fC0XqYtbh~=3Tft2`2zn@ z6(r#rsmmXMz>If%I_2|PA^A=_%s<=QJ*Dst1H}?gG~eoeqTeYUYH~7H15SR;bApPq z&6b{Yvk|TN?-14=OO{bQxmSNQFT^-69}UJ>+@aum%Q%yk9jJ60Mq^Q)8Zo;t=Z$qj zxlHDA>9yrc5$2!?(rjOoJb0!J{+;w^=f-X{hk6ioU@sW0Tbf?tPeq3** zqYQ}hHF6Bp9~sXfK~o^4^z^9ut}P;Y`!LBk%NVD?Q7jXifC=l2|4u3azp<0cF57|h z5&Z|qR^$Wd&Nd*CGAF5%kDc|ycXjHBEi5tYei$#9oAx}TzDKk)4ea7!r$2MF2nrE% zL!icdt{V6~K<_g_Ouu0chVMeLm<#v^kgK8;&80;BkSzHJH*#h|PmK~&GS}c)ITRn0 ztQZ|&?4Aeh{_SYf_@S)e6&|eGizr&=uch)_VM5cXe;t@m>$D?^CilQKfH1Z;l$LEpP1C?Q zioj)~THZV5?QAPrsNgb+E3_(jDxbk4!llVcyXs+9JgnX1%{iE>=)oHofL=E17!STKrL>`r?|oT z5m5Km-+^7sTJDM=J|m6JV1vqQM*pCIvraE9lyy<7e&!#8(d3%{jzL6C9tol>WCHV( z$^W4#0@_B%M*6Wl=A!^8gTB9|SWT&3gE#q|2~gm4?^%&zLi0XAi)zCRxDfiPi~SRp zWBwaEHbfp)dIJtpJ_}e$6doUG2>=poROWX>{C|>>re~0^x4C5pl2PmF3Y1J34V0W^ zUg5vlIO$xj%mTGf;{!o?5?>qVkY+y;1Hs-K^`Q>H*=O#W$eT4;i8A%(UMTqw6f6j8 zyuD-xQgM&c6B$lY*#28xs?9aA=EAgXgaXL-9ysDl^>C^%WbX9II8X`o2!Rp+5#ZD2 zO6V-|Zyw=n?`Bx{$B27mAuoGMGsDLTY-wO0(S0$R{x_68s-1RWuBAJWZdePp1rnzi zWUrs32?0OF37W^U|7B)Qw9~E@Q7ZAWLIgPPNqWsj%~=*h>{7(}Pwn)N@2yfTbEtcD z1)-TfaCXy3KO%zaUEzSvu#g=@Y(a8~*Q&T@6#Jht^5G3J|82Q{#=S459p^tq{Oc%* z$x~yfhbjhfHT^`2kj-3o@D)CZD-n&0li0U9yuaCP4aP{A#cpsUYf}R)3mRTs5vfoQ zXr#kRos^*RW{+~g9Qi1aR}r6Mk8*sef|TJq_VAx`yqUr-SIYX?2y3Irc!2-$gyS>Z z5IO!;_ZoaxbzbS$I06x{3P)8cTX5OZ^h3#Y0I`#u3BL`kppiO)|9c!sy)pH|E%}#s*)*FE3dyD|3NAH4k1UJUJ4bk6{gsnwz%OCGV%V3Bt~iJqu+$s(1l$qwJ}lMFv5H6y{ZN3 zTNk?EPy^z_AQV_ z=M6aXE%6jBWW$H~4<6HCTVY;<)lSSm4ud$FcMxb17f)$PS65?k3GOSFA@K*&!X3ah zTwB2KDOg?IOE>`fA1jzsmY484N&oR(#2~=JUCk;*el=Se39w)&djES)zz01gGV=P% zjKK+gE-^Q*K!4heoaxG^K*xAsRY&aB4n3^kZi<3^T~dArh#{dNlV>OgH0?pQE>mc@pHu)NQ)+c8E#^yyKSQ`!ZcSHa$Vce}vu~)#J+LiA zzT@@VmOBYio@=ONoBQT8YJkuu#Gv_GD#jt(LsQ9TIZ5bLuE7%;@oYHbDmc!9+3*qC ze5E~)anx<8Pziyqzm7@Vc3=JKlfDA}E(xc5_8;X4x~6h3mDk;5bqqre;jv^S zrXYh2IvyM7I~TNKs)F!_fx$SR5d$n7lWz?(@Y`=GY?-Ds3WnBzX7hfNaY0$&N$J|~ zh`2@l(hjOdXIY9OCF2rkI9m9dug5vUWp`*-dj|%KR|CGe^17E!ETCoF<^|dXU+k_C zm-v1e|K-bHFGav0P{TsCC^roW){8-t6^$8^UcMhBuVXnY-Sc?=A$>x+HLQT3*N^}E zn$@m@PW+2t_TI2w7RKc~X%?X47Uyw|>Kt@Cq3bE0EIe)B7D}q9%35@2{U_*q`0sf( z{D2AB%m?3KCsCh6d3aWlh)xr;m7kh)6_0*XYzyMZZW*^kOMJ+W8)WC2HlLsOySBb> zMd37KAGROAFXONq3|*!>Zp*pv3g{G8Sg6MDK!#=1c}&(Rj9lY=EL$fz*Obau(sZRU zw59Y0GA^xjkBYTARpC(;ae$LxT0s+7a-DEK3gXiK6c|z> ztd4Cd##9}GU;Xe~YUyF;?5`nQy(tX1l0AwGlEa$W8+hwihJr^wzb|h8SPIp;rBG}q6jMBV@RVe|>3Nm-n_i7!JRwa% z0%mx#c-U%Qz;*l|-$DW>EEzLuJKz z!kBrWdE6lg1E)nabO67OS*ZS%i-*cU4pG3pj=90l~V#%pQ(kEy^e?Vd7n*sN>)b2d->eM zXyBrXFgFSmA~xh__tSK{R6%f%V;3Y*M^K8L0e)hwi%7%nn)?FvNmjPAjDAyh?Fgc3 z_~ZfQ5O2^OEP=8`wv*s$Dt{Wp{*bCRy@xhGR_?t%@9f~Kqh0J(-Kh7Y&pKpR$eWzj z_W&z#Jy27Zm}gO*?jcC2Q^TZ9K2JiZRbLm7^Qqi((ZSD2eTp1YlgZ8yCYXjx_X5E8 zWV{Bd-MPk&%ULqaE|u01=bX|?0SD&a?fv21Aald3qi9GXS@ zaoJgwC;@I1MNHV(?+{Y8Vx1nDO$*k;2e|d|R<_FxBsr&!FhI$LC+nzJUhruFF{$Ps zCQU^ETCRJ3hgU&WE==g>t1qop3F&HdaP^QOy-#UzN?kFf;SQNRi~4}?kHGR#!*C*K zSp8p1r0n$r7GKnx2~ig9UBojP&+ai^zCVqnqzosBNdQvq#e_c_YGNiptswPnzm5+} zZ4CkYxZax=2A7YC{XsCd2KmNIBKW{h8gJ3~PfZXz&UH~klneuo5Xm|eO8RJ_`zGob zmND&SB)Ol%L4`vyRhG8u$L(dqNYW{c-stcK(@t0_*^_AoFRN|<05ijTO@TIMa zuo{U$oeR(ruT5Y=P*2d&BoS{NB#?DlYb<8u>qV`2@sIu{w_<7XrAqQ!SX1-K5JD4TM&h<4a!p!I@6?iX%{5pJ z6w4wY7D@SOu4+4Y6g!;P?NlD^J$TD_rHz~w>r}v$YcP+fz({9}6yQ!oEK5e=5qgGX ziTtd0KhFy@#1Na~0MwOQC$eL_>0s;4Ys{~>e+yGiVyz-jR%XzWSgIoT@neU*XJ!r^A$HeXlRPPE75=Jy+?;INEQRjwBHjg_ z>C_HG!2=q+}~fYbBd_ zvE4DYVzCzM*nA4W$;*sf0EHJ;K6SlqBY`XCxr9*}e@9DGLmkCZuo(K`9`?i|@AIS| zIYfxV#bbRI*|J`MPMMEq0Vlz^{86h8=*b8D9X}I1E%Gefi)$Y?=L}@dTM&!8^I#F3 zqj-6%a9psXs~Q%_w9jVS6O2-J=?O}i5K?~7n`<`;X~6G~;*tIH&Up9{&SWn7;3>sc z!$piV?BG;r_opF~w)rrkFY?J0F-=}6{W?7A*hkL4qRKdp*lR%RpuSlIuW(EAIFn||Q$&2>YQwi1TRRO;7?x5oz1Evg|&F^=Ei0ibz z1fM`;Ao2Q0*EU81(^iEkG`ulxa3BKtp6*2T9Zyz9m0f@Crw=^G5yI+YcmIe-d+kFE zl+dv{aH&KkuC?fM(`Bh6UPPGJmiY)lhc*Ap@IQGmH`V~9w+Yc_y+gKx8r{z1+Skw{ za5aC0rk)ucKda?XN|SPEBb{V0DImg{ky%JDG)ls z4i|Q4t|tl+DIwcoj80#5>FdvZb1ay0R_=Vigs5#e)Io1m3_vc9oqT(xnfaaUhpS*% z6GH@R=WKqHrANQw4qlRV0x{np_2!xO?P;q=hh*S2G=%RHgyMWAGEs6AC^eQde|L;o zGgOngWPEl|X?gw8p%U*?3WqVzp+Z#08kuN{-|aaF7AnDY z4in*qe@IKEGYH~HU7HJkJKVJC;lPeTdsp;Zq#@lTB*r+$++@x&qSeg|aOrUihe0Rs zuGG5a$@j5F_|t>qs$imQkp14ASQ{9D8FtDN0!%xg*d!%%q;kPQisgqiTX}PPE~^2u zd+=i!LrGZxnKz)EiW#TI3f9|XVIAj(1oxqzyVzFH z6qz8u^@W8rXNK?b<&!H0Ms|tB7Liu=I~4ss86Ah!m+E%L*mtQT33IJHq1NjG_}0$C z+Ws+fnoDP#4r;zCZ54(Q&eq;Q({vMmwKojvhH>$|7rB4G#DrFl{~_*-9ReR7vnNf>_`ERoqwXg;Ew#~p z9=C~+zr-Z%J7@H^Wl5b&^M7k;)Tl0R$e^~ zRBpE;M`6r-QE|fqxPr@{9S_G6QpJ5>#_b=Ks9~iJsr1Xkfk2R&>Y z4_$ya2IAq;YB(c#?ZQe)#h<^M3ZkrQOY~8Gof&D0;zr|u8_XtT3z3EwIbwJ@Qu83E zWpnj)_k&=lgvu7LyC>0WK6?YXZxhBg>nugy7D4`z*bL_f*@c-n+p zQ+i9q(@&3w{xDYix>5vTnPFmHQFy^g3_;3m+VUnh=8&T1CmY|9F*!T*pb|Jv;Cli` zHz1G0dQQ;B{;b6m+xK~p&O<-vXluB1qeXiKaLf#zfzR#rfCgv1`k<0=0gDUIilee4} zFzrblno)qUMkszpff)96pn|Lj0g;O|@P`Q0dJ=_w*WOF5G$8k(TN_l#tJnDShhW7E zhG=X9B~j!uPyIT~4z}{Sg3H4LpULOGlgw&4rqPnJxl%U+ZhBLU)G;Of_W&{^_U4Bz zT136_ljpxo*7Wz3#j?Y_K5zfbGftIJ_W)E@EnqlZX2Pc;(mo!#pdw%}sZ#;hLvJ0S z`oY_B`{7B|v)*SxR4Y-7cp z%8gT5eFRNkzYKHfGM(Nuq+B%}o>QVAFx7fV+upd;v~H7yq?kn1bxHni|A!U|d}Tx6 z=nx@NaCR#*;Z~*Mc>K1GTV^<=S;}#X0Nn9$zE3Fvdp3)x`J&2?@ipfJ?4uFh^Df6h zagv`1VIk-_12lTYB|(JyNk%Yt-i(2e2qfwO3MRVtMU=*rb(nhMjMYLI1zr48id5XL z&8ZoO-y8$^-Rt-qzx%2R%|lz)grOeMAw-_38Dx*1xFRiuVmMf_a#D{$2U|Pn0&^_fx z5_c2g0Pd+ZqLqn25o9m;vbV~rIYa>|>&QL>hKYG1V=p%|!W%nIR0Z)j6PB`WTQIb( z0#5FHe2d-*Xh*4{riI+Iuhkh3vRFieirGx+Xo*j`)_K->># zB#E?Fnrf}=IEUg^q?h*#x;>RC4JQW+!oZQU26FlaE`Zw5gsnc5HOz|$D-mYM?Q*Er z!Gf_hmc0z#dd3IvhR^NwUaH9@d%cneV8S65DR2Q9CYLQSgd4}?e%u)#Ol@dr2dTJJ zT^FJYZV}45YPx?=V;9?b zR@#r`whJF};zNMj1Vc>Y%88%-s(aiKjnJ6_2fY&~E;_x4XKbu&QB`!zeG@V3hP@$N}&wtmte=g6l!`X-CE!&5M)e*o6nb0Ep zX4pP|$lq5CO<3K!pm}D{!Z1hcO9aeag4tVwjAGhv;{=iq%~+{e%|Y6Or7u)u{N}M^ zROn+ZgWJ{Y$T`Vp|H!`9_6Uol7!wN1kBZ4ywx$P`6TpgF(4*E+lRm(#T0vxZX`B;2 zHgje4cC4$YdNHupkM&%{(NKcfv%)Yz8t=G_pWPG{Y=tZc0< zgi`|Q%NAd?U$_c)p?RNbpM$SQwhxlm>=$4$Eu6991 zb9wBUKx9-f>^8KvF35k=Mw*AF(c z>HQg2z7HL|i$nsQav}V{*E(^g5%@nY0JscI(bwD!Nr*yt^{{iy^h`hVC%-uh$KE>k zwU(N64+8`s#M0WI7^z#-_G9bb3HpGpYdFh%aMpUEY&AAB?j3r}sXZhr?z)4pR_9O; zMbHkiH6Eil-QOxdt#|_6W6*UOEEK^r8F}u=>$R>Sc2u9FmHosouq4UxIRG(1o^bqc zi(yD46NBqm+?K*ugvNJPc)Ylffm9MT>~2{PLZ>8i_|*6TZfM4;rfe3_{aMpvtj@Kg zhGw=D!+^^w66V2=Ne*HBqWu^K8T&<}Xi(KmJ$@pZsVKAEGvcppV@lKJ*>V5-a-(n4 zvp`JMFBLC=&Jb5Fvc{F!CphYi_y$(rb4bNRC8#~3)TDcEMTui{~Subu@obX z?vM{z)I0ouvn_a(sItP?$by?S$WPg+tnwH*i^xxYDTgiMhQ4$Kug&)Ko3K*(ISk8qaM*_s&2^T&Q?;D8>MZ{Jg{(?%oUE~`TJ zt|H{%K!S1G-wc0g_k*#>#=-RG^rk<8xJ zvb7!RR5Y`7bI$zIri^s>vL?l$>{_k@2TkB7R~hTdLtn z5oLX?8u_76tNlD6SgB@K47LB#2Wm>x@9Ynfe>r_*3Ty;5M5kc-cBAI@Oorrq%^}1H zE|2RvaH!y?uqXt>eHlwb=nC2hWjj!o+kT-{FyqL@-ko7B+bG*cYP?n``@e3cT_%2^$?<;A+D;Pn_sdvP-;2EVjZcfhdj zMzm@*Cc`oi+9@_0{8!Vl70qMawdFO&gHH{8HyU)O_Zi7VxmQq3EicX5U9wYqhQyQ=1zXIeE7^YN7#f7m zHl)V6)b$K$Pye$9$GDc4<@=h+!w8~H?s^Pt0pnal!Ln^Ji9)uJCD<+tbHUq%IH)5r zxib3c?V_9rH5<1pV^f=H&lUH%7gJ9SOZ|{~lx$(<$cj;m$U?;*$91Qr5bLxN9V&*& zoIjbDgxN*0Cez}LxJxK_8ipzK+?EfGo4HRUy+dykKJ>qadX>wMw<*?#te@zai_NPm zD>_b(5kUQ$jLkmbutTtIJo0H2vxVTNWT%Z+-;}2ThxNS3<|Vm!r9wW!f@V^ouOhJ) zm?w_XS+h{ooYEOY1~fC=FSQd5*=f*D*TN@GwYN50?KpkA1J!^db+d19jDg}8CH{%Z zp3@TTQ>Ekd@eeQSuO`?Ky6}leCJG8xZtLIy$#<348ORcG&!PI+v*3kWB_X^TEs2`_ zUG0S`0kt;oN9o4i6i$7d81+M6`$kgl@{ya-6luQWrA<8~I;@*6u;n!v9W!grCsHw; za{;w4@fP=`F-o$Qpn0K<@HzwLOKIoU)wB=Te_bH=mGBcX=M(oJi#w-c*WS=uS`PAK zlfRdt0+XAM!UXIe3)GDqr1w_DTaC;Jt?Y+7#ZvAnS=JSS{}_i3e%OvXTOld1BU=%Y zyiV501=>gM1$@dzE2=>RUognr8sc}$?pY**!*rMXdhS*6y;t0=6QGwSjZI%gt}M|3*@jQ#L+7t8_?f3roW*LsCR_;kz`b8|d|TP+uI28PV5A$`8q@&aOFs|SQA zx|jzmlBUwYr^h(M9p!3x$*odkK2lxY$&e;kvKW$!Jsl21?wN^9^$}z#{GIl2bUsb1 z#PW|QM=XS9(C>imOZ-HN|B${?!ZqI21$g<$S9<-%p9b!C4t=U%!m`ukP4E!=h;bsj zQF=1U2@mTdKi;gFI1^Q=EK&`E3y^ZGx2}2EF7)&>**Be4r*k`51eY5hG38f55M&K< zH5k)E`0#_~n8;HlszO=@vb{=Jh9=b_)eLkv9#0KN>dy!|A1p*$q4LQ)(R`n8*hk7R zS|QlhUtCd7cJk@fRy8d5F-WI|Up|YO{n-Akrj?l+v^84U{jw3Z%CyPR?s~GkQ_-8s z>@d_2>eAh#J|G3-=Xk_8LLegrltFv*JW8DZcU=1JTLYf;p5%uFeo@A8(bJLS3vEu1loT2o7JEOgF0WX`AN+;BA0HiBGf%|j*T6XN%OJqX5P4Inkd}8{E7{TGGUnaRqLgXVI(Xk!` zNwmJmO-o!g=z<+;-=GxeO0$ywFQ*3_O$iqE+945@w?*$J_GL2pE*z9k5`we%lcJ3}N<-5RR z>!%EOwX|o;;v2eU7R?w;;tj&f&qwx>q_<#n1Kh(0Cb)|e3Wa(Pq$%hjqeJfUeXSxq z1GSH~C<#~1oig`yPYU7WJ@TY2*Qzzer%F=DxF|_C+9A-77G#M&)6jLg8111qXT5Ct zDpJakD{nBrN5nq3aAQI&yG|EBlfSRalGd|Ff!uAKD#{bPtT= zm6SMf4{1QxyD>}sJ*UKOtoCF~Pq2QC3$`;gMo=R~V6R|BGBD8KbE4z}UdW}7V?yS& zEmkNGNMy$qDSYEU9QaYWlK?w$ebyT@D5kpPIkk86=Gd{mRIKB)4WZ9mVR!X3`B-Pi z#o!(VI zgM_7w@yiC#n%Gq7Zn7qLfwj z8bEoBlbos3isXn>07;QFIk9`q-CSdc30slgj?2|5P$E_Ag{`wax?o;# zALHuD7;%`$H`TT^Ztr|e;F6IXS@>%er%a?e1zdvG@l2Qo2O{-7uU8urREzns$6`Se zQWrku*KCrl$F-s$21{CMW&Q0gLDER8X>u#V&X3N6X!^XUo&{xVW$dIlL10Zk)tJK* zabaNmLZzo$+~WrFq1GsPc!(!q3UqO=sNq~B)jL^vz|#G}3JVq6!fH-brqbHw?DJfnJ#g0@Zf(?kKycc<>})ICMi1WE-$PlP|h($K^bY&dhW&fBe*l zgb{DURFL)7o{g|2?w}o%?WOnxzerrzXI?75;L9SUVnG~+FR=aObBhULIZ9tCn{qFX z3i^YOyE4;Z(8d{2ebY`R;o`p8K)h+m`RCh4v>MR?%A}$A(3OD(SGm7Wvr*pOscE{> zjiHY!ge$1@J?aKofQrDw-;;^zw6IL|HX#d^ES`R@T&j-p`pVvfsSq-jQ-0URcI|ks zv7e79mVT4+j2-^I+^Cv|Fna&H7vaPse)UzL;xZb~CTZu?z0YxYak#5sm|xw9IS16V zbgarY;ed*P`<_Rmd%6!ngDV{3HtZ=vef7a@!V8oRqHNg&c*%ks?PA@9>q}pX-tc34 z+Ycc(-3XIpPh-yLMxRUqo}ArIm`m=!Uazy{boN ziN@x;IJj~c{b;1Md6Axn7QB}o)23`7SSm$K(0=3v4#Ze?Wkbj42dq4t4z8uzZ6cCk zU^)Bh)p|Yi8VMvI4}AZtHMyYE#GW>SkI$HCHHBA{bR_uvWv{dW4;mfIT!KlN^XybL zdo~?Ej!D0-hbs5wnOJ?_Y$Lxm|+}0=o#5SvsKEn2mm?j zm^c2C%AXH4RTS{GR@;#4R*kIUnbq&c_oTSXPJ2+1|5NKFO0br)|&bFg@m zwIZ7l?mPS9$hb4Qy%ptN_MJyhw38l+jgWHr!qxasaHa8eccK!ZKM=omk4UqIAyGrg zgd1erUiI9Gx;dyAd%4yk#q-atB4E*vip(S48W}mJMAV?7!*-}XYM1s~omD>Pb&7Ef z-2K{cG-78%ba-==|N2xpF=@>`wnon2DK`-u9{}e&Jaw~aBVZDtEASeK%f3}nVf5X=rPwx{ymiXBG)cHF(L2hv;+?+!c zPgb9seHt)P`n;BWItwLgT+=>@kH^3dcD9SoL#)dA+wTQe7TtE!EuFik=~HF};lg1; zfYoPGFidhcKrMLv53h2eTnc`1mZ;58`v9CTKI&_Y5JCOsn+2?hf_vl|%~cbcA=8{9 z*eE}@LVqIv<~6e9NvDJ5((%p|rD5FvlZ#o-4}bRgka7YxRY^GvMi;kJD>suI28J@P zYT;7q&V4(+42pCRN)qIk8p_(F!xNKrGi(k$WT-roA5PqI;(7y;_bD5-%qHJlOW>@(<^Y{%!jSRGQ+p#&@_yc;~mdq%K{uFg`(0!h#o! z*h;_b*Bwt!EJ^x)3YR&T6j5EjWSMRA{^(lSz?l0X=J)yf8QU`F(;N0n} z4R16P6};Ev`1_@G%6AqgsTzcMZ+f3&%Ns~U0+JTCmH(j`dyi=xup~UVmo@U5_ zJ&b`?$fA&a>CrAVLHg^?@F^*qfifRGyYwRnu2kQQMe)wEw=~l4M`+dd6c6pqJ|}h- zUx;kYv<#3VPiOLEp{(uw3LCt)j*rEqi+-WD38$R39eLa58Ct4(zLMrMisO5|HP}S7 z)<^4?M)B5S4{&6B@%{nch%!1M(d#WpKm34_Q*|KE*)MZfx6VJ_StI&}-qy6c^@qbq zKZV#F=H@pA|7G46@B_ogYAlxNpR+4&CB*)!^3&ysKROT&c;zl=+yGeS2YbI$@lA6I zVq*?<+py}UXJ{%#8+r39&um2>LwlKX)-jITOseCx=D$GLUm3$Gnk8F>ux(|Rlq}n4 z7H%fZCtX_K`tYv(Mm#GW7g}ZsIZTHYy0dkBkn++BPxgCl8RxC1@h{cRc~Y={!eXIN zx9bbh&z28G8)mB^73D*}?&w@sXDEO>ty#W0GmwV#Jnz!M{DHb>buf^#`IbkR zUvmvus?~k0*6i^ESSpNTdEexLUY-UYPVj?mmigeg)*yRGmTdTi$bD?inF4Pn2}h*D zGi-R_obKl}nf$2CKm*QCNLs6w-wZidK?ls>F^!kCHe}VV+L1nXBy)DP&+P@UKf!GQ zsSlhp3k|}}niuq35mkLi0#exwG73R0TONJ$4UAZk;Lrxvpo@>78;Yp|V2WA}c=7Xr zIO1m{3SI3E6a3ImwNH92pV$V>UWycg-B851^RxM^y0uM`qn(n9%JW!7Uy{9!gHrO* z&Jy1Z&O$aI`1W4D>Nw#r^c=(dHCNjw_Zpx3k?QA0&6$4WXYqaUD<8HV7Evr;n0~hH zIQw3NnwvBctpgUp403^j%bQf&La~g&@~_Fu0lQ=lfmIPWVssyaTnxV`&gg?*JYIyh}6m&ym29g?`xJ z=3o@G2`hQ5xcQ!UYtLhf9!K+BoGE)zLF+n)*!NC3jI19%~x-eNcHWYbox9s|Cvt(K10P90Bnr56`z%DSZI8ohz>ozq?;C zU7hNb@&ZT~66zz8xl*GiG2Tgp;t5E=tAdvC zjoTf<t=p{_hD zkR|MPqs*$)sMFk%3ywcKJ;~b4Y?U!382T$2*^LI*HNRp`<({mw-U2VT7x-Lm8-&Hz zILAy!;`(xFE_P*?xxhIa@=KkTIy@)UW(x6lK7;ut#bn#dZj^XI4~&yx4TAbJw>p7b z&V`dp?q0fi@$-rkM}tk|&vwqHcn@k^Y)^VY3>rl4%k2YPMW@Ur#qGS#E0_b(q-H*e z6Y}wRf@d?}1aRmfL9;s%%`0=gk#~C^ROJJv2(=8zi9y6ki z9&I)JS~J1gEyeuQo+QKcC~xcof@4lNe;pOZI50 z>RtuG@F}MmN3t79w{=bDf-znzX!nzkfxO4fSBya=QdhlD>*Nuvf92726?sdP+w3)3 zG7a@+afJjc*pL{-iQ^-JdSg~4V9tIv@zo&E;K+g6-xyob-bG**)sODvgoYpSy=H2= z_tuL94b282B?k|gk78jQqM@4(*A0vcX%s8bHr~gDeYC_tQg2LtIw64v-b{j z7Lnu1;%43EBW1ZyD}?kEZ)2cVu{y27yFTBdRdEJG@@o4uPas77Hg2eJUZD1(GP;it zZyR`AOiS#0q~k$t%+#PNrm{!x%dfQ&@3p8SVi<$K{_Ru1dr$Urq!ji(0TAeXMFk#Vvl-BNz&y3SriRmf?k~eF{8C0HL8~p! zICUmu@wISb4SliFB$LPF;hGev&6iO}ox|_b$88ybjrCCLypwlXtG1dzEs~_{F<1aY zlos1yFSXoWZnCM|>$9vN*t9&Xvf5OExqs_QgyI%I&_u|MNHgY>J zTJ?JF-orrpug@=)b#{KT+B7j){ZTR}fKMFnf`X~0nNKo%))OtZK~@%-))NM9qvkiK zZU_fwynC_@Kld&o_r|9ah zlY_SozO8hbu)XipS{|GC_)ia$=cHEJOUjLAGw65pudj#mJnzvzsm{t(M_X%R-`5Ae zD58)$0n9yD%T1)*?n6h}B5QBjMde^CO~g!9K!2#={qes5JG>V9cLTSSX$jlzkJ>^n zs=Gsap(C~cpZ)zWW|>F=m2cnjF4C40Hs^m63d&}7?dqp1@{FEM4vqajX%R5vPwYUX z&zu{2`@wA(t_^#h$iMdP2E3K|I@u%ZTL?SU!K%kAwZ=fC^z#u(Jv~%!58{o6kT1(j z&IN1aH=?$87D zDTdE$tKo$Z>xkE18}CuI8yBb5HvZReRkQ!s&Yi_W;k6F{uaPD{YqFF?BYRmR6hq1` zyHsLQF?J2gIy2EoqTghv!q}HW7>2=&eHWqZ`_3Qon@b=z*?m}nsU5! zTw0niFuyIYK%Tbs_h|1=q&)ZqV22uKk|3w-1#DB&OVif6L7IirMW2d{ul=qikM}(3 z2u$q7HI&>%X1QGpJxjLF_>ugsGa9dFF| z0rayArIct&e%E>lP71{`Aqk-ygMD-8aK0aK1}&x&PI|r2-F!jossjCF8_APNJqo{m zjAGJw@2<##;Lq1m6A}i$WCp+R`e_RFB@3*@>Za-s-Knx#X1|soUbc_E&n=y=vhD4W zc;*08rZwZ4kuH-gky<}+KVbp9UVdvu!w@xY-M^Me58TVO`c~C7rsz$Ab^ZaV%lp0G zqR%S<@pv=ql_so90+br1dN^KacmN~*vO88%=kFPpg{__(Tyt_SYjlw-=Iv3{2JSa^ z=t~noi#B(EFse9xZ~CRC;QGr%9o!s=FFJ^4bZX87L+tlv5S?%D<=gh7PkP@Rx zB&Cyz(&2UFjd?7-O#@2^UU|y4`Pdx%6)22+K-BUTN?sJ+A6r+!Ux9z3dc;KY2IwKM zncTVnos7xq<+k=%>FRx~H6OY>PcN&MzI!mqgX#!g1X=txwI|I-aXueH&dbb>v7LF4 zTd8G9>fBvk&KQPw568lcG{!q)etU?$3d7jx!xZGAA9Ux{BHuVfiCrBy-%IlnVzioU zyNR)e7ByY!Xn_}xmLx(t#^=DKto?+~Ke#!&gH;3i<T}L+GZ5_9{1B}E?nkf@ZudX1syXCY2M9Ni*1;Skes zb%a((({At(|H=UK#UbU&AIeC8$DU)&VR`YS^z1DngDyd_&#WR27m~mC5qw^pl(DY+ zL(rg=XY5nfGpaoQv@^gxS2|E*f%vg{ksJ4!l;KJCOa92RIn`4NF5AAkZbAOrf}tl^3;IB)=()?Y>0qIKTDy3$-YxX zA*EBv2X4k=Lh8q-HDcYa+wCz`rB;&iU}afL1NgCuqP)WjY${vplRE{uT&wZHLSE=m zLR_4O;T?q5y=1b1@|FrHKRjd0(}*?oCD~@Mm~xT&6@?mN z%z)cVeYVc@QjbmuzI!eU8JQZo(j_gM9hdv$`(vS$c9?k9uNFaXK+{l8wc6@^!d=qmC8=s#myvE)j1r)5t;!5rjL7{`B zlSdDyTBVV@0F`}*gSFAAo)(abc6&{29{N^Qsjli|iu!Y`k^p}~=g&}kA0WD>GYBuh z>IzUfbrc?L@}1>h5a(Zb?9_zcPyGq6U+rC9wZRRXwS27D)BaPsvkFpv3Km$a)})SE zs1v@>*v}xSelq`-=^gyjd!e`ZG{3EW&III^yN|_s-;Rfr)Vp|cI>%P3Ghby#o_6X( zmr=dyY;7(T(Ixt9(jI~S6HdR6`2c>IrwDHS1eJ7aV3pct#;vpAjC!CFrIC# zVA}amxW#XK8kwi@`bQXfE_t&vO$T`OU)9tKOX$@DQSh6O7DVl8=)kU7ke_Qqn4H$d z&6T=MYL+QDT50yXgGrW=(zguq?(S5v4v^O=d5zRxC%(2=8dra6I^m%aj~FBV`Bu}q z@L{F6F{3dG8R<<)KeQBse+${$j99{OW~n04BFJeZb*aU5eO_10lwni8Yi}m zBSYl~q2h~1afP$Ye+b+9ju`*IC?_+fD;*;$oVq|W`ppfjl>gvBM z>~y4j(!1f|N<9|mHyT63Q8juV05y{n2JG^jg0kW)-+^1+3d}-FMC>r( zED=SW<$kdA8!xMf-^T@ueqsR2EMKRTrG)@lec=O7TGlPIsVhgb@Y^l#LCk3*$rBzn zB@clB`u;YN`i8x|0`QnE+MJ2Td>6L|0`jd!qW^2ODVT4q92TVzYu)?i3XaP;yN7r2Cgu{?Pzp}vENbv>>rsH|eIbc~?9hz!9r>x%YHu6^FB7|h8G9m1qa-mFiq2rbKdp~OPH(;TN^ppse?{X&Sd@r8wcQm{ z;Yhk!QQW7fl07(b`q-9To!+j*pp1GpBPa$=4kes&?j_DG=51tm*)>=`WhI?FEhVOD z4i9y&O_%2dg|i-KsYaRYi{sC%u_guge=7xR7Kh5J3GN7_vhq#d=I8!|v`2-0V9aTf zl~{Km;@%_Rjo{{#Wa9Ru8~8Zi-f4ehiV-bv+u>)$PN!c&jaH9ka}`Tql?Hb#`txR! z#w1aM
  • Articles +
    • +
  • + News
    • @@ -46,7 +49,8 @@

      8 August 2016

      Introduction

      -

      Implementation of the framework for hydrological modelling FUSE described in Clark et al. (2008) and based on the Fortran code provided by M. Clark in 2011. The package consists of two modules: Soil Moisture Accounting module (fusesma.sim) and Gamma routing module (fuserouting.sim). It also contains default parameter ranges (fusesma.ranges and fuserouting.ranges) and three data objects: fuse_hydrological_timeseries (sample input dataset), parameters (sample parameters) and modlist (list of FUSE model structures).

      +

      The fuse modelling framework takes as input rainfall and potential evapotranspiration time series (areal averages over the river catchment area) and returns a simulated time series of river discharges. It can be used to understand the variability of expected hydrological responses based on model structures.

      +

      This package is an implementation of the framework for hydrological modelling FUSE described in Clark et al. (2008) and based on the Fortran code provided by M. Clark in 2011. The package consists of two modules: Soil Moisture Accounting module (fusesma.sim) and Gamma routing module (fuserouting.sim). It also contains default parameter ranges (fusesma.ranges and fuserouting.ranges) and three data objects: fuse_hydrological_timeseries (sample input dataset), parameters (sample parameters) and modlist (list of FUSE model structures).

      Dependencies

      The fuse package, as well as the examples in this vignette, depend on a number of CRAN packages. Check for missing dependencies and install them:

      diff --git a/docs/articles/index.html b/docs/articles/index.html index ccc3942..7c80438 100644 --- a/docs/articles/index.html +++ b/docs/articles/index.html @@ -53,6 +53,9 @@
    • Articles +
      • +
    • + News
    diff --git a/docs/authors.html b/docs/authors.html index 3cc81ee..757ad07 100644 --- a/docs/authors.html +++ b/docs/authors.html @@ -53,6 +53,9 @@
  • Articles +
    • +
  • + News
    • diff --git a/docs/index.html b/docs/index.html index eb30be3..cff4284 100644 --- a/docs/index.html +++ b/docs/index.html @@ -20,6 +20,9 @@
  • Articles +
    • +
  • + News
    • diff --git a/docs/reference/fuseInfo.html b/docs/reference/fuseInfo.html index f4c9a5f..4ad820d 100644 --- a/docs/reference/fuseInfo.html +++ b/docs/reference/fuseInfo.html @@ -53,6 +53,9 @@
  • Articles +
    • +
  • + News
    • diff --git a/docs/reference/fuse_hydrological_timeseries.html b/docs/reference/fuse_hydrological_timeseries.html index 924d035..369a03e 100644 --- a/docs/reference/fuse_hydrological_timeseries.html +++ b/docs/reference/fuse_hydrological_timeseries.html @@ -53,6 +53,9 @@
  • Articles +
    • +
  • + News
    • diff --git a/docs/reference/fuse_parameters.html b/docs/reference/fuse_parameters.html index 73f4a1a..fa5919b 100644 --- a/docs/reference/fuse_parameters.html +++ b/docs/reference/fuse_parameters.html @@ -53,6 +53,9 @@
  • Articles +
    • +
  • + News
    • diff --git a/docs/reference/fuserouting.ranges.html b/docs/reference/fuserouting.ranges.html index 5f6dbb1..04e75ef 100644 --- a/docs/reference/fuserouting.ranges.html +++ b/docs/reference/fuserouting.ranges.html @@ -53,6 +53,9 @@
  • Articles +
    • +
  • + News
    • diff --git a/docs/reference/fuserouting.sim.html b/docs/reference/fuserouting.sim.html index d95f7b6..796d559 100644 --- a/docs/reference/fuserouting.sim.html +++ b/docs/reference/fuserouting.sim.html @@ -53,6 +53,9 @@
  • Articles +
    • +
  • + News
    • @@ -90,7 +93,7 @@

    Ar
    mid
    This is the model identification number (see first column of modlist). This is a numeric value in the range [1, 1248].
    deltim
    -
    This is the input time step. By default deltim = 1 for daily steps. All options: deltim = 1 (daily time step), 1/24 (hourly time step), 1/24/4 (15 min time step). This is a numeric value.
    +
    This is the input time step (numeric value). By default deltim = 1 for daily steps. All options: deltim = 1 (daily time step), 1/24 (hourly time step), 1/24/4 (15 min time step).
    timedelay
    Time delay in runoff (days). This is a numeric value.
    diff --git a/docs/reference/fusesma.ranges.html b/docs/reference/fusesma.ranges.html index 6c1d780..f19bb39 100644 --- a/docs/reference/fusesma.ranges.html +++ b/docs/reference/fusesma.ranges.html @@ -53,6 +53,9 @@
  • Articles +
    • +
  • + News
    • diff --git a/docs/reference/fusesma.sim.html b/docs/reference/fusesma.sim.html index 95c5b95..9417185 100644 --- a/docs/reference/fusesma.sim.html +++ b/docs/reference/fusesma.sim.html @@ -53,6 +53,9 @@
  • Articles +
    • +
  • + News
    • @@ -92,9 +95,9 @@

    Ar
    DATA
    This is a data.frame containing the observed time series (zoo objects). It is structured into three columns containing: precipitation (P), potential evapo-transpiration (E) and streamflow discharge (Q).
    mid
    -
    This is the model identification number (see first column of modlist).
    +
    This is the model identification number (see first column of modlist). This is a numeric value in the range [1, 1248].
    deltim
    -
    This is the input time step. By default deltim = 1 for daily steps. Other options: 1/24 (hourly time step), 1/(24*4) (15 min time step).
    +
    This is the input time step (numeric value). By default deltim = 1 for daily steps. Other options: 1/24 (hourly time step), 1/(24*4) (15 min time step).
    rferr_add
    Additive rainfall error, default is 0 (mm day-1).
    rferr_mlt
    diff --git a/docs/reference/generateParameters.html b/docs/reference/generateParameters.html index 8caa99a..6cc3cb7 100644 --- a/docs/reference/generateParameters.html +++ b/docs/reference/generateParameters.html @@ -53,6 +53,9 @@
  • Articles +
    • +
  • + News
    • diff --git a/docs/reference/index.html b/docs/reference/index.html index b669dbf..e43a5d3 100644 --- a/docs/reference/index.html +++ b/docs/reference/index.html @@ -53,6 +53,9 @@
  • Articles +
    • +
  • + News
    • diff --git a/docs/reference/model_param_table.html b/docs/reference/model_param_table.html index 2dbb61d..71db082 100644 --- a/docs/reference/model_param_table.html +++ b/docs/reference/model_param_table.html @@ -53,6 +53,9 @@
  • Articles +
    • +
  • + News
    • diff --git a/docs/reference/modlist.html b/docs/reference/modlist.html index a219539..2d01613 100644 --- a/docs/reference/modlist.html +++ b/docs/reference/modlist.html @@ -53,6 +53,9 @@
  • Articles +
    • +
  • + News
    • diff --git a/docs/reference/modliststring.html b/docs/reference/modliststring.html index 59000cc..c760837 100644 --- a/docs/reference/modliststring.html +++ b/docs/reference/modliststring.html @@ -53,6 +53,9 @@
  • Articles +
    • +
  • + News
    • diff --git a/docs/reference/readmd.html b/docs/reference/readmd.html index 74a3e39..abfabaa 100644 --- a/docs/reference/readmd.html +++ b/docs/reference/readmd.html @@ -53,6 +53,9 @@
  • Articles +
    • +
  • + News
    • diff --git a/docs/reference/readmd2var.html b/docs/reference/readmd2var.html index 8caded3..0611d7a 100644 --- a/docs/reference/readmd2var.html +++ b/docs/reference/readmd2var.html @@ -53,6 +53,9 @@
  • Articles +
    • +
  • + News
    • diff --git a/vignettes/fuse_vignette.Rmd b/vignettes/fuse_vignette.Rmd index 739bf05..6799511 100644 --- a/vignettes/fuse_vignette.Rmd +++ b/vignettes/fuse_vignette.Rmd @@ -18,7 +18,9 @@ knitr::opts_chunk$set( ``` # Introduction -Implementation of the framework for hydrological modelling FUSE described in Clark et al. (2008) and based on the Fortran code provided by M. Clark in 2011. The package consists of two modules: Soil Moisture Accounting module (fusesma.sim) and Gamma routing module (fuserouting.sim). It also contains default parameter ranges (fusesma.ranges and fuserouting.ranges) and three data objects: fuse_hydrological_timeseries (sample input dataset), parameters (sample parameters) and modlist (list of FUSE model structures). +The fuse modelling framework takes as input rainfall and potential evapotranspiration time series (areal averages over the river catchment area) and returns a simulated time series of river discharges. It can be used to understand the variability of expected hydrological responses based on model structures. + +This package is an implementation of the framework for hydrological modelling FUSE described in Clark et al. (2008) and based on the Fortran code provided by M. Clark in 2011. The package consists of two modules: Soil Moisture Accounting module (fusesma.sim) and Gamma routing module (fuserouting.sim). It also contains default parameter ranges (fusesma.ranges and fuserouting.ranges) and three data objects: fuse_hydrological_timeseries (sample input dataset), parameters (sample parameters) and modlist (list of FUSE model structures). ### Dependencies The fuse package, as well as the examples in this vignette, depend on a number of CRAN packages. Check for missing dependencies and install them: