Refine the full-field grism models with the best fit spectra from individual extractions.
refine_model_with_fits()