  We present results from a calculation of beyond the standard model (BSM) kaon
mixing including data physical with light quark masses. We simulate $N_f=2+1$
QCD with Iwasaki gauge and domain wall fermion action on 8 ensembles, spanning
3 lattice spacings and pion masses from the physical value up to 430MeV. The
ratio of the BSM to standard model (SM) matrix elements are extracted from the
correlation functions and renormalised using the RI-SMOM Rome-Southampton
method with non-exceptional kinematics. The results at the physical point
continuum limit are found by performing a simultaneous continuum chiral
extrapolation. In this work we gain consistency with our previous results and
achieve a reduction in both the statistical and systematic error.
