Finite difference time domain method