# create a 3D raster g.region s=0 n=100 w=0 e=100 b=0 t=100 -p r3.mapcalc "test_gradient = sqrt(row()*row() +col()*col()+ depth()*depth())" # compute gradient r3.gradient input=test_gradient output=grad_x,grad_y,grad_z