using last time step's cmax value when computing adsorption for inconp solver
add computeinjecprod function for incomp solver.
output watercut for incomp polymer solver.