The depletion of neutral helium atoms has been studied in an unmagnetised spherical plasma created by DC discharge in a multidipole confinement field. Knowing the neutral density profile is critical to predicting the equilibrium flow of such plasmas. A model of the emissivity due to electron-impact excitation of neutral atoms in the plasma has been derived and used to fit radiance measurements of several neutral transitions to extract the radial profile of neutral density for plasmas of varying temperature and density. We report a depletion of the core neutral density varying between negligible levels to 80 % of the edge neutral density depending on the input power and fuelling. The corresponding ionisation fraction varies between 30-80 % in the plasma core. A simple neutral diffusion model is sufficient to describe the shape of neutral density profile implied by the radiance measurements. We have used the measurements to include a drag force due to neutral charge-exchange collisions in simulations of driven plasma flow. The simulation predicts a better fit to Mach probe flow measurements when this neutral drag is accounted for. This work shows that accounting for a realistic neutral profile is important to predict the plasma flow geometry and its magnetohydrodynamics (MHD) stability.