Vn -cp.iloc_per1/(cap*rho) "Anode surface dissolution velocity"
epsl max(min(phils,1),1e-3) "Electrolyte volume fraction"
phils_init min(1,2*sphavg(h_interface,dom==1,8)) "Initial level set variable expression"