Hi Christian, thanks for fast reply! I also tried r.mapcalc with eval(), setting variables for dir_x/dir_y. But map[ , ] does not take them with the same error. Looks like I have to find another way to fasten the if-procedure >> in order to make a r.water.outlet on many many areas at the same time. Cheers, Achim