function g0 = gamma0(x, a, b, e_g0, u_g0, vt)
%GAMMA0 calculate spatially dependent mobility
%   Based on a tanh profile.
g0 = e_g0*(tanh((x+a+vt)/b)+1)/2+u_g0;
end