Even if this is a little late (sorry):
@ax3l @psychocoderHPC 
Either the old documentation $\rho = exp[ - Y * D ]$ (line 105 - red) or the new documentation (line 106 - green) contradict the old or new setting (line 109 - red and line 110 - green). GAS_D_SI was $10^{-3}$ and is now $10^+3$! 
What is correct?