[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