[GRASS-user] r.lake: understanding variable 'water_level'
Rich Shepard
rshepard at appl-ecosys.com
Thu Oct 6 15:59:35 PDT 2016
The elevation map (be_metric) has minimum and maximum elevations of 45.17
and 51.32 m. This command:
r.lake elev=be_metric water=50 lake=pop coord=2295821.12321,175257.829078
--overwrite
produces the output map with this commentary:
Lake depth from 0.000000 to 4.833561 (specified water level is taken as
zero)
Lake area 357082.807803 square meters
Lake volume 409786.922884 cubic meters
Volume is correct only if lake depth (terrain raster map) is in meters
My interest is in a specific, much smaller area. However, when I re-run
the command with much smaller water_level values there is no output:
r.lake elev=be_metric water=0.01 lake=pop coord=2295821.12321,175257.829078
--overwrite
100%
ERROR: Given water level at seed point is below earth surface. Increase
water level or move seed point.
r.lake elev=be_metric water=0.1 lake=pop coord=2295821.12321,175257.829078
--overwrite
100%
ERROR: Given water level at seed point is below earth surface. Increase
water level or move seed point.
If the entire 36 sq. km. is covered with water when water_level is set to
50 (just below maximum elevation), yet input values of 0.01 and 0.01 are
underground, how do I specify small elevations above that of the 'dam' which
is at an elevation of approximately 48.25m?
Rich
More information about the grass-user
mailing list