We present formulas for reduced Wigner phase-space functions for atoms, with an emphasis on the first-order spinless Wigner function. This function can be written as the sum of separate contributions from single orbitals (the natural orbitals). This allows a detailed study of the function. Here we display and analyze the function for the closed-shell atoms helium, beryllium, neon, argon, and zinc in the Hartree-Fock approximation. The quantum-mechanical exact results are compared with those obtained with the approximate Thomas-Fermi description of electron densities in phase space.