#run from cdos dir

echo "masses :"
massFL < ../sig.inp  #matsub sig
massFL < ../sig.inp >mstar.dat

echo "pi*N(Ef) :"
#awk '($1==0.00500000)' FeSi.gc3 >tmp  # yes ... i know
#awk '($1==0.00500000)' FeSi.cdos >tmp2  # yes ... i know


awk '(($1>=-0.0295)&&($1<=-0.028))' FeSi.gc3 >tmp  # yes ... i know
awk '(($1>=-0.0295)&&($1<=-0.028))' FeSi.cdos >tmp2  # yes ... i know


ug tmp >Nef.dat
more Nef.dat
more tmp2 >>Nef.dat
rm tmp tmp2

echo "gamma :"
gamma