Make ionization compute cH correctly from the new pH.
Sped up P matrix computation in calc_pH by using repmat and multiplying by concentration before truncating.