even when the dot-product of the Reynolds stress and the gradient of velocity is
negative.  Currently this is done by taking the magnitude, another alternative
would be to take only the positive part.

Better properties for CH4N2O
Limiting the temperature to the boiling-point in parcel