This constraint appears to assume f is always positive. Need to use the y variable to get the positive f when flow is reversed