  Phenomena such as pulsar frequency glitches are believed to be attributable
to differential rotation of a current of ``free'' superfluid neutrons at
densities above the ``drip'' threshold in the ionic crust of a neutron star.
Such relative flow is shown to be locally describable by adaption of a
canonical two fluid treatment that emphasizes the role of the momentum
covectors constructed by differentiation of action with respect to the
currents, with allowance for stratification whereby the ionic number current
may be conserved even when the ionic charge number Z is altered by beta
processes. It is demonstrated that the gauge freedom to make different choices
of the chemical basis determining which neutrons are counted as ``free'' does
not affect their ``superfluid'' momentum covector, which must locally have the
form of a gradient (though it does affect the ``normal'' momentum covector
characterising the protons and those neutrons that are considered to be
``confined'' in the nuclei). It is shown how the effect of ``entrainment''
(whereby the momentum directions deviate from those of the currents) is
controlled by the (gauge independent) mobility coefficient K, estimated in
recent microscopical quantum mechanical investigations, which suggest that the
corresponding (gauge dependent) ``effective mass'' m* of the free neutrons can
become very large in some layers. The relation between this treatment of the
crust layers and related work (using different definitions of ``effective
mass'') intended for the deeper core layers is discussed.
