Fixed calculation of n_prot in single charge model