Алгоритм программы

Рис. 7. Блок-схема алгоритма программы.

Схема алгоритма программы расчета кривых МТЗ для модели Тихонова - Каньяра приведена на рис. 7. Переменные и константы разных типов выделены цветом (целочисленные - темно-зеленым, вещественные - синим и комплексные - красным).

Вначале программа считывает из входного файла число периодов (NT), первый период (T), геометрический шаг по периодам (Q), число слоев (N), сопротивления слоев (r₁... r_N) в Ом*м и мощности слоев (h₁... h_N_-1) в метрах. Затем счетчику периодов (mT) присваивается номер 1 и начинается цикл по периодам. Определяется круговая частота (w), приведенный импеданс (R). После этого счетчику слоев (m) присваивается номер предпоследнего слоя (N-1) и начинается цикл по слоям. В нем определяется волновое число каждого слоя (k), вспомогательные переменные (A и B). Затем приведенный импеданс (R) пересчитывается с подошвы на кровлю текущего слоя. Счетчик слоев (m) уменьшается на единицу. Если он остается положительным, то цикл продолжается. В противном случае рассчитывается модуль кажущегося сопротивления (r_T) и фаза импеданса (j_T). Вместе с соответствующим корнем из периода () они записываются в выходной файл. Затем счетчик периодов (mT) увеличивается на 1 и период (T) пересчитывается в соответствии с шагом (Q). Если после этого счетчик периодов (mT) не превышает число периодов (NT), то цикл по периодам продолжается. В противном случае программа завершает свою работу.