STL2STM2BookC_StringPrint_SetupTPrintStartupSound_TabLISTBook_PrefsAct List Scene Drawing_Model Draw_Specs diagram_struct_arrayDependentsScriptPad_ListManuscriptPlayPlay_PrefsModelSubscript_Set_ListSubscript_SetArraySimulation_ModelRun_Specsentity_struct_arrayentity_structint16_arraydouble_arraytoken_type_arraytoken_type run_handle_array!Draw_Index"Poster_Info#run_handle$diagram_struct%Point_array&movie_attach'Picture_Tab(Pad)Model_Index_array*graph_type+curve_type_array,curve_type-Model_Index.AliasRecord/Panel0Sector_Specs1SET2Sense_Setup3select_struct_array4select_struct5allow_struct_array6allow_struct7Picture_array8Picture9Section_List:Window_List;Model_Window<Pad_Window=<<   D @ txxHH(FG(HH(d'`>>  ,, \ |     h    JJ  drmd   N<< T @ H PP??2x x VV^ L < Z H 0  t l p 8  ( &1 Dim_Name_1  VForest Open_Water Converted Low_Dens_Res Med_Dens_Res AgricultureHabitatb 2SLSoilb jj?@$?hTime ~*cn@N@5R@(\)??@5R@(\)@5R@5R???@+GzH@"\(?WindDayJul p p?@@@@@@@ @"@$@&@(@*@,@.@0@1@2@3@4@5@6@7@8@9@:@;@<@=@>@?@@@@@A@A@B@B@C@C@D@D@E@E@F@F@G@G@H@H@I@I@J@J@K@K@L@L@M@M@N@N@O@O@P@P@@P@P@Q@Q@@Q@Q@R@R@@R@R@S@S@@S@S@T@T@@T@T@U@U@@U@U@V@V@@V@V@W@W@@W@W@X@X@@X@X@Y@Y@@Y@Y@Z@Z@@Z@Z@[@[@@[@[@\@\@@\@\@]@]@@]@]@^@^@@^@^@_@_@@_@_@`@` @`@@``@`@`@`@`@a@a @a@@a`@a@a@a@a@b@b @b@@b`@b@b@b@b@c@c @c@@c`@c@c@c@c@d@d @d@@d`@d@d@d@d@e@e @e@@e`@e@e@e@e@f@f @f@@f`@f@f@f@f@g@g @g@@g`@g@g@g@g@h@h @h@@h`@h@h@h@h@i@i @i@@i`@i@i@i@i@j@j @j@@j`@j@j@j@j@k@k @k@@k`@k@k@k@k@l@l @l@@l`@l@l@l@l@m@m @m@@m`@m@m@m@m@n@n @n@@n`@n@n@n@n@o@o @o@@o`@o@o@o@o@p@p@p @p0@p@@pP@p`@pp@p@p@p@p@p@p@p@p@q@q@q @q0@q@@qP@q`@qp@q@q@q@q@q@q@q@q@r@r@r @r0@r@@rP@r`@rp@r@r@r@r@r@r@r@r@s@s@s @s0@s@@sP@s`@sp@s@s@s@s@s@s@s@s@t@t@t @t0@t@@tP@t`@tp@t@t@t@t@t@t@t@t@u@u@u @u0@u@@uP@u`@up@u@u@u@u@u@u@u@u@v@v@v @v0@v@@vP@v`@vp@v@v@v@v@v@v@#fffff@"\(@%@(=p =@!G{@(\)@'B\(@0zG{@.8Q@.p =q@5R@.zG{@2zG@*kQ@"8Q@&Gz@5(\@6R@05\(@*GzH@4@3J=p @00 =q@0(\@-=p =@(33333@+R@*(\)@$kQ@Gz@9\(@1=p @$ =p@(@#zG{@!p =q@ (\@'Q@0R@,p =q@*333333@1G{@0R@%\(@1Q@<(\@;nzG@-p =@!333333@&L@4# =p@-QR@4Y@2u\(@"G{@'kQ@.Q@0Gz@+GzH@-ǮzH@4J=p @1&fffff@4\(@10 =q@)@4J=p @2zG{@*p =@1J=p @6R@3Q@0W =p@1fffff@1(\@-@*=p =@( =p@7p =q@3p =@+fffff@* =p @* =q@-# =p@2\)@,\(\@(zG@8\(@-zG@5R@/\(@-p =@..zG@%8Q@*Q@)=p =@*@-(\)@)p =@0R@5 =p@5QR@.GzH@&@0zH@'\(@.R@-ǮzH@1\(@3333333@1 =p@4}p =@6R@0 =p@,333333@)Q@&\(@"zG@*\(\@+@#Gz@)p =@1\(@2\)@*(\)@)zG@)=p @+ =q@*333333@$Q@&W =p@(=p =@(\)@*fffff@-@)(\@( =p @)ǮzH@1Q@,GzH@&kQ@(\(@(=p =@$# =p@)zG@*\(@)ǮzH@.QR@.QR@,(\@)Gz@-GzH@,Q@' =p@$@3Q@1Q@&Gz@%Q@,ǮzH@-Q@-333333@/fffff@.(\@1QR@*GzH@&.zG@.# =p@+zG@2!Gz@)GzH@!33333@$\(@(aGz@,R@'zG@*fffff@&W =p@&\(@)R@.\(@/fffff@.p =q@(B\(@*QR@/zG@'Gz@&Q@*(\@.=p =@0&fffff@)zG@(aGz@' =p @.=p =@/Q@*@$zG{@,(\)@*ǮzH@)ǮzH@' =p @.\(@+Q@0Gz@.\(\@%# =p@#Q@(u\(@(@'zG@ G{@ zG@(B\(@*@%\(@+ffffff@1L@,\(\@!p =q@"Q@)(\)@+ffffff@*333333@&\(@"ffffff@'L@&\(@)(\@-333333@.Q@/ =p@'=p @)QR@% =p@$aGz@!p =@.=p =@+p =q@$zG{@%\(@"(\@$Q@/# =p@2G{@((\)@(W =p@+Q@'=p @#zG@!fffff@ ǮzH@$W =p@*=p =@,\(\@%\(@Gz@%B\(@&L@(=p =@-p =@'zG@+Q@-GzH@) =p @$# =p@$\(@%\(@%\(@+=p =@& =p@## =p@&aGz@+zG@'@(@"fffff@!\(@'.zG@0}p =@,Q@&@# =q@$@" =p@"G{@)=p @'@%W =p@.ǮzH@4GzH@-p =q@*(\@(@$L@=p =@(\)@@@QR@@,Q@* =q@- =q@0fffff@&@)W =p@)R@0\(@)W =p@!(\@=p =@1 @4޸Q@5h\)@3EQ@,Q@)R@"\(@@,GzH@,=p =@$\(@zG@-Q@2(\@-zG@)zG@1 =p@3xQ@/ffffff@(Q@ Q@'W =p@2aGz@(ǮzH@,Q@,33333@9Y@3s33333@0G{@(zG@&u\(@'Q@+@!G{@\(@,=p =@5xQ@7zG{@5L@)ǮzH@3@@3Q@'Q@1QR@.\(\@" =q@-ffffff@0QR@%Q@%@@%Q@2EQ@+GzH@+GzH !wind in km/hour" #,c @9@9??@9@9@9@9???@9@9?8C_zig25: !"Zigzaginess" of the stream network. Should be made spatially dependent when data available. The larger the actual length of the stream network within a cell relative to a straight line, the more will seep into it." #,#cT&?E??>XC??(fwi?ن?얳5L&d?E??>XC?>XC?>XC(????/7|/??8 inf_condUIF SW_to_SatW_downflow*DT=Unsat_cap THEN 0 ELSE max(Unsat_cap - Infiltration*DT, 0)2?%lP 2          ?     !)Condition for regular infiltration: > 0 - when the infiltration is according to the inf. parameter; 0 - when there is more water to be infiltrated than the available storage in UNSAT zone. In this case the UNSAT zone will become SAT and all the excess water will be channelled to SURFACE WATER " #,c m@0r@u7h S@tS!>??@u7h S@tS!>@u7h S@u7h S???@tffffg@tS!>?  SolRadAtmDayJul h? @ @@ @ @@@ @ @@@@;;@ @ @@" ʌ@$@&xx@( @*I$I"@,;8@.QN@0 @ @@1 K@2 ʌ@3 \\@4@5N@6xx@7-,@8 @9Q@:I$I @;/+@<;6@=eFTeA@>QL@?\W@@ @ @@@@A K@ANN@B ʌ@B-@C \\@Cm@D@Djj@EN@E@Fxx@F/@G--@Gp @H@H;;@IQ@I#@JI$I)@J*2/@K/5@KW5sW;@L;A@L @ G@MeFTeM@MKS@NQY@NsW5s_@O\e@O'bv'k@P @ @@PKmm@P@PǼ{@Q K@QL!!@QNN@Q{Ǽ{@R ʌ@RL\@R-@R//@S \\@SM؝@Sm@S=@T@TN==@Tjj@TΗ~@UN@UN@U@UKK@Vxy@VO_@V/@V @W-- @WPZZ@Wp@Wд @@X@XQ@X;;@Xhh!@YQ$@YQ!'@Y*@Y!-@ZI$I0@ZRv'bv3@Z*26@Z-9@[/<@[S*2*?@[W5sWB@[ӄ8CE@\;H@\S=K@\ @ N@\8C8Q@]eFTeT@]TI$W@]KZ@]N]@^Q`@^UFTeFc@^sW5sf@^ՠZi@_\l@_U_o@_'bv'r@_TeFTu@` @ @@`+W5sW=@`Kmm@`k8C?@`@`;A@`Ǽ{@`=C@a K@a, @ E@aL!!@al8C8G@aNN@aeFTeI@a{Ǽ{@aI$K@b ʌ@b,KM@bL\@blNO@b-@bQQ@b//@bFTeFS@c \\@c-sW5sU@cM؝@cmZW@cm@c\Y@c=@c_[@d@d.'bv']@dN==@dnTeFT_@djj@dha@dΗ~@djc@eN@e.me@eN@eopg@e@e5sW5i@eKK@ebv'bk@fxx@f/xm@fO_@fo{Ǽo@f/@f~q@f@fhs@g-,@g0C8Cu@gPZY@gpppw@gp@g؝y@gд @@gʌ{@h@h0x}@hQ @hq$I$@h;;@hQQ@hhh@h~~@iQ@i1@iQ!@iq؝؇@i@iZ@i! @i2*2@jI$I @j2__@jRv'bv@jrʌ@j*2@j@j-@jj@k/@k3;@kS*2*@ks@ @@kW5sW@kmm@kӄ8C@k󚹫@l;@l3Ǽ{ǝ@lS=@lsK@l @ @l!!@l8C8"@lNN@meFTe$@m4{Ǽ{@mTI$&@mtʌ@mK(@m\թ@mN*@m-@nQ,@n5//@nUFTeF.@nu\\@nsW5s0@n؝@nՠZ2@nm@o\4@o5=@oU_6@ov@o'bv'8@o==@oTeFT:@ojj@p @ @@pKK@p+W5sW @p;bv'ba@pKmm@p[xx@pk8C$@p{xe@p@p_@p;(@p{Ǽi@pǼ{Ǫ@p/@p=,@p~m@q K@q@q, @ 0@q<hq@qL!!@q\-,@ql8C84@q|C8Cu@qNN@qZY@qeFTe8@qppy@q{Ǽ{@q܇p@qI$<@q؝}@r ʌ@r @@r,K@@r<ʌʁ@rL\@r\@rlND@r|x@r-@r@rQH@r$I$@r//@r;; @rFTeFL@rQQ@s \\@shh@s-sW5sP@s=~~@sM؝@s]Q@smZT@s}@sm@s!@s\X@s؝ؙ@s=@s@s_\@sZ@t@t!@t.'bv'`@t>2*2@tN==@t^I$I#@tnTeFTd@t~__@tjj@tv'bv'@thh@tʌ@tΗ~@tޣ*2+@tjl@t@uN@u-/@u.mp@u>j@uN@u^/3@uopt@u;@u@u*2*7@u5sW5x@u@ @@uKK@uW5sW;@ubv'b|@umm@vxx@v8C?@v/x@v?@vO_@v_;C@vo{Ǽ@vǼ{@v/@v=G@v~@vK@v  h@t9@tS33333@tnfffff@t33333@t@t33333@tfffff@u33333@u9@uc33333@ufffff@u33333@u陙@v33333@vL@v@v@vfffff@w(@wc33333@w@w@x @xa@x@x陙@y0@yx@y@z @zX@zfffff@z@{D@{fffff@{陙@|<@|@|fffff@}>fffff@}fffff@}fffff@~H@~33333@~fffff@[33333@@ 33333@:fffff@i@@fffff@fffff@+33333@\@@fffff@@!@T@@@陙@@Nfffff@@fffff@@33333@I@{33333@@33333@@C33333@t@@33333@@8@i@@ə@@(@X@33333@@@fffff@@@l@@fffff@fffff@@H@t@fffff@@񙙙@fffff@C33333@k33333@fffff@@@fffff@,@P@t@@@33333@@#33333@D@d@@@@fffff@@@9@U@p@@fffff@@@@@ @733333@L@bfffff@w33333@33333@33333@fffff@@fffff@@@@@&fffff@4@Bfffff@O33333@\@h@s33333@}@33333@@@fffff@@@33333@@fffff@fffff@fffff@fffff@@33333@@ՙ@fffff@fffff@ՙ@@fffff@@@@@33333@fffff@@@@@@@fffff@x@m@a@U@H@<@-@33333@@@@33333@͙@33333@@@@l@X@Bfffff@,@@@@@33333@@~fffff@c33333@H@+33333@fffff@@33333@@@u@T@4@fffff@@͙@33333@33333@c33333@>fffff@@33333@@fffff@33333@W33333@/33333@fffff@@33333@@^fffff@333333@@@@33333@Vfffff@(@33333@@fffff@p@@@@ᙙ@@@P@ @33333@fffff@@\@*fffff@@fffff@@bfffff@0@fffff@@@g33333@4@33333@@fffff@l@;33333@@33333@fffff@t@D@33333@fffff@@@Rfffff@"fffff@fffff@@,@~љ@~vfffff@~33333@}@}i@}33333@|@|ffffff@|33333@{@{nfffff@{@z@z@z333333@yfffff@y@yT@y @x@x33333@x@@x@w@w@wFfffff@w 33333@v33333@v@vffffff@v333333@v@uљ@u@ux@uNfffff@u&fffff@u@t@t@t@t|@t`@tFfffff@t.fffff@t@t@s33333@s33333@sfffff@s33333@s@s33333@sfffff@s33333@s33333@s@s@sfffff@sfffff@s@s@s33333@s33333@tfffff@tfffff !" #,c #????????.??????8 C_Inf_Hab1 !/Infiltration rate modified by the habitat type" #,c ?=p =?=p =??333333?333333?333333?333333?=p =?=p =?=p =?=p =?????=p =?=p =?8 C_porosity0.41GJKLV   ?=p =!Porosity of the sediments. Field capacity = porosity - specific yield; ensure that alterations to porosity and specific yield are consistent in your parameterization. Units = m pore space / m sediment = proportion pore space" #,c F ? O??? O? O? O???4:QkU?8 InterceptGmin(Total_Precip, max(C_Intercep*Total_Precip,C_IntercepVeg*MacLAI))H P,  H    P     !Precipitation interception by vegetation. (m/d) If there is no vegetation a small proportion of precipitation may still get intercepted by branches and other structure. Otherwise the amount intercepted is a function of LAI. " #,c ?333333?333333???333333?333333?333333?333333????333333?333333?8 C_SatW_Out0.3=   ?33 3333!" #,c n@@+@mA&m@_1>??@mA&m@_1>@mA&m@mA&m???@deE@k ? SolRadGrdDayJul? 8ϛ@ 8ϛ@7h@ 8ϛ@@7h@kOA;@ 8ϛ@" `@$@&nu @(7h @* BL@?l@b 8@b!n@b&I@b=@c Sv@c-i'F@cMY|@cme@c=@cî@c CT@c@d@d.u @dN3,@dnJX_b@d`@dw;@d΍=@d:@e{p@e.@eNs@enY@eVH@e*ȗ~@eA:6@eW@fnu @f/V@fO@forR@f@fU.@f0d@f 8Ϛ@g!n@g08@gPN<@gpdLr@g{p@g@gШT*@gJ@h7h@h0@hQ@hqF"@h.Y@hEo@h[#@hrR@ib1@i16g@iQ@iq?@i @i~? @inu @i%߼ @jl@r|w;@r 8@r@r!n@rZz @r&I@r1?@r=@rH=u@s Sv@s^W@s-i'F@s=u @sMY|@s]@sme@s}5M@s=@suԃ@sî@ss@s CT@sY@s@sʲ%@t@t(@tN3,@t^?@tnJX_b@t~U.@t`@tl3@tw;@ttmi@t΍=@tޘ @t:@tW@u{p@uK @u.@u>:A@uNs@u^w@unY@u (@uVH@u@u*ȗ~@u6g@uA:6@uLsO@uW@ub䥅@vnu @vyVD@v/W@v?@vO@v_9(@vorR@v"^@v@v@vU/@v` @v0e @v @k@kΦ@ks@kߋr@lDoڃ@d @@`P1[@_ğס@b&@lŽ@m>@_w"sO@`z#@f2@m+^S@bzfw-@b!@fntf@o^@oKfU7@im@ku/k@p4??@p@1@pj<+4@p&k@p<1v@p @q$h@qJQ֌@jkN@l޹^[@qݮ@@r޽@rDNΚ,@ry~oq@q /@ij=-@ln?) @sVMN@sc{@sʶIӈ@t @tAwBV@ijeh@ku@t$0+@@u9D@l@n|t@nƒB@r+A@v}3Ko@vͧp@wk]F,@wF,@wӡ@w?#̍@x;Y@xY@qYSP -@whI^@n>@o 9m @tc*A*@qt w@sC@z] l@zc@{9@{g'hY@{%1@q6맹 @q.ک.b@uf@|:o@qEc5@r q@}h@}/@~54@q L@qUIwB@s^|@Ps}@u8|t@xi^M@yp8*@s,'7@s<@x p3@@7G@t62@DŗNf@y9Lj<@h8@`@$-@Y J@%u@3@#x @&r@v?@w& r@w֡a@|mgT@L@`@9@G@yɸ]tO]@}}9Nw@9恩@|@ @5*8@|R @zyyQ@|O @Ȃcz@~yL@Cä/@/g@E@P@{:?Gs\@| @uKTd@qj@~jֻ@;ՙ@'@g2@Vd@ 0 @إH$@-9 @z1@{D6@Q @l`4@z@F8@xJA @2]J@3hh@] @5a @K kn@g@[[pi@/"Y@bk@ד|@ ާ'@jP @RV@#YBu@)mqb@0'RTa@5}A5@:ҍ@?{=K@C:T@Gy,P@JŁI@@M@OBu@|8F]@~UnŰ@TK~@TC@T>t@Gg*@hK0@R<@Pp?>@N.#(@KuC|@HGTm?@DXK@@@;ud@~Q@6:8?@+Y%c@%#@IDg8@/VC@~h@ QX@LH@)}@HYT@ᅟ@&^B@|@>R@*ER@zЀt@=O@ @>Iw@re@~=׏h/@z"H8@|E4@ nTq@?@Ү@!h@T0@ߝ4@3@D@Gb@Ϸb@~5Tu@iV@T.r@yAZu@zX)@@4@~@wj@{CB@vXB@x@8OҦ+@~U)^@~_ @}]@rL46@r*7@tN?@|#XK@|g/|$=@| 12@{.z1@{.\-@{K7lH'@{Pes!`@z@oq!Ks@p.@y5@y*@y]T?@y}H@x҅Q @x`u"@p! @t k@nRJ@n?{@j9@lfOP,@v N][%@vrA {@v1@uRzcs@u9ul@upl@u1<@tǣ @t,GG@ty]>'@t=H@t k@sd$К@s‚@sTfj$@s*@rMwǂ@r(f^@rymĤ@rDڝ@r5@mK !correction for cloudy days" #,c y????I????8  SW_from_PS0\ !JSurface water from point sources (say sewage treatment plants discharges)" #,c"?PbM?PbM???PbM?PbM?PbM?PbM(????PbM?PbM?8 C_IntercepVeg0.001   ?PbM !)Interception by leaves m/d/{dimless LAI}" #,c {????!????8  C_horiz_HC1C !"Modifier for calibration purposes" #,c [@A՘P@A?@F9F@D5X@H4Z@Ex@A՘P@A @Aaw@Aawx???@Bs]O @A?8 SatW_HeadSAT_WATER/C_porosityq(  q   ! $P@A@Ai4l@A~@An0@AmzUV@A՘P@AP+@A+P8@At A:@AawyThe actual height of ground water table (m pores=m water in saturated layer, therefore m water/m water/m space=m space)" #,c @A@A?@F@F@F@F@A@A@A@A???@D@A?8 Wat_tableIC_wt_table+Datum_to_MSL|1  | 1  !" #,c LK?hr Ĝ?hr Ĝ???hr Ĝ?hr Ĝ?hr Ĝ?hr Ĝ????hr Ĝ?hr Ĝ?8 Pot_InfiltrMC_Infilt*C_Inf_Hab*C_Inf_Soil/(if C_Inf_Slope>0 then C_Inf_Soil else 1)*a?< *   a      a     !" #,c j?⓸0???⓸0?⓸0?⓸00???4Hӂ ?8 Total_PrecipRainfall+Snowmelt      !1Total precipitation as rainfall and snowmelt m/d" #,c(@,33333@,33333?@*ffffff@*ffffff@*ffffff@*ffffff@,33333@,33333@,33333@,33333???@0fffffg@,33333?8ic_SatWWat_table*C_porosity=q     !" #,cJe ????????????????8ic_UWKif Elevation > Wat_table then (Elevation-Wat_table)*IC_UW_moist_pr else 0Os4         O    !" #,c  ?]e???]e?]e?]e3???4;^Z)?8 Net_PrecipTotal_Precip-Intercept?\     !4Amount of precipitation that reaches the ground m/d" #,cm@Y@V@I??@S=@A@V@I@Q@@Q@b????@O@P? HumidityDayJul h?@@@@@@@ @"@$@&@(@*@,@.@0@1@2@3@4@5@6@7@8@9@:@;@<@=@>@?@@@@@A@A@B@B@C@C@D@D@E@E@F@F@G@G@H@H@I@I@J@J@K@K@L@L@M@M@N@N@O@O@P@P@@P@P@Q@Q@@Q@Q@R@R@@R@R@S@S@@S@S@T@T@@T@T@U@U@@U@U@V@V@@V@V@W@W@@W@W@X@X@@X@X@Y@Y@@Y@Y@Z@Z@@Z@Z@[@[@@[@[@\@\@@\@\@]@]@@]@]@^@^@@^@^@_@_@@_@_@`@` @`@@``@`@`@`@`@a@a @a@@a`@a@a@a@a@b@b @b@@b`@b@b@b@b@c@c @c@@c`@c@c@c@c@d@d @d@@d`@d@d@d@d@e@e @e@@e`@e@e@e@e@f@f @f@@f`@f@f@f@f@g@g @g@@g`@g@g@g@g@h@h @h@@h`@h@h@h@h@i@i @i@@i`@i@i@i@i@j@j @j@@j`@j@j@j@j@k@k @k@@k`@k@k@k@k@l@l @l@@l`@l@l@l@l@m@m @m@@m`@m@m@m@m@n@n @n@@n`@n@n@n@n@o@o @o@@o`@o@o@o@o@p@p@p @p0@p@@pP@p`@pp@p@p@p@p@p@p@p@p@q@q@q @q0@q@@qP@q`@qp@q@q@q@q@q@q@q@q@r@r@r @r0@r@@rP@r`@rp@r@r@r@r@r@r@r@r@s@s@s @s0@s@@sP@s`@sp@s@s@s@s@s@s@s@s@t@t@t @t0@t@@tP@t`@tp@t@t@t@t@t@t@t@t@u@u@u @u0@u@@uP@u`@up@u@u@u@u@u@u@u@u@v@v@v @v0@v@@vP@v`@vp@v@v@v@v@v@v h@P@Q@Q@@I@P@@V@V@S@@V@Q@@V@@W@N@Q@@R@@W@N@K@G@M@P@D@K@H@E@M@J@I@O@U@G@E@H@D@C@E@U@@W@R@N@Q@E@E@M@R@@F@E@G@V@W@@P@@J@D@D@O@K@O@I@J@H@S@@V@R@G@O@M@G@G@L@@@D@R@W@Q@J@I@U@L@I@N@V@@X@R@@K@P@S@C@R@Q@I@M@J@J@J@T@P@@H@I@P@H@A@@@T@V@W@M@L@R@R@U@U@O@L@Q@@J@O@R@T@@V@V@K@I@K@J@O@V@N@N@S@Q@O@N@O@Q@O@P@U@S@P@N@O@Q@@O@P@P@Q@R@P@M@N@M@Q@N@R@@M@K@O@L@H@H@G@L@P@A@H@K@M@P@W@U@R@@Q@S@U@P@M@N@P@@N@L@K@P@@U@T@T@V@O@Q@R@@H@P@M@Q@@R@M@H@K@L@L@P@O@K@O@N@J@T@@U@R@N@U@S@P@O@L@Q@O@F@G@P@R@@V@@P@@L@P@@Q@@T@T@Q@Q@R@S@T@R@P@Q@R@R@R@@T@@R@@Q@P@N@I@N@Q@S@@Q@Q@P@P@O@P@N@M@P@S@S@R@Q@U@@T@@M@L@O@S@S@V@@O@K@O@P@O@S@@R@S@R@T@O@I@I@L@S@@U@S@Q@R@@T@P@S@@N@P@@M@Q@R@U@U@V@@X@X@O@J@B@E@R@P@J@G@G@M@O@K@H@S@@T@R@R@Q@O@N@G@P@@Q@R@T@@X@@U@Q@@J@M@Q@@P@G@S@X@W@@X@L@K@I@Q@R@U@@T@@U@V@S@R@B@F@L@F@C@G@H@L@T@P@O@S@Q@S@V@P@O !crelative humidity data in % used Baltimore Airport, 1991 - Beltsville has no humidity measurements" #,c @@?@@@@@@@@????@@? Root_depth2#) !" #,c& ?⓸0??dQ?dQ?⓸0?⓸0?⓸0????8 Rainfall'IF Air_temp_C > 0 THEN Precip_m ELSE 0`$     `   ! Precipitation as rainfall (m/d)" #,c?6C-?6C-???6C-?6C-?6C-?6C-Y?????6C-?6C-?8 C_rechg0.00010   ?6 C-!Zrate of recharging of the aquifer below the base datum (loss from model system). units=m" #,f& ?hr Ĝ??dQ?6'~A)?hr Ĝ ?rN)y?hr Ĝ???4;^Z) 8UW_from_precip'if inf_cond>0 then Infiltration else 0?\s$     ?   !$P?z(?h|hh?XR@?rN)y?hr Ĝ" #,c @@?????@@@@???@@?8 C_Inf_Slope5 !Slopes are given as a map, with 1 where there is no slope and maximal inclination marked with 9. It is assumed that the larger the slope the less the infiltration rate." #,c@@??@@@@@@{@@?8C_horiz_hydr_cond5C !|The hydraulic conductivity of saturated soil, will be dependent on the soil type. Units of m/d. [NUBSAM] originally 10, 5" #,c 1?????? OutPutVar HabitatMap*06  6  !" #, fF?r$m??dڂ0#ؼ?z#&?r$m?r$m,???4Hӂ 8  SnowmeltG IF Air_temp_C>0 THEN MIN(C_SW_melt*Air_temp_C, SNOW_ICE/DT) ELSE 088D        8       !-Flux of precipitation from snowmelt {m/day}" #,g 8??????! $Hydrology!" #,c????????$??????"Soil1GC !%Determines the soil type in the cell" #,c # =>8ͤ =:V¬u?<4p(9Af?8ͤ =:V¬u==W5,/k==W5,/k????>S=:V¬u?8# UW_capil_exp$exp(-10*max(UW_depth-Root_depth,0))(),  (       !A small addition to account for the possible suction of water from saturated storage by capillary force when roots are close to the unsaturated depth. If the roots reach the sat layer there is no water taken from the unsat layer any more." #, f:?? ??? 8%SatW_from_precip;if inf_cond=0 then Unsat_cap/DT-SW_to_SatW_downflow else 02\q8        2    !$PIf there is more water to be infiltrated than the capacity of the unsaturated zone (inf_cond<0), the unsaturated zone is eliminated and all the unsat capacity becomes filled. The water in this case goes directly to the saturated storage. " #,c 3@[@[??@[@[@[@[???@[@[?& OutPutRow108 l!" #,c ???????????M???????8' C_UW_evap0.5N   ? !Nslow down of evaporation from ground relative to evap from open water dimless" #,c@@yS=?@#(@0+@9X`?i@@yS=@8dA@8dA%????@ K @?8( UW_depthMAX(Elevation-SatW_Head,0) 3#)     !&Depth (m) of the unsaturated zone. " #,-c s:?????? ?????? ??8) UW_trans_prtif Root_depth+0.10 then Water_avail*(UW_depth/(Root_depth+UW_capil_exp)) else 0(#9l  + ?  (       (     (  #      !$P??????????Fraction of transpiration that occurs from the unsaturated storage. 0.1 assumed to be the threshold depth after which the capillary action can be neglected. Once we hit the <0.1 zone an increasing amount of water is sucked from the saturated zone." #,c ?hr Ĝ?hr Ĝ???hr Ĝ?hr Ĝ?hr Ĝ?hr Ĝ????hr Ĝ?hr Ĝ?8* C_Infilt0.017   ?hr Ĝ!" #,c@B@B??@B@B @B@B???@B @B?, Elev_minElevation-C_surf_delta /K   /  !$P@B@B@B@B@B@B@B@B@B@BMinimal surface elevation in a cell. Used to calculate the possible depth of ravines, canyons, depressions that have streams in them. Ensures the proper exchange between sat.water and surface water in the stream network. Should be a map." #,c @@??@@@@R???@@?/C_surf_delta3., !SHow much lower the minimal elevation is in the cell relative to the mean elevation" #,f ?6C-?6C-??6C-?6C-?6C-?6C-%?6C-?6C-80SatW_recharge C_rechgq !&Movement to deep aquifer, units = m/d" #,c '@$@$??@$@$@$@$???@$@$?1Datum_to_MSL10 !" #,)fm?? lq????82SW_to_SatW_downflownif (HabitatMap = 1 or (UW_depth_before = 0 and Unsat_cap>0)) then MIN (Unsat_cap/DT, SURFACE_WATER/DT) else 06t%tqh 6              t       !If there was no unsat layer before, but there is one now, we assume that it appeared becasue of transpiration and horizontal fluxes. Then all the unsat storage will be filled and the saturated water will reach surface again if there is enough surface water. Assumes that vertical fluxes from surface water to saturated water are at least as fast as horizontal groundwater flow plus transpiration losses (m/day) " #,c?q9NZN%???q9NZN%BtkBtkd???UU?83 UW_deltaUW_depth_before-UW_depth(4l   (  !eChange in Unsat Water depth. Delta > 0 when water table rises, Delta < 0 when water table drops. (m)" #,f/ ?.%???.%?,Ƒ?"tk8???>UU84SatW_UW0max(0,min(-UW_delta*C_field_cap,SAT_WATER/DT))3qqs4 3    q      !9The reverse flow from saturated water to Unsat storage due to the drop of the water table, leaving moisture of field capacity. If Unsat_depth now is larger than Unsat_depth a dt ago, then the Sat water head went down, therefore the water at field capacity in the (Unsat_now-Unsat_before) entered the unsat layer." #,c  @@??@@@@???@@?5 cell_size40000C   @ !area of the unit considered m2" #,c  @@????@@@@E????@@?6 HabitatMap2ED2 !FHabitat numbers: 3 = corn 2 = forest 1 = open water other not defined" #,f/ ?xV}ݦ??zvȴ:?zvȴ:?xV}ݦ287 Snowfall0 IF Air_temp_C 0 THEN Precip_m ELSE 0 {m/day}`8$     `   !3When below zero precipitation is in form of snow. " #,s?xV}ݦ??ԁs?)_?xV}ݦ  ?88 SNOW_ICE 0 7;  !$P?DcӋN?xV}ݦ?r$mAmount of snow {m}" #,f "?$#>;?;f#Fz)7j >UXR?$#>; ?n'S}?$#>ۑ@ ?hp+089 UW_etransp#Transpiration *UW_trans_pr+Evap_UW)Ns  )  N  !$P?hp+0?Hy?;>6g?8ω>f>;>x>L>e,?n'S}?$#Transpiration loss is partitioned between the saturated and unsaturated zones; the loss from the saturated can occur even when roots do not quite reach the water table due to capillary action. Units = dimless *m/d * dimless function = m/d " #,+f? U???? 8:SatW_to_SW_upflowIF (SAT_WATER > SatWMax) THEN SAT_WATER-SatWMax ELSE IF ( SAT_WATER > SatWRunoff ) THEN C_zig*(SAT_WATER - SatWRunoff )*H_Conduct ELSE 0 qJKCqtd q J    q J    q K     q K   C    !$PVWhen the total storage capacity of the soil is exceeded by the water in the saturated tank, that overcapacity is routed to the surface water. Even before that when saturated water head exceeds the minimal surface elevation that represents the elevation in streams, creeks, a part of saturated water starts to seep into the surface storage. " #,f G? 0 ) THEN 0 ELSE -Air_temp_C*SURFACE_WATER*C_SW_melttt88         t    !Freezing of surface water and formation of ice {m/day} Since C_SW_melt is already in m/d, we assume that there is another parameter involved (1/m), but we think it is = 1" #,c ??M?????8< IC_SW_depth0t !NInitial condition of surface water height (m) above the sediment/soil surface" #,f ?[t??=p =?=p =?$/?$/?[t ?V] ?U pa 8= SatW_outflux C_SatW_Out*(SAT_WATER-ic_SatW)qq  q    !$P?:3?NQ 3?VPH?Z_uf?[t?Zv?XCf?V] ?U pbSaturated water horizontal movement. We assume that the Sat water tends to the initial condition." #,c ??????????????8> C_SW_out1.h !" #,c ?hr Ĝ??dQ?dQ?hr Ĝ?hr Ĝ?hr Ĝz????4;^Z)?8? Infiltrationmin(Pot_Infiltr, Net_Precip)    !{Potential flow of water from precipitation into unsat layer. It is a function of soils, the habitat type and slope. (m/d) " #,c??????????????@On_Map1 !" #,c??????????????AHYDRO1 !" #,c@?????s ]8?5>?s ]8?s ]8(????߭43?JDK?8LUW_mp_rUW_moist_prp/C_porosityN     !)Relative unsat water moisture proportion" #,c0( ?$#>;???$#>;?$#?$#???>ۑ@?hp+0?8NEvap_UW1C_UW_evap*Evaporation*(1-MIN(1,MacLAI))*UW_mp_r'PL90 '     P     L  !When MAcLAI is less than 1 there is always a certain fraction of land that is not covered by vegetation and can evaporate water from pores. m/d" #,c?У =p?У =p??????У =p?У =p?У =p?У =p[?????У =p?У =p?8OIC_UW_moist_pr0.26   ?У =p!\initial condition of the moisture proportion in the unsaturated zone (0-1 dimensionless) .5" #,cn@??????PMacLAIDayJulN?@@@@@@@ @"@$@&@(@*@,@.@0@1@2@3@4@5@6@7@8@9@:@;@<@=@>@?@@@@@A@A@B@B@C@C@D@D@E@E@F@F@G@G@H@H@I@I@J@J@K@K@L@L@M@M@N@N@O@O@P@P@@P@P@Q@Q@@Q@Q@R@R@@R@R@S@S@@S@S@T@T@@T@T@U@U@@U@U@V@V@@V@V@W@W@@W@W@X@X@@X@X@Y@Y@@Y@Y@Z@Z@@Z@Z@[@[@@[@[@\@\@@\@\@]@]@@]@]@^@^@@^@^@_@_@@_@_@`@` @`@@``@`@`@`@`@a@a @a@@a`@a@a@a@a@b@b @b@@b`@b@b@b@b@c@c @c@@c`@c@c@c@c@d@d @d@@d`@d@d@d@d@e@e @e@@e`@e@e@e@e@f@f @f@@f`@f@f@f@f@g@g @g@@g`@g@g@g@g@h@h @h@@h`@h@h@h@h@i@i @i@@i`@i@i@i@i@j@j @j@@j`@j@j@j@j@k@k @k@@k`@k@k@k@k@l@l @l@@l`@l@l@l@l@m@m @m@@m`@m@m@m@m@n@n @n@@n`@n@n@n@n@o@o @o@@o`@o@o@o@o@p@p@p @p0@p@@pP@p`@pp@p@p@p@p@p@p@p@p@q@q@q @q0@q@@qP@q`@qp@q@q@q@q@q@q@q@q@r@r@r @r0@r@@rP@r`@rp@r@r@r@r@r@r@r@r@s@s@s @s0@s@@sP@s`@sp@s@s@s@s@s@s@s@s@t@t@t @t0@t@@tP@t`@tp@t@t@t@t@t@t@t@t@u@u@u @u0@u@@uP@u`@up@u@u@u@u@u@u@u@u@v@v@v @v0@v@@vP@v`@vp@v@v@v@v@v@v?Q?(\)?\(?\(? =p ? =p ?QR?ə?Gz?Gz?Gz?(\)?(\)?p =q?p =q?θQ??У =p?GzH?Q?ҏ\(? =p? =p?zG{?Q?\(?ffffff? =p ?QR?\(?ٙ?=p =?Gz?(\)??zG?\(\?QR?\(?ᙙ?=p =?333333?(\)??p =q?ffffff?\(\?QR?GzH?=p =?333333?zG{?\(? =p ?(\)?\(?\(?\(?\(\?(\)??\(?Q?zG??\(?Q?zG? =p?\(?Gz? =p@zG{@333333@Q@Q@@ffffff@GzH@=p =@GzH@zG@@=p =@\(@@ =p @ Q@ GzH@ Q@ zG@ Q@ zG@ =p =@ zG{@ =p@ @ =p @ 333333@ \(\@ Q@ @ \(@ Q@ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ =p @ Q@zG@\(\@=p =?\(?zG{ !" #,(c ]#?X8,U?Vdc??Htj~?EHO?Htj~?X8,U?Vdc?X8,U?X8,U?WS庞f?8V UW_perc_rate^2*C_vert_hydr_cond*C_porosity*UW_excess^0.4/ ((C_porosity-C_field_cap)^0.4+UW_excess^0.4)nld  n     ^ ?ٙ      ^ ?ٙ   ^ ?ٙ    !Potential rate of percolation from unsaturated to saturated zone. The unsat hydraulic conductivity varies with soil moisture. The graph used limits percolation to water available to move under gravity (=above field capacity). Influence of moisture proportion on the hydraulic conductivity rate. Dimensionless multiplier of the saturated hydraulic conductivity rate. source for graph Domenico and Schwartz, Phys and Chem Hydro, p.91. Effectively the graph stops percolation when UNSAT_WAT/unsat_depth < Field_capacity F_Capacity is the proportion in the total soil volume. Poros - is the proportion of pores. Therefore F_Cap<=moist_prop<=Poros. moist_prop can become less than F_Cap due to transpiration, but not due to percolation. Multiplied by porosity, because Vert_conductivity is assumed to be measured in m of elevation, not unsaturated water content. Max rate of perc * unsat reduction m/d * dimless = m/d " #,cn@ @$b@!zG@!zG??@!zG@!zG@!zG@!zG???@!ǮzH@!zG?ZMacHDayJul?@@@@@@@ @"@$@&@(@*@,@.@0@1@2@3@4@5@6@7@8@9@:@;@<@=@>@?@@@@@A@A@B@B@C@C@D@D@E@E@F@F@G@G@H@H@I@I@J@J@K@K@L@L@M@M@N@N@O@O@P@P@@P@P@Q@Q@@Q@Q@R@R@@R@R@S@S@@S@S@T@T@@T@T@U@U@@U@U@V@V@@V@V@W@W@@W@W@X@X@@X@X@Y@Y@@Y@Y@Z@Z@@Z@Z@[@[@@[@[@\@\@@\@\@]@]@@]@]@^@^@@^@^@_@_@@_@_@`@` @`@@``@`@`@`@`@a@a @a@@a`@a@a@a@a@b@b @b@@b`@b@b@b@b@c@c @c@@c`@c@c@c@c@d@d @d@@d`@d@d@d@d@e@e @e@@e`@e@e@e@e@f@f @f@@f`@f@f@f@f@g@g @g@@g`@g@g@g@g@h@h @h@@h`@h@h@h@h@i@i @i@@i`@i@i@i@i@j@j @j@@j`@j@j@j@j@k@k @k@@k`@k@k@k@k@l@l @l@@l`@l@l@l@l@m@m @m@@m`@m@m@m@m@n@n @n@@n`@n@n@n@n@o@o @o@@o`@o@o@o@o@p@p@p @p0@p@@pP@p`@pp@p@p@p@p@p@p@p@p@q@q@q @q0@q@@qP@q`@qp@q@q@q@q@q@q@q@q@r@r@r @r0@r@@rP@r`@rp@r@r@r@r@r@r@r@r@s@s@s @s0@s@@sP@s`@sp@s@s@s@s@s@s@s@s@t@t@t @t0@t@@tP@t`@tp@t@t@t@t@t@t@t@t@u@u@u @u0@u@@uP@u`@up@u@u@u@u@u@u@u@u@v@v@v @v0@v@@vP@v`@vp@v@v@v@v@v@v@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!zG@!@!@!@!@!@!@!@!@!@!@!@!@!@!@!@!@!Q@!Q@!Q@!Q@!Q@!Q@!Q@!Q@!Q@!Q@!# =p@!# =p@!# =p@!# =p@!# =p@!# =p@!# =p@!(\)@!(\)@!(\)@!(\)@!(\)@!.zG@!.zG@!.zG@!.zG@!.zG@!333333@!333333@!333333@!333333@!333333@!8Q@!8Q@!8Q@!8Q@!=p =@!=p =@!=p =@!=p =@!=p =@!B\(@!B\(@!B\(@!B\(@!GzH@!GzH@!GzH@!GzH@!GzH@!L@!L@!L@!QR@!QR@!QR@!QR@!QR@!W =p@!W =p@!W =p@!W =p@!W =p@!\(\@!\(\@!\(\@!\(\@!aGz@!aGz@!aGz@!ffffff@!ffffff@!ffffff@!ffffff@!kQ@!kQ@!kQ@!p =q@!p =q@!p =q@!u\(@!u\(@!u\(@!zG{@!zG{@!zG{@!@!@!@!Q@!Q@!Q@!=p @!=p @!=p @!\(@!\(@!\(@!zG@!zG@!zG@!@!@!@!Q@!Q@!Q@! =p@! =p@!\)@!\)@!\)@!zG@!zG@!33333@!33333@!33333@!Q@!Q@!Q@!p =@!p =@!p =@!\(@!\(@!\(@!ǮzH@!ǮzH@!ǮzH@!@!@!@!R@!R@!R@! =p@! =p@! =p@! =p@! =p@!(\@!(\@!(\@!(\@! =p@! =p@! =p@!@!@!@!@!@!@!@!@!@!@!@!@!@!@!@!@!@!@!@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH@!ǮzH !" #,f50?R^>??dQ?R^>?R^>g8\SW_in6Net_Precip-UW_from_precip-SatW_from_precip+SW_from_PS% t    %    !hRainfall and snowmelt go into the unsat layer. Only what's in excess is made available for runoff. m/d" #,cn?zG{{=S =*=mf??'Z3=Osy/=S =*=mf=*=mf=*=mf???=O+s=LB `?_ETpmDayJul?@@@@@@@ @"@$@&@(@*@,@.@0@1@2@3@4@5@6@7@8@9@:@;@<@=@>@?@@@@@A@A@B@B@C@C@D@D@E@E@F@F@G@G@H@H@I@I@J@J@K@K@L@L@M@M@N@N@O@O@P@P@@P@P@Q@Q@@Q@Q@R@R@@R@R@S@S@@S@S@T@T@@T@T@U@U@@U@U@V@V@@V@V@W@W@@W@W@X@X@@X@X@Y@Y@@Y@Y@Z@Z@@Z@Z@[@[@@[@[@\@\@@\@\@]@]@@]@]@^@^@@^@^@_@_@@_@_@`@` @`@@``@`@`@`@`@a@a @a@@a`@a@a@a@a@b@b @b@@b`@b@b@b@b@c@c @c@@c`@c@c@c@c@d@d @d@@d`@d@d@d@d@e@e @e@@e`@e@e@e@e@f@f @f@@f`@f@f@f@f@g@g @g@@g`@g@g@g@g@h@h @h@@h`@h@h@h@h@i@i @i@@i`@i@i@i@i@j@j @j@@j`@j@j@j@j@k@k @k@@k`@k@k@k@k@l@l @l@@l`@l@l@l@l@m@m @m@@m`@m@m@m@m@n@n @n@@n`@n@n@n@n@o@o @o@@o`@o@o@o@o@p@p@p @p0@p@@pP@p`@pp@p@p@p@p@p@p@p@p@q@q@q @q0@q@@qP@q`@qp@q@q@q@q@q@q@q@q@r@r@r @r0@r@@rP@r`@rp@r@r@r@r@r@r@r@r@s@s@s @s0@s@@sP@s`@sp@s@s@s@s@s@s@s@s@t@t@t @t0@t@@tP@t`@tp@t@t@t@t@t@t@t@t@u@u@u @u0@u@@uP@u`@up@u@u@u@u@u@u@u@u@v@v@v @v0@v@@vP@v`@vp@v@v@v@v@v@v=IwKf]=LB `=S =O~k#=Ac  =5LBLx ===;=7Rte%=Ed)5@==/=*=mf=I`N=H1[=Iyx=;ONP=QYMzr=T4#˙=T:=W@ =,:=V@r=SۆHϛ=RK}"=RQAV =U.I=Q^=QΒ=T=T.O|=N֣ =Ys=U;=U]=[4$=b('3S=eem=[H,=@IC4=OT=U$р =Qj%+t=Y.=VTnt=R~t=M+^=Z4IPܩ=T?վ!=TLkJ=BMM;&=: qL=Viv=_iRZ=a~߼&|=e$Ԏ=Tl B=W'IV=U,=TCA=Tyq*=[z!=V*A=R릇=SsԮ=[!1"'=Yy =\f Y=\Ab}=[C=U5Ck=]p=[zI=PS=7Aj@=E|GU=THd=ZS"$=PX=Zv+"=^d&M=`a$=U4w,> ? ݭ?%u\?%Q>? !oҥ?3$f?1V| ^?!?$@?'qDgR?+I0?-h[Y ?.J?,*?.tQ ?;@aWD?Cy5?BOe?F=Dy?D<@=?A%*?4Yzm_?W5?{ ?1*Z?Aķ@oL?Al?8՟c?-`?&z"/?3.GjG?;F(??<#?D" ?Hbr?Ha,#'*?CٵfR?7q?7&[A5?MA?SNq?PXt?P\?Pϴ;d ?E'_Ła?P֣ ,?QB?R7{f?NUpEQ?V?[Nw6~}?`7?`ڬ?cjX d?cfN?ZЀ?Uƿ?Z?[6X?`Ӟz^&?c-wu?i,x]?lH]?n>-g?pV?o?qz b ?tc?tU?x&]?t6=?rȖԏ0?pڙ,?utv?sBň?nn ?qV`e?uܮò?y:?~@F8?y?}yan?!r?}`?|\ZZ?UG??kXCy?abnن?qbP?z_.qd$?yO*?i&cI?mܦM?xW[?yav?yT㏺{D?~uEZ?3O?x1?E?vr.w?rL?s.e?l.]e?y;;?ߒ?z+*٦?X"r? zԡF?e~B?|wiK#?}g?4?(a+Ǎ?-?س?NKx?"?x%]R?sV'*3?Fs?1v?Mbb?"/?v{?lTb?w] ?zM#?rT~z@c?hs=?dv`?4?c]?f 耡??`,o}?X?H*?c[?Ws\*0?8IT$>Q?=dg#wu=X5 '=VK&fw=Uv T=XpX =X@k=VK̒!=U>=aΓ=dr=bɐ=eΓ=_RuO=` N=btpe=_f[=Vf=R+U{=Td~Q=X1Es=W=K/=C`ʭ5=HD=LUR=Q^$=[Hj(=bDET<=cYG:=Vu=Tjm=Xyha=Y_7|g2=EW T"=M3Ȃh=X$-;C=ZD,\@=RO!=J(=K頙t=_s=\K3=C- nR==B)}`=*ɉ x=S^-p=Q}=Qgڶ =O}f=P+y~[=I=M?(8=B[ =BV+ =P=\[=T᤟0/\=c/MV=W`e=P/o =Vuٽ$=U=TKM=Sn[l|=S,=J*}4=MxPW=K P9=E* >=GH$ʅ=Fח=97=L`4=NMR=O+r !Penman-Monteith equation for evapotranspiration (Shuttleworth, 1993 in Handbook of Hydrology) m = mm/1000 (delta*Energy+1.013*ro*D/rA)/(delta+gamma*(1+rS/rA))/lambda/1000" #,cm?QJ?⓸0???⓸0?⓸0?⓸0????` Precip_mDayJul7 h? @ @@ @ @@@ @ @@@@;;@ @ @@" ʌ@$@&xx@( @*I$I"@,;8@.QN@0 @ @@1 K@2 ʌ@3 \\@4@5N@6xx@7-,@8 @9Q@:I$I @;/+@<;6@=eFTeA@>QL@?\W@@ @ @@@@A K@ANN@B ʌ@B-@C \\@Cm@D@Djj@EN@E@Fxx@F/@G--@Gp @H@H;;@IQ@I#@JI$I)@J*2/@K/5@KW5sW;@L;A@L @ G@MeFTeM@MKS@NQY@NsW5s_@O\e@O'bv'k@P @ @@PKmm@P@PǼ{@Q K@QL!!@QNN@Q{Ǽ{@R ʌ@RL\@R-@R//@S \\@SM؝@Sm@S=@T@TN==@Tjj@TΗ~@UN@UN@U@UKK@Vxy@VO_@V/@V @W-- @WPZZ@Wp@Wд @@X@XQ@X;;@Xhh!@YQ$@YQ!'@Y*@Y!-@ZI$I0@ZRv'bv3@Z*26@Z-9@[/<@[S*2*?@[W5sWB@[ӄ8CE@\;H@\S=K@\ @ N@\8C8Q@]eFTeT@]TI$W@]KZ@]N]@^Q`@^UFTeFc@^sW5sf@^ՠZi@_\l@_U_o@_'bv'r@_TeFTu@` @ @@`+W5sW=@`Kmm@`k8C?@`@`;A@`Ǽ{@`=C@a K@a, @ E@aL!!@al8C8G@aNN@aeFTeI@a{Ǽ{@aI$K@b ʌ@b,KM@bL\@blNO@b-@bQQ@b//@bFTeFS@c \\@c-sW5sU@cM؝@cmZW@cm@c\Y@c=@c_[@d@d.'bv']@dN==@dnTeFT_@djj@dha@dΗ~@djc@eN@e.me@eN@eopg@e@e5sW5i@eKK@ebv'bk@fxx@f/xm@fO_@fo{Ǽo@f/@f~q@f@fhs@g-,@g0C8Cu@gPZY@gpppw@gp@g؝y@gд @@gʌ{@h@h0x}@hQ @hq$I$@h;;@hQQ@hhh@h~~@iQ@i1@iQ!@iq؝؇@i@iZ@i! @i2*2@jI$I @j2__@jRv'bv@jrʌ@j*2@j@j-@jj@k/@k3;@kS*2*@ks@ @@kW5sW@kmm@kӄ8C@k󚹫@l;@l3Ǽ{ǝ@lS=@lsK@l @ @l!!@l8C8"@lNN@meFTe$@m4{Ǽ{@mTI$&@mtʌ@mK(@m\թ@mN*@m-@nQ,@n5//@nUFTeF.@nu\\@nsW5s0@n؝@nՠZ2@nm@o\4@o5=@oU_6@ov@o'bv'8@o==@oTeFT:@ojj@p @ @@pKK@p+W5sW @p;bv'ba@pKmm@p[xx@pk8C$@p{xe@p@p_@p;(@p{Ǽi@pǼ{Ǫ@p/@p=,@p~m@q K@q@q, @ 0@q<hq@qL!!@q\-,@ql8C84@q|C8Cu@qNN@qZY@qeFTe8@qppy@q{Ǽ{@q܇p@qI$<@q؝}@r ʌ@r @@r,K@@r<ʌʁ@rL\@r\@rlND@r|x@r-@r@rQH@r$I$@r//@r;; @rFTeFL@rQQ@s \\@shh@s-sW5sP@s=~~@sM؝@s]Q@smZT@s}@sm@s!@s\X@s؝ؙ@s=@s@s_\@sZ@t@t!@t.'bv'`@t>2*2@tN==@t^I$I#@tnTeFTd@t~__@tjj@tv'bv'@thh@tʌ@tΗ~@tޣ*2+@tjl@t@uN@u-/@u.mp@u>j@uN@u^/3@uopt@u;@u@u*2*7@u5sW5x@u@ @@uKK@uW5sW;@ubv'b|@umm@vxx@v8C?@v/x@v?@vO_@v_;C@vo{Ǽ@vǼ{@v/@v=G@v~@vK@v  h?DѢ?wzp?ۭ:`N?wzp?a?j㣨?TѢ?p >cݥ?whb?]!x?@it?@it?DѢ?DѢ?8.‹?gzp?fqT2U?0it?o6&,s.?o6&,s.?VqT2U?dѢ?e:?X.‹? it? it?[ ]R?[ ]R?X:S?), ?X.‹?`it?`it?qL?xr<~?[ ]R?2HQi?zp? it? it? it?4m?O?#3? it?cmZ?cmZ?R? ?:0?TѢ?TѢ?}!x?L?smZ?k ]R?8.‹?DѢ?cmZ?aL?TѢ?TѢ?ju%F ?6&,s.?x.‹?_6&,s.?X.‹?TѢ?xr<~? >cݥ?_6&,s.?_6&,s.?_6&,s.?*%?V6mzV?R?0it?0it?q*%?q*%?'?it?R?o6&,s.?s?B5.?M!x?|M>?|M>?8.‹?H.‹?8.‹?qL?b ?p >cݥ?DѢ?8.‹?H.‹?H.‹?TѢ?TѢ?smZ?E$?i?gzp?_6&,s.?aL?r?cmZ?9Y ?), ?tIE?0it?|"[t?{ ]R?Pit?:0?&\?fqT2U?fqT2U?zrѮ?{ ?0it?@it?@it?qL?qL?dѢ?dѢ?8.‹?Y:l?a?tIE?8.‹?8.‹?it?it?`it?`it?q*%?{ ]R?mZ?Kc?0it?0it?0it?0it?8.‹?DѢ?fqT2U?k ]R?H.‹?}!x?it??uSk]F,??b ?@(?zrѮ?@it? >cݥ?6&,s.? it?h.‹?m!x?DѢ? it? it?fqT2U?fqT2U? it?&\?u%F ?X.‹ !Precipitation in m/day" #,c 8@@?@@@@@@@@????@@?8a C_Inf_Soil4.00 ! rate of infiltration for a given soil type, read in from soils map units m/day (assumes steady state infiltration-saturated conditions) original value = 3.6 [very sensitive about .01to .05, .02 yeilded a 11.6% sf runoff. .03 {m/day} New units = mm/hour, therefore *24/1000 10*24/1000" #,cl??????8b Out_F_SF0 !" #,f,= ??xbTwn??Q&?xbTwn??Q&&8dSW_evap-MIN (Max (Evaporation, 0), SURFACE_WATER/DT)tt(    t    !'Surface water evaporated. Units=m/d " #,f A??5ƽ?;W?zG?׈X8h SW_outfluxC_SW_out*SURFACE_WATER/dt>tt > t    !YAssuming that all the available surface water is removed to downstream cells. Units m/d " #,.fxI?dYx?WS庞f??!d?D|G?cH?dYx?WS庞f ?X@6ZV?X8,U ?WS庞f8lUW_SatW_exchangey{1}IF (inf_cond = 0) THEN UNSAT_WATER/DT ELSE MIN(UNSAT_WATER/DT, {2}MAX(0,UW_delta*UW_moist_prp)/DT+ {3} UW_perc_rate) s3Vsqp     s     s   3       V   !$P?WS庞f?b~hw:?dYx?c-"?_@Y}t?Y2rV?Xt0}qk?X@Ћ?X@6ZV?X8,UThe first term {1} describes the condition when all available unsaturated capacity is to be taken by the infiltration and as a result the entire Unsat Storage is to be eliminated and moved to the Sat. tank. The second term {2} is the water added to Sat_water from the Unsat storage when water level rises. The third {3} is the percolation flux . Varies not only with the level change but also with the soil moisture in unsaturated water storage. m/d+m/d = m/d " #,cL?|1&y?|1&y??tzG{?tzG{?tzG{?tzG{?|1&y?|1&y?|1&y?|1&y/?????|1&y?|1&y?8nC_vert_hydr_cond0.007V   ?| 1&y!0The vertical hydraulic conductivity of saturated soil, will be dependent on the soil type. This may be different from the horizontal conducitvity of the sediment or aquifer. Units of m/d. [NUBSAM] Data in database is 10 sand = 3; loamy sand = 0.75; loam = 0.09; clay loam = 0.025; clay = 0.0075 0.02" #,s Q@, t@,33333?@0l:K&@%'p@-rU8@*fK]c@, t@,33333 @,wI@,wI?Y O @,33333?8q SAT_WATERic_SatW=4:l2% 0:=4 !$P@,33333@,@,tF@,@,M$d@, t@,61@,z<߄@,]@,wIWater column in the saturated zone. This is the real amount of water if squeezed out of the soil. The water head (groundwater table) will be sat_water devided by porosity. Units =m " #,s T?lyt?&E?@$k?U'ڋR?/%SfS?lyt?&E ?lyt?lyt?Y{ ??8s UNSAT_WATERic_UWl4l9 !$P??VD5?gb??&E?U?}D?5X5Z? 3?lytWater column in the unsaturated zone between the top of the saturated groundwater table and the sediment surface. This is the actual water if squeezed out of the soil. Units= m *dimless = m" #,s U?tzG{<?$   ?8tSURFACE_WATER IC_SW_depth<dh;2\:d2;h @9?8| IC_wt_table25 !fHeight [+] in meters of the saturated (ground) water table above the base datum, use a GIS map input." #,c ,q ?;Kp83;>>|y38t;`˯`_ ???? ?8Transpiration-C_transp*MIN(1,MacLAI )* ETpm *Water_avail FP_9( F  P    _    !$PPotential transpiration by macrophytes, including effect of water stress. New ETpm - Penman-Monteith formula should be taking care of all. Units = dimless * m/d * dimless + dimless * m/d * dimless * (m^2/m^2) = m/d For given daylength and canopy saturation deficit, maximum rate of transpiration of water from plant thru stomata IF water is not limiting, stressing plants to any degree. Units = hr/day * (mol m-2 sec-1 * sec/hr * kg/mol * m^3/kg) * dimless = m water/day" #,c +r?Wz?砂?I??ʅ\9?!W^??9>K7ď  ???? 8 SatW_transpTranspiration *(1-UW_trans_pr))q   )   !$P Transpiration loss from saturated zone occurs exponentially slower as the unsaturated zone lowers below the root depth (and thus rapidly reducing the capillary draw of the root suction of water in saturated zone). Units = dimless* m/d * dimless function = m/d " #,cZ@@yS=??@@yS=@_0@_0g???@LP@?8UW_depth_beforedelay(UW_depth,DT)(32 (    !hThe Unsat depth on the previous time step. Used to define whether the water table is rising or falling." #,c * ?0'"k? F"3??xbTwn??Q&?xbTwn??Q&?0'"k? F"3 ?0'"k?0'"k????? ( ?"i?8 Evaporation+C_evap*SolRadGrd/585*pan_CW*pan_CT*pan_CH BdN, B   I        !$P?"i? PG?!]U?OŴ\?po 0) then UNSAT_WATER/UW_depth else 0(s Ll, (    s (    !$P?У =p?С9?Н Hj?И?В?СD?п" G?кsH?к?ͼproportion of water in unsaturated water storage. =0 when the UNSAT reservoir is empty. Grows to Porosity as UNSAT reservoir is filled up. Used to calculate the hyd.conductivity as function of saturation level." #,0.6*C_field_cap then 1 else (UW_moist_prp-0.1*C_field_cap)/0.5/C_field_cap))#)**  #  < ?      '   > ?33 3333       '   - ?    / ?      !$P??????????Availability of water to plant root system. If the root depth is >depth of unsaturated zone, the return value is 1.0 (dimensionless). If the unsat zone is below the root zone depth, the value returned is the moisture proportion PLUS an exponentially decreasing proportion from the saturated zone. Thus water may be fully available even though the unsaturated zone has 85% water available due to the addition of water drawn by capillary action from saturated. Otherwise the water availability is defined by the moisture proportion, the wilting point (assumed to be 0.1 of field capacity) and the drying capacity (assumed to be 0.6 of field capacity) [p.4.46 in Shuttleworth]" #,-cS?-?ݡo??Ƨ?\/n`?Ƨ?\/n`?-?ݡo?1[F?1[F~?????nO??8pan_CTTif Air_temp_C > -30 then 0.463+0.425*(Air_temp_C/20)+0.112*(Air_temp_C/20)^2 else 0x       : ?ݡ o + ?33 3333      + ? 1&y          !Christian 1968 (in Haan) evap formula: temperature component. Temperature is relative to standard 20 deg C used by Christian." #,c ????333333?333333?333333?333333????1???????8 C_field_cap0.254V   ? !2proportion of total soil volume, for a given soil type, that represents water not able to be drained by gravity. Field capacity = porosity - specific yield; ensure that alterations to porosity and specific yield are consistent in your parameterization. .11 from Achmad, 1991 and Rasmussen and Andreasson" #,&c1$?,|?33332??H6U?2^?H6U?2^?,|?33332?3?3y?????fl!-?qwAT?8pan_CH21.035+0.240*(Humidity/60)^2-0.275*(Humidity/60)^3\  ?\ ( + ?θQ   <      - ?љ  <      !zChristiansen 1968 (in Haan) evap formula: humidity component. humidity is relative to standard 60% used by Christiansen." #,(c)+@t'Z?L7mtgE??Ü?艪"i?Ü?艪"i@t'Z?L7mtgE@t'Z@t'Zz?????&,0?$uD?8pan_CW*0.672+0.406*(Wind/6.7)+0.078*(Wind/6.7)^2d  ? $/ + ? lC  / @    + ? +  / @      !{Christiansen 1968 (in Haan) evap formula: wind component. Wind speed is relative to standard 6.7 km/hr used by Christian." #,c w?tzG{?tzG{??zG{?zG{?zG{?zG{?tzG{?tzG{?tzG{?tzG{????tzG{?tzG{?8 C_SW_melt.005 ;   ?tz G{!Freezing/melting rate coefficient (m/day) Changed from the previous linear dependency on the SNOW_ICE stock - makes no sense to make the snowmelt proportional to the existing amount of snow." #,+cM@!b?@!b@!b@!b?????cloudyNif Precip_m>0 then max(0,10-1.155*(vap_press/(Precip_m*1000*30))^0.5) else 0 `p `    - ?z G{  `      ^ ?        !daily cloudiness (tenths) (Nikolov and Zeller) Note: when used in eqn, constrain negative values to =0 Precipitation in this equation is total monthly precipitation in mm. Equation is: if Precipitation>0 then 10-1.155*(vap_press/Precipitation)^0.5 else 0 Zeroed for now" #,"c 9@-aoPP@s=p ?@E((?>Jzf@E((?>Jzf@-aoPP@s=p @|r@|r.????@|؊Y@}!&? vap_press:Humidity*6.1078*EXP(17.269*Air_temp_C/(Air_temp_C+237.3))L  * @nc    ( @1D /    + @m      !/Daily vapor pressure (Pa) (Nikolov and Zeller)" #,c m@A@w;?@A)88@A)88@w;?d63E?d63E????@?? Air_temp_CDayJul ;7 h? @ @@ @ @@@ @ @@@@;;@ @ @@" ʌ@$@&xx@( @*I$I"@,;8@.QN@0 @ @@1 K@2 ʌ@3 \\@4@5N@6xx@7-,@8 @9Q@:I$I @;/+@<;6@=eFTeA@>QL@?\W@@ @ @@@@A K@ANN@B ʌ@B-@C \\@Cm@D@Djj@EN@E@Fxx@F/@G--@Gp @H@H;;@IQ@I#@JI$I)@J*2/@K/5@KW5sW;@L;A@L @ G@MeFTeM@MKS@NQY@NsW5s_@O\e@O'bv'k@P @ @@PKmm@P@PǼ{@Q K@QL!!@QNN@Q{Ǽ{@R ʌ@RL\@R-@R//@S \\@SM؝@Sm@S=@T@TN==@Tjj@TΗ~@UN@UN@U@UKK@Vxy@VO_@V/@V @W-- @WPZZ@Wp@Wд @@X@XQ@X;;@Xhh!@YQ$@YQ!'@Y*@Y!-@ZI$I0@ZRv'bv3@Z*26@Z-9@[/<@[S*2*?@[W5sWB@[ӄ8CE@\;H@\S=K@\ @ N@\8C8Q@]eFTeT@]TI$W@]KZ@]N]@^Q`@^UFTeFc@^sW5sf@^ՠZi@_\l@_U_o@_'bv'r@_TeFTu@` @ @@`+W5sW=@`Kmm@`k8C?@`@`;A@`Ǽ{@`=C@a K@a, @ E@aL!!@al8C8G@aNN@aeFTeI@a{Ǽ{@aI$K@b ʌ@b,KM@bL\@blNO@b-@bQQ@b//@bFTeFS@c \\@c-sW5sU@cM؝@cmZW@cm@c\Y@c=@c_[@d@d.'bv']@dN==@dnTeFT_@djj@dha@dΗ~@djc@eN@e.me@eN@eopg@e@e5sW5i@eKK@ebv'bk@fxx@f/xm@fO_@fo{Ǽo@f/@f~q@f@fhs@g-,@g0C8Cu@gPZY@gpppw@gp@g؝y@gд @@gʌ{@h@h0x}@hQ @hq$I$@h;;@hQQ@hhh@h~~@iQ@i1@iQ!@iq؝؇@i@iZ@i! @i2*2@jI$I @j2__@jRv'bv@jrʌ@j*2@j@j-@jj@k/@k3;@kS*2*@ks@ @@kW5sW@kmm@kӄ8C@k󚹫@l;@l3Ǽ{ǝ@lS=@lsK@l @ @l!!@l8C8"@lNN@meFTe$@m4{Ǽ{@mTI$&@mtʌ@mK(@m\թ@mN*@m-@nQ,@n5//@nUFTeF.@nu\\@nsW5s0@n؝@nՠZ2@nm@o\4@o5=@oU_6@ov@o'bv'8@o==@oTeFT:@ojj@p @ @@pKK@p+W5sW @p;bv'ba@pKmm@p[xx@pk8C$@p{xe@p@p_@p;(@p{Ǽi@pǼ{Ǫ@p/@p=,@p~m@q K@q@q, @ 0@q<hq@qL!!@q\-,@ql8C84@q|C8Cu@qNN@qZY@qeFTe8@qppy@q{Ǽ{@q܇p@qI$<@q؝}@r ʌ@r @@r,K@@r<ʌʁ@rL\@r\@rlND@r|x@r-@r@rQH@r$I$@r//@r;; @rFTeFL@rQQ@s \\@shh@s-sW5sP@s=~~@sM؝@s]Q@smZT@s}@sm@s!@s\X@s؝ؙ@s=@s@s_\@sZ@t@t!@t.'bv'`@t>2*2@tN==@t^I$I#@tnTeFTd@t~__@tjj@tv'bv'@thh@tʌ@tΗ~@tޣ*2+@tjl@t@uN@u-/@u.mp@u>j@uN@u^/3@uopt@u;@u@u*2*7@u5sW5x@u@ @@uKK@uW5sW;@ubv'b|@umm@vxx@v8C?@v/x@v?@vO_@v_;C@vo{Ǽ@vǼ{@v/@v=G@v~@vK@v  h?`A7L?S@?n??Ƨ-@lD?n@XbM@Ƨ-@ r ě?S@F-@cS??9XbM@ Vt@S@TE?5?|h?ʟvȴ9@ Vt@Ƨ-@"S@)@%B\(@!\1&@E@ Vt@TE?ڰ ě@ I^@?ڰ ě@9XbM@&Q@#ql@1&@#*~"?r ě@ S@1&?vȴ?5?|h@Vu@$@-p =q@*B\(@S@Vu@#*~"@S@ 9XbM@Ƨ-@1&@r ě@Vt??1&y@9XbM@cS@"xF@$ =p@F-@$@'zG{@"S@l@%aGz@ @@#F@,W =p@0QR@!~"@ Ƨ-@lD@1&@S@+ I^@&Q@(Q@-zG@2@5.zG@5@@.\(@%\)@$@ ~"@F-@"U$/@.GzH@2R@+ =p@#NVu@!\1&@ Vu@ ~"@&@*Q@-Q@0p =@4# =p@2(\@.# =p@0Gz@2(\@-p =q@)kQ@,\(@2=p @2xQ@-Q@1zG@/B\(@1nzG@2!Gz@6\(@8aGz@8@6ٙ@6@7\(@3zH@+\(\@. =p@3\(@6 =p@7\(@9@:zH@;&fffff@;nzG@:u\(@:zH@5\(@=p =q@=zG@;nzG@9kQ@: =p @5\(@70 =q@7@6R@7\(@9zG@# =p@:.zG@6 =p@8zH@7B\(@9R@95\(@6\(@6&fffff@7\(@8aGz@8zG@7=p @8\(@9kQ@:\(@9\(@6p =@4\(@5zH@7R@8p =q@6\(\@6@7R@:@@;ǮzH@=@:(\@2\(@1R@3fffff@7(\@7 =p@3fffff@5@@6 =p@8zH@8fffff@6ٙ@3O\(@4kQ@6p =@8fffff@;\(@<(\@: =q@4# =p@,333333@'Gz@-zG@1\(@3zG@1\(@-Q@'W =p@)kQ@0Q@0(\@3Q@6&fffff@5@@6J=p @4GzH@*aGz@"Vu@$kQ@(zG@.\(@)G{@(L@#ql@$GzH@-p =q@*@ ~"@(Q@,\(@1&@Vu@&Q@0(\@0\(@0Gz@0 =q@0 =q@)\(@cS@$ =p@'\(@+8Q@&zG@TE?Ƨ-@lD@ql@O;dZ@9XbM@O;dZ@@ql@ql@ + I@*@&\(\@ I^@S@$kQ@1\(@3\(@0@*Q@XbM??n@+ I^@)\(@1\(@,Q@"1&@XbM?`A7L@Vt@+ I^@F-@&Q@E@1&@ l@) =p@%B\(@S?vȴ@9XbM@ S@r ě?5?|h?n?@r ě@F-@@ !Temperature - deg C" #,c @D@D?@I@I@I@I@D@D @D@D????@D @D? ElevationDatum_to_MSL+301(J,  1   !$P@D@D@D@D@D@D@D@D@D@Dm" #,c  @$??@v?@v?@$?@$@$????@vݙ??DayJulmod(time-DT,365)+1Z_P`(    m ,    !julian day, 1 thru 365 " #, [ 8:  V:%%%%a t@8 ~@&'(~To reduce file size the length of simulation has been set to 10. Please restore it to 365 prior to run and analyze the model.)* K%t#@#M _&'(Unit Hydrologic Model)* K%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% P!88 !9&' $Hydrology()* Kdrmd J |aX<   E% c8@Yy\&'()* V% s}T8@P?F@kX&'()* V% cC@.0 v6k.&'()* V% abq8Qa@{\Lg&'()* V% cBm W\[ 52Oz &'()* V% cVq8(!IcF!&'()* V%tD8@D H <P&'(Evapotranspiration)* V% a@ k&'()* V% a 8(s@n&'()* V% a3#8I3`@ [.8&'()* V% c A@;; :A&'()* V% a}8&@g x&'()* V% aCy8C@! >tH&'()* V% c @ )6=&'()* V% apY8r@kT&'()* V% cx'@Lxk]&'()* V% aO[8OurҐmmAApJT&'()* V% c 8;<  ]`5F,'&'()* V% cD "8 @C07 W&'()* V% aNq8ITA;@0Z4n&'()* V% c "@OK&'()* V%t'] 3e&'(Basic Hydrology Module. Should be linked to Globals module (Globals_*), Evapotranspiration module (ET_*), Solar Radiation module (PAR_*))* V% cbow&'()* V% ai8 p))GG&'()* V% a8qo [@lf&'()* V% c 8@g &'()* V% cA@yV4N[&'()* V% c{8@g@&'()* V% aqi8Vs@MnQv&'()* V% c 8@g(&'()* V% Ds@)$v0`&'( Graph 1)r* V% a8(s@ln &'()* V% a8q.U@lP h&'()* V% aR8Gozz!BW&'()* V% a7@ k"<&'()* V% cO@aTT#B\ t&'()* V% c8@gL $:/&'()* V% aP/8@PM@%K&'()* V% cCx0& : 78%&'()* V% cw1 z'j|ax&'()* V% c\ment('icoV&'()* V% a:8:@~?)5?&'()* V% c 58@C0* 2>&'()* V% c\r@&+i,l&'()* V% a+d,_&'()* V% arh b-c&'()* V% a7 P.~<&'()* V% a$8@WL /$&'()* V% aW8@00R&'()* V% c+ @?418&'()* V% aT8I?T@)z2O:Y&'()* V% c+<&@38I&'()* V% a8R@q49.&'()* V% a8r@5&'()* V% a 8  @L6&'()* V% a&8@h7&'()* V% aV8&@8Q&'()* V% cq8@9S~f&'()* V% a<&8I\cP0<<:L9&'()* V% a8=dp;aa_&'()* V%t$8$H <#P&'(#Percolation, Water Table Rise/Drop)* V% azj8 @L=u &'()* V% a*8B@j>G#&'()* V% aOU8Aq@p?lJ&'()* V% cyB@IfIg@e3N&'()* V% a 8@~A&'()* V% a8#?@w}0B:&'()* V% c+9I@K8yPCN8&'()* V% c@yVD"q&'()* V% a8@E$&'()* V% a%8#D>@PF,&'()* V% aCC8C@G>H&'()* V% a`8@wpH[&'()* V% aJ [pI&'()* V% c`@J4I&'()* V% ce@closing0K ) h &'()* V% c@L@L N&'()* V% Dp@)$vM]&'( Graph 2)* V% c|@@\xNoJ&'()* V%tL@L9 OD A&'(Precipitation, Interception)* V% a<8I'R@rP7W,&'()* V% c@L@Qz[&'()* V% c86(Rdm&'()* V%t$8@$  S (&'( Infiltration)* V% c8@T 7&'()* V% cj@~ 0U]wS&'()* V% czA(8@xVm2VW&'()* V% a{MV8j{@x WvHo&'()* V% aQ8[[@XV$ &'()* V% cN@P?wPYA[\&'()* V% cZ8@ Z}g^&'()* V% cO PO.V[? p]X&'()* V% a8AB@Y\S&'()* V% a8[["p^xx  ]V&'()* V% afZ8e@^a&'()* V%tU 8@U " _M*&'(Initial Conditions)* V% c\8`ic&'()* V% c'8N eaw&'()* V% a8e]@{!b&'()* V% ae`8(f rP6xcbw&'()* V% a3)Q8mr@xd.$wr&'()* V% c;Aelq^b&'()* V% a|98r@fw&'()* V% ab8@j g&'()* V% a#8K@Nh&'()* V% cx8@0ia+&'()* V% coA@|A>|AKjbk=&'()* V% aF8jc@o`kAh&'()* V% c+8{l2s&'()* V% c*8@CAmE.&'()* V% c4i@ A`n'A_&'()* V% a8(h@{!oc&'()* V%tY@Ys} pQ{&'(Globals module)* V% a8eq6xSSq))&'()* V% c]@m۶mrj-&'()* V% al8_+"s&'()* V% c#@0Ypt W&'()* V% cYpuu `&'()* V% sQ8@CAvD-&'()* V% aa8 ew&'()* V% av'v 0xq{&'()* V% cG @y)T#z&'()* V% avy80.@zDb&'()* V% c"Dx @{(@`|&'()* V% a@:{8@ 0|G=&'()* V% a6 `};&'()* V% d80f~{&'()* V% aR8s (qqpp&'()* V% c:u `%IkP&'()* V%% cLP ?Y]&'()* V%%% aR8Es` (``@&'()* V% f BC!~0&'(&|d|*)* V% c4@1pp'AP&'()* V% c8a?  X&'()* V% a&8qP (&'()* V% as8 6F:n&'()* V% c89:+Eu&'()* V% a*8Ky@pp >&'()* V% a0#8K^a@pp O&'()* V% ag80gq|P (bl&'()* V% cf80Y&'()* V% a?K8@4K&'()* V% d;8@+002D&'()* V% fA%@`A`A,x;1 7$& C'(& !!R)* V% a{fN e~a&'()* V% a8:8@Tpn&'()* V%% cqr8@+0)3~ &'()* V% f  ^   &'(&/)* V% d8&'()* V% a8<4 me (&'()* V% a w8 @q&'()* V% cM P@@ZD)&'()* V% ck#zmq Mx >&'()* V% aR8mc (&'()* V% at8 N7@-w&'()* V%% cf8@+0pY8s&'()* V% ag80g@`bl&'()* V% cr8@%`&'()* V% aU8y@Đ+&'()* V%% a>U8=V@ 9[&'()* V% c8@5e=&'()* V% a8@=5"&'()* V% c*@hL&'()* V%% aP8dTTKi&'()* V% c} 8"@iPp;&'()* V% a%c(8SJ@p NOh&'()* V%% a8@&'()* V% c 'QBU &'()* V%% ci,@s a resu2&'()* V%%%%% ciO( d: K;vf[&'()* V% c   "&'()* V%% ae8qe@Đ`j&'()* V%%t@J R&'(Initial Conditions)* V% c/@OC`L&'()* V% a@ 5@&'()* V%%% cl8@ z_y&'()* V% al8+l@)gq&'()* V%%%%%% az@ & u&'()* V%tK@K C&'(Rate Coefficients & Constants)* V% c8@9z~&'()* V% a8$@&'()* V%%%%%% c@ e&'()* V% c@4y`&'()* V% fC022 0Q &338'(&uu)* V% c@7gW&'()* V% c @mmm( &56&'()* V% f2  ; g &33'(&cc)* V% a8@5 !&'()* V% c:3@9d9f-WG)P&'()* V% a6]8N3@)z ?Wn&'()* V% a]Z8F@)z Ab&'()* V% ay8r@iPm~&'()* V% f:A4@yAPyAQ%F%)4&-->*>*<<'(&00--11)* V% a[8/@~?~Z&'()* V% f722 jzmg&'(&N)* V% s8  g&'()* V% fUN922 .6}B&u>6>'(& P::)* V% s:Q8@P?~*[;&'()* V% f:0 (U A-.%&0MMO'(&bb"z"zF)* V% dz8 |q&'()* V% a: 8:@~?~5?&'()* V%t"@"< D&'(Macrophytes module)* V% cE  @& 8c&'()* V% a8U8  ]]$$=&'()* V% f ;G l:m& 22'(& b  )* V% aF F@& @AK#&'()* V% a3L8#3@s .8&'()* V% aFF@&AK&'()* V% fyv=22 +(Wd& 30i0i'(& u)* V%t<8@<8V 40^&'(Nitrogen module)* V% a8#<@P9&&'()* V% cLnD $z&'()* V% cK?04 3f08@&'()* V% a41lr@>&'()* V%t: @: \% 2d-&'()* V% a8jz@pPu&'()* V%t7@7eI /mQ&'()* V% c8@0P$%&'()* V%% cy@l&'()* V% cu<A! h&'()* V% ay8y@Yt}~&'()* V% c$@Zp1c&'()* V% cj@>@NN]&'()* V% ak8Ak@Pfp&'()* V% c@&U&'()* V% c @@@&'()* V% a 80[@M:`&'()* V%% c @A@&'()* V% aC8w@jHB&'()* V% aV8@PNQ&'()* V% ceFu  P~U&'()* V% cE*&D@; "8c=&'()* V% ayj8@Lt&'()* V% cL@'!&'()* V% cGuE@0.87 4.8:e&'()* V% a}98@ px&'()* V% c 6F>B  V&'()* V%% c5@ I*&'()* V% a:88KlY-@pp e&'()* V% aw8@q &'()* V% a|@8@zC @~&'()* V% c]w38Pj\&'()* V% c8@rِ&'()* V% c 8@rِQ&'()* V% cOX\ lqD&'()* V% ckG ^U&'()* V% cNfC@0A[] >|&'()* V% c 8@O&'()* V% aJ@0&'()* V% ac>8(T@jO#hC&'()* V% c&G@W&'()* V% a>8,@'C&'()* V% c}&8@C_Q1&'()* V% c!zH@Q3.Q&'()* V% a&8&@!+&'()* V% cIK 8@YL+5Vf@&'()* V%t@a  i&'( External)* V% ck  8@3&x &'()* V% aw8@ r&'()* V% c.8@]!!;#9&'()* V% c6z @")#C&&'()* V% c3V # & B F&'()* V% a8!@\$&'()* V% a8&@%&'()* V% a."8@z& ~3&'()* V% a8@ p'&'()* V% cgN@(P'x 79&'()* V% cF} 8@y)9IS5&'()* V% cJ@*+,m&'()* V% ai8*@?+|&'()* V% a8*@?,&'()* V% azv8;@A-&'()* V% am8@A.&'()* V% cFK@ffz/9S 13%&'()* V% ff\22 0;_1 z& VgV))'(& H+b+)* V% aL8/N@|1?%&'()* V% c}8@C|2_3&'()* V% a28/:@|3?&'()* V% c:J8@P4w"LPg&'()* V% a<)8@5A&'()* V% a8qq@6l&'()* V% aT8(r@7m&'()* V% cfaq 8Ys &'()* V% a-!8(s@9n2&'()* V% ck)@:Mkx ;>?BC&'()* V% aL8:M&1@;m&'()* V% a7  W ++<22<&'()* V% fYd2 2 =}Az;&II'(& b@@;)* V% ai8:i@>dn&'()* V% a3y8:|&8@?Mm&'()* V% aUK8K@L@PFP&'()* V% fkh22 A@F?& N[N['(& bEE))* V% aNo:8@@BYQrw&'()* V% cC:8@ C6PAG&'()* V% abpF @pDU&'()* V%%%%% fT3l22 I" : 2P&##'(&nW)* V%%% c3n L&&'()* V%%%% a:48:@PP5?&'()* V% s$q Q5lXdv&'()* V%%% s-s TqC&'()* V% s;t UQ?&'()* V%%%%% a|brS@l rZNm&'()* V%%%%%%%% cLh >Z@C=`b?[YZ&'()* V%%%%% aF48@PgA&'()* V%% c|1 $K  in&'()* V%%% c!b@lY&'()* V% aR8Pom&'()* V%%%% ce qX 6 &'()* V% c%F? r 5f&'()* V%%%%%%%%% cZ_@ t{v/g&'()* V%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% cr BHE&'()* V%%%%%%% fC' (QS!G&303K'(&tt?)* V% d1S8v (J:\&'()* V%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%tM?8M?>  E7F&'( Evaporation)* V%%%%%%%%%%%%%%%%%%%%%%% c<  0;Ts&'()* V%%% cTu6 Goe H2\9&'()* V%%%%%%% aW+8.?r&'()* V%%%%%%%%%%%% c>Yz] + 6/&'()* V%% cH8   ;U?&'()* V%%%%%% a8 &'()* V%%%%%%%% dX_08 OVah&'()* V% d.<8 %37E&'()* V% d}K=8 tBT&'()* V% d/08 &'89&'()* V% daA8 Xj&'()* V% dh8 _q&'()* V%%%%%% c*f  &'()* V% c+ %I : 4YE$&'()* V%%%%% c n $0&'()* V%%%%%%% cmo +]&'()* V%%%%% au$8@0)zNN990((&'()* V%%%%%%%% gE!9L Ea|8(i&' $Hydrology()* V%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ci@@Z\[b&'()* V% ai[8Zi@nj|0[dVn&'()* V%%%%%%% c~Z8@0j~pb"g&'()* V%%%%%%%%%%%%%%%%%%%%% c -Tw &'()* V%%%%%%%%%%%%%%%%%%%%%%%%%% c//I <&'()* V%%%%%%%%%%%% cC;UP 6<"&'()* V% aBB   ((=G&'()* V%%%%%%%%%%%%%%%%%%%%%%%%%%%% cB J  5R}&'()* V%%%%%%%%% cv i?iZx&'()* V%%%%%%%%% c  8 eDZ.}-I,&'()* V%%%%%%%%%%%%%%%%%%%%%% c86(@z&'()* V% a8&'()* Vdrmd)666 İ `R d L X t l2 +ZZ t  ( ?@$?@$)rState Variables,-B 0q ,@,33333@UUUU.q,j9+7>- 08 0@UUUU.8, 9- 0t $@Xpm.t,mmmmmmmmmm- 0s 4? =p@.s,m{;HF&=e.V+ t   ?@$?@$)r Untitled,- 0  @4@"yy@4@D.,CCCCCCCCCC- 0  @4@"yy@4@D.,- 0, @4@"yy@4@D.,,0000000000=e.V+ t  ?@$?@$)r Untitled,8<8@8Dע@H#" Vp)Y paR    tM lU + |   ?@$?@$)Flows,J- 0% @Xk.%,kkkkkkkkkk- 0: @Xk.:,kkkkkkkkkk- 0 @Xk.,kkkkkkkkkk- 0= @.=,rE,&.9CJ- 09 >uMUiAEjmn.9,>f/W+ |   ?@$?@$) Untitled,- 0 @0.,a%- 0l ?PbM@.l,Q7Ir- 09 >uMUiAEjmn.9,- 0 ?ЋC%@.,/>? >f/W+ |  ?@$?@$) Untitled,J- 0 0? 6C-A-.,{~ - 0 ,@Xk.,kkkkkkkkkk- 0 ?@Xk.,kkkkkkkkkk- 0) ?@Xk.),kkkkkkkkkk- 09 <>uMUiAEjmn.9,>f/W"M Vp/ ~ L  \ <   h        drmd N 0  X  0   drmd ~   8      drmd V  p  1 2T 3bbQK      t |F L (45?fOv_خ?htj~ 0 ?fOv_خ?f?g?htj~6KK787-747E7>7@77S77777]7777717A7 777Z7T7J7g77H7[7Q77R77"7Y77/7b7777777777777r77*7'77F7B777O7|777a77W7n77 77 7 77789  , Geneva.+ $Setup #1($ 3/1/98 10:51(?Input Variables(URun #)i SatOutRate(k1)_0.0027({2)_0.0028(3)_0.0029(4)_0.0030 **  :v2:;  < *<8$= nsm5)=qro4)r