<div dir="ltr">It seems the easiest way to speed this up is to parallelize it, which can be done fairly easily with GNU parallel. <div><a href="https://grasswiki.osgeo.org/wiki/Parallelizing_Scripts#GNU_Parallel">https://grasswiki.osgeo.org/wiki/Parallelizing_Scripts#GNU_Parallel</a><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Feb 11, 2017 at 8:12 AM, 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi GRASS list,<br>
<br>
I have a watershed raster with 200 million cells. I calculate the basins and the baseline accumulation (i.e. flow=1) and have done so with r.watershed like this:<br>
<br>
r.watershed -a elevation=e drainage=d accumulation=a basins=b threshold=50000 --q<br>
<br>
That takes ~2 minutes.<br>
<br>
I now need to calculate surface flow ~1000 times over this same watershed (every day for 3 years). That would take 1.5 days which is not unreasonable. But is there a faster way? Can I use the pre-computed basins and accumulation rasters? Is it faster to loop over the basins (there are ~1000 of them)? Or do I just run the entire command as above 1000 times but with variables for flow= and accumulation=?<br>
<br>
Thanks,<br>
<br>
  -k.<br>
</blockquote></div><br></div>