import pint; u = pint.UnitRegistry()

W   = 893 * u.kW*u.hr/u.month
h   = 100 * u.m
rho = 1000 * u.kg/u.m**3
g   = 9.81 * u.m/u.s**2

print(f"W = {W.to(u.kW):.2f}")

# W = mdot * g * h
mdot = W / g / h
Vdot = mdot / rho

print(f"Vdot = {Vdot.to(u.gal/u.min):.2f}")