Hi Marcus, On 2016-04-29 at 03:17, Markus Neteler <neteler at osgeo.org> wrote: > Would the eval() function of r.mapcalc be of use? In order to hold the > interim values? Just guessing, I did not study your need in details. Yes that does solve it. Thank you for pointing out this solution. -k.