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

[Request Tracker]

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