Rather than computing the scaling factor and the χ² in one block over the entire grid of models, let's compute them band by band. The main advantage is that we avoid the creation of temporary arrays that can be as large is the model fluxes array, which can be very large. Because we loop on the bands rather than on the models, the impact on the performance should remain very limited.