We have performed a new combined set of evaluations for 90−96 Zr, including new resolved resonance parameterizations from Said Mughabghab for 90,91,92,94,96 Zr and fast region calculations made with EMPIRE-3.1. Because 90 Zr is a magic nucleus, 90 Zr is nearly spherical. Therefore, all even-even nuclei in the evaluation set were evaluated using EMPIRE's OPTMAN coupled channels code and a new soft-rotor optical model potential. This then led to improved (n,el) angular distributions, helping to resolve improper leakage in the older ENDF/B-VII.1β evaluation in KAPL proprietary, ZPR and TRIGA benchmarks. Another consequence of 90 Zr being a magic nucleus is that the level densities in both 90 Zr and 91 Zr are unusually low causing the (n,el) and (n,tot) cross sections to exhibit large fluctuations above the resolved resonance region. To accomodate these fluctuations, we performed a simultaneous constrained generalized least-square fit to (n,tot) for all isotopic and elemental Zr data in EXFOR, using EMPIRE's TOTRED scaling factor. TOTRED rescales reaction cross sections so that the optical model calculations are unaltered by the rescaling and the correct competition between channels is maintained. In this fit, all (n,tot) data in EXFOR was used for Ein > 100 keV, provided the target isotopic makeup could be correctly understood, including spectrum averaged data and data with broad energy resolution. As a result of our fitting procedure, we have full cross material and cross reaction covariance for all Zr isotopes and reactions.