enddo ! k
#else
dijCACA=dist(i,j)
- sigma2CACA=0.25d0/(restok(itype(j))**2+restok(itype(i))**2)
+ sigma2CACA=scal_rad**2*0.25d0/
+ & (restok(itype(j))**2+restok(itype(i))**2)
do k=1,nsaxs
dk = distsaxs(k)
expCACA = dexp(-0.5d0*sigma2CACA*(dijCACA-dk)**2)