4In some specific cases, it is possible to consider Dirichlet type boundary conditions, espcially at p = 0. In that case, the Maxwellian distribution function is enforced. However, in most cases, this type of condition must avoided, to keep the most general approach, and in particular to cross-check that the conservative nature of the code is well satisfied