This function produces a fermi-level dependent diffusion coefficient given by
Parameters:
is a field or function
is a variable