fixing flux integration for case with a single free parameter (i.e., a pumping well of specified head or flux). now resuts are constant in time