Simulated, steady-state mass outflow rates from the 14 planetary systems considered in this work. The values obtained by ATES (in black) are compared to those by TPCI (S16, in red). Overall, the estimated mass loss rates agree to within a factor of 2 in all cases. In those cases where ATES differs from TPCI, the differences are always upward, and are likely due to inclusion of ion advection, which ATES implements in post-processing, yielding systematically higher atmospheric density profiles. This implementation choice, however, is chiefly responsible for the dramatic gain in computational time.