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