We present updated stellar population models appropriate for old ages (> 1 Gyr) and covering a wide range in metallicities (−1.5 [Fe/H] 0.3). These models predict the full spectral variation associated with individual element abundance variation as a function of metallicity and age. The models span the optical-NIR wavelength range (0.37 − 2.4µm), include a range of initial mass functions (IMFs) and contain the flexibility to vary 18 individual elements including C, N, O, Mg, Si, Ca, Ti, and Fe. To test the fidelity of the models we fit them to integrated light optical spectra of 41 Galactic globular clusters (GCs). The value of testing models against GCs is that their ages, metallicities, and detailed abundance patterns have been derived from the HR diagram in combination with high resolution spectroscopy of individual stars. We determine stellar population parameters from fits to all wavelengths simultaneously ("full spectrum fitting"), and demonstrate explicitly with mock tests that this approach produces smaller uncertainties at fixed S/N ratio than fitting a standard set of 14 line indices. Comparison of our integrated-light results to literature values reveals good agreement in metallicity, [Fe/H]. When restricting to GCs without prominent blue horizontal branch populations we also find good agreement with literature values for ages, [Mg/Fe], [