function [Q] = bc(Q,JM,KM,jmax,kmax,M_inf)


% simple bc
% At k=1 Extrap rho, u, fix rho v
RHO = Q(JM,1,1);
Q(JM,1,1) = Q(JM,2,1);
Q(JM,1,2) = Q(JM,2,2)./Q(JM,2,1).*Q(JM,1,1);
Q(JM,1,3) = Q(JM,1,3)./RHO.*Q(JM,1,1);


% At j=1 Extrap rho
Q(1,KM,1) = Q(2,KM, 1);

% At jmax Extrap all
Q(jmax,KM,1) = Q(jmax-1,KM, 1);
Q(jmax,KM,2) = Q(jmax-1,KM, 2);
Q(jmax,KM,3) = Q(jmax-1,KM, 3);