c sigma_odl(k,ii)=hmscore(k)*rescore(k,i)*rescore(k,j)
else
sigma_odl(k,ii)=(rescore(k,i)+rescore(k,j))* ! sigma ~ rescore ~ error
- & dexp(0.5d0*(odl(k,ii)/dist_cut)**2)
+ & dexp(0.5d0*(odl(k,ii)/dist_cut)**2-0.5d0)
c Following expr replaced by a positive exp argument
c sigma_odl(k,ii)=(rescore(k,i)+rescore(k,j))*