`1pnc` disregards `useMoles` in the choice of primary variables
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?