Efficient and accurate defect level modelling in monolayer MoS$_2$ via GW+DFT with open boundary conditions