In this work we report a parametric reduced order model (ROM) based on the proper orthogonal decomposition (POD) method with Galerkin projection for solving the system of time-domain Maxwell's equations. In particular, we introduce a residual-based estimation of the error associated with the ROM. Moreover, a greedy algorithm for the snapshot selection in the parameter space is developed. We investigate the behavior of the parametric POD-Galerkin ROM by considering the scattering of a plane wave by a 2-D multi-layer heterogeneous medium.