...
- Initialise temporary parameters.
Set each owner’s upstream inflow I(o) to equal that owner’s share of the upstream component’s outflow using equation (10). (The upstream component is the next division upstream, or if it is the first division in the link, the upstream node).
- Determine the volume of the total flow based flux, fFlowLG(), for this time step to be shared according to fixed ratio, expression or time series specified by the modeller. Where relevant, also determine the volume, LossHF, that is to be shared using high flow rules (i.e. Murray-style losses), which is the volume that is in excess of the value of fFlowLG().
- If the ownership system uses Murray-style losses:
The flux based on flow up to the high flow threshold is shared according to fixed ratio. This volume is:
Equation 57 The flux based on flow above the high flow threshold is shared according to high-flow rules. This volume is:
Equation 58
Otherwise, all the flux is shared using fixed shares, expression or time series specified by the modeller:
Equation 59
- If the ownership system uses Murray-style losses:
- Determine every owner’s fixed flux total for this time step (Lossfixed(o)), where the modeller has specified which of the link’s lateral fluxes (time series, general purpose flow based, groundwater and net evaporation) are fixed at the ownership system level (noting these must satisfy equation (11)).
- Time series flux:
- Determine the owner’s share of input time series flux(es) FluxTS(o):
If the time series flux for the current link is not input per owner, and for the ownership system time series flux sharing is by fixed ratio, apply the specified owner’s ratio to the total time series flux, FluxTS, that was calculated prior to ownership processing (FluxTS has been adjusted to ensure total loss does not exceed the amount of water in the division).
Equation 60 If the time series flux for the current link is input per owner, get the owner’s flux for this time step from the input time series:
If LossGainTS < 0 the flux is a gain, this needs no adjustment:Equation 61 Otherwise, the flux is a loss that may need to be scaled down (This is to allow for the case where the time series has pumped the division dry. The overall time series loss FluxTS has already been adjusted by link processing to ensure losses do not exceed water in the link. Input owner time series losses are scaled down so their sum does not exceed the total adjusted for any owner time series gains.).
Equation 62
Set the following
Equation 63
- Determine the owner’s share of input time series flux(es) FluxTS(o):
Determine the owner’s fixed share of the flow based flux when the ownership system uses Murray-style lossesor it uses fixed ratio to share other lateral fluxes:
Equation 64 Equation 65 If the ownership system uses fixed ratio to share other lateral fluxes: Apply the owner’s configured Ratioloss(o) to the remaining lateral fluxes and add their share of these to their fixed flux total:
Owner’s net evaporation flux:
Equation 66 Owner’s groundwater flux:
Equation 67 Equation 68
- Time series flux:
- Determine the total proportional flux, Lossprop.
- Initially Lossprop = 0
If time series flux sharing is ‘Proportional’, a total time series was input, so:
Equation 69 If Murray-style losses are not being modelled and other lateral flux sharing is ‘Proportional’:
Equation 70 If other lateral flux sharing is ‘Proportional’:
Equation 71
- Determine the division’s state (‘dead’ or ‘live’).
- If Muskingum weighting factor and there is no inflow (I = 0), it is ‘dead’:
State = Dead Otherwise, compare the result of the storage function to the mass balance equation to determine division state (if they are within maxError, the division is ‘live’):
Equation 72 Equation 73 then State = Live
Otherwise State = Dead
- If Muskingum weighting factor and there is no inflow (I = 0), it is ‘dead’:
- Find the division’s storage using a method appropriate to the state of the division (i.e. whether ‘dead’ or ‘live’). These methods are summarised in the two following sections.
Calculate each owner’s outflow, O(o), using equations (15) or (16), as appropriate, rearranged so the term O(o) is on the left hand side of the equation.
Equation 74 - Adjust owner shares of outflow as necessary to ensure that none are negative. The principles are discussed in Ownership Adjustments, above. In summary, if , the O(o) < 0, the owner concerned borrows from other owners to ensure O(o)=0.
- If any owner’s outflow is negative then:
Calculate surplus/deficit for each owner.
Equation 75 Equation 76 Pass owner surpluses and deficits to the Borrow method. This will return how much owners borrowed or lent (OwnerBorrowed(o), OwnerLent(o)).
Adjust the outflows to account for borrowing:
Equation 77 Equation 78
- If any owner’s outflow is negative then:
Calculate owner mass balances to report, based on equation (15) or (16), as appropriate, using adjusted values of O(o)=0.
...
Equation 79
Dead division
The steps involved in processing a dead division are as follows:
Data
Input data
Details on data requirements are provided in the Source User Guide.
...