s.surf.tps for only 60 sites?

[M.R. Lassche]

Fred Watson wrote

> Hi,
> 
> I'm trying to perform spline interpolation over about sixty rainfall stations
> to construct a rainfall map using  s.surf.tps. I get a raster map with only
> one category. I'm generally using default parameters.
> 
> Any ideas on what the problem is? i.e. why I only get one category.
> 
> Thanks in advance,
> 
> P.S: Whilst s.surf.tps is running, I also get a stack of warnings about
> site descriptions not beginning with '#'. I assume my site file is in order
> because s.surf.idw works fine. There is also a warning about 'minimum number of
> sites'.
> 
> Fred Watson
> PhD student
> Cooperative Research Centre for Catchhment Hydrology
> Civil & Env. Eng. Dept.
> Univ. of Melb.
> Parkville 3052 Australia
> 

Fred the message about the "#" is correct for using s.surf.tps all your sites descritions should be preceded with a "#", the s.surf.tps will word fine
Use awk for inserting the "#"
awk '{print $1"|"$2"|#"$3}' <infile >outfile
 This will do the job

Ronnie Lassche