Дневник отладки

Всвязи с запарой по курсачу буду вести подобие дневника отладки.
Итак предмет ебли - процедура force. Язык программирование Fortran ebuchka's edition
успехи:

- наконец вводятся все нужные точки (достигнут за счет подгона координат)
- узлы с чего-то стали нумероваться нормально
- статическая нагрузка добавляется
неудачи:
- картинкО опять куда-то делась
- силы до сих пор не проставляются


На 23 февраля;

отображается величина статической нагрузки
- в итерации алгоритм немного упрощен в плане присваивания величин.

- найдена переменная с непонятной функцией NC

- программу клинит на операции FIPARAM = ABS(PRM3 - 180*ATAN((YI-YMAX)/(XMAX-XI))/PI)

На 25 февраля;

- решена проблема с "клином" программы; язык требует строго порядка в переносе по строкам ключей в конструкции if-then-else

- исправлена неточность в вычислении значения нагрузки; вместо центрального угла определялся угол наклонения ребра

- исправлен механизм вычисления принадлежности точек к контуру и их нагрузок; до того неправильно вычислялось значение переменной EC.

- EC перепрофилирована в счетчик принадлежности к ребру.

На 26 февраля

- исправлена ошибка в вычислении величин углов - значения стали ближе к требуемым
- все ещё не изменяется значение сил.

На момент 27 февраля

- обнаружен прискорбный момент; силы не присваиваются
- проблемный сегмент кода
IF(EC.EQ.0) FIPARAM = 180*(ATAN(XMAX-X_R)/(DH-YMAX))/PI
ELSE FIPARAM=ABS(30 - 180*(ATAN(Y_R-YMAX)/(XI))/PI)

16 марта - предыдущие ошибки исправлены
При подключении процедуры calc процесс стопориться. Выход по следующей строке из FormDD
IF ( IEN .LT. XENV(NCOL)) GO TO 500
Cообщение;
+++ ЭЛEMEHTУ 1 HET MECTA B CTPУKTУPE ENV +++++
• Массив XENV
Индексный массив профильного метода. Число элементов данного массива на один больше чем в массиве DIAG.