case="FeAs2"

def=0.0001
nef=30  #per +/-

rm response.dat_?a

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


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

cd d$a

rm response.dat_?

  echo $a >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 .


  /home/users/jtomczak/PRO/MYSTUFF/MU/mu >out.mu
#  /home/jtomczak/WORK/PRO/MYSTUFF/MU/mu >out.mu
  ../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