Model Order Reduction in Contour Integral Methods for Parametric PDEs