Suppose I have a microcanonical ensemble made of $N$ particles immersed in a (viscous) fluid at equilibrium temperature $T$. When $N$ gets very large, the particles can be treated as a continuum. Their current density vector is defined to be $\mathbf{j}(\mathbf{r}, t) = \rho(\mathbf{r},t)\,\mathbf{v}(\mathbf{r}, t)$, where $\rho$ is a pdf (giving the fraction of particles per unit volume at time $t$) and $\mathbf{v}$ is the flow velocity.
Question.What is the link between $\rho(\mathbf{r},t)$ defined above, and the Boltzmann distribution, which gives the fraction of particles at a given energy $${\displaystyle p_{i}={\frac {1}{Q}}}{e^{-{\varepsilon }_{i}/(kT)}={\frac {e^{-{\varepsilon }_{i}/(kT)}}{\sum _{j=1}^{M}{e^{-{\varepsilon }_{j}/(kT)}}}}}$$ where $p_i$ is the probability of state $i$, $\varepsilon_i$ the energy of state $i$, $k$ the Boltzmann constant, $T$ the absolute temperature of the system and $M$ is the number of all states accessible to the system at temperature $T$.