High‐order compact scheme for Boussinesq equations: implementation and numerical boundary condition issue