Accurate numerical simulations of the complex wind flows in the Mexico City Metropolitan Area (MCMA) can be an invaluable tool for interpreting the MILAGRO field campaign results. This paper uses three methods to evaluate numerical simulations of basin meteorology using the MM5 and WRF models: statistical comparisons with observations, "Concentration Field Analysis" (CFA) using measured air pollutant concentrations, and comparison of flow features using cluster analysis. CFA is shown to be a better indication of simulation quality than statistical metrics, and WRF simulations are shown to be an improvement on the MM5 ones. Comparisons with clusters identifies an under-representation of the drainage flows into the basin and an over-representation of wind shear in the boundary layer. Particle trajectories simulated with WRF-FLEXPART are then used to analyse the transport of the urban plume and show rapid venting and limited recirculation during MILAGRO. Lagrangian impacts were identified at the campaign supersites, and age spectra of the pollutants evaluated at those same sites. The evaluation presented in the paper show that mesoscale meteorological simulations are of sufficient accuracy to be useful for MILAGRO data analysis.