proc/_calir
Вычисляет массив комплексных значений коэффициента преломления для набора значений энергии и по заданной химической формуле
Входные данные: r(50) - число сортов атомов в хим.-формуле материала
r(51),r(52), ... -- Z и N атомный номер и число атомов для каждого сорта
r(1), r(2), r(3) -- значения начальной и конечной энергии в keV и число точек, одна точка возможна, r(3)<411
r(4) -- плотность вещества в гр/ксм, r(5)=CT -- двухразрядное число
T -- если = 0, то для расчета f1 и f2 используются таблицы Windt, если =1, то таблицы Chantler, если =2, то интерполяция
C -- учет вклада комптоновского поглощения, если 0, то Cумма, если 1, то НЕТ комптона, если 2, то ТОЛЬКО комптон   

Результат записывается в r(101) как E1 E2 ... D1 D2 ... B1 B2 ...
Используемая память: B C C1 C2 D E X Z b c1 c2 d i j k l m n1 n2 o1 o2 o3 o4 o5 s u v x   
r(1)-r(5000) i(101)-i(382)
На выходе из процедуры &=1, если расчет прошел успешно, и 0 если нет.

Структура файлов f1f2x.dat
первые 282 числа содержат 94 тройки чисел Z, N, N1, где Z -- атомный номер, N -- число значений энергии, N1 -- номер первого числа с таблицей для данного атома. В таблице находится N троек чисел E, f1, f2
в шапке файла могут быть указаны не все атомы и в произвольном порядке. Структура файла допускала его дополнение. При этом пустые ячейки файла могут заполняться, а таблицы наращиваться. 
Пока в обоих таблицах есть не все атомы. Какие атомы есть легко проверить, напечатав таблицу

Структура файла mmpar.dat   
этот файл содержит таблицы из 94 чисел для
Атомные веса для 94 элементов, начиная с первого  1-94
Коэффициент A0 для фотоэлектронного поглощения  E > Ek  95-188
Коэффициент A1 для фотоэлектронного поглощения  E > Ek  189-282
Коэффициент A2 для фотоэлектронного поглощения  E > Ek  284-376
Коэффициент A3 для фотоэлектронного поглощения  E > Ek  377-470
Коэффициент A0 для фотоэлектронного поглощения  E < Ek  471-564
Коэффициент A1 для фотоэлектронного поглощения  E < Ek  565-658
Коэффициент A2 для фотоэлектронного поглощения  E < Ek  659-752
Коэффициент A0 для комптоновского рассеяния  753-846
Коэффициент A1 для комптоновского рассеяния  847-940
Коэффициент A2 для комптоновского рассеяния  941-1034
Коэффициент A3 для комптоновского рассеяния  1035-1128
E K in keV  1129-1222
E L1 in keV  1223-1316
Размер файла 5264 соответствует. Для комптона использовались коэффициенты с 9 по 12, что соответствует указанной схеме.
Сами таблицы также соответствуют, что проверено.  