WIP: Stage/Level and Tidal Input mode

[dvalters]

Summary

Several users would like to be able to input river level data (Stage) and tidal levels into the model, instead of having to use discharge as an input source for reach mode style simulations.

Details

• This uses the same algorithm for levels/tidal input as was added in Caesar 1.9s and later versions.
• The user specifies a "bounding box" where the level data (water height) is applied to the model inputs
• The input file is a single column list of water heights, listed in order of time step for the stage/levels/tide data source
• The time step in minutes is set in the .params file.
• Levels are updated every timestep. It does this by setting the water depth directly from the input file, rather than having to calculate discharge from the discharge file, then applies that to the water_depths array.

Testing

This was tested on the Idealsied_stage.params test case in the test folder, and also some real data using the Beverley Ship Yard gauging station (levels) on the River Hull, UK, with a DEM covering the Hull from Beverley down to hull, run over 118 days.

Parameter File

# STAGE INPUTS
#==============

stage_mode_input_on:           yes
stage_input_timestep:          60
fromx:                         5
tox:                           10
fromy:                         272
toy:                           274
stage_inputfile:               stage_mode_input.csv

Note

Note: this still features the "bug" in #53, but the workaround is simple (treat x and y as i and j matrix indices rather than coordinates - a fix for that will be made available too shortly.