Typo in surfaceConcentration.m

[SyPark324]

Dear developers,

Hello, I found some errors in surfaceConcentration.m file.
For the calculation using three paramters model (line 37, 39), I think it should be corrected as
cs_star_p = cs_barrato(1:param.Np)+(param.Rp_p./(Dps_eff .* 35)) .* (-jflux(1:param.Np)) + 8 * param.Rp_p .* Q(1:param.Np);
cs_star_n = cs_barrato(param.Np+1 : end) + (param.Rp_n ./ (Dns_eff .* 35)) .* (-jflux(param.Np+1:end)) + 8 * param.Rp_p.* Q(param.Np+1:end);

Thank you!