A 3-D FDTD scheme for the computation of HPM propagation in atmosphere