awk '{print $1, "de"}' HvL.de | sort > jj

awk '{print $1, 2**$3}' Diets.dat | sort | grpmn | \
join -a1 - jj | awk '{print NR, $3, (NF==7?"1":"0")}' | \
sort -n +1 > sens.in

g95 sensitivity.f90

./a.out < sens.in > sens.out

rm jj