User Tools

Site Tools


start:hype_model_description:hype_routing

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
start:hype_model_description:hype_routing [2018/10/18 15:25]
cpers [Links to file reference]
start:hype_model_description:hype_routing [2024/01/25 11:37] (current)
Line 3: Line 3:
 ===== Basic assumptions ===== ===== Basic assumptions =====
  
 +The HYPE model simulates most lakes and rivers in a conceptual way. The HYPE model can contain two types of rivers, local river and main river, and two types of lakes, local lakes and outlet lakes (Figure 1). The local river represent all ditches, streams, rivers within a subbasin that recieve only runoff within the subbasin. It is modelled as one conceptual river with a length, regardless of the actual distance runoff from different land has to travel through the subbasin. Similarly a local lake represents all lakes within a subbasin receiving inflow only from within the subbasin. The modelled local lake recieves a fraction of the flow in the local river. The fraction is defined as a fraction of the subbasin area. The main river is the river part(s) of a subbasin that water from upstream subbasins flows through. Depending on subbasin delieation it may be clearly defined or represent several tributaries that comes to the subbasin outlet. Finally the modelled outlet lake represents a lake located at the outlet of a subbasin.
  
-The HYPE model can contain two types of rivers, ​local stream and main river, and two types of lakes, local lakes and outlet lakes (Figure 1). Local and main rivers are present in all subbasins and the length of each is calculated as the square root of the subbasin area. The length of the watercourses ​can be given as input. The river can be a SLC class and is then given an area, but rivers can also be one-dimensional ​(i.e. no fraction of the subbasin area are occupied by the river and no precipitation added to the river). All local runoff is entering the local river. Local lakes (ilake) receive a portion of the local runoff. The flow leaving the local river (including flow from local lake) goes to the main river of the same subbasin. If there are upstream subbasins their flow is added to the local flow when both flows flow into the main river. Outlet lakes (olake) receive the outflow from the main river, i.e. all upstream and local flows.+local stream and main rivers are present in all subbasins ​as default, ​and the length of each watercourse ​is then calculated as the square root of the subbasin area. River lengths ​can be given as input, and then their value decide if there will be a stream or river in the subbasin. If the river length is zero no river is calculated. The river can be defined as a SLC class and has then an surface ​area in addition to a length, but rivers can also be without area (i.e. no fraction of the subbasin area are occupied by the river and no precipitation added to the river). ​
  
-|{{:​start:​hype_model_description:​lakeriveroverview.png?​400|}}| +All runoff from land enters the local river. Local lakes (ilake) receive a portion of the local flow. The flow leaving the local river and the outflow from local lake, goes to the main river of the same subbasin. If there are upstream subbasins their flow is added to the main river together with the local flow. Outlet lakes (olake) receive the outflow from the main river, i.e. all upstream and local flows. 
-|Figure 1: Schematic representation of streams and lakes HYPE, and the link between them.|+ 
 +| {{:​start:​hype_model_description:​lakeriveroverview.png?​400}} ​                                                ​
 +| Figure 1: Schematic representation of streams and lakes in a subbasin of HYPE, and the links between them.  |
  
 The two lake types are separate classes. The lake classes have characteristics such as land use and soil type, which are defined together with the other classes'​ characteristics (in GeoClass.txt). Precipitation,​ atmospheric deposition and evaporation of rivers and lakes are calculated first, while river flow and inflow, transformation processes and the outflow of the lakes is calculated thereafter. Lakes and rivers are calculated in the model’s routing part after all classes are calculated for the subbasin. ​ The two lake types are separate classes. The lake classes have characteristics such as land use and soil type, which are defined together with the other classes'​ characteristics (in GeoClass.txt). Precipitation,​ atmospheric deposition and evaporation of rivers and lakes are calculated first, while river flow and inflow, transformation processes and the outflow of the lakes is calculated thereafter. Lakes and rivers are calculated in the model’s routing part after all classes are calculated for the subbasin. ​
  
-An outlet lake can be part of a larger lake. It is then called a lake basin. Lake basins are olakes in nearby subbasins. Outlet lakes that are not lake basins are referred to below as simple outlet lakes+An outlet lake can be part of a larger lake. It is then called a lake basin. Lake basins are olakes in nearby subbasins. ​
  
-A simple outlet lake has a threshold. ​The outflow ends if the water level drops below the threshold. Lake mean depth below the threshold is specified in GeoData.txt or LakeData.txt as //​lake_depth//​ in meters. Lake depth can also be set by parameters, i.e general parameter //gldepo// or olake region parameter //​olldepth//​. ​The threshold is also the the water level of the lake at the start of a simulation. The current water level is denoted ​as //wlm// in Fig. 2. For printing, the outlet lake water level (output variable //wcom//) is calculated ​in meters and you can set reference ​level (//w0ref//) in LakeData.txt to get the same height system as any observations of the lake's water level. ​The lake’s ​//w0ref// is added to the water level above the threshold. ​HYPE assumes the lake has vertical sides in the calculations,​ thus the observed variation may be larger than the simulated variation. It is therefore possible to adjust the output //​wcom// ​for the actual amplitude of the regulation volume (//wamp//). This will make the simulated and recorded water stage comparable below the threshold for a regulated lake. +A simple outlet lake has a threshold. ​If the water level drops below the threshold, the outflow ceases. Lake mean depth below the threshold is specified in GeoData.txt or LakeData.txt as //​lake_depth//​ in meters. Lake depth can also be set by parameters, i.e general parameter //gldepo// or olake region parameter //​olldepth//​. ​At the start of a simulation ​the water depth of the lake is set to be at the threshold. The current water depth, ​denoted //wlm// in Fig. 2, can be above or below the threshold 
 + 
 +For output, the outlet lake water level (output variable //wcom//) is given in a user set reference ​systemThe reference level is used to get the same height system as any observations of the lake's water level. ​It uses the unit of meter, and the reference level (//w0ref//) gives the level of the threshold in the reference system (Fig. 2). The output //wcom// is calculated as the water level above the threshold ​plus //w0ref//A regulated ​lake (dam) has two thresholds. The distance between the thresholds are determined by the regulation volume (//​regvol//​) of the dam with the assumption that the dam has constant area with depth. In reality the area may decrease with depth. This means that the the observed variation ​in water level may be larger than the simulated variation. It is therefore possible to adjust the output //​wcom// ​(and //wcav//) with the actual amplitude of the regulation volume (//wamp//). This will make the simulated and recorded water stage comparable below the upper threshold for a regulated lake. 
  
 |{{:​start:​hype_model_description:​outletlakewithvariables2.png?​400|}}| |{{:​start:​hype_model_description:​outletlakewithvariables2.png?​400|}}|
Line 19: Line 24:
  
 A local lake also has a threshold depth that is used as start value. The depth is given by general parameter //gldepi// and is then the same for all the local lakes, or by ilake region parameter //​illdepth//​. It is measured in meters. A percentage of flow from the local stream flows into the local lake. The rest of the local flow runs directly to main river watercourse. ​ A local lake also has a threshold depth that is used as start value. The depth is given by general parameter //gldepi// and is then the same for all the local lakes, or by ilake region parameter //​illdepth//​. It is measured in meters. A percentage of flow from the local stream flows into the local lake. The rest of the local flow runs directly to main river watercourse. ​
- 
-Using parameters, you can divide the lake into two parts, one with faster flows (FLP) and one with slower flows (SLP) (Figure 3). This function is used for the simulation of nutrients to simulate stratification,​ strangulation or other phenomena that may limit the mixing of a lake. With this feature, the flows through the lake follow the schedule below. The split is determined by the parameter //​deeplake//,​ which is the fraction of the lake's initial volume SLP, the remaining (varying) volume in the lake is the FLP. The parameter //​fastlake//​ determines where the outflow will be coming from. Default is that outflow is taken from the slow lake part. Increasing //​fastlake//​ will let the FLP contribute to outflow. //​Fastlake//​ equal to one gives the maximum contribution of FLP, and the outflow will be taken proportionally from the two lake parts according to their volume. 
- 
-|{{:​start:​hype_model_description:​division_of_lake.png?​400|}}| 
-|Figure 3: Division of a lake in two parts and water flows regarding lake.| 
  
 The streams may have a delay and an attenuation of the flow. The former shifts the peaks, but leaves them otherwise untouched, while the attenuation both delays and smooths out peaks. The streams have a dead volume in which the residence time of solutes is increased relative to the water and the concentration smoothed out further over time. The streams may have a delay and an attenuation of the flow. The former shifts the peaks, but leaves them otherwise untouched, while the attenuation both delays and smooths out peaks. The streams have a dead volume in which the residence time of solutes is increased relative to the water and the concentration smoothed out further over time.
Line 31: Line 31:
  
 ^Parameter/​Data ^File ^ ^Parameter/​Data ^File ^
-|//rivlen, loc_rivlen, slc_nn, icatch//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]]|+|//rivlen, loc_rivlen, slc_nn//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]]|
 |special class code 1,2,11 and 12|[[start:​hype_file_reference:​geoclass.txt|GeoClass.txt]]| |special class code 1,2,11 and 12|[[start:​hype_file_reference:​geoclass.txt|GeoClass.txt]]|
 |//​lake_depth//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]] or [[start:​hype_file_reference:​lakedata.txt|LakeData.txt]] or [[start:​hype_file_reference:​damdata.txt|DamData.txt]]| |//​lake_depth//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]] or [[start:​hype_file_reference:​lakedata.txt|LakeData.txt]] or [[start:​hype_file_reference:​damdata.txt|DamData.txt]]|
 |//gldepo, olldepth, gldepi, illdepth//​|[[start:​hype_file_reference:​par.txt|par.txt]]| |//gldepo, olldepth, gldepi, illdepth//​|[[start:​hype_file_reference:​par.txt|par.txt]]|
-|//w0ref, wamp//​|[[start:​hype_file_reference:​lakedata.txt|LakeData.txt]] or [[start:​hype_file_reference:​damdata.txt|DamData.txt]]| +|//w0ref, regvol, wamp//​|[[start:​hype_file_reference:​lakedata.txt|LakeData.txt]] or [[start:​hype_file_reference:​damdata.txt|DamData.txt]]| 
-|//​deeplake,​ fastlake//​|[[start:​hype_file_reference:​par.txt|par.txt]]|+
 ==== Links to relevant modules in the code ==== ==== Links to relevant modules in the code ====
  
Line 42: Line 42:
 | [[http://​hype.sourceforge.net/​doxy-html/​namespacemodelmodule.html|modelmodule (model_hype.f90)]]|initiate_model_state| | [[http://​hype.sourceforge.net/​doxy-html/​namespacemodelmodule.html|modelmodule (model_hype.f90)]]|initiate_model_state|
 | ::: |initiate_model | | ::: |initiate_model |
-| [[http://​hype.sourceforge.net/​doxy-html/​namespacesurfacewater__processes.html|surfacewater_processes (sw_proc.f90)]]|calculate_landarea_riverlength+| ::: |calculate_flow_from_undivided_lake| 
-| ::: |calculate_outflow_from_outlet_lake|+| ::: |calculate_flow_for_lakebasin_lake | 
 +| [[http://​hype.sourceforge.net/​doxy-html/​namespacesurfacewater__processes.html|surfacewater_processes (sw_proc.f90)]]| ​calculate_riverlength| 
 +| ::: | calculate_ilake_outflow ​
 +| ::: | calculate_outflow_from_outlet_lake ​
 +| ::: | calculate_outflow_from_lakebasin_lake ​|
 | ::: |remove_outflow_from_lake| | ::: |remove_outflow_from_lake|
-| ::: |calculate_flow_within_lake| 
 | ::: |calculate_branched_flow| | ::: |calculate_branched_flow|
-|[[http://​hype.sourceforge.net/​doxy-html/​namespacenpc__surfacewater__processes.html|npc_surfacewater_processes (npc_sw_proc.f90)]]|set_lake_slowwater_maxvolume |+
  
  
Line 60: Line 63:
 === Evaporation === === Evaporation ===
  
-If the river has an area (is a class), it evaporates and a new concentration is calculated. Normally the river area is constant over time, but with parameters a reduction of riverarea can be simulated for low volume/​flow. The reduced river area is also used for heat exchange calculations.+If the river has an area (as a class), it evaporates and a new concentration is calculated. The river is assumed to evaporate at potential evaporation rate (see [[start:​hype_model_description:​processes_above_ground#​evaporation|Processes above ground - Evaporation]]). Normally the river area is constant over time, but with parameters a reduction of riverarea can be simulated for low volume/​flow. The reduced river area is also used for heat exchange calculations.
  
 <m> frac_area = delim{lbrace}{ ​ <m> frac_area = delim{lbrace}{ ​
 matrix{3}{2}{ matrix{3}{2}{
     0 {x<=0}     0 {x<=0}
-    {(1+(x-fraxe)/​(fraxe-fraxm))*(x/​fraxe)^{x/(x-fraxm)}} ​ {0<​x<​fraxe}+    {(1+(fraxe-x)/​(fraxe-fraxm))*(x/​fraxe)^{fraxe/(fraxe-fraxm)}} ​ {0<​x<​fraxe}
     1 {x>​=fraxe} ​     1 {x>​=fraxe} ​
     }}{} </m>     }}{} </m>
        
-The parameters //fraxe// and //fraxm// are general, and //fraxe// is (minimum) mean river depth (m) where fractional river area = 1 and //fraxm// is mean river depth (m) where the slope of the fractional river area has its maximum (must be in the range between 0 and //fraxe//). //x// is the current mean river depth (based on full area extention).+The parameters //fraxe// and //fraxm// are general, and //fraxe// is (minimum) mean river depth (m) where fractional river area = 1 and //fraxm// is mean river depth (m) where the slope of the fractional river area has its maximum (must be in the range between 0 and //fraxe//). //x// is the current mean river depth (estimated as river water volume (including possible deadvolume) divided by the full area extention).
  
 === Pure delay === === Pure delay ===
  
-The delay in the watercourse (//​transtime//​) in days is determined by the length of the watercourse (//​rivlen//​) and the water’s maximum velocity (//​rivvel//​). The maximum velocity is a general parameter with unit //m/s//. The delay in the river is dependent on subbasin land area if the default river length is used. The delay is a pure translation. The delay is divided into whole days (//ttday//) and parts of the day (//​ttpart//​).+The delay in the watercourse (//​transtime//​) in days is determined by the length of the watercourse (//​rivlen//​) and the water’s maximum velocity (//​rivvel//​). The maximum velocity is a general parameter with unit //m/s//. The delay in the river is dependent on subbasin land area if the default river length is used. The delay is a pure translation ​unless the parameter //damp// is set. The pure delay is divided into whole days (//ttday//) and parts of the day (//​ttpart//​) ​(if the time step is daily, other wise the division is made into whole and fractions of that time step). In the following description daily time step is assumed.
  
 <m> transtime = rivlen/​{rivvel*8.64*10^4} </m> <m> transtime = rivlen/​{rivvel*8.64*10^4} </m>
Line 79: Line 82:
 <m> ttday = INT(transtime) </m> <m> ttday = INT(transtime) </m>
  
-<m> ttpart = transtime - REAL(ttday</m>+<m> ttpart = transtime - ttday </m>
  
 The inflow of the river is stored in two arrays (//riverq// and //riverc//) until it is time for it to flow out of the river stretch. The outflow is weighted by using the parts of the time step (//​ttpart//​) that are to flow out during the time step. The inflow of the river is stored in two arrays (//riverq// and //riverc//) until it is time for it to flow out of the river stretch. The outflow is weighted by using the parts of the time step (//​ttpart//​) that are to flow out during the time step.
Line 101: Line 104:
 <m> transtime=(1-damp)*totaltime </m> <m> transtime=(1-damp)*totaltime </m>
  
-The result from the translation of water (//transq// and //transc//) flows into the attenuation box, which is assumed to be completely mixed. The delay time in the box (//kt//) is recalculated to a corresponding recession coefficient (//​riverrc//​). The recession coefficient (//​riverrc//​) used to calculate the outflow from the box (//dampq//) is a function of volume in the box (//​riverbox//​).+The result from the translation of water (//transq// and //transc//) flows into the attenuation box, which is assumed to be completely mixed. The delay time in the box (//kt, in timesteps//) is determined form input and parameters
  
 <m> kt = damp * totaltime </m> <m> kt = damp * totaltime </m>
 +
 +The outflow equation depends on the attenuation box's volume before changes during the time step (<m> riverbox_ini </m>) and on the net inflow during the timestep (//​inflow//​). The net inflow includes the outflow from the translation (local runoff, upstream inflow, point sources etc.), removal though irrigation or regional groundwater flow, and the flow to or from a main river floodplain. The equation has two coefficients:​
 +
 +<m> riverrc1 = 1. - kt + kt * e^{-1{/}kt} </m>
 +
 +<m> riverrc2 = 1. - e^{-1{/}kt} </m>
 +
 +<m> dampq = riverrc1 * inflow + riverrc2 * riverbox_ini </m>
 +
 +An earlier version of calcuation of outflow from the box can be accessed through a model option. In this case the delay time of the attenuation box is recalculated to a recession coefficient (//​riverrc//​). The recession coefficient is used to calculate the outflow from the box (//dampq//) by the current volume in the box (//​riverbox//​).
  
 <m> riverrc = 1. - kt + kt * e^{-1{/}kt} </m> <m> riverrc = 1. - kt + kt * e^{-1{/}kt} </m>
Line 109: Line 122:
 <m> dampq = riverrc * riverbox </m> <m> dampq = riverrc * riverbox </m>
  
-The delay in the attenuation box is at the centre of gravity of the flow (not the tip) and is not calculated exactly (see Derivation of the equation of recession in the attenuation box below). Therefore there may be a difference in a peak delay between the use of pure delay and delay with damping, although the maximum velocity is set to the same value. In the case of damping the max velocity parameter ​//rivvel// may be approximated.+where <m> riverbox = riverbox_ini + inflow </m>.
  
-=== Additional smoothing of dead volume === 
  
-dead volume ​(//​deadriver//​) can be used to get extra equalization of concentrations in the watercourse. Then the outflow equation for the attenuation box (//​riverbox//​) is:+=== Additional smoothing by dead volume ​===
  
-<​m> ​dampq = riverrc*(riverbox - deadriver) ​</m> +A dead volume (//​deadriver//, ​<m>m^3</m>) can be used to get additional equalization of concentrations in the watercourse. ​In this case, the initial value of the attenuation box is equal to the dead volume. Inflows into the box are mixed with the entire volume. At present this represents a dilution of the concentration since the initial value of concentration in the box is zero. In HYPE the dead volume ​is estimated from the upstream area (including ​the area of the subbasin itself). This gives the dead volume:
- +
-In this case, the initial value of the attenuation box is equal to the dead volume. Inflows into the box are mixed with the entire volume. At present this represents a dilution of the concentration since the initial value of concentration in the box is zero. At present ​the estimated ​dead volume ​as dependent on the upstream area (inkluding ​the area of the subbasin itself. Depth is calculated from the parameter //dead// (//deadl// or //deadm//, <​m>​m^2{/​}km^2</​m>​and the upstream area (<​m>​km^2</​m>​),​ the width is assumed to be 10 times the depth and the river length (//m//) is estimated to be the root of the subbasin’s area unless it is set among the input data. This gives the dead volume:+
  
 <m> deadriver = dead * upparea * rivlen </m> <m> deadriver = dead * upparea * rivlen </m>
  
-Dead volume can be used even without suppression of water (//damp = 0//).+where the parameter //dead// (//deadl// or //deadm//, <​m>​m^2{/​}km^2</​m>​) is general parameters, the upstream area //uparea// (<​m>​km^2</​m>​) is calculated from input data, and the river length (//​rivlen//,​ //m//) is estimated to be the root of the subbasin’s area unless it is set among the input data. Dead volume can be used without suppression of water flow peak (//damp = 0//). When additional smoothing is used, the //​riverbox//​ in the attenuation equations above is replaced by the volume of the attenuation box minus the dead volume. 
 + 
 +In addition, the dead volume is used to calculate a minimum river width used for calculation of river surface area in the case the river area is not given as a class area. The width is assumed to be 10 times the depth and together with the dead volume and river length, the with is calculated as: 
 + 
 +<m> deadwidth = 10 * {(deadriver / {rivlen * 10})}^{1/2} </m>
  
 === Derivation of the equation of recession in the attenuation box === === Derivation of the equation of recession in the attenuation box ===
Line 151: Line 165:
 <m> e^{t{/}k} q(t) - q(0) = I~(e^{t{/​}k}-1) </m> <m> e^{t{/}k} q(t) - q(0) = I~(e^{t{/​}k}-1) </m>
  
-<m> q(t)=q(0)e^{-t{/​}k}+(1-e^{-t{/​}k}) </m>+<m> q(t)=q(0)e^{-t{/​}k}+I~(1-e^{-t{/​}k}) </m>
    
 Calculate //Q// for time step (the mean of //q(t)//) and replace the outflow with the  Calculate //Q// for time step (the mean of //q(t)//) and replace the outflow with the 
Line 160: Line 174:
 <m> {}= q(0) int {0} {1}{} e^{-t{/}k} dt + I~ int{0}{1}{} 1-e^{-t{/​}k} dt </m> <m> {}= q(0) int {0} {1}{} e^{-t{/}k} dt + I~ int{0}{1}{} 1-e^{-t{/​}k} dt </m>
  
-<m> {}= q(0)*(-k~(e^{-1{/​}k}-1))+I(1-0)-I~(-k~(e^{-1{/​}k}-1)) </m>+<m> {}= q(0)*({-}k~(e^{-1{/​}k}-1))+I(1-0)-I~({-k}~(e^{-1{/​}k}-1)) </m>
   ​   ​
   ​   ​
   ​   ​
-<m> {}= I~(e^{-1{/​}k}-1)~(kI - kq(0)) </m>+<m> {}= I+(e^{-1{/​}k}-1)~(kI - kq(0)) </m>
   ​   ​
   ​   ​
   ​   ​
-<m> Q=(I+(e^{-1{/​}k}-1))~(kI-S(0)) </m>+<m> Q=I+(e^{-1{/​}k}-1)~(kI-S(0)) </m>
  
 +or separated for //I// and //S//.
 +
 +<m> Q=(1-k + k*e^{-1{/​}k})~I + (1 - e^{-1{/​}k})~S(0) </m>
    
  
-Recession ​coefficient is defined as the proportion of outflow after an inflow pulse. It is therefore calculated ​here as <​m>​r=Q{/​}I</​m>​ under the assumption that the starting volume can be neglected, //s(0)=0//.+Common in HYPE is that outflows are calculated with a recession coefficient. Here the recession ​coefficient is defined as the proportion of outflow after an inflow pulse. It is therefore calculated as <​m>​r=Q{/​}I</​m>​ under the assumption that the starting volume can be neglected, //S(0)=0//​. ​
  
 <m> r=Q/​I=1+(e^{-1{/​}k}-1)*k </m> <m> r=Q/​I=1+(e^{-1{/​}k}-1)*k </m>
  
 +For a real flow the volume can of course not be neglected in this way and this relationship is only approximate.
 +
 +Alternatively the outflow may be calculated from the equation separated for //I// and //S//. For this option I is calculated as the sum of several different types of inflow and outflow to the river during the time step, e.g. precipitation and evaporation,​ in addition to the lateral inflow. In addition the initial value of the river volume (//S//(0)) is saved here to be used for the outflow calculation. ​
  
-For a real flow the volume can of course not be neglected in this way and the relationship is only approximate. 
    
 === Bank-full flow === === Bank-full flow ===
  
  
-Bank-full flow is used for erosion of particulate phosphorus in rivers (see Sedimentation). It is approximated ​by the second highest daily flow during the last year.+Bank-full flow is used for calculation of erosion of sediment and particulate phosphorus in rivers (see [[start:​hype_model_description:​hype_np_riv_lake#​sedimentation_resuspension|Sedimentation/​Resuspension]]). This flow has a probability of app. 0.5% (Nixon, 1959). It is in HYPE estimated ​by the second highest daily flow during the last 365 days, which has similar probability (Rosberg, 2003).
  
  
Line 188: Line 207:
 ==== Local river ==== ==== Local river ====
  
-The local river has a length equal to the square root of the subbasin area, if not specified as input. Runoff from the land area of the subbasin forms the inflow to the local river. The flow in the local river is delayed and attenuated as described above. Of the resulting flow from the stream a constant percentage goes to the local lake (//​icatch//​),​ and the rest directly to the main river. ​+The local river has a length equal to the square root of the subbasin area, if not specified as input. Runoff from the land area of the subbasin forms the inflow to the local river (Figure 3). The flow in the local river is delayed and attenuated as described above. Of the resulting flow from the stream a constant percentage goes to the local lake (//​icatch//​),​ and the rest directly to the main river. ​ 
 + 
 +|{{:​start:​hype_model_description:​HYPE_box_picture_v3_localriver.png?​400|}}| 
 +|Figure 3: Water flows of the local river.|
  
 ==== Main river ==== ==== Main river ====
  
-A main river is present in all subbasins. The length is equal to the square root of the subbasin area, if it is not specified as input. In areas without incoming water from upstream there is still a main river, but it receives only local river flow after the local lake. In subbasins with upstream incoming water, the flow to the main river will be the sum of outflow from the local lake, the proportion of flow in the local river not flowing into the local lake and the water from upstream.+A main river is present in all subbasins. The length is equal to the square root of the subbasin area, if it is not specified as input. In areas without incoming water from upstream there is still a main river, but it receives only local river flow after the local lake. In subbasins with upstream incoming water, the flow to the main river will be the sum of outflow from the local lake, the proportion of flow in the local river not flowing into the local lake and the water from upstream ​(Figure 4). For a subbasin with a (part of) a multi-basin lake, the main river transports the flow from land to the lakebasin of the subbasin.
  
-The return flow from an aquifer is added to the inflow of the main river.+|{{:​start:​hype_model_description:​HYPE_box_picture_v3_mainriver.png?​400|}}| 
 +|Figure 4: Water flows of main river.| 
 + 
 +The return flow from an aquifer is added to the inflow of the main river (Figure 4).
  
 The flow in the main river is delayed and attenuated as described above. ​ The flow in the main river is delayed and attenuated as described above. ​
 The resulting effluent from the river flows to the outlet lake or to the next subbasin if there is no outlet lake.  The resulting effluent from the river flows to the outlet lake or to the next subbasin if there is no outlet lake. 
 +
 +If the main river length is zero, the inflows to the main river is added to form the main river outflow.
  
 === Inflow from upstream subbasins === === Inflow from upstream subbasins ===
  
-In input files it is given to which subbasin(s) the outflow from each subbasin flows. The upstream flow enters the main river of the downstream subbasin. Inflow into the main river of a subbasin is calculated by adding outflows from upstream areas. Concentrations are flow-weighted by their relative share. ​+In input files it is given to which subbasin(s) the outflow from each subbasin flows. The upstream flow enters the main river of the downstream subbasin. Inflow into the main river of a subbasin is calculated by adding outflows from upstream areas. Concentrations are flow-weighted by their relative share. ​The upstream flow enters the main river, except if the upstream subbasin has a lakebasin that is part of the downstream subbasin'​s lake. In this case the upstream lakebasin outflow goes directly into the downstream lakebasin.  
 + 
 +=== Hydraulic geometry to calculate water depth and velocity === 
 + 
 +River diagnostic variables can be calculated for the main river. A river water depth (//d//) and river velocity (//u//) can be calculated from hydraulic geometry for the main river and the calculated main river flow (//q//). In the case of updating of flow in a subbasin with no outlet lake, the updated (main river) flow is used. The equation coefficients are general parameters (//​c,​f,​k,​m//​). 
 + 
 +<m> d=c*q^f </​m>​ 
 + 
 +<m> u=k*q^m </​m>​ 
 +  
 +=== Rating curve to calculate river water level=== 
 + 
 +Given a continous monotonous increasing equation for calculating water level from flow, the water level of a river (//wr//) can be calculated from the flow by the inverse of this rating curve. This method is used to calculate main river water level in HYPE. Three different types of river rating curve formulations can be used (//​TYPEofRRC//​).  
 + 
 +The first is a single rating curve with a threshold (//w0//) and the rating curve coefficients are input. The rating curve  
 + 
 +<m> q=k*(wr-w0)^p </​m>​ 
 + 
 +give coefficients //a// and //b// for the inverse equation: 
 + 
 +<m> a={1/​k}^{1/​p} ​ b=1/p </​m>​ 
 + 
 +<m> wr = a*q^b + w0 </​m>​ 
 + 
 +The second type uses two rating curves with the same threshold (//w0//); one for ice free condition and one for use when there is ice on the river. The transition between the rating curves is smoothed by scaling of the curves in a transition period. The scaling of rating curve parameters are determined by ice porosity and ice thickness (general model parameters). Between thresholds the influence of the ice curve increases linearly from zero to one. The minimum influence (of porosity and thickness) is used as the scaling weight of the winter curve during the transition period. If only parameters for the ice free curve are given as input, the same parameters is used for ice conditions except for the exponent (//p//) which is then set to one.  
 + 
 +The third type is for a sectorial rating curve. It has different rating curve coefficients for different sections of water level. The sections may not overlap or leave gaps in water level, even though small discrepancies is handled. The parameters and sections need to be chosen so that the whole curve is continous. Since the rating curves are used inverted in the code, the sections are turned into sectors of flow. The average of the flow given by the minimum water level of one curve and the maximum water level of another curve is taken as the actual limit between the two curves. This type does not handle ice curves. 
 + 
 +In addition to river water level described above, it is possible to calculate the water level in the scale of the local gauge. This could e.g. be in cm and with the zero of the gauge at other location than //w0// of the equation. Figure 5 below show the relation between the different parameters.  
 + 
 +| {{:​start:​hype_model_description:​river rating curves.png?​600|}} | 
 +| Figure 5: Schematic representation of river water level. ​ | 
  
  
Line 211: Line 270:
 |:::​|//​dead//​|//​deadm,​ deadl//​|:::​| |:::​|//​dead//​|//​deadm,​ deadl//​|:::​|
 |Local river| |//​icatch//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]]| |Local river| |//​icatch//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]]|
 +|:::| |//gicatch, illicatch//​|[[start:​hype_file_reference:​par.txt|par.txt]]|
 |Main river| |//​maindown//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]]| |Main river| |//​maindown//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]]|
 |:::| |//​branchid//​|[[start:​hype_file_reference:​branchdata.txt|BranchData.txt]]| |:::| |//​branchid//​|[[start:​hype_file_reference:​branchdata.txt|BranchData.txt]]|
 +|:::|//c, f, k, m// |//hygeomc, hygeomf, hygeomk, hygeomm//​|[[start:​hype_file_reference:​par.txt|par.txt]]|
 +|:::|  |//​ricew0por,​ ricew1por, ricew0ice, ricew1ice//​|:::​|
 +|:::|//k, p, w0// |//k, p, w0ref//​|[[start:​hype_file_reference:​RiverRatingCurve.txt|RiverRatingCurveData.txt]]|
 +|:::| |//wmin, wmax, gaugezero//​|:::​|
  
 ==== Links to relevant modules in the code ==== ==== Links to relevant modules in the code ====
Line 226: Line 290:
 |:::​|calc_qbank|bank-full flow| |:::​|calc_qbank|bank-full flow|
 |:::​|update_qbank |:::| |:::​|update_qbank |:::|
 +|:::​|river_water_level |main river|
 +|:::​|local_water_level|:::​|
 +|:::​|ice_on_river|:::​|
  
 +
 +==== References ====
 +
 +Nixon, M. 1959. A study of the bank-full discharges of rivers in England and Wales. Journal of the Institute of Civil Engineers 12:157-174.
 +
 +Rosberg, J. 2003. Modelining phosphorus transport and retention in river networks. Master Thesis. Dept. of Earth Sciences, Uppsala University.
  
  
Line 235: Line 308:
 === Precipitation === === Precipitation ===
  
-Daily precipitation is added to the lake and a new concentration is calculated. Lakes are assumed to be completely mixed if the parameter //​deeplake//​ is not set to a value greater than zero.+Daily precipitation is added to the lake and a new concentration is calculated. Lakes are assumed to be completely mixed.
  
 === Evaporation === === Evaporation ===
  
-The lake is assumed to evaporate and a new concentration calculated.+The lake is assumed to evaporate ​at potential evaporation rate (see [[start:​hype_model_description:​processes_above_ground#​evaporation|Processes above ground - Evaporation]]) ​and a new concentration calculated. The actual evaporating area is reduced if the lake is (partly) ice covered. Ice is calculated in case water temperature (substance T2) is simulated.
  
 === Outflow with rating curve ==== === Outflow with rating curve ====
  
-The rating curve is used for calculation of outflow above a threshold ​at several occasions+The rating curve is used for calculation of outflow ​in different circumstances. It is based on water level above a threshold. ​
  
 <m> q(t)=k*(w(t)-w_0)^p </m> <m> q(t)=k*(w(t)-w_0)^p </m>
    
-In the program the equation is solved with linearization (Lindström,​ 2016). For rating curve exponent equal to one, the solution is exact, while for other exponent an approximation is used.+In the program the equation is solved ​for the average flow during the timestep ​with linearization (Lindström,​ 2016). For rating curve exponent equal to one, the solution is exact, while for other exponent an approximation is used. The outflow of the lake is limited to the water above the threshold, in the case that the rating curve give a higher value.
 ==== Local lake (ilake) ==== ==== Local lake (ilake) ====
  
-The local lake has an initial volume which is determined by its area and threshold depth. The lake depth is determined by model parameters (//gldepi// or //​illdepth//​)). +The local lake has an initial volume which is determined by its area and threshold depth. The lake depth below the threshold ​is determined by model parameters (//gldepi// or //​illdepth//​). ​
  
-The inflow to the lake is a percentage of the flow in the local stream. The percentage is determined by the percentage of the subbasin area that drain to the lake. This percentage can be given as a fraction in GeoData.txt ​(//​icatch//​) ​for each subbasin with an internal lake, or be given by a regional parameter //​ilicatch//,​ or be given by a general parameter //​gicatch//​. If not set at all the default value is 1, i.e. the local river runs through the loacal ​lake. The flow from the local river is added to the lake. The lake water is assumed completely mixed if //deeplake = 0//+The inflow to the lake is a percentage of the flow in the local stream ​(//​icatch//​). The percentage is determined by the percentage of the subbasin area that drains ​to the lake. This percentage can be given as a fraction in GeoData.txt for each subbasin with an internal lake, or be given by a regional parameter //​ilicatch//,​ or be given by a general parameter //​gicatch//​. If not set at all the default value is 1, i.e. the local river runs through the local lake. The flow from the local river is added to the lake (Figure 6). The lake water is assumed completely mixed. ​ 
 + 
 +|{{:​start:​hype_model_description:​HYPE_box_picture_v3_ilake.png?​400|}}| 
 +|Figure 6: Water flows of internal lake.|
  
 Water outflow is calculated with the universal rating curve using general or region specific model parameters. Water outflow is calculated with the universal rating curve using general or region specific model parameters.
  
 +The above description of the functioning of an internal lake (ilake) in HYPE can be changed with a modeloption for ''​connectivity''​. The default ilake model has no connectivity model. There are two alternatives for introducing lake connectivity within a subbasin. The first is a fill and spill model for a chain of storages. The second is a the Hysteretic Depressions Storage model (HDS). You can choose to use one or the other or to have different ilake descriptions for different subbasins.
 +
 +===The fill-and-spill of ilake sections model===
 +
 +This ilake connectivity model simulates the ilake as a chain of small connected lake sections (//lks//) with different height of their outflow threshold. Depending on water volume in the ilake one or several of the lake sections can be disconnected from the chain (zero outflow with water level below the threshold). The precipitation and evaporation are divided by the lake sections depending on their fraction of the ilake area, but the runoff from land area is divided by the lake sections depending on their fraction of the catchment area. In case the water level of a lake section is above its threshold, its outflow is calculated. The outflow of each lake section is calculated using the universal rating curve using general or region specific model parameters same as a simple ilake. The outflow of a lake section is of course also added to the receiving lake section before that one’s outflow is calculated. Depending on how far upstream from the last lake section of the chain there is inflow from the above lake section, that fraction of the ilake area is the ilake’s current connectivity. The catchment area of the lake sections above this border is the currently non-contributing catchment area.
 +
 +===The Hysteretic Depressions Storage (HDS) model ===
 +
 +Another model for the local lake is the Hysteretic Depressions Storage (HDS) model (Ahmed et al. 2023, Shook et al. 2021). The HDS model does not simulate depressions as distinct (individual) depressions,​ rather, it simulates the effects of the integrated depressions system and their connection at the basin level as a lumped depression storage (ilake). The HDS model simulates the change in the contributing area, depressional storage, and outflow from the subbasin.  ​
 +
 +The model has contributing area fraction as a state. The runoff from the contributing area is outflow from the HDS model of the subbasin. The non-contributing area fraction of runoff is added to the depression storage. If the water of the depression is above the threshold of the depression, all the water above the threshold is also contributing to outflow of the HDS model of the subbasin. In this case the contributing area fraction will be one for the next timestep, otherwise the contributing fraction is calculated with a hysteretic function. The hysteretic function will keep the contributing fraction constant if the water volume in the depression is un-changed, if it is decreasing (drying up) the contributing fraction is set to zero. If the depression is filling up, the contributing fraction is increasing linearly on a curve to reach a contributing fraction of one at the same time as the depression is filled to the max.
 +
 +Input are the maximum depth (//​hds_depth//​) and maximum area (ilake slc-area) of depression in GeoData.txt (or a ilakeregion dependent parameter //​hdsdepth//​ can be used). In addition, the land area fraction connected to the ilake, //icatch//, is used. It is set equal to one if the whole subbasin may contribute to the HDS outflow. //icatch// can be lower if part of the land area of the subbasin is always contributing to subbasin outflow, e.g. if it is contributing to a main river through the subbasin. The output '​connectivity'​ of the HDS model is in the HYPE implementation defined as the contributing area fraction, aka one minus the non-contributing area fraction. This because, at the moment, the water area of the depression is constant. The output '​non-contributing area fraction'​ is a the fraction of the ilake catchment area (including the ilake area itself) not contributing to outflow from the HDS model.
  
 ==== Simple outlet lake or dam (olake) ==== ==== Simple outlet lake or dam (olake) ====
 +
 +| {{:​start:​hype_model_description:​simple_dam.png?​400}} ​ |
 +| Figure 7: Water flows of outlet lake.                 |
 +
  
 === Inflow === === Inflow ===
  
-The outlet lake receives both local runoff and inflow from upstream areas via the main river. An outlet lake may receive inflow from regional groundwater flow. +The outlet lake receives both local runoff and inflow from upstream areas via the main river. An outlet lake may receive inflow from regional groundwater flow (Figure 7)
  
 === Outflow with universal rating curve === === Outflow with universal rating curve ===
Line 267: Line 360:
 If the water level (//wlm//) is higher than the threshold (//​lake_depth//​) the outflow is calculated with the following equation: If the water level (//wlm//) is higher than the threshold (//​lake_depth//​) the outflow is calculated with the following equation:
  
-<m> outflow = delim{lbrace}{matrix{3}{3}{{0} {,~for} {wlm<​lake_depth} {gratk*uparea^{grata}*ratcorr*(wlm-lake_depth)^{gratp}} {,~for} {grata>​0} {gratk*ratcorr*(wlm-lake_depth)^{gratp}} {,~for} {grata=0} }}{} </m>+<m> outflow = delim{lbrace}{matrix{3}{3}{{0} {,~for} {wlm<​lake_depth} {gratk*uparea^{grata}*corr*(wlm-lake_depth)^{gratp}} {,~for} {grata>​0} {gratk*corr*(wlm-lake_depth)^{gratp}} {,~for} {grata=0} }}{} </m>
  
-//gratk//, //gratp// and //grata// are general parameters that apply to all lakes in the model application. //ratcorr// is parameter region dependent model parameter for adjusting //gratk// to different regions. The upstream area (//​uparea//,​ <​m>​km^2</​m>​) is included in the equation if parameter //grata// is >0. If the water level is below the threshold then the outflow is zero.+<m> corr = 1+ratcorr </​m>​ 
 + 
 +//gratk//, //gratp// and //grata// are general parameters that apply to all lakes in the model application. //ratcorr// is parameter region dependent model parameter for adjusting //gratk// to different regions. The upstream area (//​uparea//,​ <​m>​km^2</​m>​) is included in the equation if parameter //grata// is >0. If the water level is below the threshold then the outflow is zero. If the flow calculated by the equation gives a larger outflow volume during a time step than the water volume above the threshold, the outflow is limited so that the water level do not drop below the threshold. It is possible to prescribe another minimum flow (//​minflow//​) than zero. In this case the outflow is no lower than the minimum flow as long as there is water in the lake. That is the water level can drop below the threshold
  
 === Outflow with specific rating curve === === Outflow with specific rating curve ===
  
-In LakeData.txt a specific rating curve may be set instead of the general parameters. This is done by setting the //rate// and //exp// larger than zero in LakeData.txt.+In LakeData.txt a specific rating curve may be set instead of the general parameters. This is done by setting the //rate// and //exp// larger than zero in LakeData.txt. The outflow is also in this case limited by the threshold not to cause the water level to drop below it. It is possible to prescribe another minimum flow (//​minflow//​) than zero. In this case the outflow is no lower than the minimum flow as long as there is water in the lake. That is the water level can drop below the threshold
  
 <m> outflow = rate * (wlm - lake_{depth}) ^ {exp} </m> <m> outflow = rate * (wlm - lake_{depth}) ^ {exp} </m>
 +
  
 === Regulated lake === === Regulated lake ===
  
-A simple regulation routine can be given in LakeData.txt for selected outlet lakes. There are two main versions of regulations. The first determine ​a (constant) production flow between the lake threshold and a lower water stage. The second ​determine ​different rating curves for two production periods.+A simple regulation routine can be given in LakeData.txt for selected outlet lakes. There are two main versions of regulations. The first determines ​a (often constant) production flow, for water levels ​between the lake threshold and a lower water stage threshold, and a spill for higher water levels. The second ​determines the outflow by different rating curves for (two) different seasons. For both variants it is possible to prescribe another minimum flow (//​minflow//​) than zero. In this case the outflow is no lower than the minimum flow as long as there is water in the lake. That is the water level can drop below the threshold(s)
    
 == Production flow == == Production flow ==
Line 285: Line 381:
 The regulated outflow is independent of water level (//wlm// (//m//)) between the threshold (//​w0=lake_depth//​) and a minimum water level (//wmin//). For this interval the production flow (//qprod//) is used. Production flow may depend on time of year and the water level. The minimum water stage is not really a water stage, but is calculated from the regulation volume (//​regvol//​) and the lake (surface) area. Thus it does not consider the area changing with depth. The regulated outflow is independent of water level (//wlm// (//m//)) between the threshold (//​w0=lake_depth//​) and a minimum water level (//wmin//). For this interval the production flow (//qprod//) is used. Production flow may depend on time of year and the water level. The minimum water stage is not really a water stage, but is calculated from the regulation volume (//​regvol//​) and the lake (surface) area. Thus it does not consider the area changing with depth.
  
-|{{:​start:​hype_model_description:​illustraion_of_regulation_flow.png?​400|}}| +| {{:​start:​hype_model_description:​illustraion_of_regulation_flow.png?​400}} ​ 
-|Figure ​4: Illustration of regulation flow and variables.|+| Figure ​8: Illustration of regulation flow and variables. ​                 |
  
-For water levels above the threshold, the flow is calculated with a rating curve (typically flow through spillways) or all the water is discharged, ​but it is at least equal to the production flow (//​qprod//​). Compiled together as one equation, the outflow from a regulated lake is:+For water levels above the threshold, the flow is calculated with a rating curve (typically flow through spillways) or all the water above the threshold ​is discharged. In this case thoughthe calculated flow may be increased so that the outflow ​is at least equal to the production flow (//​qprod//​). Compiled together as one equation, the outflow from a regulated lake is:
  
 <m> outflow = delim{lbrace}{matrix{4}{3}{{MAX({rate * (wlm-w0)}^{exp},​ qprod)} {,~for} {wlm>​w0,​rate>​0} {MAX((wlm-w0) * 1000 / qunitfactor,​ qprod)}{,​~for}{wlm>​w0,​rate<​=0} {MIN(qprod,​(wlm-wmin) * 1000 / qunitfactor)} {,~for} {wmin<​wlm<​w0} ​ {0} {,~for} {wlm<​wmin}}}{} </m> <m> outflow = delim{lbrace}{matrix{4}{3}{{MAX({rate * (wlm-w0)}^{exp},​ qprod)} {,~for} {wlm>​w0,​rate>​0} {MAX((wlm-w0) * 1000 / qunitfactor,​ qprod)}{,​~for}{wlm>​w0,​rate<​=0} {MIN(qprod,​(wlm-wmin) * 1000 / qunitfactor)} {,~for} {wmin<​wlm<​w0} ​ {0} {,~for} {wlm<​wmin}}}{} </m>
Line 300: Line 396:
 <​m> ​ qprod = qprod*(1. + qamp * sin(2*pi*(dayno+qpha)/​365)) </m> <​m> ​ qprod = qprod*(1. + qamp * sin(2*pi*(dayno+qpha)/​365)) </m>
  
-In addition, the production flow is reduced ​from this level when there is a low water level in the dam. The flow will then be reduced linearly from full production flow for that time of year when water level is over the limit (//​limqprod//​ (percentage of volume)) to zero when no water is left in the dam over //wmin// level.+In addition, the production flow may be reduced when there is a low water level in the dam. The flow will then be reduced linearly from full production flow for that time of year when water level is over the limit (//​limqprod//​ (percentage of volume)) to zero when no water is left in the dam over //wmin// level.
  
 == Two rating curves == == Two rating curves ==
  
-Regulated outflow can also be calculated from different rating curves for the regulation periods. The rating curve parameters are the specific rate and exponent ​of is lacking the general rate and exponent. ​+Regulated outflow can also be calculated from different rating curves for the regulation periods. The rating curve parameters are the specific rate and exponent ​for the lake, or if lackingthe general rate and exponent ​parameters
  
 <m> q(t)=k*(w(t)-w_0)^p </m> <m> q(t)=k*(w(t)-w_0)^p </m>
Line 311: Line 407:
  
 <m> w0 = w0ref + deltaw0 </m> <m> w0 = w0ref + deltaw0 </m>
 +
 +The shift between regulation period 1 and 2 can be made smoother by allowing the threshold to increase/​decrease linearly over a few days. The number of days is given by the parameter //​w0adjdays//​.
  
 === Dams of specific purpose === === Dams of specific purpose ===
Line 316: Line 414:
 In addition to regulated lakes (described above), dams in HYPE can instead be regulated by a specific purpose. Dams are then governed by the main purpose of the dam, e.g. hydropower, flood control. The rules of calculating the outflow of the dam then depends on the specified purpose of each dam. The outlet lakes that are regulated by a specified purpose is given in DamData.txt. A dam without defined purpose may use the methods given in previous sections (and be defined in LakeData.txt). In addition to regulated lakes (described above), dams in HYPE can instead be regulated by a specific purpose. Dams are then governed by the main purpose of the dam, e.g. hydropower, flood control. The rules of calculating the outflow of the dam then depends on the specified purpose of each dam. The outlet lakes that are regulated by a specified purpose is given in DamData.txt. A dam without defined purpose may use the methods given in previous sections (and be defined in LakeData.txt).
  
-Dams are characterised by their regulation volume (//​regvol//​) between a minimum water stage (//wmin//) and the spill threshold (//w0//). The minimum water stage is not a real a water stage, but is calculated from the regulation volume (//​regvol//​) and the dam (surface) area. Thus it does not consider the area changing with depth. Below the minimum water level (//wmin//) no outflow occur.+Dams are characterised by their regulation volume (//​regvol//​) between a minimum water stage (//wmin//) and the spill threshold (//w0//). The minimum water stage (//​wmin//​) ​is not a real a water stage, but is calculated from the regulation volume (//​regvol//​) and the dam (surface) area. Thus it does not consider the area changing with depth. Below the minimum water level no outflow occurs. The spill threshold is defined from the lake's given depth (//lake_depth//)
  
-For dams a preferred production flow is calculated (see below). This is released if the dam in within the regulation volume, but no more water than to the minimum threshold ​is released. If instead the water level of the dam is above the threshold //w0//, the spillway equation (rating curve) is used. The outflow is the maximum of the spillway equation and the preferred production flow.+For dams a preferred production flow is calculated (see below). This is released if the dam in within the regulation volume, but no more water is released ​than to keep the water level above the minimum threshold. If instead the water level of the dam is above the threshold //w0//, the spillway equation (rating curve) is used. The outflow is the maximum of the spillway equation and the preferred production flow.
  
 Input variables //rate, exp, regvol, qamp// and //qpha// is given in DamData.txt. ​ Input variables //rate, exp, regvol, qamp// and //qpha// is given in DamData.txt. ​
Line 324: Line 422:
 == General preferred production flow == == General preferred production flow ==
  
-For all dams the general preferred production flow is calculated from input data on the dam. The flow that is preferred to be released from the dam today is then adjusted depending on purpose. The general production flow (//qprod//) may be constant or vary between two values per year. This is determined by the input variables //qprod1, qprod2, datum1// and //datum2//. Regulation period 1 between //datum1// and //datum2// has production flow //qprod1//, while the rest of the year has production flow //qprod2//. Not setting the dates gives the same production flow the whole year (//​qprod1//​). IF //qprod1// is not given //qprod// is estimated from the yearly average inflow of the dam, which is given as input data.+For all dams the general preferred production flow is calculated from input data on the dam. The flow that is preferred to be released from the dam today is then adjusted depending on purpose. The general production flow (//qprod//) may be constant or vary between two values per year. Typically it is estimated from the yearly average inflow of the dam, which is given as input data. It can alternatively be determined by input variables //qprod1, qprod2, datum1// and //datum2//, which then has to be estimated. Regulation period 1 between //datum1// and //datum2// has production flow //qprod1//, while the rest of the year has production flow //qprod2//. Not setting the dates gives the same production flow the whole year (//​qprod1//​). IF //qprod1// is not given //qprod// is estimated from the yearly average inflow of the dam, which is given as input data.
  
-The production flow may be reduced ​from this level when there is low water level in the dam. The flow will be reduced linearly from full production flow for that time of year when water level is over the limit (//​limqprod//​ (percentage of volume)) to zero when no water is left in the dam over //wmin// level.+The production flow may be reduced ​when the water level is low in the dam. The flow will be reduced linearly from full production flow for that time of year when water level is at the limit (//​limqprod//​ (percentage of volume)) to zero when no water is left in the dam over //wmin// level.
  
 == 1 - Irrigation dam and 2 - Water supply dam == == 1 - Irrigation dam and 2 - Water supply dam ==
Line 356: Line 454:
 == 4 - Hydropower dam == == 4 - Hydropower dam ==
  
-Dams for hydropower ​uses the general preferred production flow, but it may be adjusted with a seasonal factor. The factor ​is sinusoidal over the year, with a peak in December, when power output is normally high, and a minimum in June. This is done with the input variable //qamp//. If you want a different seasonal phase set //qpha// to another value (default = 102).+Dams for hydropower ​use the general preferred production flow, but the flow may be seasonally redistributed. The seasonal adjustment ​is sinusoidal ​factor ​over the year, with as default ​a peak in December, when power output is normally high (in northern Europe), and a minimum in June. For a model in another geographic area, a different seasonal phase may be needed, and is then set with the input //qpha// to another value (default = 102).  
 + 
 +The seasonal redistribution is used for dams/​subbasins with a snowfall fraction (//​snowfrac//​) above 0.35. For these dams the seasonal variation is calculated as:
  
 <m> qprod = qprod*(1 + qamp * sin(2*pi*(dayno+qpha)/​365)) </m> <m> qprod = qprod*(1 + qamp * sin(2*pi*(dayno+qpha)/​365)) </m>
  
 +If the snow fall fraction is above 0.35, a //qamp// of 0.71 is used, otherwise it is zero.  
 +The size of the seasonal adjustment can alternatively be set by the input variable //qamp// (0<//​qamp//<​1) for each dam. To turn off the seasonal redistribution the snow fraction can be set to zero for a specific dam.
  
 ==== Outlet lake with two outlets ==== ==== Outlet lake with two outlets ====
  
-An outlet lake that is not part of a lake composed of lakebasins, may have two defined outlets in [[start:​hype_file_reference:​lakedata.txt|LakeData.txt]]. The outlets ​can be divided into different outlet types depending on method for determining the outflow. HYPE will define an outlet'​s type based on variables present in [[start:​hype_file_reference:​lakedata.txt|LakeData.txt]]. Below the defining (necessary) variables are given for each outlet type.+An outlet lake that is not part of a lake composed of lakebasins, may have two defined outlets in [[start:​hype_file_reference:​lakedata.txt|LakeData.txt]]. The flow of the first outlet will go to the main downstream subbasin, and the second outlet will go to the subbasin defined in BranchData.txt. The outlets ​is divided into different outlet types depending on method for determining the outflow. HYPE will define an outlet'​s type based on variables present in [[start:​hype_file_reference:​lakedata.txt|LakeData.txt]]. Belowthe defining (necessary) variables are given for each outlet type. Note: The threshold (w0ref) is given for first outlet. For the second outlet this variable instead is given relative to the threshold of outlet 1 (i.e usually zero)
  
 === Outlet types === === Outlet types ===
Line 392: Line 493:
 Examples of common combinations of outlets: Examples of common combinations of outlets:
  
-- A hydropower plant with an estimated average production flow and a spill overflow at high water levels are simulated by combining production flow (outlet type 1) with rating curve determined overflow (outlet type 3).+- A hydropower plant with an estimated average production flow and a spill overflow at high water levels are simulated by combining production flow (outlet type 1) with rating curve determined overflow (outlet type 3). **Note** that different from a production flow and spill defined for one outlet (see Production flow above), the spill and produktion flow is here added together to form the total outflow. For one outlet the maximum of the two is used as total flow.
  
 - Same as above a hydropower plant with production flow in one branch and overflow in another (outlet type 3), but in this case the maximum allowed production is assumed higher than the average production (outlet type 2). This means that if overflow is calculated due to high water level, part of that flow is lead to the powerplant and through production to the first branch. - Same as above a hydropower plant with production flow in one branch and overflow in another (outlet type 3), but in this case the maximum allowed production is assumed higher than the average production (outlet type 2). This means that if overflow is calculated due to high water level, part of that flow is lead to the powerplant and through production to the first branch.
Line 404: Line 505:
  
  
-Many other combinations can be used, but not all outlet types are natural to combine. The following are not allowed:+Many other combinations can be used, but not all outlet types are natural to combine. The following are **not allowed**:
  
 ^ Outlet 1 ^ Outlet 2 ^ ^ Outlet 1 ^ Outlet 2 ^
Line 425: Line 526:
  
  
-==== Outlet lake (olake) ​as a lake basin ====+==== Outlet lake as a multi-basin lake composed of lake basins ​====
  
-An olake can be part of a larger lake. It is then called a lake basin of the larger lake. Lake basins ​are treated differently if they are located within the lake or if they are the last onei.e. the one with the lake outflow+An olake can be part of a larger lake. It is then called a lake basin of the larger ​lake, which is called a multi-basin ​lake. A multi-basin lake is assumed to have equal water level in all lake basins. The lake basins can have different depths below the outflow thresholdand they all use the same threshold when the water level output is calculated (same reference).
  
 === Inflow === === Inflow ===
  
-A lake basin can get local inflow from upstream areas and from the regional groundwater. Inflow works in the same way as for a simple olake (see above).+A lake basin can recieve ​local inflow, ​inflow from upstream areas (that are not part of the same multi-basin lake) and inflow ​from regional groundwater. Inflow works in the same way as for a simple olake (see above).
  
 === Flow between lake basins === === Flow between lake basins ===
  
-Flow through a lake divided into lake basins occurs ​only in one direction according to the //maindown// (and possibly ​//branchdown//), no two-way exchange between lake basins is implementedFor upstream ​lake basins ​that run into another lake basin, all the water above the threshold is drained. The threshold of the lake (same as last lake basin) is used, not any internal threshold between lake basins. If the lake water level (//​wlm//​) ​is higher than the threshold (//w0// or //wmin//) outflow is obtained.+All lake basins of a multi-basin lake are connected by the main downstream path. They can in addition have branches within the multi-basin lake, where a fraction of the flow then goes. Flow through a lake divided into lake basins ​can occurs in both direction according to the maindown (and possibly branchdown). ​The implementation prioritize flow from nearby ​lake basins ​and flow in the downstream direction when calculating ​the flow betweeen lake basins. The flow between ​the lake basins are calculate to reach an equal water level. No internal threshold between lake basins is considered when the flows are calculated
  
-=== Outflow ​of lake composed ​of lake basins ===+The flow out of a lakebasin (to another lakebasin) is determined in several steps: 
 +  - First the net outflow of each lakebasin is calculated. It is calculated as the difference between the current water level of the lakebasin and the water level of the whole lake after outflow from the multi-basin lake. First the average water level of the multi-basin ​lake is calculated, and from this the outflow(s) of the lake. Considering these outflows leaving the lake, a new average water level is calculated that the lake will have after all flows between lakebasins have been moved. This is the water level compared to the water level of the lakebasin in this step. 
 +  - The flow out of a lakebasin is the netflow plus all inflows from other lakebasins. Theses inflows are outflows of other lakebasins. Starting upstream in the lake the inflows are added to the outflow to give the flow out of every lakebasin. Note that these flows can be negative.  
 +  - The last step is to actually move the water. Preferably we want water to move from one lakebasin to the next, without mixing with water incoming from lakebasins further upstream. This we want so that substances do not travel all the way to the lake outlet on one time step. Therefore we try to remove all outflows from the lakebasins before we add the removed water into the downstream lakebasins. This might not be possible so we keep track of how much we can't move and repeat the third step if necessary.  ​
  
  
-For the last lake basin outflow is estimated by a rating curve, either a specific equation or with the general equation, or by regulation routine. Since all the water (above threshold) in the lake (from all lake basins) accumulates in the last lake basin the water level cannot be used directly, but it is recalculated to calculate the outflow from the lake. The last lake basin water level (//wlm//) is scaled by the subbasin lake area (//area//) to the whole lake area (//​lakearea//​). If the calculated lake water level (//wlake//) is higher than the threshold (//​wthresh//​) the outflow is estimated by the rating curve.+=== Outflow of lake composed of lake basins ​===
  
-  wlake = (wlm – wthresh* area / lakearea + wthresh +The main outflow of the lake is assumed to be located in the last lake basin. Here outflow is estimated by a rating curve, either a specific equation or with the general equation, or by a regulation routine. The threshold may be changing over the year as described in Section Two rating curves above (below Regulated lake). If the water level is below the threshold then the outflow ​is zero (or other minimum flow, //minflow//), or if the production flow regulation is used, there is production flow. See Section Production flow above. ​
-  ​outflow ​= rate * (wlake – wthresh)**exp+
  
-The threshold may be changing over the year as described in Section Two rating curves above. If the water level is below the threshold ​then the outflow is zero or if the production flow regulation is used there is production flow. See Section Production flow aboveProduction flow is also applied on the recalculated water stage.+In addition ​the lake can have outflows out of the multi-basin lake from other lakebasins. These additional outflows are then the flow in the branch of these lakebasins. The outflow ​in a branch going out of the multi-basin lake is specified in LakeData for that lakebasin. It can be determined by a rating curve or by production flow. If it is production flow it uses the same thresholds as those for the main outflow of the multi-basin lake, i.e//w0// and //wmin//
  
 ==== Initalisation of lake volume ==== ==== Initalisation of lake volume ====
  
-Unless a starting state is given from a file, the lakes start the simulation filled with water to their outflow threshold. That means for ilakes and simple olakes ​a water level equal to //​lake_depth//​. Dams are filled to the dam's maximum elevation (equal to //​lake_depth//​),​ except for flood control dams which are initializied with an empty regulation volume. ​+Unless a starting state is given from a file, the lakes start the simulation filled with water to their outflow threshold. That means for most lakes a water level equal to //​lake_depth//​. Dams are filled to the dam's maximum elevation (equal to //​lake_depth//​),​ except for flood control dams which are initializied with an empty regulation volume. ​ 
 + 
 +For lakes with outflow determined by a rating curve, the water level of the lake will often be higher than the outflow threshold level. The equilibrium level will depend on the size of the inflow and the outflow rating curve parameters. Depending on the residence time of water in the lake it may take time for this level to be established,​ and until then the outflow of the lake will be simulated lower than it should be. Thus a spin-up time is needed for a model simulation. For lakes in arid areas and for lakes with small catchments, the water level may be below the threshold and only occasionally have outflow. 
 + 
 +==== Constructing dams ==== 
 + 
 +It is possible to set at date from when a dam with regulation given in [[start:​hype_file_reference:​lakedata.txt|LakeData]] or [[start:​hype_file_reference:​damdata.txt|DamData]] is supposed to be in place and regulate the flow. The same way it is possible to give a date when the dam is supposed to be removed. For the period before the dam is build and/or after the dam is removed the outflow of the lake (that is not a dam anymore) is calculated from the rating curve parameters used for spill during regulation period.  
 + 
 +The removal of a dam is simulated by reducing the threshold of the lake (aka reducing the lake depth) to the ''​wmin''​-level. Similarly the lake depth is increased when a dam is built. The lake depth given in the input file is assumed to be the lake depth of the dam (as it is for outlet lakes without dates for building/​removing the dam). The lake depth is set to the lake_depth given in input file minus the the distance calculated from regulation volume divided by lake surface area for the period without a dam. When during a simulation the date of the building or removing of the dam is passed, the lake depth used in the simulation is changed. 
 + 
 +Removing of a filled dam may cause a flood wave travelling downstream. An option to have a smoother transition is to have an emptying phase with a gradually reducing threshold. This is implemented with a parameter specifying a number of days before the date the dam is removed that the threshold will start to decrease. 
 + 
 +If shorter time steps are used, the dam is assumed to be build or removed on the hour 00. 
 + 
 +==== Links to file reference ==== 
 + 
 +^Section ^Symbol ^Parameter/​Data ^File ^ 
 +|Local lake (ilake)| |//gldepi, illdepth, gratk, gratp, grata, ilratk, ilratp//​|[[start:​hype_file_reference:​par.txt|par.txt]]| 
 +|:::​|//​icatch//​|//​icatch//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]]| 
 +|:::|::: |or //gicatch, illicatch//​|[[start:​hype_file_reference:​par.txt|par.txt]]| 
 +|:::| |//lks_num, lks_dp_N, lks_fi_N, lks_fa_N//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]]| 
 +|Simple outlet lake or dam (olake)|//​lake<​sub>​depth</​sub>//​|//​lake_depth//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]] or [[start:​hype_file_reference:​lakedata.txt|LakeData.txt]] or [[start:​hype_file_reference:​damdata.txt|DamData.txt]]| 
 +|:::​|:::​|//​gldepo,​ olldepth//​|[[start:​hype_file_reference:​par.txt|par.txt]]| 
 +|:::| |//ratcorr, olratk, olratp, kthr, klow, krel//​|:::​| 
 +|:::​|//​wmin//​|calculated from //​regvol//​|[[start:​hype_file_reference:​lakedata.txt|LakeData.txt]] or [[start:​hype_file_reference:​damdata.txt|DamData.txt]]| 
 +|:::| |//rate, exp, qamp, qpha//​|:::​| 
 +|:::​|//​qprod//​|//​qprod1,​ qprod2, datum1, datum2//​|:::​| 
 +|:::| |//​limqprod//,//​w0adjdays//​|[[start:​hype_file_reference:​lakedata.txt|LakeData.txt]] or [[start:​hype_file_reference:​par.txt|par.txt]]| 
 +|:::| |//​deltaw0,​minflow//​|[[start:​hype_file_reference:​lakedata.txt|LakeData.txt]]| 
 +|:::| |//purpose, snowfrac//​|[[start:​hype_file_reference:​damdata.txt|DamData.txt]]| 
 +|:::​|//​Qinf//​|//​qinfjan,​ qinffeb, .., qinfdec//​|:::​| 
 +|Outlet lake with two outlets| |//​ldtype//​=5 or 6|[[start:​hype_file_reference:​lakedata.txt|LakeData.txt]]| 
 +|:::| |//w0ref, maxQprod, minflow, obsflow//​|:::​| 
 +|Constructing dams| |//​builddam,​ removedam// |[[start:​hype_file_reference:​lakedata.txt|LakeData.txt]] or [[start:​hype_file_reference:​damdata.txt|DamData.txt]]| 
 +|:::| |//​remdamdays//​|[[start:​hype_file_reference:​par.txt|par.txt]]|
  
-For lakes with outflow determined by a rating curve, the water level of the lake will be higher than the outflow threshold level. The equilibrium level will depend on the size of the inflow and the outflow rating curve parameters. Depending on the residence time of water in the lake it may take time for this level to be established,​ and until then the outflow of the lake will be simulated lower than it should be. Thus a spin-up time is needed for a model simulation. 
  
 ==== Links to relevant modules in the code ==== ==== Links to relevant modules in the code ====
Line 460: Line 596:
 |::: |average_flow_rating_curve|ouflow with rating curve| |::: |average_flow_rating_curve|ouflow with rating curve|
 |::: |calculate_ilake_outflow|local lake| |::: |calculate_ilake_outflow|local lake|
-|::: |calculate_outflow_from_outlet_lake|simple outlet lake| +|::: |calculate_ilakesection_outflow|:::|
-|::: |:::|outlet lake as a lakebasin| +
-|:::​|:::​|outlet lake with two outlets|+
 |:::​|set_general_rating_k|simple outlet lake| |:::​|set_general_rating_k|simple outlet lake|
-|:::|calculate_olake_waterstage|outlet lake as a lakebasin+|::: |calculate_outflow_from_outlet_lake calculate_branched_flow|simple ​outlet lake| 
-|:::|calculate_branched_flow|outlet lake with two outlets| +|:::|:::|outlet lake with two outlets| 
-|:::|calculate_branched_flow_new|:::|+|:::|calculate_lake_outlet_outflow|outlet lake with two outlets| 
 +|:::|calculate_lakebasin_average_waterstage|outlet lake as multi-basin lake| 
 +|:::​|calculate_outflow_from_lakebasin_lake|:::​| 
 +|:::​|calculate_branched_flow|:::​| 
 +|:::​|change_current_dam_status|constructing dams|
 |[[http://​hype.sourceforge.net/​doxy-html/​namespacemodelmodule.html|modelmodule (model_hype.f90)]] |initiate_model | local lake | |[[http://​hype.sourceforge.net/​doxy-html/​namespacemodelmodule.html|modelmodule (model_hype.f90)]] |initiate_model | local lake |
 +|::: | calculate_local_lake_model |::: |
 +|::: |calculate_flow_from_undivided_lake|simple outlet lake|
 +|:::​|:::​|outlet lake with two outlets|
 +|::: |calculate_flow_for_lakebasin_lake|outlet lake as multi-basin lake|
  
 ==== References ==== ==== References ====
  
-Lindström, G., 2016. Lake water level for calibration of the S-HYPE model, Hydrology Research 47.4:​672-682. ​doi: 10.2166/​nh.2016.019.+Ahmed, M.I., Shook, K., Pietroniro, A., Stadnyk, T., Pomeroy, J.W., Pers, C., Gustafsson, D., 2023. Implementing A Parsimonious Variable Contributing Area Algorithm for the Prairie Pothole Region in the HYPE Modelling Framework. Environmental Modelling & Software. https://​doi.org/​10.1016/​j.envsoft.2023.105769  
 + 
 +Lindström, G., 2016. Lake water level for calibration of the S-HYPE model, Hydrology Research 47.4:​672-682. ​https://doi.org/10.2166/​nh.2016.019 
 + 
 +Shook, K., Papalexiou, S., Pomeroy, J.W., 2021. Quantifying the effects of Prairie depressional storage complexes on drainage basin connectivity. J. Hydrol. 593, 125846. https://​doi.org/​10.1016/​j.jhydrol.2020.125846 
 ===== Floodplains ===== ===== Floodplains =====
  
Line 477: Line 624:
  
 |{{:​start:​hype_model_description:​FloodScenarios.png?​400|}}| |{{:​start:​hype_model_description:​FloodScenarios.png?​400|}}|
-| Figure ​7: Five possible scenarios of flooding.|+| Figure ​9: Five possible scenarios of flooding.|
  
 ==== Common floodplain processes ==== ==== Common floodplain processes ====
Line 490: Line 637:
 <m> a_max = fpf * classarea ​ </m> <m> a_max = fpf * classarea ​ </m>
  
-|{{:​start:​hype_model_description:​FloodVarFlow.png?​400|}}| +| {{:​start:​hype_model_description:​FloodVarFlow.png?​400}} ​    ​
-| Figure ​8: Illustration of floodplain variables of a lake.|+| Figure ​10: Illustration of floodplain variables of a lake.  |
  
 The exchange of water between lake (or river) and floodplain is determined by the respective water levels (//​wl<​sub>​body</​sub>//​ and //​wl<​sub>​plain</​sub>//​),​ the equilibrium water level (//​wl<​sub>​equil</​sub>//​) and the thresholds between the water bodies (//​fl<​sub>​body</​sub>/​fl<​sub>​plain</​sub>//​). ​ The exchange of water between lake (or river) and floodplain is determined by the respective water levels (//​wl<​sub>​body</​sub>//​ and //​wl<​sub>​plain</​sub>//​),​ the equilibrium water level (//​wl<​sub>​equil</​sub>//​) and the thresholds between the water bodies (//​fl<​sub>​body</​sub>/​fl<​sub>​plain</​sub>//​). ​
Line 514: Line 661:
  
 |{{:​start:​hype_model_description:​FloodVar2.png?​400|}}| |{{:​start:​hype_model_description:​FloodVar2.png?​400|}}|
-| Figure ​9: Illustration of floodplain soil flows.|+| Figure ​11: Illustration of floodplain soil flows.|
  
 ==== Connected floodplains ==== ==== Connected floodplains ====
  
 It is possible for flooded floodplains to limit the water flow from upstream rivers and lakes. This is determined based on the water levels in the floodplains,​ a higher water level of a downstream floodplain compared to a upstream floodplain may reduced the flow from the upstream compartment (main river or outlet lake) to the downstream compartment (main river or outlet lake). The damming flow is calculated as the flow to reach equilibrium water level between the two floodplains. The relative level of the two floodplains location is determined based om the classes average elevation and thresholds or given as indata separately. The function of connecting floodplains works within a subbasin or between connected subbasins. This floodplain features may be used together with the **Floodplain model with soil routines**. It is possible for flooded floodplains to limit the water flow from upstream rivers and lakes. This is determined based on the water levels in the floodplains,​ a higher water level of a downstream floodplain compared to a upstream floodplain may reduced the flow from the upstream compartment (main river or outlet lake) to the downstream compartment (main river or outlet lake). The damming flow is calculated as the flow to reach equilibrium water level between the two floodplains. The relative level of the two floodplains location is determined based om the classes average elevation and thresholds or given as indata separately. The function of connecting floodplains works within a subbasin or between connected subbasins. This floodplain features may be used together with the **Floodplain model with soil routines**.
 +
 +
 +==== Links to file reference ====
 +
 +^Section ^Symbol ^Parameter/​Data ^File ^
 +| |//​fpf//​|//​fpfol,​ fpfmr//​|[[start:​hype_file_reference:​flooddata.txt|FloodData.txt]]|
 +|Common floodplain processes|//​classarea//​|//​slc_nn,​ area//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]]|
 +|:::​|//​fym//​|//​fymol,​ fymmr//​|[[start:​hype_file_reference:​flooddata.txt|FloodData.txt]]|
 +|:::​|//​rc<​sub>​bodyToPlain</​sub>//​|//​rclfp,​ rcrfp//​|:::​|
 +|:::​|//​rc<​sub>​plainToBody</​sub>//​|//​rcfpl,​ rcfpr//​|:::​|
 +|:::​|//​fl<​sub>​body</​sub>//​|//​floll,​ flmrr//​|:::​|
 +|:::​|//​fl<​sub>​plain</​sub>//​|//​flolp,​ flmrp//​|:::​|
 +|Connected floodplains| |//hrefl, hrefr//​|:::​|
  
 ==== Links to relevant modules in the code ==== ==== Links to relevant modules in the code ====
Line 543: Line 703:
 ===== Bifurcations ===== ===== Bifurcations =====
  
-Outflow from a subbasin may flow in one or two directions. The main flow follows the main channel, which is the downstream subbasin ​described ​by the path given in [[start:​hype_file_reference:​geodata.txt|GeoData.txt]]. If there is a bifurcation,​ the branch flow goes to another downstream subbasin. Any of the flows may go outside the model set-ups area, they are then no longer a concern of the model. There are three ways to determine the flow in the different channels. 1) The division of the total outflow is determined in BranchData.txt. 2) The outflow is calculated for two outlets separately and then divided into the channels. For the second alternative see [[start:​hype_model_description:​hype_routing#​outlet_lake_with_two_outlets|Outlet lake with two outlets]] section above. 3) The demanded flow in the branch is prescibed. If the total flow is not enough to fulfil the need, less water goes into the branch. The main channel gets the rest of the total flow. For the first method the flow division into two channels ​by BranchData.txt ​is described ​by four parameters; //​mainpart//,​ //​maxQ<​sub>​main</​sub>//,​ //​minQ<​sub>​main</​sub>//​ and //​maxQ<​sub>​branch</​sub>//,​ which is set in [[start:​hype_file_reference:​branchdata.txt|BranchData.txt]]. Zero values of the parameters ​mean they are not used. The main flow (//​mainflow//​) is calculated from the totalflow (//q//) as: +Outflow from a subbasin may flow in one or two directions. The main flow follows the main channel, which is the downstream subbasin ​given by the path in [[start:​hype_file_reference:​geodata.txt|GeoData.txt]]. If there is a bifurcation,​ the branch flow goes to another downstream subbasin ​(with subid //​branchid//​). Any of the flows may go outside the model set-ups area, they are then no longer a concern of the model. ​ 
 + 
 +There are three ways to determine the flow in the different channels. 1) The division of the total outflow is determined in BranchData.txt. 2) The outflow is calculated for two outlets separately and then divided into the channels. For the second alternative see [[start:​hype_model_description:​hype_routing#​outlet_lake_with_two_outlets|Outlet lake with two outlets]] section above. 3) The demanded flow in the branch is prescribed (//​Qbranch//​). The prescribed flow can be constant and set in BarnchData.txt or variable and given in Xobs.txt (''​dwtr''​). If the total flow is not enough to fulfil the need, less water goes into the branch. The main channel gets the rest of the total flow.  
 + 
 +For the first method the total outflow has been calculated the ordinary way. The flows division into the two channels is the determined ​by four parameters; //​mainpart//,​ //​maxQ<​sub>​main</​sub>//,​ //​minQ<​sub>​main</​sub>//​ and //​maxQ<​sub>​branch</​sub>//,​ which are set in [[start:​hype_file_reference:​branchdata.txt|BranchData.txt]]. Zero values of the input data mean they are not used. The main flow (//​mainflow//​) is calculated from the totalflow (//q//) as: 
  
 <m> mainflow = delim{lbrace}{matrix{4}{2}{q {q<​=minQ_main} {mainpart*(q-minQ_main)+minQ_main} {minQ_main<​q<​=q_thresh} {maxQ_main} {q>​q_thresh \,\ q_thresh={{maxQ_main-minQ_main}/​{mainpart}}+minQ_main} ​ {q-maxQ_branch} ​ {q>​q_thresh\ ,\ q_thresh={{maxQ_branch}/​{1-mainpart}}+minQ_main}}}{} </m> <m> mainflow = delim{lbrace}{matrix{4}{2}{q {q<​=minQ_main} {mainpart*(q-minQ_main)+minQ_main} {minQ_main<​q<​=q_thresh} {maxQ_main} {q>​q_thresh \,\ q_thresh={{maxQ_main-minQ_main}/​{mainpart}}+minQ_main} ​ {q-maxQ_branch} ​ {q>​q_thresh\ ,\ q_thresh={{maxQ_branch}/​{1-mainpart}}+minQ_main}}}{} </m>
Line 558: Line 722:
  
 |{{:​start:​hype_model_description:​mainflow.png?​400|}}| |{{:​start:​hype_model_description:​mainflow.png?​400|}}|
-| Figure ​5: Main flow in relation to total flow for four different parameter combinations.|+| Figure ​12: Main flow in relation to total flow for four different parameter combinations.|
  
 |{{:​start:​hype_model_description:​branchflow.png?​400|}}| |{{:​start:​hype_model_description:​branchflow.png?​400|}}|
-| Figure ​6: Branch flow in relation to total flow for four different parameter combinations.|+| Figure ​13: Branch flow in relation to total flow for four different parameter combinations.| 
 + 
 + 
 + 
 +==== Links to file reference ====
  
 +^Symbol ^Parameter/​Data ^File ^
 +|//​branchid,​ Qbranch, mainpart//​|//​branchid,​ Qbranch, mainpart//​|[[start:​hype_file_reference:​branchdata.txt|BranchData.txt]]|
 +|//​maxQ<​sub>​main</​sub>,​ minQ<​sub>​main</​sub>,​ maxQ<​sub>​branch</​sub>//​|//​maxQmain,​ minQmain, maxQbranch//​|:::​|
  
 ==== Links to relevant modules in the code ==== ==== Links to relevant modules in the code ====
start/hype_model_description/hype_routing.1539869118.txt.gz · Last modified: 2023/11/16 14:28 (external edit)