Thermo moist (almost) working; very small difference with CPU version in ql and higher order moments theta, q, ..
Minor change; only call saturation adjustment if ql(Tl) > 0.