npart_dist = npart_dist + DBLE(npart_tile(bi,bj))
tmp(1) = DBLE(npart_tile(bi,bj))+npart_read
GM_BVP_rModeNumber = 1. _d 0/DBLE(GM_BVP_ModeNumber)
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))
npart = DBLE(ip)
kfloat = DBLE(ip)