Boundary condition

To simplify the problem, we assume the shock wave will not reach the boundary of our simulation region(this can be controlled by looking at the profile of gas parameters at certain times). Thus we can simply select fixed boundary conditions. f(x_b,t)=f(x_b,0).