corrected an error in the integration limit for integration in the positive radial direction - the error left one point not being set in the Xi0 array (X = i/j/k)
corrected an error in the integration limit for integration in the positive
radial direction - the error left one point not being set in the Xi0
array (X = i/j/k)

chengdi123000/CFL3D,chengdi123000/CFL3D