change procedure to fit, add force_positive_flux by fitting log(flux), allow for increased weight on bundle sides