Model Order Reduction of Large-Scale Metasurfaces Using a Hierarchical Dipole Approximation