Following a previous paper [Y. Shi, Phys. Rev. C 98, 014329(2018)], we present an extension of the density-functional theory to allow for dynamic calculations based on the obtained static Hartree-Fock results. We perform extensive benchmark calculations, by comparing the calculated results with that of an existing code Sky3D. To perform linear-response calculations using the TDDFT method, comparisons have been made with the finite-amplitude quasiparticle random-phase approximation (FAM-QRPA) method. We plan to apply the TDDFT method to a systematic description of the IVD resonances in the Zr, Mo, and Ru isotopes. The strengths of IVD resonances are calculated using two complementary methods: TDDFT and FAM-QRPA methods. For the TDDFT results, additional benchmark calculations have been performed using the well-tested code Sky3D. In these three models, the important ingredients which have major influence on the results, such as time-odd potentials, boundary conditions, smoothing procedures, spurious peaks etc., have been carefully examined. The current TDDFT and the Sky3D codes yield almost identical response functions once both codes use the same time-odd mean fields and absorbing boundary conditions. The strengths of the IVD resonances calculated using the TDDFT and FAM-QRPA methods agree reasonably well with the same position of the giant dipole resonance. Upon seeing a reasonable accuracy offered by the implemented code, we perform systematic TDDFT calculations for spherical Zr and Mo isotopes near $N=50$, where experimental data exist. For neutron-rich Zr, Mo, and Ru isotopes where shape evolution exist we predict the photoabsorption cross sections based on oblate and triaxial minima.