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

[Markus Neteler]

On Wed, Jan 28, 2004 at 01:40:06PM +0100, Geotronix wrote:
> -----BEGIN PGP SIGNED MESSAGE-----
> Hash: SHA1
> 
> Hi,
> 
> 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.

It's probably much better, if you eventually just want to extract a 2D
rainfall map, to run the computation in a single-layer volume:

g3.createwind t=2500 b=0 dres=2500

 t=2500 --> max altitude in meters that I have on my dem;
 b=0  --> min altitude on my dem ;
 dres=2500 -->vertical resolution 2500 meter.

Then s.vol.rst ...

Please report if this is faster.

 Markus Neteler