r3.mapcalc depth() and z() function upside down

[ayepri]

# computation of a 3d raster with z values from the z() function r3.mapcalc Z='z()' # export to VTK format (see attached image grid_inversed_z.png) r3.out.vtk Z output=./Z_3dg.vtk

In r3.mapcalc, z() function returns -z, and depth() function returns the elevation of the cell. On the grid_inversed_z.png (viewed in Paraview), one can note that higher Z values are at the bottom of the grid (mind the axis at the bottom left). mputed z() isosurfaces are not horizontal (see grid_inversed_deformed_z)

[ayepri]

volume filled with z() values from r3.mapcalc and visualized with Paraview