flow methods should set bounds on flow.lb and flow.ub instead of adding constraints

[pablormier]

MultiSampleFBA, MultiSampleIMAT should set the flux bounds using the flow variable and the flow.lb and flow.ub, since indicators manipulate these