В работе представлена численная реализация модели распространения КВ-радиоволн в ионосфере, основанная на решении уравнения эйконала методом характеристик при задании параметров среды (электронной концентрации) на основе экспериментальных данных низкоорбитальной (НОРТ) и высокоорбитальной (ВОРТ) радиотомографии и магнитного поля Земли согласно модели IGRF. Использование данных РТ дает возможность улучшить описание региональных и глобальных особенностей распределения концентрации ионосферной плазмы в задаче численного моделирования распространения КВ радиоволн, по сравнению с традиционно используемыми для этих целей моделями IRI и NeQuick, в особенности в областях сильных пространственных градиентов, играющих определяющую роль при формировании лучевых траекторий КВ-радиоволн. В модели также реализован метод пристрелки, использующий алгоритмы симплексной минимизации, для моделирования ионограмм наклонного зондирования. В работе приводятся примеры работы развитой модели в различных гелиогеофизических условиях, в присутствии как естественных, так и искусственных ионосферных возмущений. Рассмотрены примеры сильнейших геомагнитных бурь октября-ноября 2003г и марта 2017г, а также искусственных возмущений, вызванных воздействием нагревного стенда "Сура".