(1h) d D g d t = ( 1 − ϵ ) κ ρ H g + ω g ψ C g