`1pnc` disregards `useMoles` in the choice of primary variables
The 1pnc model uses mole fractions as primary variables regardless of the value of useMoles. So if useMoles is set to false, one has to set Dirichlet conditions still in terms of moles, while sources and fluxes expect masses.
Should we drop the useMoles switch or implement it properly?