function [pl,ql,pr,qr] = bc_Neumann(xl,ul,xr,ur,t)
%  BC_NEUMANN 
pl = 0;
ql = 1;
pr = 0;
qr = 1;
end