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.

cs_user_boundary_conditions.f90 と、スカラー質量流束に対する rcodcl 境界条件の割り当てから始めます。microgram/m2/s およびより大きな Qm 値でケースを再現し、その後、境界の扱いを体積ソースの試行と比較します。小さい流束によって質量分率がゼロになる理由を特定し、その挙動を文書化または修正できれば完了です。

索引モデルが issue の本文から書いたものです。