From fc83571742f5dc251cd87b4f103824d205ac7476 Mon Sep 17 00:00:00 2001 From: JoostBuitink <44062204+JoostBuitink@users.noreply.github.com> Date: Tue, 28 May 2024 08:52:36 +0200 Subject: [PATCH] Add support for `alpha_h1` parameter and the `GLCNMO` dataset (#272) * add support for glcnmo and alpha_h1 * Update wflow.py * fix wflow_sediment lulc_vars * update (windows) test data * update linux test data * update changelog * fix review comments * fix linux test * improve resampling; update test data * update linux test data * Update hydromt_wflow/workflows/landuse.py Co-authored-by: hboisgon <45457510+hboisgon@users.noreply.github.com> --------- Co-authored-by: hboisgon <45457510+hboisgon@users.noreply.github.com> --- docs/changelog.rst | 2 + .../linux64/wflow_piave_clip/staticmaps.nc | Bin 207815 -> 207495 bytes .../linux64/wflow_piave_clip/wflow_sbm.toml | 1 + .../wflow_piave_subbasin/staticmaps.nc | Bin 1288377 -> 1294909 bytes .../wflow_piave_subbasin/wflow_sbm.toml | 1 + .../staticmaps.nc | Bin 990934 -> 951205 bytes examples/wflow_piave_clip/staticmaps.nc | Bin 207791 -> 207471 bytes examples/wflow_piave_clip/wflow_sbm.toml | 1 + examples/wflow_piave_subbasin/staticmaps.nc | Bin 1288349 -> 1294881 bytes examples/wflow_piave_subbasin/wflow_sbm.toml | 1 + .../staticmaps.nc | Bin 990934 -> 951201 bytes hydromt_wflow/data/lulc/corine_mapping.csv | 92 +++++++++--------- .../data/lulc/esa_worldcover_mapping.csv | 26 ++--- hydromt_wflow/data/lulc/glcnmo_mapping.csv | 22 +++++ hydromt_wflow/data/lulc/globcover_mapping.csv | 48 ++++----- hydromt_wflow/data/lulc/vito_mapping.csv | 48 ++++----- hydromt_wflow/data/parameters_data.yml | 43 ++++++-- hydromt_wflow/wflow.py | 59 ++++++----- hydromt_wflow/wflow_sediment.py | 49 +++++----- hydromt_wflow/workflows/landuse.py | 4 +- tests/conftest.py | 18 +++- tests/test_model_methods.py | 2 +- 22 files changed, 247 insertions(+), 170 deletions(-) create mode 100644 hydromt_wflow/data/lulc/glcnmo_mapping.csv diff --git a/docs/changelog.rst b/docs/changelog.rst index d82ff627..35e0a014 100644 --- a/docs/changelog.rst +++ b/docs/changelog.rst @@ -20,6 +20,8 @@ Added - better support for WflowModel states with new methods: **read_states**, **write_states** and **clip_states**. PR #252 - new function **setup_cold_states** to prepare cold states for WflowModel. PR #252 - new utils method **get_grid_from_config** to get the right wflow staticmaps variable based on the TOML configuration (e.g. detects name in netcdf, value, scale and offset). Only applied now to prepare cold states (e.g. not yet in read_grid). PR #252 +- Added support for the "GLCNMO" land-use dataset, with a default parameter mapping table (similar to the existing tables). PR #272 +- Added the `alpha_h1` parameter (based on land use maps). This parameter represents whether root water uptake reduction at soil water pressure head h1 occurs or not. By default, it is set to 0.0 for all "non-natural" vegetation (crops) and to 1.0 for all "natural vegetation" PR #272 Changed ------- diff --git a/examples/linux64/wflow_piave_clip/staticmaps.nc b/examples/linux64/wflow_piave_clip/staticmaps.nc index c04951a64125345445a399b3a9dfb7edab3f1530..79e524eeb72eb1c755f8d18888455dbd40085538 100644 GIT binary patch delta 1992 zcma)74NO&496w*cQ(g~w@E#BOi1GCRMTE_a8q6L~;qjbqNXrbx7ho_c;fY(#fS49- z)e25=hpkUyE2*iwfKllIt-QWNB zKj+^686GpXpE14>pW2><0Us;%I{*)@aqik!;w(eP&r!1O!|skXeI&W=oU3E10lL-g zum=+IcMhUqcr$d0WBw@t2L`{L^4rmS{9&Rvzn3kY1_zNH)9#E!rt-^-8q3~ZQ3nRb^DRT` zHnA$faV9p0dNOTf$I@$OX>;bnf)AYOzAXi()1GA;q-gI0Lmvsz&DmeAfP3lc3Og*J zZO~3%sz?goAMt}%Bs@lsROH0Fl7-qrzg%f;FPDX=zKX-Iw8r{keSj@L1bnve#Yz4T zl4ImGA0uBPhH3Q!gJB|-xX(}Kr_os!k|!AgmNUot^!jD;KydC%#-*Jbl}|D%ET^9< zp5x|X*i);%{8G7~63T+0{%w-!qz%RikBm$E-?&y84p>s3N_bdY z?mu6w1z;>D$mYXHKS#~a6hMg?0BX-!6H*RG!zWeijnj{ZSbWc19cMq-095*1@MzL- zg!qA=1%Ti6DI7weAcXJO&Zs)4)>aB)xO5%2^sY>HT zXbbp4iTlt*npLe*O2CHSi&#O^p~{*7h)8EQYJphqDrRqc?P z#x#^oCU4EGrX41iay=IswaaN}^TRP(!mMW0o-Q8Rla&TDwvwx#v{Y$z9+VfRj6~o@ z(QXaA-9|F;)fuitfe*CljM4a$T3#L6x}H3{+o95_73&uj@k-e!(k8w@c7i(18a#&T zvC=$%%;fci!I67m0}?ifp2k|s*XX%bmN>mP>qNz@q(m^e|L8JSO~0t-^bGndN_oB) zmZ1TpG@5aOlD`C7HWc>fpM*r72Bf&Myxd?Csv&XdbO0DU^$oS2%G%scvkdWrfRFIq z4W6btPkp7gW=BCZZ`2<^eo#FD-=iG|mv<_m)G;0t9qTigIUK`JC%6URHdC<=HPK7s z`6$*$b3cB;#?J@-ZFSo?004e2WW58j6qa@i(h{^AW^BS_<)$m&DmO!6s`gzBo(So( z;&jYDA|=y`Ye?1fE{8UJKZZo1R-BvsHUu?nIl$$ysXgmToDOhzLkaVaOlSS5T}Z)G UIP$Kaxxs-W6X&!Gc07sy1uv{JvH$=8 delta 1077 zcmaJ=drVVT7(d@Rr#z$8w!(3iVECV@53o-kGNXrchSOPe^2XJ1Gr=_908j8 z%2vw4Ng249mz<(Edo!ANUM$I#3uq*f{L?g212P0@zaRaqLtYVnE5haoc6Phn%kn)% z1t<(*grX6kn)EIQsHYf(P&wNqom?Mt5yKXiX{t5&qcvuNoRj0bS%=9D*Vt2&gO(W` zMrJo}Pu@b2Qz4O_jh;hNfNg;^<6~OG$`lui>YqjoodEh_`$o7PB zeWZp)$SaJ4co|`b(i0W!@6v1I+Eyt>l_JRZN!z!JN)_*DND-uwEzCOFJ`%d@BK zm$Wh1gKk6bgH8AK^2_p`J@Ow{P7rxWss!F+hH6=QeEoBsdDB<(HY+Oqz`LAqI=QJ+ z8#lxu^iU;)OBm(p0;I$39oIyQ3&y|Tm!pXr_hH@`OQwp^PN@Ert|4eY-+LP-(Z2td z9^;L-A>n;5?6-~IY_u%qEp50)81KNne4Vwp#OG6#e+6+P6)y#t#b;Qsl`ps8Sv1aa z`5Dr^hq+d84qs;*H)V2z9l zz9!wDq~YJM-^z)rj|*EJu2@1#&R74%#BPGC3T7c{0_+Jy}tMCS%7ekD>yPnLaK#PKE z)KTg2L>M%KKm)59sHk-fRqN}j>MPgQt`yOO_yPoer{#R1nRqag#33^=5IV{CcRe$h z5~I2&XxEf4$nQ=iwe-%>Uruso3bE-=19C+urzvZ_bJ>9OKYvXiT>KnQMsyQ-{FE}8lS$R=bOU0(1n&2mnO66 zB1N3iLx41R03@JFk-_=Ln*W?m$#C}No0C{j1>=44vRMvYvat<=k< z|NN0+d4e#Qt`D_UcNClbO(*GKmvKy?JXh`GdR{1BVGm#}7fI*HaL{d=Ay_C~ zd^ErOTgL#F?coab&avu^a%(&y@1Y+?xXQeTu+O=E_YQMqKI|_zdj(!~Mf;QlR%n>f z*oNQLU-mG5lC@sa7sF(Ie4_TlVXZfaD}L3EW?j*KMJ7e`r_Z7c5_K+H9_`=Hz?m}3 z7q`~3{NGQQ%U1?Hq{ZfUw$Jfhm?b+2o8P-KbwhkIo&cs)vfsQ|OHBNW6bzK7r$;^M zdg|fE`;@lb>U8(1t!sVAKiR5eDgs&UF2YMg5$euEhd?WM+73)Ohf1`126 z=ACj)3=a}6ettVvaZGrC$1!A{7zP{z7#sr_90M2}0~j0w7#sr_9HU^5V?Y2L0~j0w z7#sr_90M2}0~j0w7#u?|s9h-`qwL_O9 zSwG)}7Lsi^rXb{|A6f1Rs985Qtobb^6ZBO&0A6{XLr`06I|QgHlbNjY=>s1OpMlF6q%gjX4=@^5I6Kvprb>1qW~$g6Y6P z0i9z_W#W}cM_w4(7Gcq{k5X#hYFNzzKQ7a)CJ9{*;QWO17i*J_wAZHLun{gQ?J~2gVq0Hz%8C)3}+8)c= z=e~axIdUC4(dG4qw^7DN#;u;(sx^jbVrob3sgC>_x;dWAKvXf3TN4oX$v81Sp*OWz zz?t>MW$ovof>E4rb6EeZpgaN@-9YD%aZltF8Y@uQD6VxY$N;S1t3;ep=J(E-T$HzH zA38gVYr}Pzhjgi2Fn6UKg{E>XxYvI`{ZqM?JYMZ+Vk#HF<7JMXkvO+zqSsToAnxF1 zbS#x?%`HBIE~RpTF|dwVp_O9BIz2bPV7fCub3#sjLDp1v!W9A)0eyeW>Q$JNiI+Gk zy`Bj{DQWtVXkQF#;Cl?0U4tX*b+oM=C!iQ7 z_itf^_}By~zQ^!OWxmG|JDL>uCiIauM0UGmddwD*%IBd;HQIk6%^0q8t7reIHOd=~ z0A(}7Ze*B|@-u1a(lggr^)|ic@tV4}M2yenR#ot4NViD8TLtreLh8)ec=7QId>@PW zbU5GQAwFyb^gI{sJ8J2MX4VMv(Wzlrc5+^2mUwv>KfW_Rhe#1cFD=B(;!mu*MHXr8 zHu3N(epG9|iWbSGQQJl(Y6ly zG~RaCUA^`fkS}&(=2el0;U^`c9b9OP>lqOJ90%SCwFN z@FHn~z3!R2HNtys+I=i$6!NRRcckzqjN%+l!+G)gRH_Gud^LW;8Cy)DM}3u9Oz5Iv zz4(!|;F7)?zd=#2I(w?u8{SBK>lM|86cJd;=4WpGzJwtx zXj=O+>b2AMkte!ZnCp096kkdl04b^xU8}){p7?Y5VR6cFOMhNW{Yl$0lt1Sddwo@3 zT?wl2)%gii+Mljeq=Eu;Sn82=X z*Tlh3>T>=0&j}AFf1HA^KZo_l9ZfT$Ucr2sSnkq&7{*t()1x^DbaVLcw`kC^Bu!w* z;?Hh$jvbmfkkQmM?>U>H4PrdMo6=t_q*;-zzb740<_T-=uj&0qor$Q@Glq zyqu!UA-_#FpsI9z3x1-v8C^-&2NrkfW0pl&p}Qi=W(}SJv>pmg*KcZOQA$qX`0nKP zg%zR{jz!@(J9p}wHF5bw;78Zxca`|>tOb4y;SXi3pdeh9I9&x2JrxaZfe6{eRalQg zuet({$^vWmqwk&iK(68d(q`xbdA|)-YK@Mg9da@gyDaT1r$3UjkG3~tpQ~AVA7R6f z{{%>kTo%ZVTulBs`o&W8_;&3^@kFkE{&KFzOqBA6KE|u$v$+;}UAJ6>F8qn<-h0X~ z({$teGTnQ_=wxr2?)4id->iA9vH$Ban$a`}&AG0R@g6ks7EjZIhVba(pZaL8F3vTo zj8a=7H}LA1+Rn=}Bj0y7yZ_=h2LmtqP+czWHRD(Z-=AHuk)d;-DE9V--f?U(F8CjPfFdQM)NprE2DD z3S-2A1+UO8dVQfVW@KGiNHg?aI%0g+AFG3D#=ba(;d}D5iDoRArZ8eP=TycAY$Ob? zw?-_Cq+2kj(Yb$c)Jz}bxJAWALs3`HGp6_VVR(UaW%`8H>yTR7grWabb~j*TRS-*d!%ykQI{RHpto z1HG*=wB-VF(P@n#t$V}-BFZo)bIxSf-BnD;BV5m2bI8r=F64q`wrhW=H=0`e&#XVXBzFY{2F`H1Z*-#jsItFYL>5H?T3Y7ONhS6}m z#*C=F1VrujmNz6#f2g|AM}BU=Fy+>l(gu8p=EiG6#b4}(N*vC(q*t=nk*4h1c^O_@9?}Jx7eb z!gwr<-&`Pb@cwecaifoLcF8XhG`ngYu#d7PAK~Qi2&0Jag7Q1_f7Afy*9enO)4qlN z7iDKW4+ow=+?kbA96Mc?I`ucCC_PI-Q@WTo{{Kbz|LsLMzVDUpTV;p3ei=$Gv?c*h3>xpZ&*lm@V_o*FM%qsGa{)VS$CYJ4utArGh+dh~iGApkCLsqv$S)cDpPYJ9LgdA|qk7mrfo zq;xelWXQOOdxS?a)dJ&#$wLMdhNr0U<>B7wZLPysoO9gLmP6VeCNIJ9UG18>Nv;Q# zA$4lpD+z>wr5EH5%h zav%Un4lpD+z>wqsLy`jwNe=cPBXzfFuVPk{n=2a)2Sp0fr<87?K=dNOF*5-AN7vAjtuSBnKFh9AHRtfFa2N zh9n2ro#cQ3Ne(b1Ilz$Q07H@k3`q_!Bssv4=3`;pjc-wDx^*wLpZ)DIvrj;0Q{@@; z-41@lJ>J1@1Zov&_A3TwvqB%mI5)rv9%qGkAm%#3?X1vOAsPob3kt`L9XBa6A7X$N z5)>vo0OEiZ5`mWMe%=X5GA)c@jO1>8YQeaC2n<%}r?fB==Ge?>V{@kBI6*a$=GNUH zKv-dbB4_G{FM|^B%Y(wqEC?4?7^rYHak&#SXW+0wQQgc!eB)}+ou4(IUIrm z2Tkr~Zo%V#L8Vd!M+(~3&D@m7A%Z^cW^U#WvB3(%l*Zi9p*TR$;504^6$@q)$Jr6P z$3}ZX+XFi+Woz5rsAUiH82->1t8@qHi@xn>%`uRda}JA6zK#^4f1G9!J4Bg3V11j( zKgyBK?8P^WEtyM!Q3LV&>^7)FjM;27_9UAJeP>iqqg{Yj#$YRZ6Gt9(6>FQI;Dff{ zrm-v|wPgA%yqR#@^ZDa5zA7SR&>Pd=(^JHC#;a#!Bm+4=9LPFK{sO4Z7@29h)73tPln$0<5P+rZqB*FAH}l23N1Zdi0Kr{p4PE@N-stkQ|!@Z07;Y?%MV8FQh+Z z2eFyJgVqCKe5ZedOxtT#A70~}jxwL))`@4XNgA$w0Nh{N+I)mn&p-br(E#_(SVP5M zAoI=$Ec3|W0=X0G#EBzqABFKH9q_NOKZqjR`NJs=&kLYg9^&&;2km+8z4$C+6sv4Zq_Iq74XIytjwTxPzr zFemqgk>1j8&mLyUrN`PfHWNQ&rI)yWuMK;)bNi$v=*n_yHb3eh%)op5@ho-eFg%4Y z##{wP@d3=p%X?-yLI-U4Jr=!U-Rs*N%raj~>wEQAo5B&h+qMSbm(mvW?ID{U{aI^U zA}&8<8?=vqxUIK%ZoVD&;mK+Gm*k*9l?&~9v}?0{i8yhg9Yu@ALR}APEq8@%Q!7MapF$d z*|2(vN=uQk(G;k=C%fT^1ohGcL6iDM2}f5?Z;Z1x;1gRc7jiE6C!D2KXaNA(()AKjpodRB0@Bn^%z= zrP8Bs$fP^C$MC&RCibSVX-d9~ZgA?WTjd~~z8+cH|9=!oKRqLN-1l*L= zOP~81PlM8Gy(;9%e^My@=cdS9O0H3vskiS&@<#cN@-SDmEG~)jU zHQI0&z^pQA&nltEfA1!VBq${%7>3frtxo|SzJvQujMf{CXIzceD4U-3PZ8^^R_JZr zk{%abSIYsXVf~GaU#&)sKxzH3A{-CDVY6=PoyV^>*bwpj z%}+yPxf_jWmFBb|o$mbK5EGi-%|v}Zv{Y-QkUGuL?>h0g zMViZc@$WI@sns?4Aed;Y=fK9s{~%`QHZTB)87P4+m378XH-C^it@>W-0S(9h{g}QN zGin)HyF31A;@#_Nh4+H)p!0slU)M}S{`AO5p>_wO)Gp{Q@-J6hrK~aIL6y~SN2$jX1bh_Umh-`j5Agq1 zNn0m2sHnJ~;Jd#Lg4xDfb!?oi=5Cy{$IjvZ;lL~o`^SErb-*hK`KP?Ew`%Kn*Xg*# Sd&p8iz8}~OiB7?Bp8p@Bv)jc0 delta 11042 zcmeHNc~lg~x9;lh*bg~5Ro$KgZw=aeI%P?#HhVW)DE<>J{HLftFZoprO6Wh7*8dwFB-ohjrF}nPug9sq zm)xL(V)(DM)3(+)PpN0v#_ux;l#wzD zSBg4i8hpAwUt`wV_DQ!HD1*d(w>0*gbAiW3d+-sBU_w-eM0HHtGPYk@LVrq4ktg(^ z5LM7_=K%2`X_r^S*wfz&P1+iz7YiWeJQkXS4%1X`!GGG&1$FgHO+LubvAtM2>O)hX zJC~~4u$AYNs_USX%TTn=dEXM&>XwH{ye_p5uF)Y` z^PPRHa{{)is4T+S{=Bu1vhhyQqG$ferkh97ZNs_f6iqgA_OVsCORi_4yS3-7vrv7X z76Wy3CaPZ$s78$Yce?szL9eJullHW8c;UkyLP-o&GoaV(xKz?7R5fW+JlnfiA5YY}4o&gr3Mc8Q}y5}q|l-@Cpk)R(R@c!?tm#@sEN z62I40i6@R!;*2aM_8Fzb=W~>JN{|pX5@_L5MMbA64Ai5ymG}^+#B*;dak4Kflmri8 zO;+NIGn9B;og7zDRM9oLfuaJe;_!n4e~DG%nnWeuy{}e+$jnLF=W_~##wNm zxwIM=N1-fcjPOuXab1O~AI_aevv_pPXsb!S#S)GB(ofqf7en}^FP>Cf=h82t|XS%rS9?B7f;dNSv} zke(jG{PxnVEGEGTz{Pv%RX(2>m#*Pj!wgu8)!fJ(cAc0 z1MSZcF~NTOXdi~C2=1|uFv-4nd&g~cV4 za*FNa9K|L16I^=ElDG*UZ4#lP%N=&2Ab92~(-<7^DXV6B{>68s7e@wNM})T@ARWH5 zYrF>z+^`TTIHi%cd4KZT=2w~n zyfZ~y5N2cl7|rW=S0f$X`a%0TM!Ct_L{aI&%Jx(85I211x}oMDqe+5`xrblS%l9*# z9?9E(Y7S##XFcjVfI|b>SOeE-DdYDFQziO{6Ki{ssRx;*s<-gmUBtP1>bO|0dKWV` zoM}=NmHZ1~*-e>CxE&AGCJ%1XLnb4T34AQC&g#S3*o~dTmhwKH9$ZvAHa(CjZ}y2s z?ENWLom7M>gC#{AZl+0(ubW1SFschJlFx|+ypQ#XyqcnI#3%Gt`F*LHR!;R`14sU^ z@#LDOv056N`f1cS_g!rk<51igewEzmd)yu&m7Z~%FMYSd=Pyr{ja~I3Eu8mRA(Ysp z$H|>lzw-aITv3F|&6CkQ6#Mj^Hj8%87X>8qpGr! z)UZ;2!O>j9X+8JIsX}I^R^i!elgc*8x54hAoP=a=@?uF@Kg0$~gIO;O=82<|91|VI zummDDNTS^&dQfFWO(I^qN#$Jy>mXv=NKR_rDaQc|B4S^a7oJ*CQXW@1J|EUc#I}`a z5B@$#2iOTpYp zGNTmV!y4>St7^eqf26@3YE`Y6EgiHtVlHjLS+%MxdXF79)T%bHi*8Mh!+Bd&p|+h# z;*|C2&|6%-{bM@PJvh(lngfW-+p5yj^tLKox=mHe`rJRhNLbl8va+qtWu8b;z3QPR zQ@}j4_Vwn_H|e=WVQMl+20oJy@-zL^Hg<0szT511T91Eu-+YVGhb%5t5REX?cXr*0A#<{|u=pQWu3w>Y!$+C-Lfsc-}t74NjksT2;6 zI{xs;pEe8HCI0>L$#h!whEBL~f+h#w*{IUdZiur_Wup=P)I(!qH+82Lb6pkGI7UBhAP8zE~j0g`fVV&dFe|iv8k}ivxIL45?}O z=RJSt>?NiohI5bC_=hqPf4gL5Gp?)A%)!N##P2W7>gkf-naRB0?3esjvX_iqUzc1~ zZqZmnnX}CNpRVwbRLTrT8csjMXEK6)PVmnJlRqCXYl7jGD-0GmVWs8y9$X(=v*lrTfBt<8FLgSQ! zTJKi-++f@On#8Ir#yBeMLrtHk@vMVdKjzd{HC}sA8&K7MBe?dITv2}Rd9d}8Y?^@+ z>=nZuml(fFegBvtm)At&F!z%-Aatx^mk|IT8Es zB4PGb|GGm?=>71lbX}Br_G|CD$?a!TM0?XmX+CoMs+pqwtJDYk<#zi{(f;7hj|1fP z7B>XDCHl|=Z@K-7mr8tSbm!SAIZ;CdN61>%bX^=FCn|SHMD&d*7v)5+>k?6X{jc}s z1mmR^i%Ki#_?DdbDo!Fc{yN4>PAr`)5xU>cC{m^)cS=OQZpnH19UeC%qH$f=O*wJS zOCwgEKAn9eCsxOcg!##$NUfY0{3c#HiVV*7U3^pIc8gOgX_5cINI7wlh_ujEMJEr| z%ZYVQBx1%N!*4BeLUgP&^MYcJ^+#=@8 zH8Hm`&Dj`g!G19sEe@YSav4#DI#a8E6UCVF67x#puTl>mE?VKbODz z$D*uKB!7)NXInCj;}3}26hy6T`+Zk@`iilo=q@Opu>yqcv0h`fJsUnKSvV~IE}gx_ z#5U*IhV{rB(3RhP4p_e4qqE<`WwDSWD7EqBsOKh4e@`(C$+ zi~__i#8L5h@ML`jlW|eUeK}cw3$epPj?ck+rs!u_Va1aApPH3$={a9aHTT<8{f$Em zUhT$9{7%F!3K}s`WRgk^C&HNhrv!$G|FRN;mGvr~pCPAzM9@=v*An{2JKboJE?dXn zt?m!>4`&V5mMX7j?S;-+3CEED(xUZOoTbj&osx9Jq~;C7_HPWE{*@a4|4ogP3^1R8 z6=DLpcye4+U~+0yU~j8*`V*`auC67Sma zN!Yc+ldx+)Ct=t2O?-Zex@}iJk-k8Q$@WdS6+$$yyAtO}32@;3dnx&#CMCX!_=6nT zIg8IvQFUD;gJ7npNe(6Mgq8TGV@mwd6D7_M5#L!tgF1{-;=7AQTm=q$8kG)J_mw#A zRdGQBf5QkR{${=sf7qbJ83&cv=Q9C^z&kSPe52A~%3H8_h+PuB+vE&1TX1P9n#%=D z*X%_2;YEX)ed=w@eOcHp$%IYn(0{#Ja=r6|9)Y|X;Q!NCOLAd@3P9KZgRlVxVFL`p z1{j15FbEr95H=nwgbf@(*Z_mD0R~|M48jH&gbgME2pjMqYy>NX4IDt&0E4hWfgo(a zgRlVxVFL`p1{j15>URknIDoJL24MpX!Uhn@6Y^X% zE~!LoxaBttWUm;b)=TLC-x{NpE$IC1*W#A%Y>km@`P};IIs5m>%SG46F zWjco}U-%qpb>pir84GQaU*YB8%8KFc#)Oco7>f8Qqd|O3I5y7mc?UL#z>- zAdx5;;*Houuob%2GIoGq3!^AEF{?eHq^KAM0AdG9cj!fOw7qh)V*(iu`0ROOU=J7w zh)tIK-1#I;0!iDK+w);qAT~u}RWXHc*e8?WfkQ7Cy<5S!K8oxhIHNvlZ$J6*DB;&TKY0ntWm zhQ#Q#@j0o9g+zOZfaBAoMS(=Q1ro9RujTlU`-D<^MUlNYr`%CU?J@DcNb8UAdmAxKZbc0&|}+hGygV`cQti((i7HY=QfD;@-)d_l27b0 z?&Chk^YYCk^U_tyyx5+kY~+8^*i#>cs75~da^#%q->yklO- zRXNOq4lzCEnz++#ETs1Lx_0*E`@0%{8f?+x+9br~3^C4<2D&hU*Tn?s$5cMax|n zZ^;egdhfBQXirzge4_5;DCz~5+F-fOG~4+GrT$#;DVKWC@*9I6UN=;8=?V0hF#abB b)trZo|B8mVSI=8LK}Y`FxYZNfKa2VgY{`gr diff --git a/examples/linux64/wflow_piave_subbasin/wflow_sbm.toml b/examples/linux64/wflow_piave_subbasin/wflow_sbm.toml index 73aba1d1..74d658f6 100644 --- a/examples/linux64/wflow_piave_subbasin/wflow_sbm.toml +++ b/examples/linux64/wflow_piave_subbasin/wflow_sbm.toml @@ -65,6 +65,7 @@ canopystorage = "canopystorage" glacierstore = "glacierstore" [input.vertical] +alpha_h1 = "alpha_h1" altitude = "wflow_dem" c = "c" cf_soil = "cf_soil" diff --git a/examples/linux64/wflow_sediment_piave_subbasin/staticmaps.nc b/examples/linux64/wflow_sediment_piave_subbasin/staticmaps.nc index f3b9d3b7f18c0c6e0206f66baba714c1f23242a1..3d378d7a117f94b9d955d7235066292fc770ef9f 100644 GIT binary patch delta 9704 zcmeHNX;>7;x~}T(S$bv|mVsdr5fD&80S$;NJBUjJfw(J-G6;%*1A;qxjWH2|kcpOC z6c6IALD94t6PF~SF^QT)5@QgTMB|S7M9uN2w|c3eC+GgX&-ruad8Fs9_p7g#ufL^g zZ24?Z!}0+&-T@6u^??6`3x5GHJZ&gUO-W3rbO$JZp^Z+gjornJ-FaVC`>+@$H)X?m z$gExe;I>sokA=zBZ2Dv7ar?N#sIXBv%4%SqM8;>PBmy;{0C2_{a=O==SwWKy=%i8B zo@qG^L8ubLYf(< zHJXIcGRpO6;TG02?Z9pae_;&@ehpdx@fuENs)b$<#0yUV4!T@>B}gfGDMvvn;j*Z^ zL+WNgf~U2YV}bSk4;NaEPJzswx6a;Zig8*F;ZbX|vsu4}vQV>mvqWS zZ>@F-rLQ}wtjAp3m@SK(4c4DsmNJ7FkYvp=creKq{9~H-v6CR0oT$@Te{;>#@*LI@ z|6cnv%q)S{*A9xNR%2h*P$Q@Z*dXF(Og3wc_BDr~zqLP<@S4B1t7y>liEk36k=j8& z29FTRpm}jRAI;GLV)IwstmmC48vSjn8FQp6Bsl}W2*62PPC}bPSk9fJ-S1l6yUUg{ ztYzmtyXAOU2M3QNOH!2Q)Qd}U@9DxEcKT$a^@3wBREJonrs+L9Q59~TY@JuYx0ZyB$)0)RK-hZu)oVjO~raR?^HA($A4U}7ACiE$*{VH_er zj6*On4#C7Y1QX*BOpHS?F%H4RI3kuAhX@ek5KN3iFfk6n#5e>K;}A@YLohLpfO$J@ zxTv%)f$mM6BZ{->90nlg5saI)@o8?8c2yS~%Ikh>Zd$oA7B&^4;4HMPJ6wczZe#mf zQ=k377;e5?i>?L2x~5alkA=}nMMrP+V<{YnqIdy5Ih7md{*B|XG zQ!U`N6rtGJa3pGof<|U&#Q7xaKDLg_>IXjv?_rcCOoF^EK#@QtJ>$`dI9)JG%!dZ! zaxqz9PkT4;ldr_^Vg#soJ=+<*TFmkmnae91;+h52L7iO4ObY zgPq*)CV7HGtxNzogN6N2wqO&s zhf$ZBkTTrKiJu zsgPQmh|U&4Z|X)Fx?KeQsj9U3qI5QC%Rl8)eOchnNEiSu2o! z3^j!sC+8x@QohJk0@pF2zWaNi*ClWawNnsSz}xyfbIaB31M5n~5EOfxj%53p>ZrU2DCb zh{|4#v?-fLp=;@MZ zu8HR2(mWCzA>$^IuH`7Ak}^;$+R@NTDuB{9p`uEvGj;MeR9#8=&^Q{Rk1Hup8izr2 zmS^1Op*xk7KXqjn(pOP_R6{H3T}AmulQ;(%A4#`foLhluhPj|)y8o?X99_QrCn%{2Yakr5M!qKDCO8-52p43UkG;XrCZqe>uGUBOr>d z!KQuv_3a4Mdz*3!wKPcD*_XAH(P_chYol+KXu?M{dnj zdb}SM*i9Ng@*Ud9Z2Tdq4;L51(^;!{x{T2O;xPT~txN}=c6Q{tHPhq#sAfWcaU0Wj zm*s1v-uIvLy$i=m%2#w^jEwiS6@?OgAeh1QF+z`e6VskQe)OlP9Qj7*a=7#VpzrBy zUlR;yaS!zr%4wXW(xU?RQLO5*OHvtM&8svKx~&n@djkWMJiW%i6`xa{*PzG?%pcs) z^Nfo#wSOEUAzs*>; z-U}M;Y99r3p)!e3Kr$A3wZtRBPH0;WFVsMU?v^Wz;eOvvoMi7xxKQE_?l}q>KjIC- zJ6AR>b$02T7Swx=!o7kdRLFFdgp|BcmYHN$h%O%A2#zCM4~f$z6qx2Fi|iu{>M~d1O&!{T2J_E~7Bpe5!iP>)tI>kFiq2CU z28an2GNDp;X1IAy>Cp_*940lZ(@dp#$rjVp^!ErwPKQ)Rf#;{Bkz`awMtMy8j|{i_aX943xxJV7Sl$xi3`ISGWT*YtV%Ib+Vmfae4;84Ujnk zF+KMYRO3Hl&L;^~=;aO~2+JQ#jUbZ1L=3xtz64Ii&<8~E5SMpAstNk>U-{1_OnZa= zA^_l13_XB!N&?^q3_Za>!rsHM3m8P87F##~3?@*lL81v9j#+OIL*O(O$OT>?mSC|= zi6d|emiDooR>4UOJA(wmwh8RRqFv+|D1^wbx`{ljrhr8-2;579bX0U&?bUa{{@h{0 zftTeEynK0;)4gG?Sdr<0RGER(XAT|%fL9O7eS$yfWGfgLzFCJaq>MVw_HAH%eL1|=~U2?ydOWUH}KuP>I z5l>0N;@|gPbXqRXot>)0k!8R{X$FqvSjcb6`LE?d=*WJx5LAiIN+RB}9gCM4bC02k zTiJB_daMJQwUzD69O_q}hYq)DaZ`%ek37o4?4VGZRQO_cI8*HD<56QyHXwW`4`A9oZgJC-cYK!>>J7xVN|b4(mv<7s*DUb zS|f7(La9ZhTSzkQmOc5yWxVD(QRIHI>SP#oP7oOb8D+zt z+d5c;*maq4m@f77StLEnKR2>SVW=%#yX506R>Vl$Yz|T0H5}7c!?Nu;ur`;i*Rq9==YyS z=Kln===d|$B<7ud1_xa3R?5yJCnxiz(?pd0m&z3fs9diyUtyoPwh;B#VU!_D+Qq;S zwUG%enSG0AH*kEgX~`MDX8At7D)R>p{=tb}bo;;kNJNHnKM#0}kMBCRB- z>k>?bMa6Sd%=y#uO1l;kngn@mn%OedTv*B{*w|v0&{4>ALT#vUGl&xp@TwhzAXz~t7Shv9v!|I$4qhfS$r8@!$qBXdqyh)= zgbNvxFW9MB@uq2oxX|btLpYK=STmkPmU&KbQHq9-&3p|)GKosyipo?v!o(ueNwr7L zDF@O(vWgaKAcHJ(By7B(5CBcOxg^tKDo)DHt%?()8e~Yy(UugOO7l$xSr&6{Mu{x| zT9t}X+<|x5G^nX{#d!sP)ZwUc1smmctGXbNujkc5G@{|0lo@K`#|LxcpM&9WQS`azJcPeEVdjynV?{ae(|&s;=_( z6c0w2kzvBj9vkI^<1mp(zj*cRhQBB0+G*l@+i~ySJX1G8kdocbH_eGtNy;e6&X*8gmkw zJvkyXMm#_OkM1G{qbZzj$5G6a`n5O2UoEx?aRz`GGb&ooPSkguE{THBhwCVGX+7K5 z`|M+JfKHMuc}CRV{EG|OfG35v*0T$PNp6wEfT`Oh7a>=X`3OX>Q;F&P)i^hAJNuY| zHRzqa%0TNUYy)cFtyEIZ!Yx-OO53Sa(@XbjxxG4#D@0p1vx88YMKcm*HY!~i{hI!p z`O0n_^38)RU3p#07b6eQomrY(oDdg=)V_&t7K6~y7>ALl`id@Xxv@=$|2c5RHvRD! z{yJf9Rv-fC+8(OqvhjM}ke-_}Yl{ehQ#qd7`HS8T@67SJ|+Koeja4Cn@7i%Wd zecf~jPg58^2o>#8R-*;GaCN%QvhIc0qwQDrn=pF95FNjqjS!9BffrHP%Kby68AlDK z-=5o{<^cMY7nzj4=(H%(D~=c0vr8`Wau*g^dE{tUQDpQ|4Vjke8Z!)K_l1m>I9pDu zj?*=;{U92l*SIqE5A^BWg?7y=KkjO~W}^<_D?Ve>adtEm%z*w~K_^ECNx|$`tibW! ztCPHfpBZlBElqd`-%YB~ucwhJe4|C<>J_`Od5G{DywOeEOyG;3FVjzDqHipk5iCi} z(a5hvv+ch>V*O`-#NrR~!tD%xiuF1Gw!0YI#(!-0FCqkwEK>a=;Y0k18-GauiWYl} zZ({6t+l35n<;MrMs~ABH0NeeGgl$(Z62@09!Z$LYt@f385}Cv}hJV8#Wv;~I!j|+G z#s{{#RuZt)vJ$q{uoAYFt`fG@suH%9s1ml7r4qK4qU5^7R)I5k_nxK1KwHj}c6Kj9~0zxO_TFbLU*Nj@bzR z6i~+vo~^|>7k?G6b`ZYWxq{AZ*LtDYXUc4ZZ-iWt?_;IY>PN~%{!0|9Yt~NU+Zi5u zp!2^gqfzQ^Z8DY^$))bp`t;PxPI_c&R?F2ONG``U45&qW$;WUCizNS^9W`(Yb z%etohF_gbekh1Ua1if@awyad^%J>ZW@&MmuGg@rWskvXBbXmgPgR*D}4f%SJzry5; zMjShxa1<2-X)I9a`4&2u$5Pl_Olpa3QjI9e;yYxi^3>bz^P;7Sb77lQfV zr^jmx_{uYg=KqsXeVlHBUQdqUp!rMda-@sb_4Pg~or(kb3E~fVu1N_h;wH!I6pXX& g>X2VVewUuTxYA_ZL%O4#@8@~xg@G#_rn^J`FP(N9tN;K2 delta 25402 zcmeIbXH*ki)F?a?0wfT6uL2^VqF|#5XUlGKih`(M57-bD0a3?Zu{RVo z*u~y^!`{2teeB+SCK(=h-uvTQ>sxo-KbN&)vQOP-_q``MuFem*zCWHpw>T305!%Gvi zwn|ln6RS&;q@VbKF7mEz+C_&2Gjc}4Fb1$fl1I~}a$zL}qtNW~WW)AqhG7d>CPxRF zTI_p&zI3K`56kLB>XGbPhKbXngA?`7_pbP+tt?{eB#nR901TAJFAAeulR{5im&V{L_ngOudM>?Hh zizQ5^Ftp71%oE{I+VoKR3=4;v_gNSOpF`9WnUbI-SQw>4M;QsDOvY*py^`8dsx6)5 zOwZ1rv%U&sHegt`G|;H8*3fbtzT~CEqO|;La;c?pBkAE=i#Jtp#`!Eev-E^Xvf=K# z`odCLjCIuMeWd}W-6W3;N@tli*Xn*!R(jCXLi+pe%buk-O(#jyrA%09yqUGMaK*Y} z)v(*zec95R<|FiJJFF4Y>BJ7&n#We$v}`MHI!ki^ndbeS=Dp}Y?Yde_shd->*@FMH zf0!q)*yR)=Q6BhDhq%ze|FpwG!e7;Y+QspQ|I;pxKl-0`ar`lLyB#B*K%Af%&#co@ zvE~L#wMtvLJL=60P&XfQFTH5e$Fk2FwakGuq#ZvN+8tz?Cr1rlC&{hlDkbzmh^-Azq?5rN{#_2JlE^s z_2H*B(z**Wn>y8phmWrh2W+hiE8*l1>ozbXhcORuvAdGZRV|7>%#co9o0%k1dnvR}>Xt zwqhvXBwi~?$_Uo7qWo(XD_z}1a*e;_#xAco{_UVw>HMz~_{toKKOfkbJufYB%tDsfYR|kp-`Jly4Idjas}f$eWlqBYM`r&l8^T!V@XA4K1GbP~Hi&gn zG=Uy#Cc7k|D^t#&ffk8{jDI_b^_JSmzFnr2nW0h*W}mTcE*j*RmY~8mjCJVcR+^T` z_FMVSqCE@%|z*lv?0kK=BdUE zWg}UMmpOkXgSD3!S@JJ4SQm*~BilrlFEVh%o%tPb*XVe{Bg*GTEj)4v^{a_Bj$q zt(A2q3{-(2`>;+7^Ii35(7$MgPygScSrsa=bR6(C@6o>l7yfSsCP`mcoDAS00{A&k zVPHwcTs3%%BvN1UdisB;d-0uQIq#Scr{{y_Yw;gev4>Q6x1^bts?iL|cuz^-jM@on zfgP;0p~|#p?J!R*mGOSb1W!rf3>9!KmTfk@a#nbKnE>Z>rQ6uatW*c{fL{I1sPR<;WI%*al{aVsU(zBQ!`&!~EsqV(Fcr9^|6i?<4 zy#`u>x#S4o#cPho~Jgv{TqpkWK$;|zmYgM5)C3NbP!HxfsOTE8W2pcyKQfI$QXRhsfuZ(wz0F#q63xBo{zBwjOXNvM`+p!qIoY{QiG{I)Xf{*v68ly-#Wy*|PU?-XaQ@>}m~ z=a$~pcHmn+&|b|)jnsXj%J?k#dO&hXSZf?LZT`Id4y}dKQe(bA<#J<}TGEg3F$G`E&oN)KGKx~Ot4)iyCw#kou4oFzxa?l0lc z-TGqTCpr(7Z4f7gWi4iDmfbfNy4Mt`ermd(O{d)h#qLR)p*t;evO>K$e=*ol%D)+A zcvMw9P&!RfAHHwOM}Y|Z$EGzxkV~-78x7b9L zdtZ7}!Kbv)exVxct!-)|i7TRH{+Be9`IGImEu~pSFE&&9l$oi{@6bA`r|OufEisjB zDAF9KXao#Wntgrx@4D3rm1(IqYbyCE4(hQ623c+>QPvw&eUuK0DAugxa8uh-nc6Gy zPJQ_YJNk;7ewvj!i_$H~RCY?_$PuD%Wu>Cs)uEm7{Y57W`lgws?noBBEepOu{imQa zx{t`n3LQCgfaHz|zvHI1lVnwU{`yUAN6Da({MVb>4FX%N1?vTWtkgq2jgbep@m4G1 zw+CWhV@;oa`}D-oq3Yg)rhM0GZF|YoU_OJsch2SutF;}q%-8AiTdK7kl5baQ5+!MK zaq0FM+3LyDKZ$)_nm%c#<6(Rw|~YY?ltb#kLljp|YsN5vieBi31D^WxV$- z?KCYDr8YISxms4UC47V1+D;v8U*6i>E;_smqazZEq0y=h(^|6#qjsj?-Gl~blksY| zZ0vhmoo<=#2s5}t1N!}}PP^!cP)3LEc~^T?S9?Z8nW=;+R<#*4a{)keM#PeN44@s; zQryo1h-txyUXxL~-VS=d)94;Ltm0qOy`-3r`6$Vfrk)wyiK=`@m(j*)Eg3OR`1_x~ zh@MWWajF&%v_~aKYDA%XLv$>q3+J~Ct|yrp2|6H|T!}_9D~7|M>+w_U>JF-}i8qYZ znpiB(8M7MZGR|Xob1-)ww~s{LJ-= zBrWcaoZqUaYr)^oQT*kL%@wowxf>Pcffs|DRD_0w*N;a0hpO4Wd%8MD(jKMe{+vT< zVQYF*5Q5TeUu~?8ef0DEmpQ(gYI)h6t%cSD ztd0f#Ew4Olua(LlsFase^7stNc{`#bSPSgryOblOYOs8n;1zH6ZOM^vY4LN8yQ zqbpN2SgU*3OOh=PalHvcZdbXD`*#QrL+xgt-k^Sn{4%-*UMh>io(I5BA+Qwj^&GLPNLcxDdYE>_mG;R|4-J%bM~ z*{rb|@1v$|VQ7-z$jJKpRWMRrzBpE9&gSr|VZ2ooEau1@n95ya7K3QoGke>l-C&~m zJF&99(n<55KH-P?u_LK<7{}w0Y!}`;PG&B>@jf_PAd~K?sBT!X%bp*Wpp&OR?$dit za^NYqf8jI!L!7Luozs)9WaG$AhHzPa(UPwj-2W zPN_?Ly5VyXtMGx*T`_ZL^wEoHa$GoBPCj&gjsK8zqDmM+T>2Z6MEx$LPkA>owP_^r zKOabR$LNrNw^wn4DHn0bQ3o>ORx=Xj+=i^Ul10KFG49V_~%Y;L~0+aP#^0HcIn{#H zL^_jUK}O`!4+WWj+l<^=(TeQo--(=yn1O?uuE7oaBl4wSB>DLylU$#iO1A3yw;)Cp z1*CjU2vPmj1;O$sQU`D1ZOlan3BGcEy(&!?FBjoDsST>+z4AQ z!Z6416Tb}nbzL7a!8C#B-*6>CO`l_F-!S6yp&!1g(+8VhFv7PFmEnapCvo^8d$P`; z68H9QM&6hh!cjk99M*t@=XNDVW*l*MRN+c42*lfs3@CtO^!6vstywZV*O4@7YeXLQ zZcTDQbCvGW7*(j^lO zJ;?c{tw2bckOzxgiTQuP^Y}=ybAh!2;nm-Mt5~s4*1K_0an2 z@r&9H8|B6ghq!_HO|ZXVF}@!B6n_byCy4)-a0%|RFO0bTH6a_$yu;c48*$d+&3JZh z1TmgH39Ll~*%ZzYo!(Ed^zIPi-NTqvDkI^%gNb5SI*4a`kQQT-GB5{3b{N@7;>p!D zc_5bKi1x1_veP-23~Ha|KqfW2i6_sziyJwGld#0OIPBp`{NY>%X)wh{z|~&)yZ?Gl z6v`uChQo+p5I6~Gpp8h67)NZeB?&KY-I_G-c@KXcWQUhoya9W36ST;cyf0M14o`4Y zKpN5K8A{LJh975 z%>7{KcN%Hh_6x3wj3M6FIkK;9D>9** zKiMP?CgsoCkWLeuhD&LQI%HTD^5lWnM@!Bnq@s5?Hc-*HbJZMrQ zqVVZJT8y+N)1OC^r#GbJ;`CRzdLrzY5P;i8+~-${=nOl_T!Gr-Us7^ zH+JOVqHe@BuYl}gGsxytJ8&mg!2JBgm&4kFPooAk!1d&p`_9BSTb>zw*%>(Go5HnNF`2B0?3>{dZgtXdjae4Lm2|GBb2^T zvq3tw$-TZ;v2WB>?0n=4rmXTc4kuamI^L?sJPqUB`PDyOCwrT7n&n0{dtJKIt9U z(2Katd$_BUVhaJaFy*iyso*iTbjZMAi8nyH5=j3h14wZ4 z2ojT>K9^2}!-8QwM$(EEnlT(&0j?ScuL`x{XD!u9rEjcK&O2d0sy zR)3(MC8f(+2`mnbXh>`}n3KN=J;9d55ZMBT+-dLx+sSR6$Upf`B*w>&Sa)hqj(r*q z{#ynHbQivs+=$4(I zEJ>T$O~}g;0VKBdPy9IUG+tL|LSjFUCUejO%rLo}_q>T**5bQ3RXYj@+lut|wkK=1 zv>;D;bE3Pa1zEDofvk>vg{}1iflekQ9NUwgSL_^!+50}gfqp=emZbCD0J1o?FX>wS z7q4t-CJ-`EN%b@@;|x0p}|S$y~t{v!7$XJ)w(`xkBk zy1_~+ii3`{7x)tMRY3;#e}Lb<*CVA3+JW#kBOc?wVvFV-VOAV5Er}yJqo$K{ioOJg z1;9>^ag(-PaKxW4xOid!8KmVvZZ>=bAqygxPdboY3DzW_nK?N(Q%+Wu_z{nwcEoB^ z1AfK&0729}+sJnyWuXg&!y*;|Pv4cd@| zH+rPoC^<2h>qv~Q_<}14BES2Gl9{V*!GnFr_Nk2o>q4DQ*h5d8oL|Ll(|*Y98Lv-* zpH|_C=aY$Lp&y9~Xe;PJ_t@5iVt#G48<}-9l$d;&Ml3sLlXqL%kiJE)FhB0GH8ziI z0MhvbpMBeb=-gHk*Q)m9akYwK7o+UlHwI64_aOPxqsdY81YpBh5`Kt~@!fJs+qYdo z{i2C+OE-whF>yOM2P=p966(fjWKP!l;hdW;gvF6WIZR6kZyX0ElXuM>iQpy{8WC60 z4J@}0*jC`)snxi}7fUi}$PALODW2q56AN;&V>Y?XK+JP9ix|YTCW>r3fwj>}N`w`q zaxhqDK#Lj>_kR5e%ZwAekJv}P z=}u(TEEzu5gM1H>k&vUENbHUaxc?v-ga^&Y#trG@lfFA>Xgv zVnWtL4kLc=3`mBuIb;f)fL4O?TCxqqnoYp*O(mK4x`1>HCggb|un=t^5PXNtGF?cS zcXRT>Vu0Wf^(xzw+c=htzMcoUL=;eB5l$;Zq{Wvn_;zj+2(l;PEc=hRe-lTd^`|L> zo!La%w*?sxZZG6O8=Jhq{kHfJD^CxSvZj@cto+!H?Cmp?@-I>8zO*IYx3?`sgyF<< zVLPxK!(nDe0Tb#`cDVE=ty?t&%kd3c)MSFynM0b5=HTEl5G!R7OM@@?_oG5kb09{O zXJC%ybJw0u<*E~lxi!t*G1s^yx&KQWA`gH3 z+ocD@M$>?E-|^hO#$Yep!NYn%G!jj!CqBnC;&xx^NFLkhlLMpUG5;$HJDlKgS`vbQ zuLh63>P4Q<4I)!t*MO@(fv3DS#OpmH$o8aL*l%kTE^#(90Cr!()K)sTZ!8Ek2CQt4uI6! z66ARke#$=v0vHL&lgIDP`1tG&V7l}`HTPl@lLO$!Tp);WC%sIpiOVCFXku1=+B2NSK+~Alm=+uNi6d;H?lEQ%TLA9z-hIufj*8KY)(s7Kp5Z)F0d!z_q7LKn?{yi zlhs8cO}Tu;RKYO!8)^+Pk|&A#xCBpn8zGqcg>hEoRMa;-ZM`m`d6eF$K%!&Uh(zBA zB+;>>LGEti38py^o6f<35wYa_;7(9JFaft$j@_5e5au19-3Hu;1(|)w2V&VV#5~fL zIGhL-Q%inyf=-&UeMlaV$cyZNSt8c0aPM`IwtwkOw~29tL+3n4fO zf|$%1d{G8O2M)xQWe8=iT>(eRGDFFf*E{eCTbAs#w{ zapOLX$+%W+NYt3_g0Z2hYuVYH+zXmPvii7!h|i!%q$!tKq+!DmSAd+GF|F_zw~Pcs zc?06yNRpfy2Cm)*N)Cv0Iv@iX-;cHTbtV(L!~$9O;43Kx#NLA?V|Jc}2&FwFJij=a zWBL|%CA0)US>kysM>si^kv67UK)G;YY+3|GlQBel+#e};sL>RYWToSpNHQDu0XrB5 z{!Its{_%K9L>jp+^&}@gC`iDNX%K|$6plex{wKc?>GDe|)B~suj{1=VS;2LVCN5WX zJn=)jEWvH~ZEYoB2D=H;yEwKjrqO&_QVhgh?%;r|NCW~uQKfnWbkn8S^i__%8cK>L z%@!Px$M2!!_evYmt+K1&srwk^;%1}u3C%+&S1FN?RXCG3mpYLQt{tJ3DX5g;vdqRf zLwOKrdQ-4^l=v-oCJ9xKkl>G;r)|1`*~nY>apH`**y+hNoVYQMjA@gvfNZcUl~OZW3c1y z=1{G;hhKk;hGIk>;eT}`Cx@M(v~0wD{>J^>h3-&CXz1W zW|1Q!yAz}C^RUmbJbb?CQ>dFo31*D4k`jz2S(LpL73!?*Uu+YEO0nx3N7cm3*OII- z3nh75Q@~wHNXEe^p~7KmqJ)IZ0P0I-5O<6gOkanI9uTA0+K|C6cA$*6utSGdLYi%s zIfN*;#Rx7UN|`f%yx?W0ueA|)O0!!k*Ho4%KWX+#+uRNw7Nm}%MX5%$lA=s`8#U$? zrgb7pHOc*1q?`2=LQ&m%q6NrIvJla>F?A9xa&4UrvS$lWIwe@(4ZSWyE&g+I3S8f^N$A)f0{hN*(nS|YVV)L`_8!KB}o zu_R;J49M?Xz?D~F$AK+KL25n}pCTaazXBDLhH(6-|AAU=afF4yd0IrG@z*DPOBm=P zP%Szo{Ez`GVGZtM0|Wr0fFj>GN)^hu?)*wY_pTp3joD>-#Bj;(KDk|vCQpx=+VcUd+wnLr>jG;`3 z@BEFUr59>csFCUIeH8yVh6H3s{cI9Nz%=mE|}@DPSgf%23)Y~PLa8dgLc2abaP zUk@_nyV&GXO9&brgk7d3QRP_w1}aT6p|}A?f~JeK{Au=nB-D{xfq7%WPAH%#XDLWq z!L)f$+?q=Yri~|OoB0uQBWpnkr=@iwCtGAf?(GLV-^W3N5xH`rgk=0p5q8Q=8dqTi zCN7DrHZmm7WA0*VYzOAA!|sj+xY2`LT)t*I{ys$x35=ZlnBoh&&4rj_KGaWpkj_gw zLh$QB{I|>!NJ`Jy`{Eq(k9(n@x{-EWp}^iny;H*cMGYaa@`BeVeD169h*O$S13`6t zauO{bQN=hke-spg<49MV1So(&xoDLhgf-o8<-tn)I?EdBZ2iDT1`--IWPfi5dix!Z z`64A2Thd77{i%X03qNEnNV3#nmSDa&rR59fRPRLjvD9ZUdF$N}3Q5nw{ri!ID?>#ODUGRdMY{Ay+Gg-r4)^a5<4lOkT*g|QmM0O7uvLzlGH2XWhlc&$Qf~ zRFDdh@h|R`*%l$VpaBHcHt7-@V*0)if7tX@yg zG%N0&x0M&poS-NqB`aDB0fF_I{-6n?a5~AyrSb9LMy3i+rQlf>=WyaRS1|3CP?-vh zP~z6DdO^AQEk}!QTis3z&y)b|xjjujW|2jejfj2f7JPV~L|`2)k%~rNEEd!E{ZRR& zXi|qlk3=^=zfoWCd(=kJ7j0i#qc*sQ+;3wYLL({aEf?(3gW}0a!L-y}b>V`6l!C{M zDY+@-s9B>gDlK!=PsR7l(L6`}XF5FcoP0|@o^~!*y^52;Q&#=Cn%)EZp!ZMr#iN8? z-7n!6(1q#8`E9zgJpB!4rXYOM8?7}6=7&}5cC{_bEJWuwo_UXRY5O-I|b$D*#)d2pUB0+yxmi&29q8|k9qc1}q>(Od6h z$SG+)a(!?HNe-MsP1@c;I=wUL0Oyn+Wf{nE%~BN7a77U1lKfjXtYk>ZhUTC=baP@e zw5p^SE%Av$5UHYDi5*e1@5z8o9{M?{45_#m=+T%R=vC}DX)QmO&PLEUjRq*06s28z0evI-fSmu z#7MyQiQg=Ar$G`jF&QRsiM>A)4bobI^lnT7EEmBV^H8PZF(Aucx}m@FAKh@kY&EcV z6cG468q=l(*g6$yk6Q&?D?>{y#>3~~sID;6&yoYXQ43igbZt}#@}4&n^)%ZGtY3gu zT$u;LwjIcD9)xTT`Z!|>O0OD-^4G6|G0V`8GA}fAc`Fdc+33TjN^jt?j85+7skXW0d0(_|6FI!M{A(h(RoK{$$w%Stq{OK(&t z?E?#S1126ty)KjOOPi6&g`;DHQB1fVB=hN#8J??E+2 zrf5^aTC^*nDQME!Ao{ZI)DN_-h+^AV%kAKN`oTu|Vh0KGQW+usb`g4&^%!*X4C+?- z2esR@6}cUJft1B7;FSonIR!o4&uqckUel~88mDKWynFv`Rmb; z<*NYgPE>m+h7QUucBz8?Qh}WE`Vy>o7kvrW1qWcjHMOlq<3626*Pb2+aaa!LUn9uZ zQj5*#*;gJ2U5qmOodc~H3?i~dU`C;ICT!FgSe*;P8X1%_b4JiF%Y4Byb*r3!9>tZQ z>O?&dvCo0k>ptS3gr~;n?Au6`J%1q_I~fGG0FYP$9GwC(=Z^xLZUE6(2jVgt;5ZQA z-oVLqD9zOg$!w+z1Xq68gg?D#fG4*ZU=%UEewIp#(?=N z7ig>ErUgCmYYux|4D#O*hJ*BXMS)Ry2u!*Om>Y1j9$Rzg&MpS&Ujmbg(T6!k zX#JZ*0s^-3b?Do(oxlryVEzz#ibhJs-2EVodx1tQw{gQO!2tIzKPHf7-`*LZ4X=Qa zx}1tz2kN&GEYNM(>3IU< zrk>W}iVv@bQ{6(te?CSzIlQn5Wj;JGgA?rp6S_i}``2j|U_Jr7vk&_7!y(A^!Ox(v z=Cgu^EFTWKGy~NbO@vckMkj`CLY*eCT(3(`+^_HHQqU}$7k>A#V&^Z&-cZ`6w z%%%y&u(NVfvlk#{=fHKG6L$N4>@aB8LNHn5(8;J%C^&cr>Js=9eDfw)WDU^qKH^?J zhJ{SwvC#!!_BM3bauXPc7qIgRAVmy#VV9t9&td{QO}G@4oWC3R@)fjw0?geCGzMco zl7%Oki_n8bzUC;mn1u=LKoMGlEa!p<=2FbUmD`@U3$B7{oQkVNx3LNtb$18TX3o7o zcN^%L3z!TCPMd;HY7cTc6jomkG%P|t%ao|mgSP0!j=Z2NolSu@S*ZAM8+0nF05o`= zV8q>*#-McupM$uU-BFL2FVKGsmB#M_X72`i+B$KW{Wha=Hw5w-0}DBW+?JwoHkoM9 zJTEXxv(#9YWg37}8vtif!TFYf5HA5?9|E*10Ik^!_Av#7Jr|7aWT5d9h&C1h#vkD% zYk_9qc>`U8^lmf^4AP4YymxmB$nsdQVAE02<56I&vw-9C;FMQ@#)}17&-*2TFF8E| z7gPlHsmLA-Zv=D@a23OmD;p174g^266wRz86vH58Ow3&1%rNkdrD)rewSsx~pS1{P zZw2%BS#WSvWp7T@7F@?w}NC6fFMj?+4%np4_}&kAOyj;LhiOEJvcx zgMI+5E<>PL4o=ij;Kj^F4Nyq1E1ImB0@Vu0j>5U%fiWU z=B0uuoAsFC>|Y!OdwLZqw_OF(mWv!tl%Re?m0(>q3L19q>^LwMML|~SCSY@aLo9U; z#3~V4Y{?g7c5|vdw>t6%YHIru#%sejW74SI23OTeh{I$yr%{6EOKVUcId@?6Y7p?= z=>9%i!FPLgJS>FBN*aJ)alZ>P*cn7&oDkhr9Elc`ehx|l4s`>MdlMEJ4esP5ycq8I z?E4@>E5KLzfmxGrE4SNm!NF2)bLvrG+G%i}y4=*$hFoaFOXzjhZThPsuH&*1utou3 z`yEVWAn3ClSF+-fIYcAgpbR-cvz8E-u7%JMe)B==n3NM#F#nJ>Aa`Y86)d??Ilf@Q z2Lo45fVIwn7-%5qWHbtGxJ`(REen%C2D2g9of{O}dL?**kr2r(6C%6V*0Y6RjasM< zlBM8n3y@+L22DvvsUeeKpW|rE%0K8~l{vSiLy0l?A8&{u=L*cZ;=Y;Bd#9VIY%pan z_-89f03U)G+78^%O;3pd zZ@Ct*nE|p`4))*((C9tLSUzwrn1*>Mpu&RD9L2?-LS-Fapq5pizz(cL&ev}US=r}7mtk!K87IGd z37to0DK>qT5JUz!u0e3d$3Z*uwgio7Q-$XJa^Yy`p*^lU$Wj(ad>QD?a-h>y2(S)< zM4Lc5ep|3A)Q4DXNfA6#x%L$ZT2G+FrQZZ!J#+SA^y|WP;HDAxCg(TAS9?%T{~e&i zpM*76mAqG=$I55O&TBJjVsS%Fp9KpR1kQDE3+lLRPEd=H;X=@E)IAn#*)ke!0X;BB z{cZ8#oq)AIe0M_1QLc2Yt^&t$5iG}k8cRb=eHIz&?nFLTkI+$K$F2IaAL%BnMBAFa zN9Wc&2W);o(0>X-W<74Kc@0>H%fLyg0v;`%M)`eLAl0TX6r*9vo@V(Fap^%U@(pzD z1>j!_61^Ac?DD>nmhvfp#U!lKFpTLg22HUk5y{lOY5!D_c$)UcN z+Nnp!51>4c<0v}&1x;*FvzaFmjYF$^|DwHfpP|EnRbmpw5Be^1ps9D5a;w{Q&|Z7^ z_Fn)uNXvjrIJiiNoK8cY_J<5ab}z~xtDFZiG8&D7Si@oq^+1)t(d|NP(XVqLoOw5R zu>Ek>G6)-P2wD0e*p;SpG%a&GDC2Yk8gaSB4?$ObV4h(uCI#xE4>3)M#SaMW~C{-37G~~SQD!8Vf&0w=9 zXukU%aIF<+dya%Vx&JFr@R>k0pAYFsHFY^UVs{>G>Cl+lnfC_{`5Xo5ouGm-L}@r; zAypPyDFDIW3I_NvS{z#;yx!HgaYL2H+=E5#U}Rl68oHUj9|B_L!liYz;7%>p=Pv!y z=AIRpa+ku&!IiE?)B3*#S8@b>ynGj^(U2>WJ91D^f$vm*238q?fSGgSKDFZ-OnD7R zz-d6mq7k0&k!A>3}c)2-f=)ctblbZ>|=6<4p(3UyVWDn{(f~ zUqK&i3^}_OQZB)h#IFFw?725xy}5?BO}G)a7*2a!TTZ99j7!x23O3jho(6ex z?`lFht8{<(ZoqsgHH#FJWF&cyJ!#o*xkN zn{Z_9In+7Oj_Xro%k?mHfOx=yo2>YaZslhP(i!&f7KA;mxqu-aAnDrNu%KIrlJj6s zQ*Q5_KQP9EJ7d_2Yt_mQzLsbP7T~oE%tIQw6qZ1->8m{Z)?L7J`;&hJKi#cz3plC< zf~}=7`i>44{J4*tU#INY48bZiDOl7uVEybJW%0?hagA+{mVgn_mYSks3q zI&8qX9V|!L+7IAVHDCdBpiP@|N;2lQKjDC0ehTi$7i(^M$`tg)BSP)BfA_BzM(W+r zhx0cB2Y(GRgVSI~7_M@LEoXM&2ePWx5oow@!Atah?LFk$SzD0NqQi{@Zp{0o$90-e zg1&mZfhd**P1Od>8*{zO8*vO1#<{pO;kwQK4weAwOGoFTHfP>JiJ}H1!rUCbIPnZZ z87FSxf)-r2${OfcqmFQCz38vg5x_=|>mWf~cm51$(-g2NmO%NIoP!x0XWCB`Qv3(9 zEf(nW3ItRisDFiKw-6Qog&gHISh_!8Ov=|_bd`>6!6s!pWiSemw-z;oQVm` zy*+XQst!LO)^yUH^5ck;B1a!#=@)v#X*;idII{O>sEj8eywL!jg zpzS%td!IlS??aeU4)&G-Z>rD5cm9QL_Mx9!(U=LaYR9El`EVzO*+ZJrh%-QQto7w z9O@p;pdfAxR>Br`)#A*zJcrePfz=NM+E~FiQq4h#Jpdbbh!q2&3?Ixz{V*45w~#Dx zwsOp6D6K*2*9lP@)Q0{CW+o;!=CDbaK-`z5KHNGmxKJ|_G$JMTtsKti$WaQv^?r&L zzj_Mg!8Zc7|LDrOfdehM8><+wGcpibf8bU)U=j>+*IroA#Ka%og8{=1T-%5skjnrF zw%&uZAB9Z!5=3!}V6*d-49;51KKoNSN{n|rtHoW_xdk$5!0J zZY&heoWaHWaIITCMqMsAa2KYB3cI?r8pj+GzKi&^HdSRgkCg#*j! z2F<`~7;r8-72K>SNAR>wxV#J>_(hHl*9(0GwQ>W|YR_$7B4|b4fJi$huKz1rj^=t@ zt5-lOwi*3=UyUNK-bb(Qzd#!%KZWuj3-CF+3)@>IEyVE z_>w%?l=!%hb`bsQandC|Tncc1+2NAW){KH}{5xn*gGYkQm(- zl2uAR`trRpnDqx z&QsC4<6BX4*SA6bKYM~pKR`Dt;crcrnYZ7!w-XTBlM8z|MyQu{nve&ntQRmS2vH>p z`MnIaev^eV&P;@)VIu_04N*;-6g1@d?V#(gMnO>KM-`={dxI&R1azObL!5KDAy^#+ zczq{$Jpj>K4-r8ritt^6Bu$q=(S0gBLMTKF9>|gLvN*70bHFIiL5ED%pv6ZM096Ob z;Pz8A=f)nN-IRMhpW)ix&j-PO1nHVL5Kk^Rr2jsjxIAYbV3-eMFfgjWP~kgPodOQ+ zJS3;bA^JK2FCMhdl4tW}RqmbFUxXvbz~V23K=-Qf1g6q)bv^rcXZsof6B@gF8{Nk4 z!a1pA!;3C6SN$ZfkNB)^N)0@#yWFY$aQdun1>bKL{JwF|3|W}>DLaZf1S!o#h(#tchAYlNSbq4Y!`pstP|5cA|$^3 zd(AB*th9}K5k?-HI3hWV{&;%p7T5y*Sg9`q*T2v|tq2lM=Bs~H%1Zszjf|pw;?RuD z(P7DHDXAlxE*HD~HQoAA$yo!EhmN2>4S%*m?A2Vg#E8A!)6z$LUN7$0Rx?+wnX8P1 zpZcmVQ(>innojMItc>BE(vjF6q-j^5t<&60v4fU2=$}G;e?v;*$du&an>Dwpuu_Pk zgg=wI-SH_|NdqoxZh&Gl;HL@dKly<8c8MuNGt!f;9TkZRzbF&_3ZEwXz^=hbEUa8Hai|N4z2PW?;ZKN?wR zOfTf3KSFI~rJh_Qy!rNDv$5_kppQS&<4sHqd-Gcibj_vB)^}J-e>mKSS8@h2t@Eb( zD)<>RX89pT(5pQ1>TI3f>+QtebQ42!C#Nq5PlwW5(rPb3tG)W{-&@k8m(Ok(!h1i^ zUaK0lUzX}cg@3t+Wdh+nLR&ivuqsi{aNBc#ovRURVg#! z@AI-8zRPS`EMJ@qm(d8Ha{br{t&JNt@a1!4gZLsZU2FaxFPj4k$GQh znmMxG68bTm_C(NeIecL$B-LAH0eE1(tUc}ZlhKlVmi>t!Ah-=bYlAF93N3;3x3kgk z>2aXYPIrf0?4@)h?HfW@4X(Glum{CO%Aa2%dqO{`;6sJ1cELm*t&*kEGukQODOeD# z%ZZ!yrF}_s!D;+~b+VN(-9~1A;zJO9{>xXbk&UK(DT-kVJlK}C-Dz>S)49_6isUPCelF(wQ~gi zB>jH`&FS3K+QF_A>{%a-t_1~A0m-4F5D9VA)US)W!}O}L^~$;(t~UHeS2<^<)99|O zZo7^<+}GmFbJwx*faWE{du+n2}&!v!!lF)nQFXw%+BtbvsUd zR4--fzWJ!MZpS4TzQj$=nI2A-EUMc+En3}vy_6TOG7|4dh@7f!uj&75P2KU9mFo7Y zz}OCT+wZ;<+RbtgDdlzB%j~4;BkcY>YDC?R!JRc7pO^mnT(_g8Qqyrf>hAHn9pB3} z9n;1*Kd#$R{z%y6smV3hicbVI@23wl6Z%Q6Ku=jTGVH2a&4%Lo^+Ks{`s$qj8R=qgPUv za95K!<2~H%D2RZ8?(TZBp?ak02*}zvx9;W>%yrc>s*~l3bvr^rH61N)OgUJ$LorIz z(Xs#4Yjr!Wtk-m0URF|XkEyqH`P_aWnA2Sw9jZHH)OYm^(`$Q*QtNhbO?l^dV8w^n z7N6?2>-1H(pJ1nTsoQ>Tmb(2#GsdTG`<#>N_HyfxRdw4Je-qlxT0d~uP`ACzL$0Pq zM}=EX-HuH?)g97ywu5K!u?g~y^3|)@T;n>!7)5^kBI<7>QM1{J}6iB=o9>IEmEnZ1y?6D%;Ref34aJGjKa8KO~KH;6EkF zk)iK8^+pUc?_-w?o=cYd*x%KxVyBrQzFYZMjZ~E+%UP)oy<62_tNQ-1_v>_otHIQ_ znT0AA@Jln~fq+o*V*M-BUB~fGsjOTpYkl2u!{oQ6Gxodx;7iiv6O>}0rrD9XcdAq( z_;T@+jb#SD?m11Ett22kf> zWQ;bS=*L=H3bC@tN{S?N+K{lTW$URp_p0;U%HvZ3<7{O(Z2bUxxP+$eKwaa?v_q}m4sI;J2gj@_N zN$I_2YyaQfp#OL82BrBATt6xBiRuUatGyXZ4KlU2VQJulrs`W_#PL&u>%*Njw?>Nn za4(iHQ-o`;$AVTN|Ej$cOS3`kWmp=x_983|47XnipTq@gvqMdPZC0p(Ym-3@T$=!D z_Gd(k0(L|AtIhQ^aBX&{fotWrLc9iQru9~N)Z;dQiMgV6k$;-MOf5I5f-&lghj2Sa6PRQ z2Z&lJ!lG7+u&9+HENZ0)i&`ndqE?EqsFeT<``2ovDAb}>iW5Yw6k$;-MOf5I5f-&l zghj0sVNolY|7xW;K-5YR7PV4@MXeNJQ7c7Q)JhQst)!9HOf9ve?z~9>jym$BzM;y< z+~5h^Bx-K(`2WYhXSH@I4ANQkGJNiB?JufNW(M9(`Jh!YC~JI{wBnmb8^~1CyBeg% z)AF6Bkg{ZizNV1Us=cAPbp4}A_QDtLd9(vYFWhF3;8kaI(i+WZ>jQ8u$}6#ka}Yi^ zQmyTvt^n4*V@(BjtX*x+=U&(Ct*+iQeV|!JS@oY~4g-8E(kiZgl@L6D?Y$;$tTY(9 z>7PE7r^E`DNN83XqV8n$x+M>f&PvP-88{HAE;O>zP)(=)@XW*!X^H9aS;+&vB=_{` zy@Wi=``*&-FCU~SBzEK98|jXat1ne6gu-WXXC;9n8xICsdTsC_c zQ2ACJF4wTElx`^DIUm*;!}MGWdxkb?88Tyt`r^TT8h-ty{}BWHi3#aiP3j`7KszS# zBv<080=4MHN~c`XuI8OgbXW2h%?)KLPg6sXaj5pKpm##b_!p^$GgZ04hF6;V2p1}f zo6&WcejU|2_X}0Oi+O_hKD?txIh80mOYh?6?`Ipzv>?~uUu7HiHf)R4lUSzj`UiS^ z_Z&kX`;nRgZX;oN&>_uxr<*TFs#fP1YD-OO5unHE7*w!cwRf!HJFS1)>6FutT3lA` zRT{q1;#C#;UHSK+l1SCq4iYdjaG^0AZKuArs&R#We@%AK!h&{fW%ZWlAIXx#OHXwq JOaBn}zW}K;Wk>)3 diff --git a/examples/wflow_piave_clip/staticmaps.nc b/examples/wflow_piave_clip/staticmaps.nc index 46a1e6e45b81ed78d74673e7fc70174a2f28c990..ccb69299a0cc2e7fb961c5cab7e47b94764eea22 100644 GIT binary patch delta 2056 zcma)7YfM~46rRUjmdCQdve4yCp|DV3aS@~sZI$vWG*+M%T9?OGwxM+!V5ORB6pKxN zs4GN^fEs3vpZ((k zi081U6zucY#(!_D@H`a9;+rA$uEhpY6W^rlUhZ}%qn6a0YH2?9P_{oy442gxzovl1geIfm&~ zfx$4vmGWMhOr?=oRtryK2uRMH;Ft5)oCl0cXPRA;ORd^UvnrC)FAbjKmSWhdS3Z6% za)nBXEC{NWNpq*kV4R}S?3(;l)MD)){>>zC)hls|zH}>e_cIHGhvOjGU(!uR;1bd?}b=1HF zp?UJ#I_XFayuX+BcA&%3owKlhgV2J(M|ATon3$gQg1D}-p`ilgAT!gof8%5>2MV8 z&X)UiDDT(&nrP3lrz-mco!&BG+A1X;K+QXZmV3vGiO;qpY#YytjcwBSm&hvMS4^1n ze;6mWP3^{VgQ5pVQKZnc=h3<4_)!(l+Nm4UFV;-Cd@wVhN{$nPVe) zEY}MX%N3ew7dvl?qH(eC{m8yCx+xm3T_5UncXc>hJM1g~c_$S36~RZui(Ss115S6V z$F)DpjM)L7hWti*j_7aKm=pY=CXT{zXqrm^9v~W3h@Ig+#nPwVK^<1jVivv7F};Ga)qJa*yA zhAQdvG>ioLay+(YoP@x|@1Utk-{H#Y@-nbcMS*E&Zw?%hu4G~#j#6^m)jz$^QLm};!!}K+#)i~a+9qa`O8SVhYD`R{4ZVBsq6t5Ge%#5N z`R0D#`OcZ>b&GP#ayaKDrG%hZEYZIRWnNzyUe{3Dge_-qfu1AYp_MGB9Hqr}OQlyO zPibd+N5u+Rl9x*5b!&eZ8;bUVR|FsowD4b#k0p{8@Uj3==*@BW``npMWKm0#NS3nL zAhRNCoQ_?4t)aYM1QTzVJ@xTi*VP5fYML8rQ3ApU*$|-i;&46arx?wlh9i~JTrMCB z#X#z540~8vM$oF*3<83>w;l+w;f$uFpHqzAQa)Y>8O&-grAO@9GOM-^YHE*36X zJE!P{kCx2pig=HUn8b%=XC(NHE%O&aJ=+IGY`;H0`EVFl7!q!G#Xs^j|S z)`~&jLGjdo>SNIqQHYpwYyy5ih4@+E%LHMaxwgqIM5CJk9&Utlq-O-U-wWCF zWKZ*X^<6JK-pt2rq)GjD1a`eb-y-N!UUeH5<)3}LadfrMSBFa61}PlED5Uoko}m(c z-FSDv*UaY)m+QEVd=T-P_jJ z)wW|>DA@5{I??R80I^rLg_qdLdVD8S-D@XSLbF^gejZs4d3++5nDu5&=8_o2 z)A3l&0rHwP`30ds6}#kks%lcACSE2%naDr!JE(FV0m^$_gr+@nmTDIc4u3tX4i=Co Sd7h8iuv|9(p$@vp-{e0AoIRHS diff --git a/examples/wflow_piave_clip/wflow_sbm.toml b/examples/wflow_piave_clip/wflow_sbm.toml index 0c615ac6..8b90f558 100644 --- a/examples/wflow_piave_clip/wflow_sbm.toml +++ b/examples/wflow_piave_clip/wflow_sbm.toml @@ -65,6 +65,7 @@ canopystorage = "canopystorage" glacierstore = "glacierstore" [input.vertical] +alpha_h1 = "alpha_h1" altitude = "wflow_dem" c = "c" cf_soil = "cf_soil" diff --git a/examples/wflow_piave_subbasin/staticmaps.nc b/examples/wflow_piave_subbasin/staticmaps.nc index 9365dc5b269a806e67ac134e7764daf382909c28..4d97e6b6b8a895f2d7e0bd7ef6daccea78e4c12a 100644 GIT binary patch delta 16792 zcmeHOdt6n;_TRJjJ})?j!{H60f}&~i5JSZbQ4mBEG(c2T;3OcX0wO9t5=*akvuoI~ z#=5>A`Mmjn+E)0wn2*w*saF#-(`#j#J#ML4nf+$=>~+rGfcUt-`}zHI4 zF_AqYyM!0z7j}({j*O0p>{3(f({j=Gr6{L553NCBNmQq^4udcZoh7r-Da3VhBh2+#Wj}DV z9Rm61IHsQ!5c&@5z6a;QYI4&C$8P-PWOaeT)95$4&?HMOoJC>#j?*o2EU>P#A8X6n zkGK&EHRH({c?p>*eK|dKVdW5)-ywSw>n-opD)YR!LEBW~_(ub^ffn*y8KPGm5_^AB zy~r~V`CntrV;yQIYgO@FsbillNUU0LczpHewp4_ka71_)cyMmA8PYlT8})X_O7C>^ zisL8mD2KgcfSB7SD@w@)*`bcx`XhPxz^wLPYi6h)#Uv)Fd|bH%4*AcS(^Z<&6{YgR9;*RQOckf&)L0bl_QLJY z1|_Gb4CV}u6&(UlE4(TeJ%UD%Dv#hFVv(P3rG@1jo~tRIxrwd!6x^Z*iVD7{q=J5m zb0|+OeLPqy?18Q#>xi{2rg2}99P7)uZytVK%I%2)c=6Pao~`kWIP_=>!Orq85+PgY z#@nYAl#!$(?Ttfeh@YB-J>yX@zi%pzj7Mz*GG$2Pyz%T`EqXjY9*q<{Jq-?YYX;ty z$0y_B1hg7Q_CbE`fjh@_l<#wv(>OX2jS%b&^md9r)z_ZI#l!fgg&k^ca9{MEIPq}l zT-@5z)K40cgkIPy2sU}xX5bNe{$ue}VZw38qF*{lLk;`|chm*jMDq}AbvNmSOdD;E zFCOX77vkPa-CoC=eE3oLLT@A@y(1d;$Ay8spBVFx0Xn>@yKSX()MzV;!cpDKdelv2 zXI&#Z^ZasumdqYdvt~6 ztZY7j=WFqjY(A9Vw+p|Y&9~rZ@4;uX`F4Dp=kXued{~Y%&x!%coHM5vOiHy^WR;Ac zJfQ$qHxdIC+Cu1zu@g$i7xXDBu)}IcVvxd{2;aA0g1tDOtbF|32)-qMV>(_rf^Wmu zZo&IUkhaV@k1vnlLjo`-VlFVm+mDEE)zVh(ye?UJ5;C`2DM+EPios=Ge=7Z>j}T- zHOe1$@+03Rd|3GGXDM&9O7jm1!$XCj;ZVQrH-vw8?36o7{nTG+{c%qIYOv@hj{Sbp zQwr}BE)u?$@mWaRHE;f0#q_yRA~9{r2GbR@NmYH<%7seR)i2YkIZi%)2jS<>n<%@k zvPmH=#WQC7>Iv~{Y2$Ix%O-5iZZHya@y$!r=ADnb(cxI>{Z5u}$W5x?b-b6rjR42H ze06?eT>b2p)G;$$ZC}DL{Nyy9&sC{IQ?rh9WUknLB2A?o1go zuwY8T_?Oq=BmHzOy24!miER}%Jypwn8g2(j4AT$=y4FG;qX92Z(X|nb2?kt8F8Do@@MT)gTZi5I>spOT{e)$DBt|H<#U~3Y z%5VPw;z&h|((2kMdt)(5kr@Zu%ZrlB?4uV90}38E%;_a1r4s4YG;Tt zwv#}P`di_Lu0WnKOQ(*}aR21g3c5tJu5e4=acPloFGG3~?!F_*LiW~^iJIkI|As~( ze4Mz#hCffzJuSRSwAf3A?Do!;&*hXhMR~qKWXAm>r!-W@b8+%qC_43S6-%o- z@}T!Owvyg|9evYVvToC_LPCizkHsDOUfyk=TDVpH9QkfPMu_u*t}MqF-?c86@^!#T#5LV>p#U&%3s_=hK%WmP2h zRH~rSeJjj3B$ueC^|au}bBX$-2^O49FP5yZ;OV)>R>H317W`Lwam~$uYsm$V_u+H7 z#t^=y4&NcA@aq~6Hs9^eaX4xx!<6ZDdmon9>eWP8tJg*weZ=9P+?>UTA147}DkMbO z)@$EMp@qgyRea2JoOaU~?|s3z>^t=W@q2;Cw(a^}wzf0VLevD07yWD`yX|fI2rv9~ z@N_lj_9O1|Gci^&YID6>?$t(M$2re*|M(>}5tJ^-UcB4%ty-)mF3wSEw2Vq>uO^nf ztq{e3cz0G4xtA4U<3|HtP!sI~MdrMw%J&g9@$CSGSP+;nO-QTC=njLsjhM6 z1BGZkp62M`ytw_f^hX^Y-mTnbDoMsxh>z<` zBV#&^V^uk>X4{ti4e!sRix}%V-5j`cd>OXuP2GaeD$Qx3ScLLU(x*D9T2Gpet?_1e z9Bej~MMGj#bh)bMX89(k^iSZfj8u*|At%JoaXO*tW- z5#fk>lP8Ytz?<=pzI>e2%Vx?K=!5-m9P&54NqT0pJ`v&fHtBuHb|dC%>u2n7{tNmz zKLqFH8;!!KOYWF#SA6lF7UBk6MTQDTZqYmN7q^Vjt*5n2ZlirvBEvUV-dDVj4*Zkc zKD}4Y|E$PH{&`8C6?yNCXAjZ|zNU$}xN|r22>FSg|4ipVHa|Q#+2C?fiGwDa{KQ53 z_vh0}tE8ptO%tN9`D1gITuEBK&a^dJSUTOje()z`FqaN|xr{k}Vl`{a`A>=DMrS(IJb=w)s`V%NmB>5T6EvD=X- zp{fMv7mg8nY`oiKHU0J~p#O&4eZ;K`%QrB(fOjXg&t3xh#AEJ0K}DmM80Db7-a%-VCJ90`?Ni- zsncSzm(v(%eOimp*K6^zPqp}F3A;7G2<%#=V8uT-_oPX&3|ab%7bf14}g zC!9UpV%I99g|V9v81&83;(rYF#@hpUU#ax6zCDkx_H_3WUzpWA$4(2B$Ip8755zH1jBsoxkBnKFh92kTo2RtM>z>wqsLy`jwNe=uwlN=~O zk^>A$4lpD+z>wqsLy`jwNe(b1Ii>$(>!u_*P=F){7?K=dNOFK7$pMBW2N;qZU}us8 z0VFxVkmLYEk^>A$4lpD+z>wqsLy`mROmZNABnKFh9AHRtfJu_0PY_S=xMVU~A-&kw zoopuk-;_@Wb|CRc=?8qKqds<5Ni$=$V|WoW6hxgnF9EmSaXXO5PKx1DfM)Tg%j)* zWDela9%lba2tE=AD+>Cg67qg39cBR%pHc`*|Kh^2W5>fJ zL1Ko&yO9@G0}Jfq>|+XG#vm;9F#EQE34_ETiY}j-$z}~b3(CX2pEkPlUs-WgJL_&$ zl&--|7q@VuzWh7Dyg(<^#yL zh5UOF+2L+->eq_MJ>tz~x9{TF-on@!NBSD(;fi?TED2jq;D6Owdj=&S zj@IIm=HE{F!}U3PQaOq?q{h$3pC`I+mKOFlpL>~qd!oAa$ITA$xaIYBY#N4(mLpO8 z%r$C_E&57aE?9iTE#kTpYPG~u4b^&fC2C)0WZqd*X@MQiF>NElUJH z-OG62V#_MLu*hP>2c}svajP+wa2z?;VwQRpSvp4v+haW?Qehwx4^%`Y&?8%aom6_# zVtkVRl0E7@==s<^C-~YMR@p#zDH^bp7^#ES?&)x$`!+*$>4%jM z;GFnQQmcazwXErD>PqDnqn5n)0}Lfpr)AR)s8p@RdZ(7awR5QnH{~;P`BcRAiS96-uRo7zau233*)QybM*P#aMjchyi!wbvRq zjV3JUlJf-kV@6>*f%HMj`*W$WSm|`H5jQJU_r58Y&dDYBx_P^*{@-1ZI)ri2I%E@y z@w7gSAzaD+!o0_Me2kK9g+FH9B#J?fh?ku+=PFH}8_N&3B^7`QLOEHW%*gkU510HnatD zW({fk?;0X?wfl2PwZOQ#HXCmE)-uh+|6!dr-2`x}^xCyb)cwDImrN6kvJxCa>0%p8 zfFnn3|Bll}y>a!c*(;S#*ZNlp8?07nZNpOCFZH!n516LyH`9Nu87%^%ji(i3Tzu1R zozlDaU+b_ba}wsVxo@?5&~x@51kJh7a!{jTsq{<^ThoKCs=@lNHQY41+nk`r z)GnLTQ#n>fx0(FMnxU@Zx>h_YeRRpySB)v(6MYZ|6Pt;R}Ec03})6E1+cmPKa3l7 z9q0he87M(c%Le^dSAUpmwdwo04s3eb4hJ^Jr5*Z9iip1PtRmU{2)RHYH4YKuxNOSLXt?hrmc5P7FJRDW3P za)a%y`H;MawbQts%QC$J{$}|X6Tx9;k$*&a6@GgaU(_+5d?#101*VmQ3fR<5JiG38jYIh{Sq^p zQa?4gqQtlaV;d1rV~lZ$7(bUplek3Va^sRjqr~W~UJgu~d+)#d+&|wu5A?kCo;qDs zr~0i^)w6E?&^4!1-}CFRW{ZUpKHbIt3{%q_y?EOnD@vHE!Yx!1%-y%(b=m*|k-kcFM@)h;#qAh&K^Bv&Y! zWez!(sypDR7gDNfe>4q2D8&(Ku^ZOwrJit1!@Iwna#UH;(0a#jmTu|O;lRsUd8HLEHj#`x?-df@++kX?+ ztBXL*MYYl1ujtrxRJ%A(jX3VSu69{ado^v>-e9}ujL+RgcN~NF9=@5FPDg>Mq755_ zBQ0+SIA(SE-y1N!Sie@;fI!EFZmH0SjW_y4NR4pJ3N_=mdwFdprbITh_CRPiw$ja0 zno>yW5c7*wi5GNI;^HhNPRmtd(`Y4b%2VRWLE;pm@fAOnOjL%VfVuyc67#$g&%UL^ z3BFkL1Rc;1QQ|YRm3Y-oXI#lJwyVws4CCz~%{p-SV}cT|O;+OVJCu0p6(xSzqE9KOhq=+C~{Au}gIWQ4jS_322rT-u>fFUq|Auxa;Fn}R2fFUq|Auuv-4Gb7S zU;sm407GB^Ltp?yU;sm407GCT>=YOHJFqR6_ws52T#t8^y! z`y#97-6*~2i)AiyndhYEdK~eMLg9W7+^{PVgt=0WJSgxRH`D|H2Y*KEBt}}C>618p2Qe~{McWtBrFC6adfiK zW?4w-Dx;bd#h_en%uclfJ)S|v9AFbk@c^`*gvBAN>pLkg8N`*#?Hq}SM>*WC7wq=p zuHXCHi2q)8EXVG0@-b{6+J~=tJ9_bk1a#^Em%d9Js~IFm%{JlduEX{@mOcsQM`+l~ zuBb_Lu6yP!iNl*Du!5aRy1(qUlGJ#xqv+r-W4kRRa~j#Px;Z|7H@ct5+UTOP|Nh>= zQfpNNZ8U_CIp;NM(x_s+jEqrSn%FjSWuD&QBNV_ecg8%86V^Y8-?&S97+;m+Kwcfh zrjp^FtkqwWRXo4d@Zo=BnC~BV;$e2HP3$9SIaWGrGo7{KQ^vS*LRW%6`S2NRp)0ra zkrsyk26?)WT3iiFJVb2e?0(Yy8~2?r{QF>l{&k*P!p;ieuC~@pBy`XNNd9N+s_wVj z_1Tb^G;jd(W?xyjXvc?qfAAr~Vg^9_=aN=fZl9c9R8nd$&LUy^SZ_8hmn7|D16k8( zGI1a4%QoeaMf=!}?BpP_V;|dLbc;A+Z+U?9Y4#~YY!#U$MN=o)K@t%5k!d$drx#8t zEwaZ=w%b4=5N;#$I?4~UPqG!~Q@J1!b!=NU!IuoHV}n@z5K>l0+j3?$A$4p3NA-a0 zt7HBBTLc;h$X>K5wu#~ms3eefB|1Af;DvQu|u>+~3I{pCE)in|^;dV#+(j~m;ad*iWI(2N=6@&HW~Y2MAP!nr5)Tlaj! z8Qu9Uhq;k~+|icG2|v@yhVwaI{H{l8%9lJ6c%LWWTJGb|Cj! z3x8%XzhUYnFxQh@XuFe*GuXJ*P|qg zj$H6%Y~_OEBZz6aX+7WFp)N7-Y4?Qc+gohznN(JMG1BWR!N+?hzG)`DpeG4ER5$tT z2$k8yJ#T(tFy@oZ@wCp4$8jtjktz$n^Hc7q3m-gI#hSTce?q?}O{GSGx|&(7j=$CO z_Z=!^=GL`X^{iG|dmKrKx~_Z>=<+_`D{ef(f%5dxi_@oke2Tr?Ui4!EnOUv!?FFj` z!tG@u}kg?0BJM7W+yAa&Go{F*o9JQ#<{zVR#ein5fYP#h)~qg2c1P!=wiX@}X$ zCL~U_jl1|5DA<7L(@ILp3X8`Nw3n7mSl^YzZdH})cG@H@oJ$^T(v*?V8dVp(@;_ER zB_nH8q1KudDQ)pF^p{pG{|{0&HVcVky_DGUdD~PPmfc)QCTv%g;+kJ#F`-gFg#=`o z4+s-to*3|sdY(V2(hcB#ZHWtKQd{wjpC=@?M(={2e=A%~V;uh!Vti}>wK9`;dbe89 z7!R?jX^fXqt5>I}=?q=nwZN`0xO5DEXd>@!QjO)VxA@B1OntcqKmAf+dhxK}tM^}U z$95Nlb`{MV?`%bO_AvVi-MMb!7Pgq)&}poe-{v)?2m%3ns`VZ^Qn+1}Lz+r;9%NaT zE<*Qg@jrCFH=EHdQF!?ODAC0uvUIY}gRWNV7o}HAtJSN*@|UYEla}!}g?R8J_A?d} z|JU6kLKx=PsEPsNZb&EWhJ@`?esJis7dKlz_x9fJx)7!pjTgszzq-Tlga69^IG0$)|_qIG6rv($h}%_ zF!^DV#)Z4^SL?%j2G59D>&L%dVNBt8Zji=5laxUSt%=~37%&T*}@$cGFiHSM*5-M;G6d+>mawzG<eWcckk9WTD zyFjla?{zb4d1E(oKg@^aX%}>3-5EKWe=MekS%00|BKg_1g!*B}0NYTnwT{Zs{`G-U z^JrFcKHfhG;Rv~HY9k{>(<$|L1+-3VKmeX<}e!?vqQM1$7i zyMcwAv$t%dE}LL4v!&O%kh})17pMJ5O{O(y+j3uSQ!Wy^CH(zQ|R&z+2LC9eoyj=Jec5&gpcY$Hn=wE-;7C)??j zbHRFYx=~xm?t)>Mb&w9ys`QE8B(`EAIv;kp=0^>9Z4@Ft@KWkv3K$%vGMc zMb5;?_sN4yU5xoz_3S&&dJcc zpCF!xGGk(;%*4}=5zb;bH)g9!(hj?HGRE|i&k~tTR*>dw77zHc^G1ik;_pI_CB=_T zP3GElxjok#rK$`w-hQC%i>qj@R$^NZB}csRuQ*~V(DnA)E$syvU)dvE3qR`cF+Z+= zHk}-pqER7n`OeV|tC8$$Us%6ssb|gbax(5U-E0rgP1z;0?cCQdec+-=<8JqJC3D7^ z+v3c?8+AfQT-UltM#9;qd6Xr`J+)q1@ey+*>7OGQawP`&h~|>W{W=B2?15h6VLdXn zd%3t#`A_LEnvo7dc~b}MK6Y?VAy*sJnq;(HIO!!pmc?ZQMq9l9rUT@CK5nA!Hb;u5 zqALDjt?TP~98G)3L5to>%D1ahiHoP+iZwbkUD&a%8fMj&?8(pX5UOQutHzPMPJnMbrt&W5gz~K`E0U!x_-6?T@7Gd#A_~JRw8(;;=i7u z|EZCivBpJ~@$~*~9yoxE@|NiCrEVv}xvEn@SCmlt$jZ@$&h!@u#%IksN>A*&+%D0> z*9*08#{+%njGITB1@-HAaWGZ}%WVpv*PVA<-Y56ZxTWZZ%ca9f!oS?A|9x)y|9fsq zj|YY0IWZ>kc?UzcTMA5fT?*{(AqgqipG#2Umd%)0i(y)JVKQ#nfl(~npJ~~5$p-7Y z%88VWTM{lAw}?$7*eQ}!6=D=|&f6b{0BD{=fwO3c2j#Ofv`K8b|A z5W^t({>-vSxj-ytm>jzjn~4&CbySJxG%N9-5a~@847p*n5}$on!j)jK`k>Na+;2+U zwWG9xf&GI_CEmM8i7WRhadd+c|J5ks5O{&Y>_4b97(E{r4QWU^?lC$7&GhKE6fNK{ zS7<&#fP&UAz zY=A-80E4mt24xel_+){yfdMERU{E%|plpCa*#LvG!BYUr20SPm#nO`v$_56YY=A-8 zz#%9b@StpfLD>Lg2z@Th$|LF@5S=qn< zlnpQ_8(>g2z@Th^LD>LvlPwrxU(CE9HCqW&_Rl?$We+?iZ%|mK?Z} z964|BsLW3I?Y{Exgsx%Je|nUT@Id)iu0!%lCKrw`EVdQpL7)*%kXd#}9z+}AL?|rC zgODRUNG^0~^#k!oK{!cf(GZvf2q!~X0i6eA#sxzgj!pow;)21Cqmd{41w$Z5LryMU zF!=gHtPxI?>#5=kCfO=z%!%n+gI6U48{uKHfjXgtz6?r(s3Qz7t=y7}3kr+M;%%kE z@iI&Tgwy1Tt$cb(VG&FUgwti-g}xn1vD+rv#@k_r(DyOl8a#twf*=gJG@Fs3m?h$` zT-yInzjC?TM@j~w_qz)_Cr^Cm^eVb3sb_PE@vv6QJu3NIfZw`k@W4La{gX-RR#i4> zN?>&^h8W4sn+(~dGEQ-q&IsGdyT&MckQhE;c(7n zyc@CgG43X#)y5Tg@0w@3NJCF!6`8FyYRMiSV+Jv3j3LCogHg|iX^q{(xq~a*cv>R` zVdH{@`NDzh^FV%1g>j9Z&^MM=A#=z?();7a3C2vW;+SaP7m;DF89i{Vf&3^`x`6rE z663ONT+CB9o<42##gW%)%8wkpXQ6n|$v3Sw9uMc#VFsE0;UuO1k)Kf}9&Zx=Q${Vn z_G=^U#K}Cc+4P*+d~^KiEz#!FQ6@TM%a#IL6k+nfcXHO}NgkSdnzVdK4-@r}^Euc= zj-xiW{OXp8jUY|=CM~(S$h3r?lW*!1&ILUL!{%kw@MzrNZIWT$3hJzGx#{z8?p7yv zp~)HaJScJ2nbVO1-6G2Nu+#O^78tXql6!+_0z{7=Af`!IegClCb11S6B66YZ01kro3Au-Z5M!}haD{* zc=yv0n}va{;y>wC>*CO2)1$1m^&6b5+0xviPWG-FMZ2dn|0*#H8c_3{;?As2?nya)mM`Nw-G7cddc!O4h7%tn`kR zc3jl#fV<#}cNgBLeMtJ3ShRe7v1Lp+AGp_|VofbBg%h&QX-n4fefC-IajoUTt-?GU w8~HvB7J7EW9-_x6?pAg}IJq3Jui~?lgs)c)xKu4%K1E0WT=?ZvTo1(j7nMueKL7v# diff --git a/examples/wflow_piave_subbasin/wflow_sbm.toml b/examples/wflow_piave_subbasin/wflow_sbm.toml index 73aba1d1..74d658f6 100644 --- a/examples/wflow_piave_subbasin/wflow_sbm.toml +++ b/examples/wflow_piave_subbasin/wflow_sbm.toml @@ -65,6 +65,7 @@ canopystorage = "canopystorage" glacierstore = "glacierstore" [input.vertical] +alpha_h1 = "alpha_h1" altitude = "wflow_dem" c = "c" cf_soil = "cf_soil" diff --git a/examples/wflow_sediment_piave_subbasin/staticmaps.nc b/examples/wflow_sediment_piave_subbasin/staticmaps.nc index f3b9d3b7f18c0c6e0206f66baba714c1f23242a1..4a1cfcb2155df2ba4e5c643741de7af4c76cea10 100644 GIT binary patch delta 9702 zcmeHNcUV-%+MhG$Y&pBjQWh3OL_olT0VE(+I*28LNbCx$EJYD;L9nBrF(yI~veB_c z;UZ#*21Rq!7+VrcOrj={L=9p|GXr|v_VTY=st5chRQL;k`AA-KN^(LP)a`|SLK~e>6JzxX1I|1O5AGiv);l*V zJ7RF3;6CBOVV(<%me)KiqLP|3sd|c8yZ*s#tBM{+C0R4+PnpN<;}4*MCgmurfq4=U zm!6ye)BplOS5}W2Gb@as13KXez}nL-rzjAt#PC|A(0yj0<@PyxUT48RNg1yES`GhU+b${((Np*!bV!}aNdl845D zA9KZym;?O_#P4_S;922M~p#HN{~|WQjUUD{AE#hhty5~_)gY09P_Ooe6-MNbP8bZynW_I zbF|ZPiaKOH>1@`oVflptK!Y8t-ac#{;F1Pi^wuhuVEVd~%6i1bjoG}&*d&H#?a^QxbGK>u;{PTAssN;@)r1!ORlM`r1L!+-mI08ma~L09!@; zjLl@NQ9kA%^tbkh5?=GSb`=epKl4e18mS%hWAF&U44N0K^VS?1AU1#9-FnV>lF`q$ znz08fgOWz$w*feb%ZYDO2+O&1sK;Hadk@)ihP7)>7^$&wUhS@q(Q+Hv&yZ7wz_6PQOMY|m? z2ij)-mBjC(xE0{|Z}#|oBV)&>Cfj58dG`3q4tw0tYLD|B$y>&&EeGJuSa|~gcP#dJ z_hx%Mr^OyeUb4qSVtDuLSdiXgk1y2AFyD2<1*7p(VjO~raR?^H5ioD34HuNwC6s$}m++!YI*S3wc?9ETZF-vB ztXJU8~DK-(FO>4rcqI-KDApyp}=~Gn*QT8X_qpGc^2MqIEZ0 z%VqSVJ`C+?lqO7qye&YHKqWon(6LxuFO-l+8H~%tOolz}y@1cY7Q>4XprZ9`7xZd9 zs|xhNN<}w(?I7y!BMn2IUk-7IM_e8i3GXRUdmh!x$qjE(C$PVj@du}|urG?3PK9*~ z!C~18-2T32f3O69?FQ;IZw~}Zu&g)uqN-CY_+EilOsBe1326It%13tt`+)~2t;rk( za!?Cyftyw6`E)9j(LOxlkJr7L<7QAdDQ3Z${7u*%MqNsL@^B|7etMiA`f?^^>GV zbul#-ZWja=@U}k3+;VmM(7IAF1jXE@BiR1yUuOvuLRICu$td*`Hif%bN-bUw!yt++ zq9(u>lJ-SLn2m~-(Q&9cTGIzDs$&Duw+1+we>WI()B_GfPNnz-=^kG*3q7i*YB~L< zy4(-?|hS34C?ZsPo zkU8<({T?rJkLIpSf{S{?>HrQ~AeKU#RzPEU`m?b<_D06(q-s>w+Z5v>8CW$kBQ3Ag zY$3sqVhj@FRF<7zSe9Yo*|NN>l3WrGDaKV2(Fk1=%ta-+Bsfxxn?$;nq3{Z5fGgV3 z&J`dfgfPV1G4y3PyzObPc^{IqD zQ6$b$jJKp)FU~FBG{c-IxY%f(&RE6M zrG)+$hw1Nbr91Gnvm@92Bt702o+R{_w=sQpS)Nwveg8S%yI`E8d_`x*$aq~_Q6SNK zdoh?kLg-O%VY<^#AO9&TM?T@YEbiPt=zBWb*LVY3+*3UjI*pf9dX(cnid8*#NGjv1 zc$FqXw>4pUSAf5gr`H&`qO;0#8WeG!`GXsJj&X5@d&VOY;)T6=2d%!2yZldh(FuNo zY7Ja%I( zW|CPcoC(B_jweQ!6k3vI7fd6GmSO@V+L5QRX5J@x36gp#CQy=a;3d*?E#~6f!u%|p zi2)dgDT?VX3Ft@U70fVa@nR#e(@az0Tq6kdkObI-h)^*a^ok{1Pl?mU=bPpxiR>c_ z>N;291@~`8gZXE93z|4r;Y}y0)o8(7MVF}#1H^<%F~L%IX1IAy$)VAtIYeq!rfMv_q#qss~lv##D#Bi~Efe0u0(h^lK@H*~aAVP?eyw{;4OrTU`2xB%mb5 z0>gz4$bF%rzT6#uhHIiDm{p_e-dCoF$DHG&8N z6EN%w`Vu$|LvIktLtNegsV3;hf8{@$Fzp5UivWPnG4ue^DG7idG3*2e680X3UBMs% zwb;V`U@(DV4H8A*aLjsvXaaLoARBms7=pzzC6>U=SlZimS_Q{3>;mEm+a|CN3wMxX zpb#Rz?k@6F6$BQ+Ak;1zq@lu7YR|p{_GAwe4!kUX;N{D!tR8i9#fnT1q{2BncX>5~0}T-p}(21??$iFk4%7XP8oqEm8l?#whLjw}NvNi%RP z!$Q7O&wV2oLI?M#g`i4wRuJ))tysL&n0*9I+QO#M*JB*etSxL8W`Do>Ty&sSi<^?g zCf~NSF>fXD*AriqlHJ0)R97y7y7_v(rx2% zVGHZZoQoa5PR2fNBJ3Amsa=`h2mjDT#{ToFz@}4L|3jy{2^T*OqY(8zoJd1wz38eF zinXoSA#Y~492VY#ccn^%s)nHGZ|?zJ6ou%eE(D?cY)`!YMcvwa-_nIgL?uB^WLu z-GY*Ex6COYE#o!UiX!)uRL4W$SwUnh#V8y8Z0ukWW`E_X=c;Z7L!T+q;0Ly0E#K@s z-ARb#fB=YEM(vYGB1(J;9B@%?J!fou8?K`ib!+=J>}u!3G@T+(CQld%J)bl+;uE zHA@j^(#rOO0H5}7cnKeS><9=(i*Rq9==YyT;mH$EWb6 z(S z%_gyL_-Qq+>4h?69}#J;!XhtkcKRbhB zNJv+{J$gi9Mtq~Uc&kS-B%0GBVuv=GNGl2Ix_DDTVbR=Vb6!qvNw)$*lOV56HCv{c z3rhF|8&kv*Iue;T-c8-0nntGi_;=F=C>QL0Qf{paN zRh1vW*Yj#28qw6{Bs<=vdR_8tIX`AD+PZ>`c6pT-2;}YkmZQOyxV<{9I9=ZEe?VwA z?C)Q2LEe5BS1^4(V53~3f4YBJPCW5LOFm$u4f~=Vu9LUFNTX3yCGN1tXgWSBU!=SDfxSQ%7<*3FU>_eGuBN`e7sX4 znzG`X7{IJs$HgS=fCC$>LS4YF(hu2^hez3t^_x%@{BV*C0( z1|5~RPyAVEH&}LmAe>6cTkPko5*MBE(LlSk-Hwxpvgpqb$Tc;UNW>4pB}sCkY>Py! z+Uhr5PDK485uNWoNRbl?janGobAfZX-75AE6tffKft;(^NW*Z?lkdno4jV6aG_+>d z9+MN!6%x_buxX;4_~np9oHZ8t%ZW`7BqG-|H$hHh_h5y|F-~juoJNmUv%_FbUgKcI zQq23^)hQ@;4OJbH*2jDa}aj)%;X`n5O2 zUo9GiI0Hbm85OQ)C+WLQmqbDE!*vi{T+jCPI`dc@ptB@Pp5gU3|Kfr+;7Oq^_3Xl4 zB)3Rnz|{Ssi;%0xd<5c+6A5Yj)i^hAEBhEy)#%+_$^h$UYy)clM5%<%!Yx-iO5Lth z(@Xbgxm`MqD}^@Kv4c>mMKcnmHz{2i{hI!Dd}TKt`Q%b8U2$E@7b6eQomrY}oDdfV z)x3po76Z|tXor!g>WVIPxv@=$|2c5RHvQ3P{yMCVQ^<7yh88W^rkv=(p9T18 zl;*MhZ>Xh8I!M%Kby68AlDK-`O6P;%r&1I!@QX_M^}Uy~dTPf1pp}&bMn``Epm= zH5+vZU-226kFsMZ!3^j>6?ArVkQB_0!3rGjy*kb-_=Vv%-qM5y@!g~v{dy{?!Z%tp zuAVU)PYw}YgEza2n+bgI^JUtJbo8x7GlC_FISTm}Ya0LiBi8@#k68RcUbvmXPu>Ig zkL@l-gur(HBH>3Css54he*VObKcs(oi#^6SF?PJ|LI$@A7YVkj7(omG+x?4#ZC5W6 z##b)FH!`5D_LX=NnZ!7Tf5RYUuEgWQmh>0K2e!Ia60p^>61LT_61J7D61LT<61J77 z61J7461J718ExC`) zm-(2~`MhiH{LXe2ie?V+0c)x`}!{IcEd3!LJPI;&gqz4oSz}0BLi{b4_w! pAvYyXr(m3IcZ2*a7CSc0P!vQ3d%%XM2#7lNioIbC zcCq)~u=g%@AG>#-NrnfX_x||S`qo|d&&^tx>{IvIeeX$*%L;?8?F%locfICl#L$DBxupW$@=GdGK_C4 z_IBSDf3#))W9>en>NO1+L(STO*_#`}+GY&%Mi-i98VvL(Z)q57*R^XHbJSfuh@Bx} za#;GA9jC>N=M(Q~+bX{|(m6lbR5Jip`#`5NY^j9F7KTb3O@! z;B%OIB2yN!2n(Zh=qMv$l*wpqp;uBjO0}tjoax^2Q_dHG2?h+ymIoX4)*4)?!?ONiDZDZX`W)bK!<ugYp@sURqsG z$jT3xT1bE0dC|T6hUo-prj&^&Pc*ZZ7B5>md1 zpLTKlq5rgt1aLkA2m*Y$Ccrb*H> zDR1~JLj&_|gMBRiWsWc({t=lY!Xk4-SY(a}i_8&WkvSqPGDn0(=4jx5nIjGmnIpm? zb3|BVjtGm)5n+)zA}lgTghl44;d;yw2Z+oOVUal^EHX!gMdpaG$Q%(CnIpm?a|D?3 zX~1xu@&eXsn?uX&r1E)2W4Gzec4DPRmq_%tm6Tw1+hD#)qE?EO5v=95$}c&rbY&;W zRsNzIyL8*JuLr!#=YAQ-SLI0p`QXOvIcb?=UUGTzHxHFx6Sk#-pX9+>^9N(FfM$9(&Q)7{CN)DNsm7W^ z=>7d32Kc-4B!jhVWU)aldJ^U$KCn$~Ty|(?tA&dKsi!OqODPI6bKC$!r0O$M+n8jb)&?U@(n8V53` z;A2B(MbeAb%qbY)$n2Y8Ll_GkUO9kmz!vjM2Cz{HgwMS~pE0#w+Bu@1lFt7(aD8_~g7AjtM9{NX`tj*X+I>4)%$ z&^FBFbWyVz=C`Cmi|;g;bujYNOq34G7?k>cj%w6kHky@qoAal$SbK?)CI2Feb&2YYYohIj+I%c&YDVku|UL>O*Kz=FAvt2Kd@HOnBQzF(b00UmI#A* z81#oX4VNTCgb}cIN`>10YfE^N%{Jw6kW7cS&yzT6Ew4ABzX}A|k9A_0Z>om_{zWr# z>i-VS>TrprW4|v25C0uF|9>+uMfz&ni69;#fFE-d29{LJRRc#!qV*-Orv8Vz7v4&i z@{WaYdOl>X7XN+)dr*aUNt$V?8cmao^^^optDCSI*uhF0s!Y4r4fE7e8Sj&f^OO`% zQvuiF*=AF#W<>ri6X2XKcNsmAm1^QIMz>R9myF!hobOY3^H&n*2BHaNrREx*=n4Iy zsoBHQN6g@3Ur8ECx|i~!UrAggHC^~+uOtqV(uw@RS3pZJm;8fQ5*sbBll(8*^TdX? ze=Tv5Z0Nw_*AnMOqCsS(4w~hS^7?0HW|*e)B!;{PH{64t$II+AH~(;ku7iS)U|d z_De1bYmK3%Es(e0t~Fm;Zp;^{94~4a8>vnm*6QIbaVgd;m-bZ4Qu=z=yx*#Vr&EIzE3E-oukT`d?-_zM_)_ebdZRcchBmmIdFSektgTuEVl&!iNv;C%J9HZ@;1KBw5jx zzjj00Q8HjS|K)~ugW%??!Fs`qm3pYBG4jw>K59k$dSC2otm)Hlo0&W^T-|%Xla7}0q9pYamu{PttDZddquA%I>63Oy9WgAw zu{w}vrDD>`c52sCY-_F=DvL=TmL8sy+|RIB#{1mV&d@SZYEx62uVpn;!Z*03?bOcp z#m$XvVk0{-IwGMM8m-zeEj5cUYG(@GO=xg78LM{7#{Rd|>6YpC2!q=+px?*pw26%h zXLR`PceGb@wr50?nJSoKRhM5g=K(ZlL@b#{0NODv#QiLQm==uaH5s+*ZKwAugYKcj zD*iRyi%R*p50V^d`stA!sLHo@8EKr+f)V3{KmU0}b$3#YQ8mA>Jt9d`D+=9fqGKtY zKet8bUy_-Xqyv)4m1!h%+fW#EEpd`v{Xu_g;&o%SCYFkGMy-T7UaoP=>(5dBI!qUS zP`2}$4ydEmqYK9MyZPUnWGPh7Tbrm?YcW%UflRJzX*M#85l$}?Jb@P5f1HA!|4_dt z-?p8?+)O?qdo`tv*z3PRuI|kosW6wmJUc*_pT1U+qQ%{o^PBZ_E%X^&8If7U^@ur<9Q2toOlFE&<3KLmLG%N&1A zHa?zr^k)M>Oy#d~ehS&YjAs9|l#=4nr_Fy$6h~HE6$F8Q@LJN3&mJX}OQPp!h>$)= zZYe!yy!y4;?k^SFJ+kHIAT({6rI*Uhq@CUVyGzrB>R904BFodZTIu}$D%lxI>_RhY z44w7<2lXe2)IZ-w*HYT$@$Pm!SoCGQ@}Yr@pP#0iPd_Pm?}r9mRN0}rkGv#T#2Isj z!Hl6gmlbrzQWbchQ{IwB<(h-^UIo41V>(kSxktJ-bM zze9K!a+^oDbd7mOsZ6H2S)+T-OLAPC@k%ZO)89qUNIe6_$y6KlWL4f0>sh>D*&hYT zE=sfhDaZWMF0vSvELdh~Ey)l%M;oa_(>Y(&_U0kKE>>pSH04w97Qs4bzj~%UkQoiE zv1TCD7Fgsw=Zku9!@OJHuhaTL4a=YCCUch9>+-L<$r>1|HQRFevI4cv^G5M9dn;cM z9S5dya>^X092T-=9=TQU%q+f5ysT-%=fFsN1|M3qQDZgU#Z1}6&?Lc;k@XGO#z=Mf z(s-FUo5!z&@m4Xgm?N`)GIy0(2%>4v>}i#8or&dd$IE(4C(M2Fm>&|r4yV>(43CGi zop_%FnYr}(yU<*LOu8pyx?sgldwxigPJ#ZIkMB6i{wLhN`A_-x39`<1PLDg2OH-PW zPcQf5iHdNXw85Ns+V011e9Y;L;YRYWCYb)s`T2z#x5bn;hP8b>!i0?Q_aY0v4kkgx zQ^<(C3%J`=Ez)=05c2bI3>ooHHE#AT0e1;lg!R5Y!KH@`NI_N(P7gbRR~uO1^rMAz zp*X93wh!?a4<|zLT5NHh`*!pa=dyefZ}Ql{zvD*yFsKvRe0&Va9hpwrH|Rq)6!jw` zyJh06KbUMd*_G^II}kUn2*2HAL{6V;LF^(C`PsHPG0WzN;r1A^;$k--n{}o8 zMcPkjW`UJk((tauO-cX4U3l~?50cq8jqF<7hEQ@jr7!a9g3m^+!23sb#>~NyM=oT@ zaq&bsdEfCB{!KE8Drp#T>1#}q^*fPXm0if>rqLwuTrklcr9*<=T)_<{UBF>S9LTVn z%}9iEE3)iz4vD;QGj7sr7|BR+Cr4&ABnx>P(r$baxxY7)goXAXODuwj-xg02Hqnw)tD3{fYw^PC z@9?~!$@s&D^VrGNkfa~IM9CGX>|FB?_rY@^MpKu=Ccp8{r9LFz*qAsUUW}XWxrsL) zk0hCHrBuX4ev}cCB37YleHV#3Umrq-ol5uVYc3cVUFU*0a^IVnqFj_X%f-D?n*+M zKEu-95ybC(AACip7dAg{gl`?J!1HZR;K+maWQ{=;?&;Hvyf!g}qrS&Dq5+A_?@Wx$ zIO6W8!c|-dh_@T*R|LoC8A!aWSu!)}~PNd_;~ZA?ZDZ$f5-29V;4x#Z9N zBoGclOvgBlZ?0Zqu_oqdV*v}UArlnzPA@gR5aT*NxJf{64EM^sBYHYQE%2NB1oMdZ@+ zVshv1B=Ypjb38gWlgyttPxzba-i=t)8VcxoXnpbcNo|LXa{ao4T>ru*IMA>ZUkiPL zKS$0H#Q$@o1b5mSLEQeBkaefu;@rUXIA`HTJTpIv7|)yl)*^~*h-8RP&&ODLXAtq} zW=yJ-(Qw{@M6ojy#Ir3(i!n*-p9dm4gzO-RRIWnB*tmYOnm&cP%Fh<>AjmVMHhhoP;#cMx#h(V);UyNY!5-ZJEpjFAiWRWKV;mEdLG<~CB)g9(pk_q8ZrtS>->$&> zZ2!R?`xwyCPxynUBj{};Q0_Bc*7Y#9e{PG%cY1-j->d`4aFr3GavxIhE(9!?BYD?R zpM2N07Oa~4;%j)`?qxVx-cQ@Pjmau6 z1$jGoE;ju13@`m;16p(s@0(wabLVd7cK>?91izQF4OUZ?)FLBLy z*fA*xw~o5UwTOMrt()dYF0cNMx2EjF?U%d@#YwO2$b$u4h-*O+*~w;+jp=saPA-G_ z`GGG*v<1i48H6s9oIR66yyg~?y`vq-;}OPWc#m3qy`L@V*t=3~7b#{{+&XObt{oY> ztp|zjiO8Ccaip8YAd<;?lbS;VNT0zUu<_PP?A-P<2P(t)UwF^$-Q$M`Co!byo{?nx za#w<BG za_n`;&Ey|=;u#Cl__h^3r8eh2kA;8eVgq= zTM|cP^B8iw!DDPEw{;@_6grVOzW`$0p)EQ3aVYq285q!=_-bk+BL7+-gbK=0t&H)y z)7!9h2P1MX*oI^;wjoPDwIluHo@8r<0^*G%Vw_`1TFq=iUJMH&@hyMgM+v9!nqm_Y z|7j$dh3;d9$>)6LjOTI|-ofeGF+f;f($mMDtlrd|JmJlW?(XJf(M|`lGWsR9)(-|c znUF|qPr6^Wb0B8#dI1Oe07+Vqj(38{!uZ~#bIl*TyoH%S$Y3Sa(}JwmT<^gTxZ4eE zK{?DJiom$H(E~wwKA*fK7~;#%?1;;RKBVeT7gBpZ96V<(__xVqbGyNyAm#YZ{vPC* z={@X9dVu2Hf;Dbqqpt7p`;NY#zLw;C(-9=R;Y3pUFrO&C3?<*g8h|@9Z%j5mUyX|* zBEYi6l6HrziGyw_p$@Y2&`11P9!O5la3l85-2`+)l~fc59BwP{CGLxY4D5R!zj>!e z${Vx+;cZ4d#(u#TUhQF40x>O1AbBIElCz551V;qHPLFVt)}3(F@6Wh&d=MF+6B9~4$kex}^B&eA=IXhiWR+I%0kB~ORYC{8pPDPU6B!paC+rXZfm;^zj=}Xwe z2?8tM)<)nCOMU^%PEnqBR9bjVz)kj4VM&8lB`jpBf+{*rxWqO6Q>qdb6d0@aJ$Fqlh7yCc>K9kVp$wOVuD%=deAk#C83yK zUFk+<90?~T@23#Uj=ALRrdFhP$xFF0{`jw-|_dZsZVyxRykbYbUTaR!NDltXvKT>ojOl1LEGNFJYN6g7-0dw;e~~ z_GH|p?gFcw4{MQmyE8%f6G5(duo!Q+PXpHC*Ct(wtcE2+M|+TOVKNeSqyvfHejfK7 zAcOFr8CkzBlYG>72MsO6mrisb1#g4MwNPKMK4~UoRrC-N@YaB2DZL<5-~_ailvh)2 zAl7UGj&Cx_zS{|;V;~{V8i9pq1%cpOY?kdpDtx@ibBlh0L)5EkOK#zKGV)phL0+HsQKjT~ZO(4jgfOG6W;J!^9iPrC?5O(GgY47HwU!=W|1FdiJ9QWDeN31+O zNZKl28Cm|J4cXIcI^|!o(tU9YymwD)hzKKz=lnKcIflZ__5vob>kU-H4T9qIeB8F6&$1Hy6x zC*D3N(F4UL!_9Yz{75K!yAg)K@_r5 zO`j2Y`?#x5CUZ5(rQ9kncg!_zLGJz3hR7oj|8nUDvC$OZ+&4VCw=viYckr;@5RJr= zn(@ytjkw(xJCa8>`egseM9lw;!4AiHoRNYc;H$->uXvMZvqQ+_SGC~kkK;+N4Dnje zD6%!>CJxvfgUg)F41nDiF}0P>Z5s;$l&!2uHpKmP5wKwpWCs4^MjM6%c-uhi*qD6r z_>2{~{(xR@!O1}eMxM8{1P^Ng!sAT*OM@V_wgh?JfS>S>fB;59^5pSrBR)2>9hfdX zP|ZEq#AH9XF&79T+({1;YvS^dC7PI(pYl{IRrx5v88EFIPjvLqj}JV%Et zyxxgiI28j1;|^@llqj3~koHTQ!Bhv42^*h6G+{|P?QKRH-G3v5##B-Zr-qPiZCBtU zvF}003-Qn&i}1{Y4Im8g2N&3#@O#^WugxHfugdBpk)d2VY_edO`wX^*7|D|)d{~4h zyonOb{rm(gax&&Cp0ZY#&^$_SL@?2@YeZtN2b0+Nksxvh)rkV;HY?VZeRx} zADDpKtHkb$rwQ|p&1?nk!-C8_=m)XvC}JLMOB{}ei>W0)GD#;x**2^INaRhnL-PB( z`ZdU84J5%oAPpp>;lojoVB3?cPeRGt+W8Qigg{K@48AA}q5}uw$})s9*RF^o71`lr z(yQ%wm@P~8*jtl9gdoX;W+WTaJxShNDJ;~VM2#CoEWCc= zdt>CJdskOd(ysxT7M=`V)&MsjJ^_2H+;OMi@wjoX#$=3dD-tuRt6*%X>RNU*CwD`p zk(^$xAmY;~5*f-R78%%Z*kvH+MocR_#x0`3P+o^PH=3lTM}VvMgOURx9rnvW#`j_E zy&cK;PVqq2-S~2v0kQXB$*3J?AVO&i3C~ZC=9vDaoe3=gP?mV!%o9#dWu%p<7Emsd z7@L+r(PR|S9`joY9%>}TBvt9SDw@p1y}%AefPd40xPL636qP})Nj=H&_X-j;XbJ=& zJA`AR)*JasRleB5lL zKB0LCXV#FY4Ft}NApy}v zB(HmffYE{~dvfCa5vY~@5i&qpGb=oG5HjZPkY4g&j3R|y^GkOLJ5r8qD7qkU$YWm` z*}Er;H2zr%7O)>EWYw zJFg*$z6JJU-bH+SiI-s5X+ z4C0QFg6V5F-UDJ3TN^Ua#SWD5CU$7&E2PS`N-r!>2z za!qBK@{?w-w9W0nAwlXWT9j&3D=Es9w=tt$Vp=DnRFmABLAqE^A{5ol$D4!9qzVyj zD^n-YBG=W~AbYj|rBi}s&L&G+Mnjb*0TL8-VJSw*CN_f9F#u|iUj)WbhS5lkiWM~g zRQTgtuEM788sgaw6__eGttC<`L=8r_I86F%8cnj6OoRN+1zdSGcI@Ar6r~qJ@hJ+z z{>xA?X$Z%U`5&m|mqu9#oTo)38h?G%w}gQ%0@Y&EA`cqC64u~8)&h<|bnn`cQCkgh5W*qc$=Y(+($2yWB`@Z7l zt*prGJC0EKt)_zRES22L5w=ZAYCY(&z!=Je#ExG%T6&>Ig&LWjK1cBPqewt@L>`32 zk3xlnh92N-6c1t8Bq&e0!}eWBk0B+*vHutd@bw^5zJpCJwt%3~LD*$V3RRA^uc6X3 z9f})ZBxt%w%b#ZNhC?0M7t9+Ac0vI~IZHv>il)qg;?`_ZG-WI~(=33P8(9lVI3=S4 zIng{Da_<1x`5q1#h{)yRWhCoQny^!T%9w3NVB%89N+Uz^Ebb1b#^gd5$@ z$CazL;%}4Wkif{v_euV++iZwA=0g3n8|k>HJp{iVByiIVfu!`DJul26|F{xdzsj+OI$Iy`k->yU4Y}XifZl$?qdrTC#ik5Wb#JoZ$|4V13z95#m?4<& z4H<>PIn_H+ek}GINZ$A~ghJ9YaQ^|M;qoxj(IHrvm8Bd63g2ER;ZlmFubfPdb}bV4 zL@7n%q2vxqDC7+jl2qy}+JrZ)qa-!U!cP5S`Bsp}SwLtOAkeosHA~Rdi)(Ypsm<5G ztoFkTMs&nuhYqNlO?RMw6lcm45l)^qi%gMECp6-c**wB6u2kY31r4c~7o>HrX#(}% zX;3{661;*qiN0UZ%Ovj&LgDF12Kn$hSnzo}?;Atdc15sC)W_4K6i%H;?BqY7cv%aA z)EolUZ%|tJ!sXbn7aq_+9gc2R*1i7X;%8d!PAN)<$oMDs(rl9uT+jf5YMXSC4KaOJ zjNfngAk->e<>Z3bzRDf4e8(kxxDJa;@Tl1n1e%q0E!fNpXHHTSld@$kg@C~NbYIYf z5jc|+;_}2qa3hn2r&91Ni*q>siYuCOQ>aV@M=5bj-yTqIe#6n?+h(^@!ZRg6dv@E(RU=}bz6l?iBN13fOQfRF7mLO8eIHalDVo%w&?C{!&#l)N{2sMY^pCc$u2mb{ zgYLJm4xy10^_KH?>p}74gkW0guR4ETe@ell#+2NYa@4HRKPoM=)KA5C&(b_c{bwpX z@|<{6K9+VaRlQ7*!BbZK*_z({d!hGt*M%d5Ufs`;=h6A8$M`L}vI6~erzatN!UwH3 z2<3;==ytZP$Sy|b)}KW2Egzy@B@a|w;c*oNMnU5_qg zG(>mr{Nl^7JgjZc%8e+;J`Ksgx*=~>17s^dfrjroi{{TuK#hBqpn(rhqNSHtp_DbZ zQB-do&h@?)%7TJ-8ie4IX}?G^zEo;WgIenw=v`<)Ic9Ajp)_^hO*iTXqS}%+Esp7a{?*)^MN=fW>*V^wU~&_Dm+) zqCFaQt|@@?Y!a|6Ph5x^Oj=JDjkI%0>5ks`97RqkbCK))(@3)aBx=(7Hqz;tO$Ruq zeXqzuj;j`tzY8^&^+H!ilp&uv(@}S`&A|G3XxZgCAZ%NK4Cg?|W}y$$CZWvg{-|*6 z3K+8leXsCFgO~b(FwR8pH&pomclCkW%aL)5$H?o(Aq4scJZKEkkOOqwjIOVcaf|$n zxgQnhkxi2Y6zdRWmnuhK*f`-RE-uT__)a}hv9uQ~)CHJ$1ogOh0{x@A7HKUCN3$LL zKv?I&&gW6@!3WUc&o|H$D}C-l#tIOpc_6$qkgso76!D-3nijqs1miBc^Lr_pw|hET zce)s@UpF7+%^ihiJerFpzZ!ztratq4frkLIr6{fc(~x_6gV5{F-O*+DUm@3y428uG z!_h0ykdSq1Y;WS#!Blh_vPP)n(D0D-u%W=L`zUD8edK-T2--i=fioQD!Wp-?A}s2y zN(80N04*tl6Xl@qM?atm&1CSzrc4mWK}*dLdoLJhJQz+s8fonh1agf>OFk8$)HRui ztlo%D4Cw_t?T*G>YUhCrhXtX>0fwmg@NXftMy6;((Q33asVQjEnGpJC-JuU?T?xgu zv6kC`x%7jL^7(cW;;k}5{H+r7GUpNK=4sTW>Nje$VKZ_&@Ej>im%(2$$mS&UWL+~f zFq#KqH33Br1FJCyI9Y|>77(1)MdN?A3*&+78ttHT~8I&Opr53J5gO;uUv^!Akp%^+SJK3cR`b!0J(yNQG z;vMukQWqS60oT;F295c63SE7248&n6oPU)dUyCg^qNiVYAap6p?sFEjVjzgfDuEfr z(&?~KV_uik@9{1r`?zf~>V@GGJbct{o{y z(RNkB$(&NApvk8;ffXBsHoNT)nJ_;LHX8-zuTr3`iklMhIKT_`x)2h$Jpza5?~DPX z@&K50128w>W<0Xy&YoEa(!U5Mm7?v-hoIAjm(kkS2L%Ld8212;SZ@8gmx2NAS$R|-&E7rJKpS2HBXv0yw+7U2Jy@Vyu-Rwe&rYy%f8bm1 z7;I?{TIF>d>~bd8jXijUP8x*Y(O2x zv0RUfPTbE4KkY#f2BWQy{ZYW?{?Iuc)OVDCwalgo#jvAtLbK-}W@o{5oE3KacJvTv z*L*NpW6+71lPENFJL(kt1AOxaSY#E@@gCw{Jc5Nx;jz(qVD=Vt$Z`W1i082LG9X19 zcwv{2uTSHGJB+&+l3KV6`0@p`eH_f)3^WE~Kb(apnG4YU1^(tJzm$auZ9oxPf-L8Q z27hO@exX;jk6qP6L z1!nI8dfGa1*?l&mN;d@Z83zkFgWQ&*F*ey~z#MNdN;A}0R%9E1Q|kw3QNj6^fDkVN zVIKsvD*~@BTn+|5 zv=~jVA{4_AWnA2B;LH&4j^${}O5q6H(q05+are`CF|{^sxh-2b^@*TA4aOTB=DVy<#;q0Fu0egA{DYsk!)0U4M zj+dc6gOy-i)(aYT_RJVC79}B8nI>R!e?csD7Q`wUS!^m4WOie^J-0IYJ8EkC0>*2@ zH)ArW-3C|XE5uC-8+aD4_WF-y2FT39X8SDt6 zFh+>(wjGWYlztY<01kBlk9z|a842#>1pG1Fv6=TkgqDG?3IMYvo zfDGnBusb^>zU6Z81j8YcTOveu@hxWx!5Xzt?IerA+ZG|kP7IooiPFO+z&^*&sO7)W zgKBebQ@b)_@IO8fL(UeMbJ=|(U+`8pS=nIH9`MgrkN`daGqe@Bq07;LpZbzukQ*8T zW>$ZJAZQFsUJ0E4CD_B6`S&2soDP<{Et;AZ2i|fuU^5M5u@daTVW81Fkg-DGTqq6m zP*A%cLObeq5lEB|7QY-cOHT%0Z!_i~zEor5+%^m{qkRGa=}J43?6?e92C|`lz>i-5 zZP4aGFu5kJe*$(bz`aHQCluV+OFY;#ImkyK%smgEPlAt`0^;!<<}7z=^08$ zJc%mWKSwR9KY|@tj-0Pu7qYTX11`bZ1~N{5=^{FZ&QNT6DD-Kc9k*s>)w+5&oDj{4iuLpuO#efaK#l%rhfTvH8>dmqwGT8_3feTUAjdIs2hhoJu?gv@%}X!BaI4wrzF zQUyF(JcSB-FGH#gpD9K|l-z`5SmZ0{+H=6a93*-V(%)GL_I<7(M``_apit<^ zU5dyDp1(wYoIZjbdj+;@A$nW87$U0OLXty$FSSz-kL^bV9>-8@?sJ;hpk~ugAR32O z`~N|EWL2<2wCYoNXM@a?}KZh)2nmvmr(5ILQKJnc8>kL;dT zKvp>iWMm{70kMX~ChCEzfTLT5*rHFzU^w$G@L>DktQ8P8ToSidr#f5MlVmCT!Mo{j5gs_;JrkzRC z(6}uX5aI5E!1F8Q)HUEhVo+koF=*uB=G;){8noey6wJ_NQ2Zyz^V<^$t#3gtsl!n* z+n9a@49F27(Vu*36$tRAzme!%Bj-kTI|8xX9Y|zqz!%R0`t1OFe+&iKyK-sMf1v*T zCjPGQMWDWX)Ff54FQzN6rpelmm29;3PLyTP??LtFDC+=+c(fPzm2s`V0ZaL8vUMDI8ij6q7nVe_f7&`JRa{$?=1htR_KZNlH%S~qU6 z(wMuyz#WXND@Q{&(|3bF%v`vP_7>d9h5FpZpW58hB2(^SL?yV=m1s)eSKvwxqYszv z05uwNC2~g&3M%lO>QBHbBM>liZp_CvT!Trk0Et41M22#J);{n<7b#ZY)Q&*TbVcCp z;{YA-4dN$h#N!wd-Z{-o}u#doJaYJXsDx6hQnkP|Tit z-Pwn0c*}$vc9Y?>$F$~jddj$D{V!mHE#Ya9C-=5CoU_UdgzpB-l~S`vF-cXz^LY4< zsx4QZI2~fabKtvI(ToBz!p-1!s(=UQAny4NF~12%R-Z*3gYCFpCAM5QLkEZlEVzk^ zU+89GmLQ!G4{k!((~=7srBI?$%gIi(nL zTOV^kuK)#i__H-PHEk05>=C8*+rRqO2qX2b>%;k*frGyanZYTrBMetH&6YDe{~cM? z=m<2NKko&4xB4z}?WirtXvv{Q0ypOT)Z;pgD??vAUPBbif~IN%=8d_Ym5n%tiQrsZ zns8laegjJY^`#@TQLEE$p+r#&5@Bu*Uz~Ufp^OtZe_nI0OI0m&tWrm~v|jYb=`diU z$F-9nt}A~Uv}qF96ic9d3(mm|jx*&43M>5$*%k}*c?kll57fU*vs;LY|3HrN3M}1k zuyVD~_XB8l7gW0umOBjD8U=i{ zLnzoSjGq8mSQ(@!3!z}=PNR#CO8*Ns+{}D$SS64%^>^bo{`i7+T|a}`b%I!>#dXkG z8CNXzcuW1{Z2(f<~mpzmdZk9XU$jH$G3$!k15=Jos9`_8(n2*T26dcYOr|c18w5 z8wlKr1WZCf?%E0qnwSK_dnjPoj%yti0&*Dy!PYyF_9Kw#UW6!a0c>`TlEGPP$)`Xn zN6CqfXSBF0Iya%la~EW!2s!<86U?|N@Ysr*--U(3nKQU}KdzR$R9@ZGg=IBLXtj95VbJs3DXU4|Ly3^^9*mNmkVFK=%(26k2heYivM)6+#Ty zd$1AM?$%uYjgEr;ZD0N#J)dC?iv@GDP$aONZqN*@h5_fYL&436aRg7>ge%DMgJ0y> za6Ql$P%AeOt+w3OMS@lo^ozE0;`+X{!pVjIzqcQq*b%02Y*-gC5W;u9zj zvS7cT(hLE_BLc*v17|SXfrBib8=Bsj>ye_SC)I$hkJ~{76@y=B$*oY?aCh&>DcK4V zi>3C0C;K#DIc%~BP4<6{N=LszjdV=7c{hx(x#sdLm9zMQ24{28Sc_x3vR2P z0T|Xs+%Cs3Fy~=hZL9MrJnN&7p8e|EoU_=}jxQ^aO-fAoUX*_1zRAUAs6v@Jhbi`G4#{Zq2}6y1?TibIw?y%Vov}g8vENZec499sUia z8hwz8HXt?jaH;^lxBz(bTF)B3hUp7*JWWg6XoBuXNY@E#0pHvvQ)Au+loB&(Es^w0kS5l+)DV6!9G)4LGh{(=&A zd+y`3X54(823%ibFYcuLDta|{I{bA9(m4qe#Rcg5x1&K|FVvDp;eIE>xWeA8Veu#E z{6E{EI`#-@8|4a5a?H>nn$)m`e) zb|`&Tw~X&I1AgDQdz!4;Te49c66GYfR1Uc}Gw1K0hr=%ebH6XPVWlk~ai=dy*UYQ? z)xS<aY$z1i1N?B=ux{*<|O&*+;Ju)ISBQ1Sc)1_i}pr%_tCN-yD>fmAY zr{PbRiM^VumKd?OYewd^4^K-Sx>0kh3M+*uN_d&n?@mn1N$GbF$x{KL#MBYlS+kCbeQh*-xqug6(wTTPF#&@(hU zGchc#OLW(;jw>tl;hq?4{?%(qg8G-hzcsSZm|n<5e}w95rJh_Oy!rQCwZ8r@ppQM& z<4sHqd-9tMbj_vB*0x(se>mKWS8@h2t#hXOD)<>RX8A!z(5pQB@=U$nYi-2dOcO(M zC#TN`PKDE3(&{cjtGoK_-&@k8m(Hvk#QWUWUacCjPnPaYg@37vdt3KWr?Fh!`mcbu zmX?>Y+Wh`vnLR&aqihWC^GV`C#ovRUQ6)3s@A0xczSB%uJYSj$m(d8H@&ni?t@Z2H z@s+b=1NaheU2Fa>FPjAm$#drUv5;KMdq-Ga$HS|LlPXS7qmQ?L+PmlHSZP5V;lf>ZeYYh=q|x{b{K`1=t0 z{Fkp@B^ycm(iCxV7W8R07YCJXzPeP_oR5h$45ZY`l^q`v#o0r3Qr!;p45o92@u)=R zszMuOEm^vELpbuI*7T3>%*#5{CToiENS%Dx3|b3{WxeSy!n63X0c;$_E0Z$Eo-Wo? z2ILPGhAm6eoxtDcWeK_fQZ%XE8&3x%)y)z3lk)!&G^caZ>jt}0u>0R&Y#k_w3P=tW zmBbmcrs|H;o`!K0w4&@$2tRA*E3LFr=pl{k=DPEPyUP`-!!GhrR%`Gk1>e(Ej?BKK zORv<|p+23|9i~@|t(VvDaJAvryUIB;okn+L_1ksa;l3U>SSdR1W^w(F58X8#DPNW) z)$iB^cTrR)z>FmG+D-L4st;*8a`i6Vso!z(gL)}b_l-xS^*b)Q@MUgt&h${KWI_G* zDY5GIYvsIfm63QyLiA*Hdu`vJtLl%ptWvjE2gkRo-+t$X&~BE0P${q9USTIyA7R(0 z5yR?t4D6`s__X-vr}`Z&l$wrXF?WvD@Ay`!>6kLg`BD9j%7^L>Q^OCt|2kZrId9Sf z4sB4e??n9xC|uqC=|qb*_1m@c)$Osv-Yl!%ZcwRi|Mm)E9R(3E(B1WyY^WY;I)ZZ6 z&#u4uICEX~jG9z=a{Z35a7{;x>yr-D?@)};bhPh#+jV-Y+mEwTI@NDKJ44<6ycy$H zzkSvTb$g|C*oyk?3%?5OW-aeKtgGK%;UQO3qrJi{uYSjd?&=O{8{2_1`1mAwd-=*0 z>eyo@zeOf%!e2>}pO?P;9deOhoh%=wm3vpgJNA`Z*|pQ$H+rB_9gM`QFaBVdZW4Oc zLxMyaa3=RGzcEGbrUiGC@`Ey12mWJ<92xqrQE$L7b3Sy+;<;40pZy)pY&*5ThWlRB zFRQ@n;i|G!IV;tnH?A6NR^J=;ZjFv`C7Ajavv9>cesPvO7!XNasDGKd=>*;>ot10l ztgSz8i2Rmx+CKO1d|8HkoKg(XG&?bOPnJsrA1>Zb$HFxvzbRfU{PTle;m88w|IoOi z?ut})Cd&+_F--qg-KPk%6`CtcSh)M*8=b_khb2siaK88JUym0rozd5PxNqC%9&hZ& z#_|UZAne$`Rgd1jGFpAP-dlYLJqFGwDdJt^x?Oa4<_Mu0y-h`*&Qo82Mkq}@>HTMN zK8t9Z^GSZJo}}v$-YHgQ$^S9Z6>jRw0_t3dj??Co16XTIAx0KiNRebt88Ub1U5={$ z1LPkh>^$DcLeY?)n98opR%qeVO zZdmIgjxs=ggY96veqZ@8o8wn(ub?!yvhig4YvSkNluU3GV1X*Q_41WNk>cB{)}i*z-|byx*Sgf*JX7YxGs;=z;&sc z2ChrlG;m#}rh)5HGbTgbQkREm2Gpfq8n`aq(!g~|mIj85O86vFB=g$eaMPRmeDzp~ zi^|yCAVaI%%Eyp@7hxck9$Ee^`ERAbvN~j6!|}T8Ek#66gL#L zQiMgV6k$;-MOf5I5f-&lghj0sVNok7{FhdW14OM9VNokZSky`p7PV4@MXeNJQ7c7Q z)JlMb{p+++6lzf`#R;NTim<4aA}ng92#Z=N!lG7+u&9;Hf3;E^AZn!ui&`ndqE?Eq zsFfltYNZH+R?>)Tx|Z5ecif-=M;(1a-%w>_Ztxgx5H&Y=^#9}EGg>tJlrug_e0i@}1DR@SXM^-aTDH>^P!Z(oC z`P|+lOX|?roaF4V z{{4aKLL(~;*L3O+%}ySck(`;BliJ@~a#x?;N653h|4r?_@&TFxVi*3Mk?t_L`ZBe0 zxNUKBI_sq523JO?uSOKhi|WgW8S#o=;ci4L#$U6OaP_|>BNi5!V(`2T7_pG3xeZOM z9BQsD7E6VU`UW+zI0*Npfu2}1Vo{K>2QO>Gh^0WL6@U#GPTV*aKvxEfp!%LAAyWw;Fu4Vb2Lk&HHE$kU74A(MbMiKQb zgL^eh`a}Qk1$c?s=W0#9LWc{l@jS_wc&b1bda%+-7qx46CllS}`~`DEnab1D5Cj`a zd&}qrkTU*xy5V$HeyHJPFF)b-L~%1=p+4=^JNF4!zl+I%_&&6~M_IZH3z*I%rz z)H0A`@K19Mdm6UJ>Uk{Fd+mKaK043P&wjY3VB1K8toYvH#*5*qWqF3$Qj`nR115%) broadleaved evergreen or semi-deciduous forest (>5m),40,0.6,0.3,0.0,308.0,0.039,0.5,0.0,0.0065,5.40 -50,Closed (>40%) broadleaved deciduous forest (>5m),50,0.8,0.8,0.0,430.0,0.036,0.5,0.0,0.001,19.20 -60,Open (15-40%) broadleaved deciduous forest/woodland (>5m),60,0.8,0.4,0.0,430.0,0.036,0.5,0.0,0.01,5.40 -70,Closed (>40%) needleleaved evergreen forest (>5m),70,0.8,0.1,0.0,382.0,0.045,0.5,0.0,0.001,19.20 -90,Open (15-40%) needleleaved deciduous or evergreen forest (>5m),90,0.8,0.4,0.0,382.0,0.045,0.5,0.0,0.01,5.40 -100,Mosaic forest or shrubland (50-70%) / grassland (20-50%),100,0.8,0.3,0.0,406.0,0.039,0.5,0.0,0.02,5.40 -110,Mosaic forest or shrubland (50-70%) / grassland (20-50%),110,0.6,0.456,0.0,286.0,0.07,0.2,0.0,0.015,5.40 -120,"Mosaic grassland (50-70%) / forest or shrubland (20-50%) ",120,0.6,0.5,0.0,179.0,0.127,0.05,0.0,0.03,1.97 -130,Closed to open (>15%) (broadleaved or needleleaved evergreen or deciduous) shrubland (<5m),130,0.6,0.5,0.0,432.0,0.07,0.1,0.0,0.035,1.97 -140,Closed to open (>15%) herbaceous vegetation (grassland savannas or lichens/mosses),140,0.6,0.24,0.0,457.0,0.09,0.0,0.0,0.05,1.97 -150,Closed to open (>15%) herbaceous vegetation (grassland savannas or lichens/mosses),150,0.6,0.02,0.0,137.0,0.04,0.04,0.0,0.35,1.97 -160,Closed to open (>15%) broadleaved forest regularly flooded (semi-permanently or temporarily) - Fresh or brackish water,160,0.6,0.3,0.0,308.0,0.04,0.1,0.0,0.001,5.40 -170,Closed (>40%) broadleaved forest or shrubland permanently flooded - Saline or brackish water,170,0.8,0.8,0.0,308.0,0.036,0.2,0.0,0.0005,19.20 -180,Closed to open (>15%) grassland or woody vegetation on regularly flooded or waterlogged soil - Fresh brackish or saline water,180,0.6,0.15,0.0,107.0,0.127,0.01,0.0,0.04,1.97 -190,Artificial surfaces and associated areas (Urban areas >50%),190,0.6,0.011,0.75,179.0,0.04,0.01,0.0,0.0,1.0 -200,Bare areas,200,0.6,0.01,0.0,0.0,0.04,0.0,0.0,0.0,1.0 -210,Water bodies,210,0.7,0.08,0.0,0.0,0.04,0.0,1.0,0.0,1.0 -220,Permanent snow and ice,220,0.6,0.01,0.0,0.0,0.04,0.0,0.0,0.0,1.0 -230,No Data,230,-999.0,-999.0,-999.0,-999.0,-999.0,-999.0,-999.0,-999.0,-999.0 +globcover,description,landuse,Kext,N,PathFrac,RootingDepth,Sl,Swood,WaterFrac,alpha_h1,USLE_C,Cov_River +11,Post-flooding or irrigated croplands (or aquatic),11,0.6,0.2,0,390,0.127,0.01,0,0,0.2,0.97 +14,Rainfed croplands,14,0.6,0.2,0,390,0.127,0,0,0,0.35,1.97 +20,Mosaic cropland (50-70%) / vegetation (grassland/shrubland/forest) (20-50%),20,0.6,0.44,0,397,0.127,0.01,0,0,0.27,1.97 +30,Mosaic vegetation (grassland/shrubland/forest) (50-70%) / cropland (20-50%) ,30,0.6,0.56,0,400,0.127,0.01,0,1,0.25,1.97 +40,Closed to open (>15%) broadleaved evergreen or semi-deciduous forest (>5m),40,0.6,0.3,0,308,0.039,0.5,0,1,0.0065,5.4 +50,Closed (>40%) broadleaved deciduous forest (>5m),50,0.8,0.8,0,430,0.036,0.5,0,1,0.001,19.2 +60,Open (15-40%) broadleaved deciduous forest/woodland (>5m),60,0.8,0.4,0,430,0.036,0.5,0,1,0.01,5.4 +70,Closed (>40%) needleleaved evergreen forest (>5m),70,0.8,0.1,0,382,0.045,0.5,0,1,0.001,19.2 +90,Open (15-40%) needleleaved deciduous or evergreen forest (>5m),90,0.8,0.4,0,382,0.045,0.5,0,1,0.01,5.4 +100,Mosaic forest or shrubland (50-70%) / grassland (20-50%),100,0.8,0.3,0,406,0.039,0.5,0,1,0.02,5.4 +110,Mosaic forest or shrubland (50-70%) / grassland (20-50%),110,0.6,0.456,0,286,0.07,0.2,0,1,0.015,5.4 +120,Mosaic grassland (50-70%) / forest or shrubland (20-50%) ,120,0.6,0.5,0,179,0.127,0.05,0,1,0.03,1.97 +130,Closed to open (>15%) (broadleaved or needleleaved evergreen or deciduous) shrubland (<5m),130,0.6,0.5,0,432,0.07,0.1,0,1,0.035,1.97 +140,Closed to open (>15%) herbaceous vegetation (grassland savannas or lichens/mosses),140,0.6,0.24,0,457,0.09,0,0,1,0.05,1.97 +150,Closed to open (>15%) herbaceous vegetation (grassland savannas or lichens/mosses),150,0.6,0.02,0,137,0.04,0.04,0,1,0.35,1.97 +160,Closed to open (>15%) broadleaved forest regularly flooded (semi-permanently or temporarily) - Fresh or brackish water,160,0.6,0.3,0,308,0.04,0.1,0,1,0.001,5.4 +170,Closed (>40%) broadleaved forest or shrubland permanently flooded - Saline or brackish water,170,0.8,0.8,0,308,0.036,0.2,0,1,0.0005,19.2 +180,Closed to open (>15%) grassland or woody vegetation on regularly flooded or waterlogged soil - Fresh brackish or saline water,180,0.6,0.15,0,107,0.127,0.01,0,1,0.04,1.97 +190,Artificial surfaces and associated areas (Urban areas >50%),190,0.6,0.011,0.75,179,0.04,0.01,0,1,0,1 +200,Bare areas,200,0.6,0.01,0,0,0.04,0,0,1,0,1 +210,Water bodies,210,0.7,0.08,0,0,0.04,0,1,1,0,1 +220,Permanent snow and ice,220,0.6,0.01,0,0,0.04,0,0,1,0,1 +230,No Data,230,-999,-999,-999,-999,-999,-999,-999,-999,-999,-999 diff --git a/hydromt_wflow/data/lulc/vito_mapping.csv b/hydromt_wflow/data/lulc/vito_mapping.csv index b2edf877..3e00c0cd 100644 --- a/hydromt_wflow/data/lulc/vito_mapping.csv +++ b/hydromt_wflow/data/lulc/vito_mapping.csv @@ -1,24 +1,24 @@ -vito,description,landuse,Kext,N,PathFrac,RootingDepth,Sl,Swood,WaterFrac,USLE_C,Cov_River -20,Shrubs,20,0.07,0.5,0.0,410,0.07,0.1,0.0,0.05,1.97 -30,Herbaceous vegetation,30,0.6,0.15,0.0,106.8,0.1272,0.0,0.0,0.045,1.97 -40,Cultivated and managed vegetation/agriculture (cropland),40,0.6,0.2,0.0,390.4,0.1272,0.0,0.0,0.3,1.97 -50,Urban / built up,50,0.7,0.011,0.9,257.4,0.04,0.01,0.0,0.0,1 -60,Bare / sparse vegetation,60,0.6,0.02,0.0,10.7,0.04,0.04,0.0,0.25,1.97 -70,Snow and Ice,70,0.6,0.01,0.0,0.0,0.0,0.0,0.0,0.0,1 -80,Permanent water bodies,80,0.7,0.01,0.0,0.0,0.0,0.0,1.0,0.0,1 -90,Herbaceous wetland,90,0.6,0.15,0.0,106.8,0.1272,0.0,0.0,0.05,1.97 -100,Moss and lichen,100,0.6,0.085,0.0,136.9,0.04,0.0,0.0,0.04,1.97 -111,Closed forest evergreen needle leaf,111,0.8,0.4,0.0,382.2,0.045,0.5,0.0,0.001,19.2 -112,Closed forest evergreen broad leaf,112,0.8,0.6,0.0,294.1,0.036,0.5,0.0,0.0065,19.2 -113,Closed forest deciduous needle leaf,113,0.8,0.4,0.0,215.3,0.045,0.5,0.0,0.001,19.2 -114,Closed forest deciduous broad leaf,114,0.8,0.6,0.0,429.8,0.036,0.5,0.0,0.001,19.2 -115,Closed forest mixed,115,0.8,0.5,0.0,406,0.0393,0.5,0.0,0.0065,19.2 -116,Closed forest unknown,116,0.8,0.5,0.0,406,0.0393,0.5,0.0,0.0065,19.2 -121,Open forest evergreen needle leaf,121,0.8,0.45,0.0,396.1,0.0575,0.5,0.0,0.01,5.4 -122,Open forest evergreen broad leaf,122,0.8,0.55,0.0,352,0.053,0.5,0.0,0.01,5.4 -123,Open forest deciduous needle leaf,123,0.8,0.45,0.0,312.6,0.0575,0.5,0.0,0.01,5.4 -124,Open forest deciduous broad leaf,124,0.8,0.55,0.0,419.9,0.053,0.5,0.0,0.01,5.4 -125,Open forest mixed,125,0.8,0.5,0.0,408,0.0553,0.5,0.0,0.01,5.4 -126,Open forest unknown,126,0.8,0.5,0.0,408,0.0553,0.5,0.0,0.01,5.4 -200,Open sea,200,0.7,0.01,0.0,0.0,0.0,0.0,1.0,0.0,1 -0,No data,0.0,-999.0,-999.0,-999.0,-999.0,-999.0,-999.0,-999.0,-999.0,-999.0 +vito,description,landuse,Kext,N,PathFrac,RootingDepth,Sl,Swood,WaterFrac,alpha_h1,USLE_C,Cov_River +20,Shrubs,20,0.07,0.5,0,410,0.07,0.1,0,1,0.05,1.97 +30,Herbaceous vegetation,30,0.6,0.15,0,106.8,0.1272,0,0,1,0.045,1.97 +40,Cultivated and managed vegetation/agriculture (cropland),40,0.6,0.2,0,390.4,0.1272,0,0,0,0.3,1.97 +50,Urban / built up,50,0.7,0.011,0.9,257.4,0.04,0.01,0,1,0,1 +60,Bare / sparse vegetation,60,0.6,0.02,0,10.7,0.04,0.04,0,1,0.25,1.97 +70,Snow and Ice,70,0.6,0.01,0,0,0,0,0,1,0,1 +80,Permanent water bodies,80,0.7,0.01,0,0,0,0,1,1,0,1 +90,Herbaceous wetland,90,0.6,0.15,0,106.8,0.1272,0,0,1,0.05,1.97 +100,Moss and lichen,100,0.6,0.085,0,136.9,0.04,0,0,1,0.04,1.97 +111,Closed forest evergreen needle leaf,111,0.8,0.4,0,382.2,0.045,0.5,0,1,0.001,19.2 +112,Closed forest evergreen broad leaf,112,0.8,0.6,0,294.1,0.036,0.5,0,1,0.0065,19.2 +113,Closed forest deciduous needle leaf,113,0.8,0.4,0,215.3,0.045,0.5,0,1,0.001,19.2 +114,Closed forest deciduous broad leaf,114,0.8,0.6,0,429.8,0.036,0.5,0,1,0.001,19.2 +115,Closed forest mixed,115,0.8,0.5,0,406,0.0393,0.5,0,1,0.0065,19.2 +116,Closed forest unknown,116,0.8,0.5,0,406,0.0393,0.5,0,1,0.0065,19.2 +121,Open forest evergreen needle leaf,121,0.8,0.45,0,396.1,0.0575,0.5,0,1,0.01,5.4 +122,Open forest evergreen broad leaf,122,0.8,0.55,0,352,0.053,0.5,0,1,0.01,5.4 +123,Open forest deciduous needle leaf,123,0.8,0.45,0,312.6,0.0575,0.5,0,1,0.01,5.4 +124,Open forest deciduous broad leaf,124,0.8,0.55,0,419.9,0.053,0.5,0,1,0.01,5.4 +125,Open forest mixed,125,0.8,0.5,0,408,0.0553,0.5,0,1,0.01,5.4 +126,Open forest unknown,126,0.8,0.5,0,408,0.0553,0.5,0,1,0.01,5.4 +200,Open sea,200,0.7,0.01,0,0,0,0,1,1,0,1 +0,No data,0,-999,-999,-999,-999,-999,-999,-999,-999,-999,-999 diff --git a/hydromt_wflow/data/parameters_data.yml b/hydromt_wflow/data/parameters_data.yml index 99231a74..53004431 100644 --- a/hydromt_wflow/data/parameters_data.yml +++ b/hydromt_wflow/data/parameters_data.yml @@ -1,7 +1,7 @@ corine_mapping_default: data_type: DataFrame driver: csv - kwargs: + driver_kwargs: index_col: 0 dtype: landuse: int16 @@ -12,6 +12,7 @@ corine_mapping_default: Sl: float32 Swood: float32 WaterFrac: float32 + alpha_h1: int16 USLE_C: float32 Cov_River: float32 meta: @@ -23,7 +24,7 @@ corine_mapping_default: esa_worldcover_mapping_default: data_type: DataFrame driver: csv - kwargs: + driver_kwargs: index_col: 0 dtype: landuse: int16 @@ -34,6 +35,7 @@ esa_worldcover_mapping_default: Sl: float32 Swood: float32 WaterFrac: float32 + alpha_h1: int16 USLE_C: float32 Cov_River: float32 meta: @@ -45,7 +47,7 @@ esa_worldcover_mapping_default: globcover_mapping_default: data_type: DataFrame driver: csv - kwargs: + driver_kwargs: index_col: 0 dtype: landuse: int16 @@ -56,6 +58,7 @@ globcover_mapping_default: Sl: float32 Swood: float32 WaterFrac: float32 + alpha_h1: int16 USLE_C: float32 Cov_River: float32 meta: @@ -67,7 +70,7 @@ globcover_mapping_default: vito_mapping_default: data_type: DataFrame driver: csv - kwargs: + driver_kwargs: index_col: 0 dtype: landuse: int16 @@ -78,6 +81,7 @@ vito_mapping_default: Sl: float32 Swood: float32 WaterFrac: float32 + alpha_h1: int16 USLE_C: float32 Cov_River: float32 meta: @@ -86,10 +90,33 @@ vito_mapping_default: source_version: 1.0 path: lulc/vito_mapping.csv +glcnmo_mapping_default: + data_type: DataFrame + driver: csv + driver_kwargs: + index_col: 0 + dtype: + landuse: int16 + Kext: float32 + N: float32 + PathFrac: float32 + RootingDepth: float32 + Sl: float32 + Swood: float32 + WaterFrac: float32 + alpha_h1: int16 + USLE_C: float32 + Cov_River: float32 + meta: + category: landuse + source_info: landuse parameters based on corine classification (https://land.copernicus.eu/pan-european/corine-land-cover/clc2018) + source_version: 1.0 + path: lulc/glcnmo_mapping.csv + koppen_geiger_mapping_default: data_type: DataFrame driver: csv - kwargs: + driver_kwargs: index_col: class meta: category: landuse @@ -100,7 +127,7 @@ koppen_geiger_mapping_default: regr_chelsa_mapping_default: data_type: DataFrame driver: csv - kwargs: + driver_kwargs: index_col: source meta: category: hydrography @@ -111,7 +138,7 @@ regr_chelsa_mapping_default: roughness_river_mapping_default: data_type: DataFrame driver: csv - kwargs: + driver_kwargs: index_col: 0 dtype: N_River: float32 @@ -124,7 +151,7 @@ roughness_river_mapping_default: riverbedsed_mapping_default: data_type: DataFrame driver: csv - kwargs: + driver_kwargs: index_col: 0 dtype: D50_River: float32 diff --git a/hydromt_wflow/wflow.py b/hydromt_wflow/wflow.py index 903d9575..e0d2b99e 100644 --- a/hydromt_wflow/wflow.py +++ b/hydromt_wflow/wflow.py @@ -867,27 +867,28 @@ def setup_lulcmaps( self, lulc_fn: Union[str, xr.DataArray], lulc_mapping_fn: Union[str, Path, pd.DataFrame] = None, - lulc_vars: List = [ - "landuse", - "Kext", - "N", - "PathFrac", - "RootingDepth", - "Sl", - "Swood", - "WaterFrac", - ], + lulc_vars: Dict = { + "landuse": None, + "Kext": "input.vertical.kext", + "N": "input.lateral.land.n", + "PathFrac": "input.vertical.pathfrac", + "RootingDepth": "input.vertical.rootingdepth", + "Sl": "input.vertical.specific_leaf", + "Swood": "input.vertical.storage_wood", + "WaterFrac": "input.vertical.waterfrac", + "alpha_h1": "input.vertical.alpha_h1", + }, ): """ Derive several wflow maps based on landuse-landcover (LULC) data. - Lookup table `lulc_mapping_fn` columns are converted to lulc classes - model parameters based on literature. The data is remapped at its original - resolution and then resampled to the model resolution using the average - value, unless noted differently. + Lookup table `lulc_mapping_fn` columns are converted to lulc classes model + parameters based on literature. The data is remapped at its original resolution + and then resampled to the model resolution using the average value, unless noted + differently. Currently, if `lulc_fn` is set to the "vito", "globcover", "esa_worldcover" - or "corine", default lookup tables are available and will be used if + "corine" or "glmnco", default lookup tables are available and will be used if `lulc_mapping_fn` is not provided. Adds model layers: @@ -900,20 +901,27 @@ def setup_lulcmaps( * **PathFrac** map: The fraction of compacted or urban area per grid cell [-] * **WaterFrac** map: The fraction of open water per grid cell [-] * **N** map: Manning Roughness [-] + * **alpha_h1** map: Root water uptake reduction at soil water pressure head h1 + (0 or 1) [-] Parameters ---------- lulc_fn : str, xarray.DataArray Name of RasterDataset source in data_sources.yml file. lulc_mapping_fn : str, Path, pd.DataFrame - Path to a mapping csv file from landuse in source name to - parameter values in lulc_vars. If lulc_fn is one of {"globcover", "vito", - "corine", "esa_worldcover"}, a default mapping is used and this argument + Path to a mapping csv file from landuse in source name to parameter values + in lulc_vars. If lulc_fn is one of {"globcover", "vito", "corine", + "esa_worldcover", "glmnco"}, a default mapping is used and this argument becomes optional. - lulc_vars : list - List of landuse parameters to keep. - By default \ -["landuse","Kext","N","PathFrac","RootingDepth","Sl","Swood","WaterFrac"] + lulc_vars : dict + Dictionary of landuse parameters in ``lulc_mapping_fn`` columns to prepare + and their internal wflow name (or None to skip adding to the toml). By + default \ +{"landuse": None, "Kext": "input.vertical.kext", "N": "input.lateral.land.n", + "PathFrac": "input.vertical.pathfrac", "RootingDepth": + "input.vertical.rootingdepth", "Sl": "input.vertical.specific_leaf", "Swood": + "input.vertical.storage_wood", "WaterFrac": "input.vertical.waterfrac", + "alpha_h1": "input.vertical.alpha_h1"} """ self.logger.info("Preparing LULC parameter maps.") if lulc_mapping_fn is None: @@ -934,12 +942,17 @@ def setup_lulcmaps( da=da, ds_like=self.grid, df=df_map, - params=lulc_vars, + params=list(lulc_vars.keys()), logger=self.logger, ) rmdict = {k: v for k, v in self._MAPS.items() if k in ds_lulc_maps.data_vars} self.set_grid(ds_lulc_maps.rename(rmdict)) + # Add entries to the config + for name, wflow_param in lulc_vars.items(): + if wflow_param is not None: + self.set_config(wflow_param, name) + def setup_laimaps(self, lai_fn: Union[str, xr.DataArray]): """ Set leaf area index (LAI) climatology maps per month [1,2,3,...,12]. diff --git a/hydromt_wflow/wflow_sediment.py b/hydromt_wflow/wflow_sediment.py index 0a643978..747a75e3 100644 --- a/hydromt_wflow/wflow_sediment.py +++ b/hydromt_wflow/wflow_sediment.py @@ -2,7 +2,7 @@ import logging from pathlib import Path -from typing import List, Optional, Union +from typing import Dict, List, Optional, Union import geopandas as gpd import numpy as np @@ -309,47 +309,41 @@ def setup_lulcmaps( lulc_fn: Union[str, Path, xr.DataArray], lulc_mapping_fn: Union[str, Path, pd.DataFrame] = None, planted_forest_fn: Union[str, Path, gpd.GeoDataFrame] = None, - lulc_vars: List[str] = [ - "landuse", - "Cov_River", - "Kext", - "N", - "PathFrac", - "Sl", - "Swood", - "USLE_C", - "WaterFrac", - ], + lulc_vars: Dict = { + "landuse": None, + "Kext": "input.vertical.kext", + "PathFrac": "input.vertical.pathfrac", + "Sl": "input.vertical.specific_leaf", + "Swood": "input.vertical.storage_wood", + "USLE_C": "input.vertical.usleC", + }, planted_forest_c: float = 0.0881, orchard_name: str = "Orchard", orchard_c: float = 0.2188, ): """Derive several wflow maps based on landuse-landcover (LULC) data. - Currently, ``lulc_fn`` can be set to the "vito", "globcover" - or "corine", fo which lookup tables are constructed to convert lulc classses to + Currently, ``lulc_fn`` can be set to the "vito", "globcover", "corine" or + "glmnco", of which lookup tables are constructed to convert lulc classses to model parameters based on literature. The data is remapped at its original - resolution and then resampled to the model resolution using the average - value, unless noted differently. + resolution and then resampled to the model resolution using the average value, + unless noted differently. The USLE C factor map can be refined for planted forests using the planted forest data source. The planted forest data source is a polygon layer with planted forest polygons and optionnally a column with the forest type to - identify orchards. The default value for orchards is 0.2188, the default - value for other planted forests is 0.0881. + identify orchards. The default value for orchards is 0.2188, the default value + for other planted forests is 0.0881. Adds model layers: * **landuse** map: Landuse class [-] Original source dependent LULC class, resampled using nearest neighbour. - * **Cov_river** map: vegetation coefficent reducing stream bank erosion [-]. * **Kext** map: Extinction coefficient in the canopy gap fraction equation [-] * **Sl** map: Specific leaf storage [mm] * **Swood** map: Fraction of wood in the vegetation/plant [-] * **USLE_C** map: Cover management factor from the USLE equation [-] * **PathFrac** map: The fraction of compacted or urban area per grid cell [-] - * **WaterFrac** map: The fraction of open water per grid cell [-] - * **N** map: Manning Roughness [-] Parameters ---------- @@ -363,10 +357,13 @@ def setup_lulcmaps( * Optional variable: ["forest_type"] - lulc_vars : list - List of landuse parameters to keep. - By default: \ -["landuse","Cov_river","Kext","N","PathFrac","USLE_C","Sl","Swood","WaterFrac"] + lulc_vars : dict + Dictionary of landuse parameters in ``lulc_mapping_fn`` columns to prepare + and their internal wflow name (or None to skip adding to the toml). By + default: \ +{"landuse": None, "Kext": "input.vertical.kext", "PathFrac": "input.vertical.pathfrac", + "Sl": "input.vertical.specific_leaf", "Swood": "input.vertical.storage_wood", + "USLE_C": "input.vertical.usleC"} planted_forest_c : float, optional Value of USLE C factor for planted forest, by default 0.0881. orchard_name : str, optional @@ -381,7 +378,7 @@ def setup_lulcmaps( ) # If available, improve USLE C map with planted forest data - if "USLE_C" in lulc_vars and planted_forest_fn is not None: + if "USLE_C" in list(lulc_vars.keys()) and planted_forest_fn is not None: # Add a USLE_C column with default value self.logger.info( "Correcting USLE_C with planted forest and orchards" diff --git a/hydromt_wflow/workflows/landuse.py b/hydromt_wflow/workflows/landuse.py index 70acc00b..3cbaf905 100644 --- a/hydromt_wflow/workflows/landuse.py +++ b/hydromt_wflow/workflows/landuse.py @@ -11,8 +11,8 @@ __all__ = ["landuse", "lai"] -RESAMPLING = {"landuse": "nearest", "lai": "average"} -DTYPES = {"landuse": np.int16} +RESAMPLING = {"landuse": "nearest", "lai": "average", "alpha_h1": "mode"} +DTYPES = {"landuse": np.int16, "alpha_h1": np.int16} # def landuse(da, ds_like, data_catalog, fn_map, logger=logger, params=None): diff --git a/tests/conftest.py b/tests/conftest.py index 0158c0ef..5001735b 100644 --- a/tests/conftest.py +++ b/tests/conftest.py @@ -25,7 +25,12 @@ def example_wflow_model(): logger = logging.getLogger(__name__) root = join(EXAMPLEDIR, "wflow_piave_subbasin") - mod = WflowModel(root=root, mode="r", data_libs="artifact_data", logger=logger) + mod = WflowModel( + root=root, + mode="r", + data_libs=["artifact_data"], + logger=logger, + ) return mod @@ -34,7 +39,10 @@ def example_sediment_model(): logger = logging.getLogger(__name__) root = join(EXAMPLEDIR, "wflow_sediment_piave_subbasin") mod = WflowSedimentModel( - root=root, mode="r", data_libs="artifact_data", logger=logger + root=root, + mode="r", + data_libs=["artifact_data"], + logger=logger, ) return mod @@ -76,7 +84,11 @@ def example_wflow_results(): @pytest.fixture() def clipped_wflow_model(): root = join(EXAMPLEDIR, "wflow_piave_clip") - mod = WflowModel(root=root, mode="r", data_libs="artifact_data") + mod = WflowModel( + root=root, + mode="r", + data_libs=["artifact_data"], + ) return mod diff --git a/tests/test_model_methods.py b/tests/test_model_methods.py index af3f77e7..72fedbd1 100644 --- a/tests/test_model_methods.py +++ b/tests/test_model_methods.py @@ -775,7 +775,7 @@ def test_setup_lulc_sed(example_sediment_model, planted_forest_testdata): example_sediment_model.setup_lulcmaps( lulc_fn="globcover", planted_forest_fn=planted_forest_testdata, - lulc_vars=["USLE_C"], + lulc_vars={"USLE_C": "input.vertical.usleC"}, planted_forest_c=0.0881, orchard_name="Orchard", orchard_c=0.2188,