docs: def soma_surface_areas(nrn_pop, neurite_type=NeuriteType.soma):
      '''Get the surface areas of the somata in a population of neurons
  
      Note:
          The surface area is calculated by assuming the soma is spherical.
      Note:
          If a single neuron is passed, a single element list with the surface
          area of its soma member is returned.
      '''
      nrns = neuron_population(nrn_pop)
      assert neurite_type == NeuriteType.soma, 'Neurite type must be soma'
      return [soma_surface_area(n) for n in nrns]