Support time-varying diffusivity coefficient

[edwardlavender]

It would be great to be able to have one diffusivity value for resting behaviour and another one for active behaviour.

[edwardlavender]

probably just deal with D==0 and D != 0 is enough