flux_fortify_quadratic function

Fortify fluxes with a quadratic fit