[GRASS-dev] bug: g.region vect=3Dmap breaks region

[Hamish]

[trac is non-responsive right now so posting here]


Hi,

In a lat/lon location if you import a lon,lat,z file:
{{{
G65> awk '{print $1 "|" $2 "|" $3}' fort.100 | \
       v.in.ascii out=fort_100 x=1 y=2 z=3 -t -z
}}}

then set the region to match the vector map's bounds:
{{{
G65> g.region v=fort_100
}}}

the region will be broken:
{{{
G65> g.region -p
ERROR: region for current mapset is invalid
       line 18: <t-b resol: 1:00:00.0196>
       run "g.region"
}}}

removing this from the end of $MAPSET/WIND gets you operational again:
{{{
top:        4560.63
bottom:     -7.39548
cols3:      8
rows3:      9
depths:     4568
e-w resol3: 1:02:32.898867
n-s resol3: 1:03:00.119988
t-b resol:  1:00:00.0196
}}}


t-b resolution (elevation) obviously shouldn't be in degrees..
?


thanks,
Hamish