We present a proof of the Riemannian Penrose inequality with charge in the context of asymptotically flat initial data sets for the Einstein-Maxwell equations, having possibly multiple black holes with no charged matter outside the horizon, and satisfying the relevant dominant energy condition. The proof is based on a generalization of Hubert Bray's conformal flow of metrics adapted to this setting.