case="FeSb2"

def=0.0001
nef=3  #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 .
  cp -f ../inp.imp .

  ~/WORK/PRO/MYSTUFF/IMP/Aimp/Aimp >out.Aimp
#  ~/WORK/PRO/MYSTUFF/IMP/MUimp/mu >out.muimp

##  ~/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