#!/bin/sh -f
# shell skript

libradtran_path='/local/libRadtran-1.5-beta/bin/'
tmp_path='/home/users/merk/radtran_2010/rte/tmp/'
script_path='/home/users/merk/radtran_2010/rte/'


# create files for different g and l
echo 'create files for different g and l'

cd $script_path


# vary g and l to find a graphical solution

for g in 0.9 #0.7 0.8 0.85 0.9
   do echo 'g: ' $g
   for moments in 70 71 72 73 74 75 #10 20 50 70
	do python rte_heneygreenstein.py $moments $g 		# call python-program (param1=l, param2=g)
	echo '\n\n'
   done
done

###echo 'done. files created. \n \n'

echo 'run libradtran -> phase'
cd $libradtran_path

for g in 0.9 #0.7 0.8 0.85 0.9
   do echo 'g: ' $g
   for moments in 70 71 72 73 74 75 #10 20 50 70
	do ./phase -c $tmp_path"rte_g"$g"_moments"$moments"_.dat" > $tmp_path"phg_g"$g"_moments"$moments"_.dat"
	echo '\n\n'
   done
done


cd $script_path
echo 'make plots'

for g in 0.9 #0.7 0.8 0.85 0.9
   do echo 'g: ' $g
   for moments in 70 71 72 73 74 75 # 10 20 50 70
	do python rte_heneygreenstein_plot.py $moments $g 
	echo '\n\n'
   done
done