I am trying to evaluate the impact of raising the downstream water level
has on the sediment transport in a reservoir. Essentially I have two identical runs with the only difference being the downstream water level.
Based on the MASS-BALANCE output, the two runs have similar BEDLOAD FLUX entering from the upstream boundary. When I integrated SOLID DISCHARGE along a segment near the upstream boundary using BLUE KENUE, however, I got two different discharge rates. Also, these two rates (from BLUE KENUE) are not quite the same as that from the MASS-BALANCE output. Does anyone know why this is the case?
Other information on the model setup:
1) conlit.f - I have implemented a sediment discharge rating curve in this subroutine. It calculates QBOR based on the total river inflow. I also skip the call to disimp.f to avoid QBOR being overridden by Q2BOR in the boundary file.
2) The boundary file has LIEBOR = KENT and LIQBOR = KSORT for the upstream boundary.
3) I have armored the upstream boundary nodes to prevent erosion.
Edwin