Use proprietary dispersion correction technique to remove dispersion caused by time domain finite-difference (left figure). Our correction technique is amplitude and spectrum persevering. The numerical errors are around 0.01% (right figure).