Продраться через атмосферу.

Подождите, там нормальная сканеная книжка, в djvu формате и 2005 года, вы о чем?
Весит 10 мегов

Вот картинка из этой программы и модели моего двигателя. Тут взята виртуальная поверхность в сечении сопла на выходе из него. Само тело, создающее эту плоскость исключено из модели как твёрдое тело, поэтому программа расчёта не видит его и оно нужно тока для снятия параметров газа на выходе из сопла. Слева то самое окошко, в котором записаны параметры газа по сечению, их на самом деле ищё больше, в настройках их можно выбрать, так чтобы не путаться лишними строками... Там есть и вкладка с интегрированными значениями, но тут можно обойтись и средним значением параметра по поверхности. Цвет в данной картинке показывает распределение температуры газа.
Слева находится дерево построения модели.
Inlet Mass Flow1 - задаётся секундный расход газа по горящей поверхности, там же определяется куча его параметров, часть которых берётся из PROPEPа
Real Wall1 - задаются стенки сопла как реальная поверхность
Environment Pressure1 - Окружающее давление, в нём выбраны те стенки области, которые не являются частью двигателя а ограничивают расчётную область, определяя давление окружения как нормальное атмосферное.
VS Heat Generation Rate1 - Подвод тепловой мощности в тонкую область, которая покрывает всю поверхность горения в модели и тем самым моделирует процесс передачи тепла в продукты сгорания, само значение мощности вычисляется из данных хим. равновесия PROPEP (тепловой ээфект реакции * на расход).
Ниже идут уже инструменты для снятия параметров.
Например XY Plot может строить график в экселе по заданному эскизу в модели для выбранных параметров.

Картинки эт конечно здорово, но все-таки хотелось бы увидеть цифры
У нее есть скриптовый язык?

Люди, давайте проверять в лоб. Если у нас лобовое сопротивление равно весу, значит, совершаемая им работа равна работе, совершаемой силой тяжести.
Поскольку мы считаем, что большую часть пути проделываем с постоянной скоростью, то, значит, вес и тяга падают экспоненциально от времени. Очень удобно интегрировать. ДО этой скорости мы разгоняемся равноускоренно. Значит, тут работа растёт квадратично от времени. Наконец, когда полезный груз оказывается в апогее, у него есть потенциальная энергия, соответствующая mgH. Просуммировав работу на двух участках и помножив её на два, и добавив последнюю компоненту, мы должны получить искомую механическую работу, которая должна соответствовать (с домножением на кпд цикла Карно) энергии той карамели, которая была потрачена на подъём.

Другой способ проверки - посчитать характеристическую скорость, вычесть из неё гравитационные потери, которые исчисляются, как произведение времени работы двигателя на "же", и, считая, что аэродинамические - такие же, вычесть ещё раз. А потом посчитать высоту баллистического подъёма, считая разгон одномоментным. Правда, этот вариант грубоват.

Ну, прикинем. Число Циолковского 16,5. при УИ 1200 м/с это даёт ХС 3310 м/с. Очень немало, заметим в скобках, значительно больше, чему Фау-2 (2900 м/с). Так что удивляться подъёму на 100 км не приходится. Если мы летим до окончания плотных слоёв (20 000 м) со скоростью 250 м/с, это 80 сек, т.е. гравитационные потери 800 м/с. Плюс аэродинамические - те же 800 м/с. Остаётся 3310-1600=1710 м/с. h=v²/2g=146 км. Плюс-минус лапоть проходит.

Где подвох? в исходных данных модели. Они нереализуемы.