Solving time domain Helmholtz wave equation with MOD-FDM