Numerical implementation of the harmonic modified mild-slope equation