<div dir="ltr">Hi,<br><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Jan 3, 2016 at 2:23 PM, Ken Mankoff <span dir="ltr"><<a href="mailto:mankoff@gmail.com" target="_blank">mankoff@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hi List,<br>
<br>
I'm a grass newbie. I'm trying to route some water through a DEM and think r.terraflow is the right tool. It is crashing. Can someone point out what I'm doing wrong? I have not found much when searching for this problem on google or list archives. I think there are several issues occurring here. I'm having the same problems on grass64 (installs easily on OS X) and grass71 (on a linux machine).<br>
<br>
<br></blockquote><div>I don't have any experience with r.terraflow, but try using r.watershed. </div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
I'm working in the following region:<br>
<br>
> g.region -p<br>
> projection: 99 (Stereographic)<br>
> zone:       0<br>
> datum:      wgs84<br>
> ellipsoid:  wgs84<br>
> north:      -657600<br>
> south:      -3349500<br>
> west:       -638000<br>
> east:       864700<br>
> nsres:      150<br>
> ewres:      150<br>
> rows:       17946<br>
> cols:       10018<br>
> cells:      179783028<br>
><br>
<br>
1) r.mapcalc produces file of 0 size<br></blockquote><div><br></div><div>no idea, have you tried testing it on a smaller region? </div><div><br></div><div>Anna</div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<br>
When I run<br>
<br>
> r.mapcalc 'phi = (1000*9.8*bed + 917*9.8*(surf-bed))'<br>
<br>
The file has 0 size. I now do the same calculation using nctools and then load the netcdf file, but I'd rather do this in grass if I could. FYI, bed is 589 MB and surf is 190 MB. When I calculate phi with nctools (ncap2 -O -v -s "phi=1E3*9.8*bed+917*9.8*thickness" <a href="http://in.nc" rel="noreferrer" target="_blank">in.nc</a> <a href="http://phi.nc" rel="noreferrer" target="_blank">phi.nc</a>) the resulting file is 1.4 GB. Clearly the 0 size grass-produced phi is due to overflow int/float/double issues? But a) there is no warning/error and b) wrapping the variables above in float()'s doesn't seem to help. The file still has 0 size if I do a similar but simpler computation:<br>
<br>
> r.mapcalc 'phi = (bed + 0.917(surf-bed))'<br>
<br>
<br>
2) r.terraflow crashes with 'Bad address' error.<br>
<br>
> r.terraflow --o --v elevation=phi fill=f direction=d sw=s acc=a tci=tci<br>
<br>
The results of that command:<br>
<br>
> cell phi header compatible with region header<br>
> Elevation stored as FLOAT (4B)<br>
> WARNING: raster phi is of type CELL_TYPE --you should use r.terraflow.short<br>
> Region size is 17946 x 10018<br>
> STREAM temporary files in /home/user/data/grass/GR/PERMANENT/.tmp/host/26900.0  (THESE INTERMEDIATE STREAMS WILL NOT BE DELETED IN CASE OF ABNORMAL TERMINATION OF THE PROGRAM. TO SAVE SPACE PLEASE DELETE THESE FILES MANUALLY!)<br>
> Memory manager registering memory in MM_WARN_ON_MEMORY_EXCEEDED mode.<br>
> Reading data from <phi> to stream<br>
> </home/user/data/grass/GR/PERMANENT/.tmp/host/26900.0/STREAM_baaaaX><br>
>  100%<br>
> total elements=179783028, nodata elements=18914292<br>
> largest temporary files:<br>
> FILL: 9.38G (10067849568) [179783028 elements, 56B each]<br>
> FLOW: 11.99G (12869498880) [160868736 elements, 80B each]<br>
> Will need at least 23.97G (25738997760) space available in<br>
> /home/user/data/grass/GR/PERMANENT/.tmp/host/26900.0<br>
> MM warning: limit=314572800B. allocating 156610272B.  limit exceeded by 3706B.<br>
> EMPQUEUEADAPTIVE: starting in-memory pqueue<br>
> EMPQUEUEADAPTIVE: available memory: 297.937MB<br>
> EMPQUEUEADAPTIVE: desired memory: 297.937MB<br>
> sz_stream: 270400 buf_arity: 200 mm_overhead: 8666496 mm_avail: 312409214.<br>
> EMPQUEUEADAPTIVE: memory overhead set to 8.26501MB<br>
> EMPQUEUEADAPTIVE: pqsize set to 37967839<br>
> EMPQUEUEADAPTIVE: starting in-memory pqueue<br>
> EMPQUEUEADAPTIVE: available memory: 297.163MB<br>
> EMPQUEUEADAPTIVE: desired memory: 297.163MB<br>
> sz_stream: 270400 buf_arity: 200 mm_overhead: 8666496 mm_avail: 311597910.<br>
> EMPQUEUEADAPTIVE: memory overhead set to 8.26501MB<br>
> EMPQUEUEADAPTIVE: pqsize set to 37866426<br>
> EMPQUEUEADAPTIVE: starting in-memory pqueue<br>
> EMPQUEUEADAPTIVE: available memory: 296.389MB<br>
> EMPQUEUEADAPTIVE: desired memory: 296.389MB<br>
> sz_stream: 270400 buf_arity: 200 mm_overhead: 8666496 mm_avail: 310786606.<br>
> EMPQUEUEADAPTIVE: memory overhead set to 8.26501MB<br>
> EMPQUEUEADAPTIVE: pqsize set to 37765013<br>
> MM warning: limit=314572800B. allocating 155798940B.  limit exceeded by 3706B.<br>
> rusage/gettimeofday: Bad address<br>
<br>
<br>
<br>
Thanks for any help,<br>
<br>
  -k.<br>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a></blockquote></div><br></div></div>