#!/bin/csh -f


set measured = ba_104_24_04_97.refl.dat

gawk -f filt.awk < $measured > $measured.plot

\rm -f lut.dat

foreach c (10 30 60 90)
\rm -f results.$c.dat
foreach soil (0.05 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8)
foreach lai (0.001 0.25 0.5 0.75 1.0 1.5 2.0 3.0 4.0 6.0 8.0 10.0 12.0 15.0 20.0)
	echo $soil $lai $c
	if( ! -e xx.$lai.$soil.$c ) then
		kuusk < $measured -cAb $c -l $lai -rsl1 $soil | gawk -f filt.awk > xx.$lai.$soil.$c
	endif

	set rmse = `paste xx.$lai.$soil.$c $measured.plot | gawk '{r1=$2;r2=$4;e=(r1-r2)/r1;sum += e*e;} END {print sqrt(sum/(NR*1.0))}'`

	set lai2 = `echo "e(-$lai / 2)" | bc -l`

	echo $lai2 $soil $rmse >> results.$c.dat
	echo $lai $soil $rmse $c >> lut.dat
end
echo "" >> results.$c.dat
end
end