Model:TopoFlow-Channels-Diffusive Wave: Difference between revisions
No edit summary |
No edit summary |
||
Line 14: | Line 14: | ||
{{Model identity | {{Model identity | ||
|Model type=Single | |Model type=Single | ||
|Categories=Hydrology | |Categories=Hydrology, Terrestrial | ||
|Spatial dimensions=2D | |Spatial dimensions=2D | ||
|Spatialscale=Landscape-Scale, Watershed-Scale | |Spatialscale=Landscape-Scale, Watershed-Scale | ||
|One-line model description=Diffusive Wave process component for a D8-based, spatial hydrologic model | |One-line model description=Diffusive Wave process component for a D8-based, spatial hydrologic model | ||
|Extended model description=This process component is part of a spatially-distributed hydrologic model called TopoFlow, but it can now be used as a stand-alone model. It uses the "diffusive wave" method to compute flow velocities for all of the channels in a D8-based river network. This method includes a pressure gradient term that is induced by a water-depth gradient in the downstream direction. This means that instead of using bed slope in Manning's equation or the law of the wall, the water-surface slope is used. | |Extended model description=This process component is part of a spatially-distributed hydrologic model called TopoFlow, but it can now be used as a stand-alone model. It uses the "diffusive wave" method to compute flow velocities for all of the channels in a D8-based river network. This method includes a pressure gradient term that is induced by a water-depth gradient in the downstream direction. This means that instead of using bed slope in Manning's equation or the law of the wall, the water-surface slope is used. | ||
}} | }} | ||
{{Model technical information | {{Model technical information | ||
Line 27: | Line 27: | ||
|Start year development=2001 | |Start year development=2001 | ||
|Does model development still take place?=Yes | |Does model development still take place?=Yes | ||
|Model availability=As code | |Model availability=As code, As teaching tool | ||
|Source code availability=Through CSDMS repository | |Source code availability=Through CSDMS repository | ||
|Program license type=Apache public license | |Program license type=Apache public license | ||
Line 37: | Line 37: | ||
}} | }} | ||
{{Input - Output description | {{Input - Output description | ||
|Describe input parameters= | |Describe input parameters=These inputs must be provided as grids: | ||
These inputs must be provided as grids: | |||
*flow_codes = D8 flow codes (Jenson 1984 convention), | *flow_codes = D8 flow codes (Jenson 1984 convention), | ||
* [NE,E,SE,S,SW,W,NW,N] → [1,2,4,8,16,32,64,128] | * [NE,E,SE,S,SW,W,NW,N] → [1,2,4,8,16,32,64,128] | ||
Line 66: | Line 65: | ||
*Sinuosity: Scalar 1.00000000 [m/m] | *Sinuosity: Scalar 1.00000000 [m/m] | ||
|Input format=ASCII, Binary | |Input format=ASCII, Binary | ||
|Describe output parameters= | |Describe output parameters=This component computes the following variables, as grids: | ||
This component computes the following variables, as grids: | *discharge, Q, [m^3/s]; | ||
discharge, Q, [m^3/s]; | *flow velocity, u, [m/s]; | ||
flow velocity, u, [m/s]; | *flow depth, d, [m]; | ||
flow depth, d, [m]; | *friction factor, f, [none]; | ||
friction factor, f, [none]; | *free-surface slope, S_free, [m/m]. | ||
free-surface slope, S_free, [m/m]. | |||
The user can choose which, if any, of these to save. Each may be saved as a grid sequence, indexed by time, in a netCDF file, at a specified sampling rate. Each may also be saved for a set of "monitored" grid cells, each specified as a (row,column) pair in a file with the name "<case_prefix>_outlets.txt". With this option, computed values are saved in a multi-column text file at a specified sampling rate. Each column in this file corresponds to a time series of values for a particular grid cell. For both options the sampling rate must no smaller than the process timestep. | The user can choose which, if any, of these to save. Each may be saved as a grid sequence, indexed by time, in a netCDF file, at a specified sampling rate. Each may also be saved for a set of "monitored" grid cells, each specified as a (row,column) pair in a file with the name "<case_prefix>_outlets.txt". With this option, computed values are saved in a multi-column text file at a specified sampling rate. Each column in this file corresponds to a time series of values for a particular grid cell. For both options the sampling rate must no smaller than the process timestep. | ||
|Output format=ASCII, Binary | |Output format=ASCII, Binary | ||
|Pre-processing software needed?=Yes | |Pre-processing software needed?=Yes | ||
Line 81: | Line 78: | ||
|Describe post-processing software=None, except visualization software. Grid sequences saved in netCDF files can be viewed as animations and saved as movies using VisIt. | |Describe post-processing software=None, except visualization software. Grid sequences saved in netCDF files can be viewed as animations and saved as movies using VisIt. | ||
|Visualization software needed?=Yes | |Visualization software needed?=Yes | ||
|Other visualization software=VisIt | |||
}} | }} | ||
{{Process description model | {{Process description model | ||
|Describe processes represented by the model=The diffusive wave method for flow routing in the channels of a D8-based river network. | |Describe processes represented by the model=The diffusive wave method for flow routing in the channels of a D8-based river network. | ||
|Describe key physical parameters and equations= (Source: TopoFlow HTML Help System) | |Describe key physical parameters and equations=(Source: TopoFlow HTML Help System) | ||
------------------------------------------------------------ | ------------------------------------------------------------ | ||
Line 97: | Line 95: | ||
Pw = w + [2 * d / cos(θ)] = wetted perimeter of a trapezoid [m] | Pw = w + [2 * d / cos(θ)] = wetted perimeter of a trapezoid [m] | ||
Vw = d2 * [ L * tan(θ) ] + d * [L * w] = wetted volume of a trapezoidal channel [m] | Vw = d2 * [ L * tan(θ) ] + d * [L * w] = wetted volume of a trapezoidal channel [m] | ||
|Describe length scale and resolution constraints=Recommended grid cell size is around 100 meters, but can be parameterized to run with a wide range of grid cell sizes. DEM grid dimensions are typically less than 1000 columns by 1000 rows. | |Describe length scale and resolution constraints=Recommended grid cell size is around 100 meters, but can be parameterized to run with a wide range of grid cell sizes. DEM grid dimensions are typically less than 1000 columns by 1000 rows. | ||
|Describe time scale and resolution constraints=The basic stability condition is: dt < (dx / u_min), where dt is the timestp, dx is the grid cell size and u_min is the smallest velocity in the grid. This ensures that flow cannot cross a grid cell in less than one time step. Typical timesteps are on the order of seconds to minutes. Model can be run for a full year or longer, if necessary. | |Describe time scale and resolution constraints=The basic stability condition is: dt < (dx / u_min), where dt is the timestp, dx is the grid cell size and u_min is the smallest velocity in the grid. This ensures that flow cannot cross a grid cell in less than one time step. Typical timesteps are on the order of seconds to minutes. Model can be run for a full year or longer, if necessary. | ||
|Describe any numerical limitations and issues=This model/component needs more rigorous testing. | |Describe any numerical limitations and issues=This model/component needs more rigorous testing. | ||
}} | }} | ||
{{Model testing | {{Model testing | ||
|Describe available calibration data sets=This model/component is typically not calibrated to fit data, but is run with a best guess or measured value for each input parameter. | |Describe available calibration data sets=This model/component is typically not calibrated to fit data, but is run with a best guess or measured value for each input parameter. | ||
|Describe available test data sets=Available test data sets: | |Describe available test data sets=Available test data sets: | ||
Line 109: | Line 107: | ||
* Small basin in Kentucky. | * Small basin in Kentucky. | ||
* Inclined plane for testing. | * Inclined plane for testing. | ||
* Arctic watershed data from Larry Hinzman (UAF). | * Arctic watershed data from Larry Hinzman (UAF). | ||
|Describe ideal data for testing=Several test datasets are stored on the CSDMS cluster at: /data/progs/topoflow/3.0/data. | |Describe ideal data for testing=Several test datasets are stored on the CSDMS cluster at: /data/progs/topoflow/3.0/data. | ||
}} | }} | ||
{{Users groups model | {{Users groups model | ||
|Do you have current or future plans for collaborating with other researchers?=Collaborators include: Larry Hinzman (UAF), Bob Bolton, Anna Liljedahl (UAF), Stefan Pohl, Tom Over and others | |Do you have current or future plans for collaborating with other researchers?=Collaborators include: Larry Hinzman (UAF), Bob Bolton, Anna Liljedahl (UAF), Stefan Pohl, Tom Over and others | ||
}} | }} | ||
{{Documentation model | {{Documentation model | ||
Line 121: | Line 119: | ||
}} | }} | ||
{{Additional comments model | {{Additional comments model | ||
|Comments=This component was developed as part of the TopoFlow hydrologic model | |Comments= This component was developed as part of the TopoFlow hydrologic model. For more information on TopoFlow, please goto http://csdms.colorado.edu/wiki/Model:TopoFlow. | ||
When used from within the CSDMS Modeling Tool (CMT), this component has a "config" button which launches a graphical user interface (GUI) for changing input parameters. The GUI is a tabbed dialog with a Help button at the bottom that displays HTML help in a browser window. | |||
This component also has a configuration (CFG) file, with a name of the form: <case_prefix>_channels_diff_wave.cfg. This file can be edited with a text editor. | |||
The Numerical Python module (numpy) is used for fast, array-based processing. | |||
}} | }} | ||
<!-- PLEASE USE THE "EDIT WITH FORM" BUTTON TO EDIT ABOVE CONTENTS; CONTINUE TO EDIT BELOW THIS LINE --> | <!-- PLEASE USE THE "EDIT WITH FORM" BUTTON TO EDIT ABOVE CONTENTS; CONTINUE TO EDIT BELOW THIS LINE --> |
Revision as of 00:30, 16 February 2010
Contact
Name | Scott Peckham |
Type of contact | Model developer |
Institute / Organization | CSDMS, INSTAAR, University of Colorado |
Postal address 1 | 1560 30th street |
Postal address 2 | |
Town / City | Boulder |
Postal code | 80305 |
State | Colorado |
Country | USA"USA" is not in the list (Afghanistan, Albania, Algeria, Andorra, Angola, Antigua and Barbuda, Argentina, Armenia, Australia, Austria, ...) of allowed values for the "Country" property. |
Email address | Scott.Peckham@colorado.edu |
Phone | 303-492-6752 |
Fax |
TopoFlow-Channels-Diffusive Wave
Metadata
Summary
Technical specs
In/Output
Process
Testing
Other
IntroductionHistoryPapersIssuesHelpInput FilesOutput FilesDownloadSource |