(March 2016, J. Skowron)
INPUT: (mag, error) AND COEFFICIENTS FOR THE GIVEN FIELD.CHIP: (GAMMA, EPSILON, UNIT, BACKG, MAG_01) if error > 0.003: new_error = sqrt( (GAMMA*error)^2 + EPSILON^2 ) else: F = 10^(0.4*(UNIT - mag)) B = 10^(0.4*(UNIT - BACKG)) FLUX_NOISE = sqrt(F + B) MAG_NOISE = 2.5/ln(10) * FLUX_NOISE/F new_error = max( error, sqrt( (GAMMA*MAG_NOISE)^2 + EPSILON^2 ) , 0.01 * 10^( 0.6*(MAG_01 - mag) ) ) end RETURN (mag, new_error)