<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p
{mso-style-priority:99;
margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
span.E-MailFormatvorlage18
{mso-style-type:personal-compose;
font-family:"Tahoma","sans-serif";
color:black;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 2.0cm 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black">Dear GRASS Gurus,<br>
<br>
During the last days I've discovered a strange behaviour of r.watershed and D8. The used DEM is the result from a landscape evolution model and there are no depressions. I cannot use MFD as the D8-FlowDirections input data of a variety of tools I’ve written
during the last years.<br>
<br>
I got obviously wrong values for the flow accumulation of r.watershed and D8. The number of cells in the Flow Accumulation grid at the outlet point was significantly smaller than the area of catchment analysed catchment (spatial resolution in north south and
east west direction is 1) also computed with GRASS using r.water.outlet and a conversion to an Area. I checked Flow Accumulation at the the outlet with an alternative code of a colleague and the contributing drainage area was exactly the value
<br>
of the catchment size computed with r.water.outlet. <br>
<br>
I further found some deviations in the flow direction grid. In some cases the flow routing leads to the lowest neighbouring cell which however might not be the steepest decent. r.watershed chooses sometimes a diagonal flow direction while the the local channel
gradient would be greater in horizontal direction. <br>
<br>
I started some test runs with very simple geometries. I've defined a ramp with r.mapcalc that dips in north direction and I expected that the this will also be shown by the flow direction of r.watershed. However, this is not the case. The synthetic DEM features
a watershed!!! See attachment. It seems that r.watershed fills the region (without need). This strange behavior disappears in by applying the multiple flow direction algorithm.
<o:p></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black"><br>
<br>
I have tested with the latest versions grass70 , grass71 and grass645.<o:p></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black">Here is a test case (grass71). I would be great if somebody could reproduce the problem and show me a workaround.<br>
<br>
<o:p></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black">---------------------------------------------------------<o:p></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black">g.region w=0 e=1000 s=0 n=1000 res=1 -p<br>
r.mapcalc "DEM = 1000 - y()"<br>
r.watershed -s elevation=DEM@PERMANENT accumulation=ACC drainage=DIR<br>
<br>
<br>
Please help<br>
<br>
Jörg<br>
<br>
<br>
<br>
<br>
<br>
<br>
<o:p></o:p></span></p>
</div>
</body>
</html>