[GRASS5] [bug #3163] (grass) g.region bug(s) with grid 3D

this bug's URL: http://intevation.de/rt/webrt?serial_num=3163
-------------------------------------------------------------------------

Subject: g.region bug(s) with grid 3D

Platform: GNU/Linux/i386
grass obtained from: CVS
grass binary for platform: Compiled from Sources

Hi,

there are two bugs within g.region or the underlying libgis code:

I want to set a region:

GRASS 6.1.cvs (spearfish60):~ > g.region n=4926000.00 s=4921000.00 e=597000.00 w=590000.00 res=100.00
GRASS 6.1.cvs (spearfish60):~ > g.region -p3
projection: 1 (UTM)
zone: 13
datum: nad27
ellipsoid: clark66
north: 4926000
south: 4921000
west: 590000
east: 597000
top: 1.00000000
bottom: 0.00000000
nsres: 100 <- ok
nsres3: 29.94011976
ewres: 100 <- ok
ewres3: 30.04291845
tbres: 1
rows: 50
rows3: 167
cols: 70
cols3: 233
depths: 1

Then I want to set the resolution of the 3rd dimension:
GRASS 6.1.cvs (spearfish60):~ > g.region res3=100.00
GRASS 6.1.cvs (spearfish60):~ > g.region -p3
projection: 1 (UTM)
zone: 13
datum: nad27
ellipsoid: clark66
north: 4928010 <- oops, was setting to default?!
south: 4913700 <- oops, was setting to default?!
west: 589980 <- oops, was setting to default?!
east: 609000 <- oops, was setting to default?!
top: 1.00000000
bottom: 0.00000000
nsres: 30 <- oops, was setting to default?!
nsres3: 30 <- oops, was setting to default?!
ewres: 30 <- oops, was setting to default?!
ewres3: 30 <- oops, was setting to default?!
tbres: 1
rows: 477
rows3: 477
cols: 634
cols3: 634
depths: 1

Combining print flag and 3d resolution setting:

GRASS 6.1.cvs (spearfish60):~ > g.region res3=100.00 -p3
ERROR: region for current mapset is invalid
       line 15: <depths: 0>
       run "g.region"

The WIND file contains:
cat grassdata/spearfish60/PERMANENT/WIND
proj: 1
zone: 13
north: 4928010
south: 4913700
east: 609000
west: 589980
cols: 634
rows: 477
e-w resol: 30
n-s resol: 30
top: 1
bottom: 0
cols3: 190
rows3: 143
depths: 0
e-w resol3: 100.10526316
n-s resol3: 100.06993007
t-b resol: inf <- no good

The bug is probably in lib/gis/rd_cellhd.c (?).

Markus

-------------------------------------------- Managed by Request Tracker