[GRASSLIST:2398] Re: The g3.region settings ....

Geotronix mail-list at geotux.it
Thu Jan 29 11:39:40 EST 2004


-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

* mercoledì 28 gennaio 2004, alle 15:43, Glynn Clements scrive:
> 
> Geotronix wrote:
> 
> > I'm trying to get s.vol.rst to work for the computation of rainfall map
> > using a 50m resolution cell (1000x1000 cells grid).
> > Now, the site file that's input for s.vol.rst has 4 columns
> > x|y|z|%rainfall(mm/year)
> > 
> > When I tried to launch s.vol.rst the first time it told me that I had to
> > set the 3D region. SO I launched g3.createregion which reead the
> > defaults from the current 2D region. I gave him the following options:
> > 
> > g3.createregion t=2500 b=0 dres=1
> > 
> > as to resemble:
> > 
> > t=2500 --> max altitude in meters that I have on my dem;
> > b=0  --> min altitude on my dem ;
> > dres=1 -->vertical resolution on 1 meter.
> > 
> > Now, I used the maskmap option for s.vol.rst so as to narrow the
> > computation only to a subarea of the 1000x1000 cell area while all
> > rainfall station (25) fall in the 1000x1000 cell area.
> > 
> > I launched the program and it is takeing a hell of a long time to
> > compute. After 4 hours it's still at 0%.
> > 
> > Did I miss something in the settings of the 3D region or does it simply
> > take so much time to compute? Yet, 1000x1000 cells don't seem to be much
> 
> It's a volume, not an area: 1000x1000x2500 = 2.5x10^9 cells.
> 

Do you mean that 2500 I set is in cells? I thought it was in meters.....
The 1000x1000 is a grid with 50 meter resolution. So to set z=2500 must
I give it 2500/50=50? So that I get 1000x1000x50?
Or maybe the vertical scale resolution is indipendant from the x,y ones?

Cheers
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.2.2 (GNU/Linux)

iD8DBQFAGTdMIRaQQEznKxURAhRrAKCD3eyyP34TLBvf71Rcgsis/lylTQCg0DGf
Dy/0a76lLHsvyoFZc93HxK4=
=Cd3m
-----END PGP SIGNATURE-----




More information about the grass-user mailing list