Corrected index.

The counter of Q_flow_fixed was incorrect.
Changed source.

This avoids a numerical Jacobian for the initialization problem.