The authors report on the application of the FDTD method to the simulation of radiation from vertical dipole antennas near a interface of air and lossy soil. Maxwell's equations are solved in cylindrical coordinates taking the rotational symmetry of antennas into account. The truncation of the lossy soil at the boundary is performed by some special PML technique. By using non-uniform mesh generation and restricting the computation to the computaional window, where the radiated pulse is passing, the authors can solve large-scale problems and compute far-field values directly by the FDTD method. The authors provide amplitude and phase distributions of radiated fields for different types of soil and antenna positions.