@jngaravitoc I'm not sure what EXP gives you or assumes for units, but let's say you have to set a mass scale, time scale, and G=1. If you have output from an EXP sim as positions exp_xyz and exp_vxyz, you could do this:

import astropy.units as u
usys = SimulationUnitSystem(mass=1e6 * u.Msun, time=100*u.Myr)
(exp_xyz * usys["length"]).to(u.kpc)
(exp_vxyz * usys["velocity"]).to(u.km/u.s)

Thanks @adrn! I think for EXP when analyzing external simulations what we will do is: usys = SimulationUnitSystem(mass=1 * u.Msun, length=1*u.kpc, G=1)

exp assumes nothing about units, other than G=1. In other words, the user provides positions and velocities and that determines the time units. This seemed like a clean choice at the time (EXP actually predates stuff like Gadget).

Anyway, Adrian's solution is tré cool!