function B = b_depz_eps(vdx_fn,f_n,sqrt_M)
global delv; global delt; global eps
[~,Nv,Nz] = size(f_n);
B =delv^2* repmat(eps,[1,Nv,Nz]).*(vdx_fn - f_n/delt)./sqrt_M;
