Расчет характеристик рассеяния агрегатных частиц проводится с помощью 
метода суперпозиции Т-матриц, и соответствующий код описан в статье
Mackowski D.W., Mishchenko M.I. 
A multiple sphere T-matrix Fortran code for use on parallel computer clusters // 
J. Quant. Spectrosc. Radiat. Transfer. 2011. V. 112. P. 2182–2192. http://dx.doi.org/10.1016/j.jqsrt. 

Здесь выложен соответствующий exe-файл (mack.exe). 
Для этих вычислений также необходимы два файла -
- файл конфигурации с координатами составляющих агрегат частиц (см. Agr_structures)
- файл mstm.inp с остальными параметрами агрегата и процедуры расчета (дан пример для спектра ледяных частиц).

Описание всех параметров, приводимых в этом файле, дано в указанной статье.
Практически для данных расчетов надо корректировать следующие строки файла:
  параметр			пример
number_spheres :  		50
sphere_position_file : 		Dmack_80_30_50.dat
output_file : 			ice-L100-050r11_C.dat
length_scale_factor : 		0.69120E+00  (параметр размера частиц в агрегате)
real_ref_index_scale_factor : 	0.13015E+01  (действительная часть показателя преломления на данной длине волны)
imag_ref_index_scale_factor : 	0.00000E+00  (мнимая часть показателя преломления)

Значение параметра length_scale_factor определяется радиусом r составляющих агрегат частиц
и длиной волны, т.е., 2*pi*r/wavelength.

Если необходимо сделать расчеты для нескольких длин волн, надо после параметра 
" fixed_or_random_orientation
  1 "
объявить начало нового цикла : new run
и ввести те параметры, которые соответствуют новой длине волны.
Пример дан в прилагаемом файле mstm.inp

В результате вычислений формируется файл, который содержит как входные данные,
так и полученные значения эффективности экстинкции, поглощения, рассеяния, 
элементы матрицы рассеяния в зависимости от угла рассеяния и
коэффициенты разложения этих элементов по сферическим функциям, 
которые на следующем этапе могут быть использованы для расчета переноса излучения
в среде, состоящей из данных агрегатных частиц.
Один из таких файлов приведен здесь в качестве примера. 

=========================================================================================

Характеристики рассеяния одиночных сферических частиц вычисляются по теории Ми
с помощью программы, разработанной М.Мищенко и доступной на 
https://www.giss.nasa.gov/staff/mmishchenko/Lorenz-Mie.html.
Программа снабжена подробными комментариями и не требует дополнительных пояснений. 

Здесь выложена версия, sphera_wl_mr.for, адаптированная под выдачу характеристик для нескольких длин волн.
Она формирует два файла 
- *_ff.dat для значений фазовой функции F11 и др. элементов  матрицы рассеяния в зависимости от угла рассеяния,
      	он также содержит эффективности экстинкции, поглощения, рассеяния частиц;
- *_co.dat для коэффициентов разложения этих элементов по сферическим функциям,
	он также содержит альбедо однократного рассеяния и число коэффициентов разложения.
	Этот файл необходим для расчетов переноса излучения в слое таких частиц.
Примеры этих файлов приведены.


=========================================================================================

Результаты расчетов характеристик однократного рассеяния запакованы в папках, 
названия которых соответсвуют названиям материалов частиц и их количеству в агрегатах.
Характеристики одиночных частиц находятся в папках, содержащих 1 в названии.