Reimplement energy exchange between gas and dust

Now using my own exact formula.
Different dust species are treated separately.