CC BY
18О ДРИФТЕРНЫХ НАБЛЮДЕНИЯХ ЗА ПОВЕРХНОСТНЫМИ ТЕЧЕНИЯМИ МИРОВОГО ОКЕАНА
Касьянов С.Ю., Музыка Г.В., Никитин О.П. (Nikitin@sirena3.ru)
Государственный океанографический институт Росгидромета
1. Введение
Большие возможности для изучения циркуляции в навигационном слое Мирового океана открывает широко практикующийся в настоящее время за рубежом метод слежения за дрейфующими с поверхностными течениями буйками - дрифтерами. Эти свободно плавающие буйки оснащены якорем-парусом на определенной глубине и системой радиосвязи с искусственными спутниками Земли (ИСЗ). Путем определения местоположения буйков через определенные промежутки времени с помощью системы радиосвязи (в настоящее время повсеместно используется система ARGOS) оказывается возможным с большой точностью трассировать траектории буйков в океане. Указанный метод позволяет изучать течения в лагранжевой системе координат, следить за их характером и изменчивостью на протяжении длительного времени (месяцы и годы). На рис. 1 схематически изображен один из применяющихся типов буйка - дрифтера, ставший с 1991 года стандартом.
Рис. 1. Схема дрейфующего по поверхности буйка с заглубленным парусом.
В бывшем СССР ответственной организацией за разработку системы наблюдений за течениями с помощью дрифтеров и ИСЗ был Морской гидрофизический институт (УССР, Севастополь). С распадом СССР ни на Украине, ни в России подобные работы практически не ведутся, зарубежные данные дрифтерных наблюдений за течениями Мирового океана надлежащим образом не собираются.
В то же время история таких наблюдений за поверхностными течениями насчитывает уже более 20 лет, поскольку океанографические буйки (дрифтеры) массово используются в рамках национальных и международных проектов исследования океанской циркуляции, начиная с 1979 г. (запуски дрифтеров до 1979 г. носили в значительной мере экспериментальный характер, в них отрабатывалась конструкция дрифтера и методика измерений). Дрейфующие в поверхностном слое океанографические буйки запускаются в самых различных частях океанов до нескольких сотен экземпляров в год и в последние годы в Мировом океане единовременно дрейфуют и передают информацию до 1000 буйков. На рис. 2 точками указано местоположение 920 дрейфующих в Мировом океане буйков по состоянию на 1 мая 2000 г.
Рис. 2. Местоположение дрифтеров в Мировом океане по состоянию на 1 мая 2000 г.
Большинство дрифтеров оснащено датчиком температуры, с 1993 г. стали широко эксплуатироваться дрифтеры с датчиками давления, а с 1 997 г. - также и с измерителями скорости и направления ветра. По спутниковым каналам помимо данных о местоположении дрифтера передается также информация и об измеряемых параметрах среды вблизи поверхности океана. В соответствии с ранее заключенными международными соглашениями оперативная информация с дрейфующих буев поступает по каналам глобальной системы телесвязи (ГСТ) Всемирной метеорологической организации (ВМО) в метеоцентры многих стран мира, в том числе в Гидрометеоцентр (ГМЦ) Роскомгидромета, где она декодируется и архивируется. Однако, поскольку ГМЦ не является профилирующей организацией по изучению
океана, на практике используются (для объективного анализа, составления прогнозов и климатических исследований) в основном только данные о температуре в точках измерений.
2. Создание архива данных наблюдений за поверхностными течениями
В Государственном океанографическом институте (ГОИН) Росгидромета было принято решение о сборе исторических данных дрифтерных наблюдений в рамках подпрограммы 10 "Создание единой системы информации об обстановке в Мировом океане" федеральной целевой программы "Мировой океан" и о создании на их основе информационно-справочной подсистемы о фактических течениях в поверхностном слое Мирового океана (ИСП "Течения").
2.1. Поступающие по каналам ГСТ данные
Первоначально предполагалось использование накопленного в ГМЦ архива с целью его применения в качестве основы создания ИСП "Течения". Однако от этого отказались по следующим причинам.
Во-первых, из-за организационных трудностей, необходимости модификации программного обеспечения в связи со сменой кодировки передаваемой по ГСТ информации в массиве ГМЦ содержатся многомесячные пропуски, когда поступающая по ГСТ информация не декодировалась и следовательно не архивировалась (до 1.11.1991 г. использовался код DRIBU, затем до 2.11.1994 г. - код DRIFTER, в настоящее время используется код BUOY).
Во-вторых, эта информация контролируется только по формальным признакам и содержит немалое число ошибок.
В-третьих, по каналам ГСТ передается информация далеко не со всех дрейфующих буев, некоторые руководители океанографических программ по разным причинам передают дрифтерную информацию в национальные центры сбора позже и поэтому в оперативном режиме поступает информация только с 60 - 80 % буев.
В-четвертых, информация передается с буев разного типа, в том числе с массивных метеорологических буев, неподходящих для адекватного отслеживания течений, а также с заякоренных стационарных буев.
В-пятых, из-за ограничения на объем единовременно передаваемой информации по каналам ГСТ некоторая часть данных вовсе не передается (например, не передается информация о наличии/отсутствии у буя подводного паруса).
В-шестых, оказались доступными более качественные архивные дрифтерные данные.
2.2. Данные, подготовленные в рамках международных программ
Такие данные были подготовлены в рамках международной программы исследования глобальной океанской циркуляции WOCE путем соглашения между Marine Environmental Data Service - MEDS (Оттава, Канада) и Atlantic Oceanographic and Meteorological Laboratory - AOML (Майами, США). В соответствии с указанным соглашением поступающие по спутниковой системе радиосвязи ARGOS данные дрифтерных наблюдений (являющиеся более полными, чем передаваемые по каналам ГСТ) проходят первичную обработку и контроль качества в Майами. Данные о местоположении дрифтера и об измеренной температуре воды, не соответствующие некоторым простым, но тщательно подобранным критериям непрерывности, отбраковываются из временных рядов различных параметров, образуемых на основе поступающих телеграмм. Данные с тех типов дрейфующих буев, которые в результате многочисленных проведенных полевых тестовых экспериментов признаны подходящими для отслеживания поверхностных
течений, интерполируются с шестичасовым интервалом и раз в полгода передаются в MEDS для архивации и публичного распространения.
На основе указанных архивных данных в 1998 г. в рамках программы WOCE был подготовлен компакт-диск (именуемый далее как компакт-диск WOCE-1998), который содержит дрифтерную информацию с 1979 по 1996 г. Летом 1999 г. упомянутый диск был прислан авторам настоящей работы по их запросу. Содержащиеся на компакт-диске WOCE-1998 данные плюс данные по декабрь 1998 г. включительно, которые в настоящее время находятся в открытом для копирования доступе на WWW-сервере MEDS, и составили исходный исторический массив данных дрифтерных наблюдений, формируемый в ГОИНе. К этому массиву были добавлены поступившие в ГОИН ранее на компакт-дисках дрифтерные данные, полученные в рамках проекта TOGA за 1985-1990 г. (имеющиеся на них данные, подготовленные в MEDS, не полностью пересекаются с данными, подготовленными в 1998 г. в AOML).
Одним из недостатков программы контроля и архивации дрифтерных данных в рамках соглашения между MEDS и AOML является то, что подготовленные данные в дальнейшем пополняются только новыми данными, поступающими по системе ARGOS, в чем можно убедиться, сравнив дрифтерные данные 1979-1996 гг. на компакт-диске WOCE-1998 с данными, находящимися в настоящее время на WWW-сервере MEDS. Однако, как уже отмечалось, ряд собственников дрифтерных данных могут передавать свои данные в центры сбора только по истечении значительного времени (годы). Даже в рамках международного проекта WOCE руководитель того или иного подпроекта имеет право в течение двух лет не предоставлять полученные данные научной общественности. В настоящее время в Интернете доступны не имеющиеся на компакт-диске WOCE-1998 данные дрифтерных наблюдений по ряду экспериментов, например, в рамках программ IABP (исследования Арктики), FASTEX (Северная Атлантика), SEQUAL и FOCAL (Тропическая Атлантика). Эти данные были переписаны авторами настоящей работы для пополнения базового массива. В будущем предполагается пополнение базового массива и за счет некоторых других архивных данных. С конца 1991 г. Национальный институт океанографии (НИО) Индии осуществляет собственную дрифтерную программу исследования поверхностных течений Индийского океана. В рамках этой программы в период с 1991 по 1997 г. было запущено 39 дрифтеров. Первоначально дрифтеры закупались в США, однако в настоящее время налажено собственное их изготовление. По договоренности с индийской стороной ожидается передача в ГОИН полученных по этой программе данных в рамках совместного (ГОИН - НИО) проекта исследования циркуляции Индийского океана.
Помимо содержащихся в базовом массиве данных о местоположении дрифтера и об измеряемых параметрах среды представляют интерес также и дрифтерные мета-данные. Такая, заимствованная с серверов AOML и MEDS дополнительная информация о всех дрифтерах, данные которых содержатся в массивах AOML и MEDS, включает следующие сведения. Это сведения о моменте потери подводного паруса - драги согласно показаниям специального датчика наличия/отсутствия драги, о моменте падения напряжения на батареях, после которого показания термистора становятся ненадежными, о типе драги, о причине, по которой буй перестал функционировать (например, был выброшен на берег). Данные о наличии драги (стандартом в проекте WOCE является подводный парус с центром на глубине 1 5 м) важны по следующей причине. Ряд сопоставлений оценок скоростей поверхностных течений, полученных по данным дрифтеров с заглубленным на 1 5 м парусом и не имеющим такового, показал, что дрифтеры без паруса могут перемещаться быстрее дрифтеров с парусом на величину, составляющую до 10 см/с (Reverdin G. et al. Seasonal variability in the surface currents of the equatorial Pacific. - J. Geoph. Res. 1994, vol.99, No.C10, pp.20323-20344). Это происходит скорее не из-за возможного резкого вертикального градиента скорости в дрейфовом слое, а за счет прямого действия поверхностных волн и ветра на поплавок дрифтера - искажающих
факторов, влияние которых минимизируется в случае наличия подводного паруса. Поэтому скорости, определенные по дрифтерным данным без паруса, могут оказаться завышенными.
3. Приведение дрифтерных данных к единому формату и дискретности
Следующим этапом после сбора глобальных архивных дрифтерных данных стал этап их приведения к единому формату и дискретности. При этом за основу была принята форма представления данных на компакт-диске WOCE -1998. На этом диске массив дрифтерных данных разделен по годам и по четырем океанам на 40 файлов, представленных в зазипованном виде. Ниже в качестве примера формата записи дрифтерных данных на компакт-диске WOCE-1998 приведена одна из строк файла 1995_arc.lis, в котором содержатся данные, полученные с буев, дрейфовавших в Арктическом бассейне в 1995 г.
0000467,1995,01,01,1200,69.512,-15.659,6.5060,0.0644,-0.0042,0.0062,0.0080,7129,64931
Вначале указан номер буя (467 в рассматриваемом примере) или идентификатор платформы по классификации ARGOS, который может быть семизначным. Далее следуют соответственно год, месяц, число, время (по Гринвичу). Приведенные выше строки соответствуют 12 часам 00 минутам 1 января 1995 г. Далее приведены широта (положительная на север) и долгота (положительная на запад от Гринвичского меридиана) местоположения буя на указанный впереди момент времени. Как широта, так и долгота даны с точностью до тысячных долей географического градуса. Далее следуют значения температуры, которые приведены с точностью до тысячных долей градуса Цельсия. Далее приведены компоненты скорости буя соответственно на параллель и меридиан, вычисленные по последовательным координатам и времени. Эти компоненты соответствуют компонентам лагранжевой скорости поверхностного течения. Знак минус соответствует направлениям на запад и юг. Размерность значений компонент скорости - м/с. Далее располагаются оценки ошибки определения местоположения буя в долях географического градуса соответственно по широте и долготе (описание расчета упомянутых ошибок представлено в статье Hansen, D.V., P.-Marie Poulain, 1996. Quality Control and Interpolations of WOCE/TOGA Drifter Data. J. Atmos. Oceanic Tec., 13,900-909). В тринадцатом поле располагается номер эксперимента и в четырнадцатом - номер буя по классификации ВМО (отсутствие этого номера в записи означает, что данные с рассматриваемого буя не передавались по каналам ГСТ ВМО). Следует отметить отсутствие важного индикатора наличия драги в формате данных WOCE (в формате данных TOGA присутствует колонка с указанием глубины подводного паруса, при этом 0 означает его отсутствие).
Поскольку все полученные из различных источников наборы данных имели различные форматы и за исключением данных TOGA и WOCE (которые также существенно отличались друг от друга по формату записи) были неэквидистантными, они были приведены к описанному выше формату WOCE. Для приведения данных к единой шестичасовой дискретности использовалась линейная интерполяция. При этом пропуски, составляющие более двух суток, не интерполировались. Следует заметить, что для приведения данных к шестичасовой дискретности в AOML используется более сложный алгоритм, описанный в цитированной выше статье (Hansen, D.V., P.-Marie Poulain, 1996).
4. Краткое описание экспериментальной версии базы данных системы
После образования единого массива данных была произведена попытка его препарирования, поиска и отбора данных при помощи MS Access. Однако, как быстро выяснилось, с полученным объемом данных MS Access не справляется. Поэтому далее была осуществлена экспериментальная загрузка массива дрифтерных данных в базу данных (БД), реализованную под управлением СУБД SQL Server (версия 7) с базированием на ОС Windows NT.
Помимо хранилища данных, в котором данные представлены без разбивки по годам, была образована связанная с хранилищем таблица регионов (или океанов), а также таблица, содержащая дополнительную информацию о дрифтерах, включая сведения о моменте потери подводного паруса - драги. При этом в основной формат записи были добавлены два новых поля - код региона и признак качества; в результате число полей увеличилось до 16. Объем таблицы для хранения дрифтерных данных в символьном формате - 986 Mb, а в числовом - 692 Mb. Хранение данных в числовом формате приводит не только к экономии памяти, но также позволяет оптимизировать организацию индексов и поиск данных.
Следует заметить, что рассматриваемая БД была создана в технологических целях исследования эксплуатационных характеристик и поведения на различных типах запросов, выбора стратегии индексирования, а также редактирования исходных дрифтерных данных. В дальнейшем запланировано проектирование, детальная проработка и создание производственной БД с предварительным определением необходимого числа таблиц, полей в каждой таблице и связей между таблицами. Такая БД будет постоянно пополняться за счет новых данных и будет допускать многопользовательский режим работы, а также работу в одном сеансе не с одним, а с различными регионами Мирового океана. Первоначально предполагается организация удаленного доступа к производственной БД по корпоративной сети с клиентских станций, а в дальнейшем организация удаленного доступа к исходным данным и результатам информационных расчетов через Web-сайт в Интернете.
5. Контроль и редакция дрифтерных данных
5.1. Уже при заполнении хранилища выяснилось, что полученные из различных источников записи требуют значительного редактирования. Помимо ожидавшихся при объединении различных массивов дублей были обнаружены дубли иного рода. На компакт-диске WOCE-1998 в файлах 1988_pac.lis, 1989_pac.lis, 1990_pac.lis, 1996_pac.lis и 1995_atl.lis было найдено соответственно 4134, 6560, 794, 72 и 15104 двойных записей, относящихся к буям с номерами 8010, 10037, 10042, 25134, 20657, 20659-20664, 20666, 20669, 20672-20678 и 20680. Каждая пара таких записей соответствует одной и той же точке пространства и одному и тому же моменту времени, что недопустимо. Обе записи в паре могут полностью совпадать или буя, ошибкам определения местоположения буя и номеру буя по классификации отличаться (по неизвестным причинам) в полях, соответствующих температуре, компонентам скорости ВМО. Ниже приведены примеры таких двойных записей:
0020657,1995,07,12,1200,29.326,32.158,26.1060,-0.0176,-0.0141,0.0041,0.0087,9044,13567 0020657,1995,07,12,1200,29.326,32.158,26.1060,-0.0177,-0.0141,0.0048,0.0066,9044,13567
0025134,1996,05,10,0600,22.158,-121.197,25.9460,0.5256,0.6174,0.0000,0.0000,1325,, 0025134,1996,05,10,0600,22.158,-121.197,25.9500,0.5260,0.6176,0.0000,0.0000,1325,21906.
5.2. На том же компакт-диске WOCE-1998, а также в списанных с сервера MEDS файлах с данными за 1 997-98 гг. были обнаружены грубые ошибки, выражающиеся в потере знака при превышении компонентами скорости дрифтера значения -1 м/с по абсолютной величине. Ниже показана последовательность значений (с шестичасовым шагом по времени) зональной компоненты скорости течения, относящаяся к дрифтеру №20118 (файл 1994_pac.lis), в период с 00-00 4.02.1994 по 18-00 12.02.1994:
-0.9523 -0.9611 -0.9693 -0.9789 -1.0000 1.0205 1.0177 1.0093 1.0153 1.0257 1.0314 1.0346 1.0355 1.0345 1.0315 1.0263 1.0183 1.0057 -0.9924 -0.9879 -0.9826 -0.9785 -0.9807 -0.9800 -0.9771 -0.9718 -0.9634 -0.9504 -0.9260 -0.8899 -0.9040 -0.9974 1.0320 1.0036 -0.9998 -0.9988.
Можно видеть 4 скачка скорости: с -1.0000 до 1.0205 м/с 5.02.1994, с 1.0057 до -0.9924 м/с 8.02.1994, с -0.9974 до 1.0320 м/с 11.02.1994 и с 1.0036 до -0.9998 м/с 12.02.1994. Судя по соответствующей последовательности изменений долготы при перемещении рассматриваемого дрифтера, очевидно, что указанные скачки являются ошибочными. Выполненные с использованием экспериментальной БД запросы выявили свыше полутора тысяч случаев подобных потерь знака. Отсутствие в массиве AOML/MEDS отрицательных значений зональной и меридиональной компонент скорости менее -1 м/с означает, что описываемая грубая погрешность является систематической.
5.3. Наибольшее число из обнаруженных ошибок оказалось связанным с наличием в записях ошибочных и в, частности, постоянных значений температуры. Например, в файле 1994_pac.lis на компакт-диске WOCE-1998 зарегистрированная дрифтером №22270 температура поверхности воды от нормальных для приэкваториальной зоны значений 27-29 градусов Цельсия в заключительный период функционирования дрифтера с 22 по 30.10.1994 г. постепенно упала до отрицательного значения -4.959 градуса Цельсия. Ниже приведена соответствующая последовательность значений температуры (с шестичасовой дискретностью) в период с 00-00 21.10.1994 по 12-00 30.10.1994:
29.1390 27.1790 26.7000 28.1190 27.9510 25.5290 23.5370 23.1650 22.2810 19.7750 17.7820 17.4110 16.5270 14.0210 12.0280 11.6570 10.7720 8.2670 6.2740 5.9030 5.0180 2.5120 0.5200 0.1480 -0.7360 -3.3550 -4.5770 -4.3300 -4.9590 -4.9590 -4.9590 -4.9590 -4.9590 -4.9590 -4.9590 -4.9590 -4.9590 -4.9590 -4.9590.
Подобные указанным ошибочные постоянные значения температуры обычно встречались в конце записей, но иногда и в середине.
5.4. Контроль и определение качества записей производились путем последовательной выборки данных из БД и просмотра файлов записей, относящихся к различным регионам Мирового океана. В результате просмотра ошибочные записи отбраковывались или исправлялись с присвоением им соответствующего признака качества данных. Записям с сомнительными значениями температуры присваивался признак -1 . Для зональной компоненты скорости признак низкого или сомнительного качества был определен как -2, а для меридиональной - как -3. При этом нулевой признак качества отвечает записям, в которых ошибки не были найдены. Дублирующие записи были исключены, а потерянные знаки у компонент скорости - восстановлены.
В качестве средства визуального контроля данных была использована созданная авторами, работающая на базе ОС MS-DOS программа визуализации дрифтерных траекторий, а также графиков компонент, модуля скорости и температуры воды вдоль выбранной траектории.
На рис. 3 показан пример записи, полученной с дрифтера #19161, с ошибочными постоянными значениями температуры в период с 18.08. по 29.10.1997 г. Значительные всплески скорости дрифтера, достигавшие 2,09 м/с (31.07.1997 г.), вероятно, реальны и обусловлены перемежаемыми инерционными движениями воды.
Рис. 3. Пример визуализации графиков модуля скорости M и температуры T вдоль дрифтерной
траектории.
5.5. В результате работы по массовому контролю качества собранного массива данных глобальных дрифтерных наблюдений за поверхностными течениями и систематизации обнаруженных ошибок авторами был подготовлен соответствующий документ, который был отослан в Англию, в штаб-квартиру программы WOCE (WOCE IPO). Рассматриваемый документ был помещен дирекцией программы на официальный сервер штаб-квартиры, а также обсуждался на заседании Комитета по подготовке данных программы WOCE. Одним из итогов обсуждения стало принятие решения о дополнительной редакции и замене находящихся в открытом доступе на WWW-сервере MEDS дрифтерных данных с учетом сделанных авторами настоящей работы замечаний.
6. Краткое описание разработанной программы визуализации
При работе с массивом дрифтерных данных упомянутая в разделе 5.4 программа визуализации предоставляет следующие возможности просмотра:
- всей совокупности дрифтерных траекторий (лагранжевых траекторий течений) в заданном районе Мирового океана;
- всех траекторий буев, дрейфовавших с поверхностными течениями в заданном районе Мирового океана в каком-либо задаваемом месяце года;
- индивидуальной траектории буя, перемещавшегося в заданном районе Мирового океана;
- графиков компонент скорости на параллель U и меридиан V вдоль выбранной индивидуальной дрифтерной траектории;
- графиков модуля скорости M и температуры T вдоль выбранной дрифтерной траектории.
Выбор перечисленных опций и переходы между ними осуществляются при помощи появляющихся на экране ПК функциональных кнопок, управляемых мышью. После "нажатия" кнопки меняют свой первоначальный цвет.
Для просмотра индивидуальной дрифтерной траектории, а также графиков компонент U, V и модуля M скорости и температуры T вдоль этой траектории следует выделить ее графическим курсором на сводной карте траекторий, после чего программа визуализации открывает новый экран с разрисовкой выбранной траектории; далее можно по выбору перейти к экранам с графиками U, V или M, T. Для выделения следует подвести курсор к нужному месту экрана и нажать на клавишу "мыши".
Для просмотра всех траекторий буев, дрейфовавших в каком-либо месяце года, следует "нажать" на одну из 12-ти графических кнопок, расположенных в правой части экрана. Конечное положение каждого из дрифтеров на "месячной" карте траекторий отмечено крупной точкой. Находясь в опции "месячная карта", можно также просмотреть графики компонент U, V и модуля M скорости и температуры T уже вдоль выбранного "месячного" участка какой-либо траектории, для чего на "месячной" карте следует выделить эту траекторию графическим курсором.
Программа визуализации работает на персональном компьютере типа IBM PC с процессором 80486 или выше в среде операционной системы DOS 3.30+.
Описанная программа является прообразом разрабатываемых авторами графических средств визуализации хранящихся в БД данных о течениях.
7. Некоторые параметры подготовленного массива данных
Сформированный и отредактированный с помощью образованной в экспериментальном порядке БД массив данных глобальных дрифтерных наблюдений за поверхностными течениями Мирового океана охватывает период с 1979 по 1998 г. и насчитывает 6328307 записей скоростей течений. Из них 94% записей содержат также данные о температуре воды в поверхностном слое; 81% записей имеют номера ВМО. Массив содержит записи, полученные с 4238 дрифтеров, из них с номерами ВМО - 2739. Последнее означает, что только с 65% дрейфовавших буев информация поступала в оперативном режиме по каналам ГСТ. Число записей, соответствующих одному номеру дрифтера, варьирует в широком диапазоне от 2 до 9382 и в среднем составляет 1493. Вне БД без учета метаданных массив записей дрифтерных наблюдений, представленных в описанном в разделе 3 формате, имеет объем 540 МЬ. Наблюдения посредством дрифтеров неравномерно распределены по акватории Мирового океана. Наибольшая плотность наблюдений соответствует Тихому океану. Некоторые районы океана вовсе не покрыты дрифтерными наблюдениями. На рис. 4 приведена карта Мирового океана, на которую нанесены с помощью описанной в разделе 6 программы визуализации траектории всех буйков, дрейфовавших в океане в период с 1979 по 1998 г. (траектории изображены желтым цветом).
Рис. 4. Освещенность Мирового океана дрифтерными наблюдениями (построено по данным за период с 1 979 по 1 998 г.).
Во временных рядах измерений имеются пропуски, в периоды которых информация с дрифтера по разным причинам не поступала. В рассматриваемом массиве данных пропуски не маркируются. Продолжительность наблюдений с помощью буйка-дрифтера в среднем составляет 1 год, хотя по отдельным экземплярам превышает 6 лет.
Настоящая работа выполнена при финансовой поддержке РФФИ. Код проекта #99-07-90160.
