distancek(k)=0.5d0*distance(k)**2*sigma_odl(k,ii) ! waga_dist rmvd from Gaussian argument
c write (iout,*) "sigma_odl(",k,ii,") =",sigma_odl(k,ii)
c write (iout,*) "distancek(",k,") =",distancek(k)
c distancek(k)=0.5d0*waga_dist*distance(k)**2*sigma_odl(k,ii)
distancek(k)=0.5d0*distance(k)**2*sigma_odl(k,ii) ! waga_dist rmvd from Gaussian argument
c write (iout,*) "sigma_odl(",k,ii,") =",sigma_odl(k,ii)
c write (iout,*) "distancek(",k,") =",distancek(k)
c distancek(k)=0.5d0*waga_dist*distance(k)**2*sigma_odl(k,ii)