#!/bin/sh
#
# Written by Jaro Hofierka, for v.surf.rst
#
# this is a script for a cross-validation analysis of RST parameters
# OUTPUT: CSV table
#

# j - smoothing - the number after decimal point, e.g. sm=0.1 is defined as j = 10
# i - tension

INMAP=precip3d
ZCOL=precip

OUTFILESTATS=/tmp/data_cv.csv

######### nothing to change below

rm -f $OUTFILESTATS
echo "tension;smoothing;mean;population_stddev" > $OUTFILESTATS

j=10
while  [ $j -le 90 ]
do

  i=10
  while  [ $i -le 150 ]
  do
  
  TNS=`echo $i`
  SMTH=`echo $j`
  TNSSMTH=`echo "t"$i"s0"$j `
  echo "Computing tension/smoothing $TNSSMTH..."

   #interpolate sites CV differences:
   v.surf.rst -c input=$INMAP cvdev=data_cv_$TNSSMTH zcolumn=$ZCOL tension="$i" smooth=0."$j" segmax=700 --o

   #calculate univariate statistics for sites: 
   eval `v.univar -g data_cv_$TNSSMTH col=flt1 type=point | grep 'tension\|smoothing\|mean\|population_stddev'`
   echo "$TNS;$SMTH;$mean;$population_stddev" >> $OUTFILESTATS
   i=`expr $i + 10`
 done

 j=`expr $j + 10`
done
echo  "Finished. Written $OUTFILESTATS"