GRASS logo

NAME

r.topmodel - Simulates TOPMODEL which is a physically based hydrologic model.

KEYWORDS

raster, hydrology

SYNOPSIS

r.topmodel
r.topmodel help
r.topmodel [-i] [elevation=name] [basin=name] [depressionless=name] [direction=name] [basin_elevation=name] [topidx=name] [nidxclass=integer] idxstats=name parameters=name input=name output=name [qobs=name] [timestep=integer] [idxclass=integer] [--overwrite] [--verbose] [--quiet]

Flags:

-i
Input data given for (o/i)
--overwrite
Allow output files to overwrite existing files
--verbose
Verbose module output
--quiet
Quiet module output

Parameters:

elevation=name
Name of input elevation raster map
basin=name
Name of input basin raster map
Created by r.water.outlet (MASK)
depressionless=name
Name for output depressionless elevation raster map
direction=name
Name for output flow direction map for depressionless elevation raster map
basin_elevation=name
Name for output basin elevation raster map (o/i)
MASK applied
topidx=name
Name for output opographic index ln(a/tanB) raster map
MASK applied
nidxclass=integer
Number of topographic index classes
Default: 30
idxstats=name
Name of topographic index statistics file (o/i)
parameters=name
Name of TOPMODEL parameters file
input=name
Name of rainfall and potential evapotranspiration data file
output=name
Name for output file
qobs=name
Name for observed flow file
timestep=integer
Time step
idxclass=integer
Topographic index class

DESCRIPTION

r.topmodel simulates TOPMODEL which is a physically based hydrologic model.

Note: (i) means input; (o) means output; (o/i) means input or output

The -i flag indicates that input data are given for (o/i). Without this flag, all inputs (i) and intermediate outputs (o/i) should be given. For example, [belevation] map will be created from [elevation] and [basin] in every run. However, given the same [elevation] and [basin], [belevation] output will be the same all the time, so r.topmodel can directly take [belevation] as an input with this flag to save time.

Selected Parameters:

depressionless map is created as follows:
r.fill.dir input=elevation elev=depressionless dir=direction type=grass
This option can be omitted if [elevation] map is already depressionless.
belevation map is created from [elevation] with [basin] mask applied:
r.mapcalc "belevation = if(basin == 0 || isnull(basin), null(), elevation)"
topidx map is created as follows:
r.topidx input=elevation output=topidx
qobs
Compare simulated flows with observed flows and calculate model efficiency.

REFERENCES

K. Beven, R. Lamb, P. Quinn, R. Romanowicz, and J. Freer. TOPMODEL, in V.P. Singh (Ed.). Computer Models of Watershed Hydrology. Water Resources Publications, 1995.

S.C. Liaw, Streamflow simulation using a physically based hydrologic model in humid forested watersheds (Dissertation). Colorado State University, CO. p163, 1988.

SEE ALSO

r.fill.dir, r.mapcalc, r.topidx
How to run r.topmodel

AUTHORS

Main algorithm sources are rewritten in C based on TMOD9502.FOR.
Thanks to Keith Beven.

GRASS port by Huidae Cho
Hydro Laboratory, Kyungpook National University, South Korea

Last changed: $Date: 2009-08-16 02:28:19 -0700 (Sun, 16 Aug 2009) $


Main index - raster index - Full index

© 2003-2011 GRASS Development Team