n_Li_s should be calculated with Area

The old init_cs doesn't contain p.Area in the calculation of theta_n.
Also in the calculation of p.n_Li_s, the calculation should contain
p.Area.