from pint import quantity
import numpy as np
import pint
ureg = pint.UnitRegistry()
#shape of the galaxy
diameter = 1e5 * ureg.lightyear
thickness = 1000. * ureg.lightyear
volume = np.pi * pow(diameter/2,2) * thickness
# number of terrestrial like planets
planets = 1e6
planets_per_vol = planets/volume
planets_per_vol.to_compact()127.323954474 1/kilolight_year3