Deriving values from a raster map

Is there anyway to determine the value of a number of raster layers at a
position defined by a sites file ?

The example our Hydrology want to do is :-
they have a sites file containing their tipping bucket raingauges and we
have various raster data sets, such as potential evaporation, standard
average annual rainfall, etc. They would like to be able to 'simply' get
the value of each raster at each raingauge site

I've shown them 'd.what.rast map=x,y,....' and they think that's great -
I just feel that there should be some way to do it in a more elegant
style.

Thanks for any help,

Steve Culshaw
NRA North West
e-mail : nra@cix.compulink.co.uk (use this for general usage)
or sculshaw@cix.compulink.co.uk (Private - only checked at best
biweekly)

N R A (nra@cix.compulink.co.uk ) writes on 1 Sep 94:

Is there anyway to determine the value of a number of raster layers at a
position defined by a sites file ?

yes, s.sample, available from
ftp://pasture.ecn.purdue.edu/pub/mccauley/grass/s.sample.tar.gz

does nearest neighbor, bilinear interpolation, and cubic convolution
sampling.

--Darrell