#!/bin/csh

alias rt2 /data/leiden/scat_software/rt2

set M = "0.4 0.2 0.05"				# vegetation moisture
set K = "26.2 111 209 9.0"			# wave numbers = 2*pi / wavelength
set G = "0.9 0.7 0.4"				# soil moisture
set H = "0.05 0.2 0.7 1"			# height
set D = "200 300 400"				# nden 
set R = "0.04 0.09 0.14"			# rad1 
set W = "0.03 0.07 0.1"				# rad2 

foreach k ($K)
\rm -f output.$k logfile.$k
 foreach m ($M)
  foreach h ($H) 
   foreach g ($G) 
    foreach d ($D) 
     foreach r ($R) 
      foreach w ($W) 
 	 echo $k $m $h $g $d $r $w

	gawk -v wave=$k -v mv=$m -v mg=$g -v he=$h -v n=$d -v r1=$r -v r2=$w -f SUGBEET.awk > xx.rt
	rt2 xx.rt
	echo $k $m $g $h $d $r $w >> logfile.$k
       	tail -1 xx.out >> output.$k
        \rm xx*
     
      end
     end
    end   
   end
  end
 end
 paste logfile.$k output.$k > log-and-output.$k
end