function [pu,qu,po,qo]=BoundaryCond(yu,uu,yo,uo,t)
pu=[uu(1);uu(2);uu(3);uu(4);uu(5);uu(6);0;0;0;0;uu(11);0;0;0];
qu=[0;0;0;0;0;0;1;1;1;1;0;1;1;1];
po=[uo(1)-1;uo(2);uo(3);uo(4)-1;uo(5);uo(6);0;0;0;0;uo(11);0;0;0];
qo=[0;0;0;0;0;0;1;1;1;1;0;1;1;1];
end