An advanced ray-tracing numerical method is coupled to an ionospheric plasma transport code to model the self-consistent propagation of an obliquely incident high-frequency (HF) beam in the ionosphere. An HF radar propagation scenario is presented as an application. During the daytime, a density increase in the lower ionosphere is predicted causing the defocusing of the incident HF beam.