Целью данной работы является моделирование донного таяния для всего ледяного щита Антарктиды. Для этого была разработана трехфазная численная модель на основе одномерной однофазной задачи Стефана. Для дискретизации уравнений используется неявная разностная схема на неоднородной сетке первого порядка точности. Реализован параллельный алгоритм и выполнено моделирование донного таяния Антарктиды. Получены распределение скоростей донного таяния. Средняя скорость донного таяния составила 29 Гт/год.