case="FeAs2"

def=0.00001
nef=20  #per +/-

rm response.dat_?a

for i in $(seq 0 1 $nef)
do


a=$(echo "$i*$def" | bc)
echo $a
mkdir d$a

cd d$a

rm response.dat_?

  echo $a " 91.37  ">inp.doping
#  ln -s ../$case.weight .
#  ln -s ../$case.energy .
  ln -s ../$case.kgen .
  ln -s ../$case.outputkgen .
  ln -s ../$case.struct .
  ln -s ../phi.dat .
  ln -s ../inp.temp .
  ln -s ../inp.optic .
  ln -s ../inp.ef .
  ln -s ../inp.trans .
  ln -s ../awtrace.dat .
  ln -s ../inp.imp .
  ln -s ../inp.gap .
  ln -s ../inp .

#  /home/users/jtomczak/PRO/MYSTUFF/MU/mu >out.mu
  /home/jtomczak/WORK/PRO/MYSTUFF/IMP/MUimp/mu >out.muimp

#~/WORK/PRO/MODELS/SPECIAL/SEEBECK_coh/a.out
#mv muT.dat inp.efshift_list

  ../wientrans >out.trans

#./wientrans >out

  cd ..

less d$a/response.dat_1 >>response.dat_1a
less d$a/response.dat_2 >>response.dat_2a
less d$a/response.dat_3 >>response.dat_3a

echo "">>response.dat_1a
echo "">>response.dat_2a
echo "">>response.dat_3a

done