npart_dist = DBLE(npart_tile(bi,bj))
tmp(1) = DBLE(npart_tile(bi,bj))
tmp(2) = DBLE(myCurrentIter)
tmp(1) = DBLE(npart_tile(bi,bj))+npart_read
npart = DBLE(ip)
kfloat = DBLE(ip)
X(L) = SNGL(1.0D0/DSQRT(DBLE(FLOAT(4*L*L-1))))
Y(L) = SNGL(DSQRT(DBLE(FLOAT(L))))
Y(L+LMAX) = SNGL(DSQRT(DBLE(FLOAT(L+LMAX))))
C = SNGL(DSQRT(1.0D0-DBLE(S)*DBLE(S)))
ZZ = ZZ * DBLE(FLOAT(L+L+1))/DBLE(FLOAT(L+L))