This work develops a dynamic homogenization approach for metamaterials. It finds an approximate macroscopic homogenized equation with constant coefficients posed in space and time; however, the resulting homogenized equation is higher order in space and time. The homogenized equation can be used to solve initial-boundary-value problems posed on arbitrary non-periodic macroscale geometries with macroscopic heterogeneity, such as bodies composed of several different metamaterials or with external boundaries. First, considering a single band, the higher-order space derivatives lead to additional continuity conditions at the boundary between a homogeneous material and a metamaterial. These provide predictions of wave scattering in 1-d and 2-d that match well with the exact fine-scale solution; compared to alternative approaches, they provide a single equation that is valid over a broad range of frequencies, are easy to apply, and are much faster to compute. Next, the setting of two bands with a bandgap is considered. The homogenized equation has also higher-order time derivatives. Notably, the homogenized model provides a single equation that is valid over both bands and the bandgap. The continuity conditions are applied to wave scattering at a boundary, and show good agreement with the exact fine-scale solution. Using that the order of the highest time derivative is proportional to the number of bands considered, a nonlocal-in-time structure is conjectured for the homogenized equation in the limit of infinite bands. This suggests that homogenizing over finer length and time scales -- with the temporal homogenization being carried out through the consideration of higher bands in the dispersion relation -- is a mechanism for the emergence of macroscopic spatial and temporal nonlocality, with the extent of temporal nonlocality being related to the number of bands considered.