cs_surface cs[mol/m^3] "Insertion particle concentration, surface"
i0 i0_ref_distr*(2*(csmax-cs_surface)/csmax)^0.5*(2*cs_surface/csmax)^0.5*(cl/1[M])^0.5 "Exchange current density"
Eeq mat2.def.Eeq(cs_surface/csmax) "Equilibrium potential"
eta phis-phil-Eeq "Overpotential"
iloc i0*(exp(eta*0.5*F_const/(R_const*T))-exp(-eta*0.5*F_const/(R_const*T))) "Local current density on particle surface"