Explicit finite dif. method seems to work now.
I get evolving concentration profiles. Start modifying for the vacancy
diffusion in an electric field.
Calculation runs, need to fix units.