Small boundary mass flux of scalar leads to null values of concentrations in CS v6.0 and v7.0

Putting a small mass flux Qm of the order of microgram/m2/s ( rcodcl (ifac, isca(iscal), 3) = -Qm in cs_user_boundary_conditions.f90 ) leads to zero mass fraction values, when larger values of Qm give finite values. I tried different turbulence models (k-eps, k-omega-sst) with different wall functions (2 scales + scalable with/without log laws). My y+ are quite small (<1).
Following one advice from Y.Fournier, I rewrote the routine as a volumic source term close to the emiting surface, but the problem persists.

Start with cs_user_boundary_conditions.f90 and the rcodcl boundary-condition assignment for scalar mass flux. Reproduce the case with microgram/m2/s and larger Qm values, then compare the boundary treatment with the volumic-source attempt. Done means identifying why the small flux produces zero mass fractions and documenting or correcting the behavior.

Written by the indexing model from the issue text.