The function computes a total flux at each perforation (connection), not just one total flux per well. Update documentation accordingly.