Example table of saturating flux. The tabulated values correspond to the function values of 'saturationAnalytic' and 'der_saturationAnalytic', psi(psi0) and d_psi/d_psi0(psi0).
For use with 'Base.Math.interpolateTable':
y = interpolateTable(psi0, saturationTab) psi = y[1] der_psi = y[2]