Негидростатическое моделирование атмосферной динамики на основе полулагранжевого метода

Изменение климата является глобальным процессом, тогда как его последствия — подъем уровня моря, уменьшение размеров ледников, изменение ареалов распространения растений и животных, оттаивание вечной мерзлоты — носят региональный и местный характер. В последние годы некоторые из этих последствий стали вполне ощутимы. Научное мнение, выраженное Межгосударственной группой экспертов по изменению климата 1 (МГЭИК) ООН, заключается в том, что средняя интегральная температура по Земле поднялась на 0,7° С по сравнению со временем начала промышленной революции (т. е. со второй половины XVIII века) [129]. В частности, значительная доля потепления, наблюдавшегося в последние 50 лет, вызвана деятельностью человека и, в первую очередь, выбросом газов, вызывающих парниковый эффект, таких как углекислый газ (СО2) и метан (СН4).

Научные исследования, мониторинг и предсказание климата, по мнению Всемирной Метеорологической Организации (ВМО), являются ключевыми элементами смягчения влияния последствий изменения климата на жизнедеятельность человека. Согласно резолюции международного саммита по моделированию климата в XXI веке [128], состоявшемся в 2008 году, основная задача научного сообщества заключается в консолидации усилий по созданию системы численного прогнозирования климата с высоким пространственным разрешением. Развитие высокопроизводительных компьютеров с параллельной.

Фоль группы экспертов МГЭИК, учрежденной Всемирной Метеорологической Организацией (ВМО) и Программой Организации Объединенных Наций по окружающей среде (ЮНЕП), состоит в сбалансированном представлении существующих точек зрения на основе научно-технической и социально-экономической информации об изменении климата, полученной со всего мира. архитектурой позволяет уже сейчас разрабатывать модели климата с горизонтальным разрешением в несколько километров.

Уравнения гидротермодинамики модельной атмосферы должны быть согласованы с выбранным разрешением. До недавнего времени в большинстве моделей атмосферы использовалось гидростатическое приближение («приближение длинных волн») [9, 4], которое предполагает малость вертикальных ускорений по сравнению с ускорением свободного падения. Подобное допущение может нарушаться при шаге сетки по горизонтали 10 км и менее. Более детализированное представление орографии в численной модели (за счет увеличения разрешения) может приводить в существенному росту значения вертикального ускорения, развиваемого в процессе обтекания гор с крутыми склонами. Гидростатическое приближение в таком случае может нарушаться.

Другим важным процессом в атмосфере является конвекция. В тропиках нагретая поверхность воды приводит к значительному вертикальному перемешиванию воздуха сопровождающимся формированием кучево-дождевых облаков с высоким содержанием влаги. Такие образования с характерным горизонтальным масштабом порядка 10−20 км могут перемещаться под воздействием атмосферных потоков на большие расстояния вплоть до средних широт. Учитывая важную роль кучево-дождевых облаков в глобальной циркуляции влаги в резолюции, принятой на WMP (WCRP Modelling Panel) в 2008 году, рекомендуется описывать осадки конвективного характера на основе явного алгоритма, не прибегая к параметризациям.

Первым примером климатической модели с описанием отдельных кучево-дождевых образований является разработанная в Японии модель NIC AM [104]. Негидростатические уравнения гидротермодинамики атмосферы в этой модели разрешаются на основе явной схемы по времени с расщеплением по пространству и физическим процессам. Высокое пространственное разрешение NIC AM (3,5 км по горизонтали) практически стирает границу между блоком климатической модели, ответственным за воспроизведение динамики атмосферы, и динамическим блоком 2 современной системы численного прогноза погоды. Поэтому создание блока, воспроизводящего динамику атмосферы, является актуальной задачей для моделирования климата и прогноза погоды. Подробный обзор некоторых современных моделей численного прогноза погоды приводится в разделе 1.2.

Впервые решение задачи воспроизведения атмосферной циркуляции на основе двухслойной, квазигеострофической, полусферной модели [86] было предложено в 1956 году, а уже в начале 1960;х годов появилась первая 9-уровенная модель, основанная на полных (без приближения квазигеострофичности) уравнениях атмосферы [109]. Результаты экспериментов с первой совместной моделью общей циркуляции атмосферы и океана были представлены в 1969 году [73].

В России создание гидродинамических климатических моделей началось в 1970;х годах. По инициативе Г. И. Марчука Отделением океанологии, физики атмосферы и географии АН СССР в 1973 году было принято решение о разработке математических моделей климата, основанных на моделях общей циркуляции атмосферы и океана. Такая модель, базирующаяся на использовании законов сохранения и методов расщепления, была построена в Вычислительном центре СО АН СССР [12]. С тех пор модель (ныне модель климатической системы ИВМ РАН [5]) значительно усложнилась и в настоящее время участвует в международном проекте по сравнению климатических моделей CMIP (Coupled Model Intercomparison Project). Отметим, что в этом.

23десь и далее под динамическим блоком мы будем понимать блок модели, ответственный за воспроизведение динамики атмосферы. проекте принимают участие около 30 климатических моделей, созданных в разных странах.

Одним из основных направлений развития климатической модели ИВМ РАН [5] является существенное повышение разрешения. Один из возможных вариантов нового блока решения уравнений динамики атмосферы в новой модели — это усовершенствованный динамический блок глобальной модели прогноза ПЛАВ [18], использующий полу-лагранжево представление адвекции. Модель ПЛАВ разрабатывается в Институте вычислительной математики (ИВМ РАН) совместно с ГУ «Гидрометцентр России». Решением Центральной методической комиссии по гидрометеорологическим и гелиогеофизическим прогнозам Росгидромета от 29 ноября 2007 года модель ПЛАВ рекомендована к внедрению в оперативную практику.

Увеличение пространственного разрешения является одним из способов повышения качества прогноза. К 2012 году предполагается, что горизонтальное разрешение модели ПЛАВ над территорией России возрастет с 28 до 10 км, а число уровней по вертикали — с 28 до 50. Основу модели составляют уравнения гидротермодинамики атмосферы, записанные в гидростатическом приближении. Увеличение пространственного разрешения модели ПЛАВ приводит к необходимости создания нового негидростатического блока, ответственного за воспроизведение динамики атмосферы.

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

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

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

Вторая проблема, решение которой обсуждается в диссертации, является математической. Это:

• разработка алгоритма построения покрывающей земной шар сетки с квазиоднородным разрешением. Исследование точности решения уравнения переноса и уравнений мелкой воды на сфере для таких сеток.

В качестве сетки для глобальной полулагранжевой модели динамики атмосферы в диссертации предлагается использовать редуцированную сетку, построенную по разработанному автором алгоритму. Число узлов по долготе в такой сетке зависит от широты, что позволяет избежать сгущения узлов в окрестности полюсов, вследствие сходимости меридианов. Предполагается, что редуцированная сетка будет использоваться в модели прогноза погоды ПЛАВ и, возможно, в новой версии модели климатической системы ИВМ РАН.

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

Внутренние гравитационные волны являются одним из погодооб-разующих факторов. Горы на пути воздушного потока со значительным содержанием водяного пара способствуют образованию облачности. Изменение температурного градиента атмосферы по вертикали может приводить к отражению 65−100% энергии волны (в зависимости от ее длины) обратно к поверхности Земли и усилению приземного ветра [59, 62].

Гравитационные волны являются дивергентными и могут оказывать локальное сопротивление тропосферным воздушным течениям, влияя на общую циркуляцию атмосферы [65, 3].

Интерес к внутренним гравитационным волнам усилился с распространением планеризма в начале 20-го века. Оказалось, что наибольшие высоты полетов на планерах (до 15 км) достигаются с подветренной стороны. Отметим, что в некоторых случаях такие колебания способствуют развитию сильной турбулентности, представляющей серьезную опасность для авиации [60].

Первые эксперименты по моделированию обтекания гор воздушным потоком проводились уже в 1949 году [10] в приближении несжимаемой атмосферы. Развитие идеи Н. Е. Кочина разделения атмосферы на слои, каждый из которых характеризуется индивидуальными значениями скорости, температуры, плотности и др., позволяет получать важные научные результаты и в настоящее время [13].

Первый аналитический спектральный анализ характеристик внутренних гравитационных волн на основе уравнений атмосферы в линейном приближении был выполнен в 1941 году [88, 72, 106]. Затем появились модели, в которых нелинейная система уравнений атмосферы для случая однородного по скорости набегающего потока сводится к линейному уравнению [67, 68]. В 1959 году были представлены [105] первые результаты исследований возмущений атмосферы, возникающих при больших числах Рейнольдса.

Появление компьютерной техники ознаменовало новый подход к описанию процесса обтекания на основе моделей циркуляции атмосферы. Первое численное моделирование нелинейных гравитационных волн было выполнено в 1967 году [40]. Исследование было продолжено [44, 41, 75, 35] с использованием двумерных версий трехмерной мезомасштабной модели прогноза погоды [123]. Амплитуда гравитационных волн расходилась с данными наблюдений, что являлось следствием упрощений уравнений гидротермодинамики атмосферы трехмерной модели. Позже, в [60], было указано на важную роль верхнего граничного условия в численном моделировании гравитационных волн.

Первые конечно-разностные модели общей циркуляции атмосферы использовали явные методы. В таком случае максимально допустимый шаг интегрирования по времени определяется частотой самой быстрой волны в системе. Поэтому с начала 1960;х годов широкое распространение получили модели квазисжимаемой атмосферы, в которых малозначимые с метеорологической точки зрения звуковые волны исключались из рассмотрения [84, 27, 66]. С помощью гидростатического приближения удалось существенно упростить численный алгоритм и отфильтровать распространяющиеся в вертикальном направлении коротковолновые колебания. Отметим, что гидростатические уравнения атмосферы включают также горизонтальные гравитационные волны, фазовая скорость которых сравнима со скоростью звуковых колебаний.

В 1976 году для негидростатической сжимаемой модели атмосферы был предложен метод численной фильтрации звуковых волн [118]. Ответственные за звуковые волны линейные слагаемые уравнений модели интегрировались по неявной схеме, а шаг по времени ограничивался быстрыми гравитационными волнами разрешаемыми явно. Работа предшествовала появлению полунеявного подхода [97], основанного на методе расщепления Г. И. Марчука [11].

Полунеявный поход предполагает интегрирование линейных слагаемых по неявной схеме, а нелинейных слагаемых и правых частей уравнений, описывающих стоки и источники соответствующих величин, — по явному алгоритму. Границы разрешаемых частот в таком случае определяются выбором линейного оператора модели. В сочетании с дискретизацией по времени на основе полулагранжева алгоритма [111] метод позволяет использовать достаточно большие шаги по времени по сравнению с явными конечно-разностными схемами. Эта особенность объясняет широкую распространенность полунеявного полулагранжевого подхода в задачах динамики атмосферы.

Применительно к задаче мезомасштабного моделирования в 1978 году в [108] было показано, что полунеявный метод может быть неустойчив при интегрировании высокочастотных гравитационных волн орографической природы с соответственно высоким разрешением по вертикали. Решение этой проблемы было предложено позже в 1981 году в работе Робера [97], после которой полунеявный метод получил широкое распространение. В [31] была продемонстрирована возможность одновременного описания гравитационных и звуковых волн на основе неявной схемы интегрирования по времени. Метод получил развитие в [117], где был успешно реализован полунеявный полулагранжев подход (SISL).

Эффективный метод решения двумерной негидростатической задачи обтекания на основе полунеявного подхода был впервые предложен в [110]. Здесь также указывалось на важность осреднения слагаемых, ответственных за стоки и источники энергии, вдоль лагран-жевых траекторий для случая трехслойной схемы по времени. Трехмерная полулагранжева полунеявная модель сжимаемой негидростатической атмосферы для численного моделирования орографических волн была представлена в [85]. В настоящее время полунеявный по-лулагранжев алгоритм успешно применяется при решении широкого круга физических задач, в том числе при моделировании климата и численного воспроизведения мезомасштабных процессов.

В представленной в диссертации двумерной (в вертикальной плоскости) негидростатической модели для воспроизведения динамики атмосферы применяется полунеявный конечно-разностный подход с по-лулагранжевым представлением адвекции. Алгоритм решения уравнений гидротермодинамики квазисжимаемой атмосферы [124] представляет собой существенно переработанный метод, предложенный в [101]. В качестве граничных условий здесь используются условия непротекания на верхней границе расчетной области и скольжения на поверхности Земли.

Преимущество конечно-разностного подхода, используемого в разработанной модели, заключается в локальности и линейном росте числа операций по мере увеличения разрешения по одной координате. Однако, конечно-разностная трехмерная модель атмосферы на регулярной широтно-долготной сетке, помимо вносимой фазовой ошибки, имеет существенный недостаток: неоднородность шага сетки по долготе. Например, при разрешении по широте 10 км шаг сетки по долготе в окрестности полюса составляет около 150 м. Этот недостаток приводит к существенному ограничению на шаг по времени в эйлеровых моделях (определяемый критерием Куранта), неоправданным затратам на расчет «лишних» узлов сетки (до 20%). Для решения упомянутых проблем в настоящей диссертации предлагается использовать редуцированную сетку, в которой число узлов по долготе зависит от широты. Шаг сетки по долготе при этом постоянен для каждого значения широты.

Редуцированные сетки широко применяются в спектральных моделях атмосферы [51, 82, 127]. В конечно-разностных моделях в настоящее время используются регулярные широтно-долготные сетки, икосаэдральные сетки, сетки типа «кубическая сфера» и др. В данной работе представлен метод численного построения редуцированной сетки для полулагранжевой конечно-разностной модели атмосферы. Достоинство предлагаемого подхода заключается в возможности видоизменения критериев, используемых при построении сетки для учета особенностей модели атмосферы.

Редуцированная сетка характеризуется квазиоднородным разрешением, что позволяет использовать ее в задаче моделирования атмосферной циркуляции и численном прогнозе погоды.

Цель диссертационной работы состоит в решении двух взаимосвязанных проблем:

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

• Разработка алгоритма построения редуцированной широтно-дол-готной сетки для полулагранжевой конечно-разностной модели атмосферы. Исследование точности решения уравнения переноса на сфере в зависимости от относительного уменьшения числа узлов сетки по сравнению с регулярной широтно-долготной сеткой. Сравнение точности решений уравнений мелкой воды на сфере на редуцированной и регулярной сетках.

Научная новизна результатов диссертационной работы.

• Предложен усовершенствованный метод численного решения негидростатических уравнений гидротермодинамики квазисжима-емой атмосферы, основанный на работе [101]. Обсуждаемый в диссертации алгоритм обладает большим запасом устойчивости по сравнению с [101] за счет использования двуслойной по времени полулагранжевой схемы SETTLS [50]. Метод расчета геопотенциала, предложенный в диссертации, имеет меньшую численную ошибку.

• Насколько известно автору, в работе выполнено наиболее полное сравнение характеристик орографически возбуждаемых волн, рассчитанных моделью квазисжимаемой атмосферы и моделями, основанными на уравнениях сжимаемой атмосферы. Показано, что процесс обтекания изолированной горы, отношение высоты которой к ее полуширине не превышает 0.25, с хорошей степенью точности может быть описан на основе негидростатических уравнений гидротермодинамики квазисжимаемой атмосферы.

• Разработан алгоритм построения оптимальной редуцированной широтно-долготной сетки для глобальной полулагранжевой конечно-разностной модели атмосферы.

• На основе численных экспериментов показано, что точность решения уравнения переноса (адвекции) на сфере, рассчитанного с помощью полулагранжевого метода, является монотонной функцией относительного уменьшения числа узлов оптимальной редуцированной сетки по сравнению с регулярной сеткой.

Научная и практическая значимость.

• Повышение пространственного разрешения является одним из важных направлений развития климатических моделей. Трехмерная версия представленной в настоящей диссертации негидростатической модели высокого разрешения может быть использована в качестве блока, ответственного за воспроизведение динамики атмосферы, в новой версии климатической модели ИВМ РАН и модели прогноза погоды ПЛАВ.

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

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

Личный вклад автора.

• Разработка и реализация алгоритма численного решения негидростатических уравнений квазисжимаемой атмосферы на основе метода, предложенного в [101].

• Воспроизведение процесса обтекания изолированной горы потоком невязкого газа на основе реализованной модели. Сопоставление результатов численных расчетов с аналитическим решением и результатами других авторов.

• Разработка и реализация алгоритма построения редуцированной широтно-долготной сетки для полулагранжевой конечно-разностной модели атмосферы.

• Исследование точности решения уравнения переноса в зависимости от относительного уменьшения числа узлов оптимальной редуцированной сетки по сравнению с регулярной широтно-долготной сеткой.

• Расчет редуцированных сеток для проведения экспериментов на основе модели мелкой воды на сфере.

Апробация работы. Материалы, вошедшие в диссертационную работу, докладывались и обсуждались на научных семинарах Института вычислительной математики, ГУ «Гидрометцентр России». Они докладывались на Весеннем коллоквиуме по моделированию погоды и климата (Триест, Италия, 2005), Международной конференции по вычислительно-информационным технологиям для окружающей среды «CITES-2007» (Томск, 2007), Международной конференции по решению уравнений в частных прозводных на сфере (Экзетер, Англия, 2007), Семнадцатой международной конференции по негидростатическому моделированию (Бад-Орб, Германия, 2007), Международной конференции по наблюдению, моделированию и информационным системам для окружающей среды «ENVIROMIS-2008» (Томск, 2008), Научной школе-семинар «Современные технологии прогнозирования погоды» (Москва, 2008), 49, 50 и 51-й научной конференциях МФТИ (Москва, 2006;2008), Научной конференции посвященной 175-летию гидрометслужбы России (Москва, 2009), Международной конференции по вычислительно-информационным технологиям для окружающей среды «CITES-2009» (Красноярск, 2009).

Полностью диссертация докладывалась на семинарах ИВМ РАН и ГУ «Гидрометцентр России». По теме диссертации опубликовано 2 статьи в реферируемых журналах ([20, 21]), рекомендованых ВАК РФ, 1 статья в реферируемом журнале [19], б работ в сборниках тезисов.

Диссертация состоит из Введения, трех глав, Заключения, Списка литературы и Приложения. Диссертация содержит 15 рисунков и 4 таблицы. Каждая глава разбита на разделы.

3.7. Выводы.

В главе 3 представлен разработанный автором метод построения оптимальной редуцированной широтно-долготной сетки для полулагранжевой конечно-разностной модели атмосферы. Такая сетка характеризуется квазиоднородным разрешением на сфере, что позволяет использовать ее в задаче моделирования атмосферной циркуляции и численном прогнозе погоды.

Метод построения редуцированной сетки использует критерии, основанные на среднеквадратической интегральной ошибке интерполяции некоторой наперед заданной функции. Сетку с минимальным числом узлов удовлетворяющая критериям мы называем оптимальной.

Сокращение числа узлов сетки сопровождается ухудшением точности численного решения уравнения переноса на сфере вследствие увеличения ошибки интерполяции. Приведённые в разделе 3.5 результаты расчётов показывают, что при использовании оптимальных сеток точность численного решения является монотонной функцией параметра вф критерия 1, однозначно определяющего редуцированную сетку. Таким образом, приведённые на рис. 3.1а-б зависимости позволяют определить сетку, отвечающую заданному соотношению между числом узлов и точностью приближённого решения уравнения переноса на сфере.

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

Наличие дополнительных источников погрешностей таких, как неадвективные слагаемые, операторы осреднения некоторых величин по долготе и прочее, наряду с ошибкой интерполяции может привести к тому, что сетка не будет удовлетворять заданному соотношению между точностью решения уравнений гидротермодинамики атмосферы и числом узлов. Однако, основной целью редукции является разрежение приполярных областей, необходимой для корректной работ параметризаций процессов подсеточного масштаба. В современных моделях прогноза погоды расчёт параметризаций занимает 60% — 75% всего времени вычислений. В климатических моделях — еще больше. Уменьшение общего числа ячеек при переходе на редуцированную сетку приведёт к сокращению вычислительных затрат и в блоке переноса и в блоке параметризаций, что скажется на общем времени интегрирования модели.

Полученные сетки предполагается использовать в климатической модели ИВМ РАН нового поколения и глобальной полулагранжевой модели прогноза погоды ПЛАВ [18], допускающих применение редуцированных сеток.

Заключение

.

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

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

Трехмерный вариант модели может быть использован в качестве блока, ответственного за воспроизведение динамики атмосферы, в рамках модели климатической системы ИВМ РАН [5] нового поколения и модели краткосрочного прогноза погоды ПЛАВ [18].

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

Проведенные эксперименты по моделированию динамики двумерных течений с сильно неоднородными полями скорости демонстрируют возможность сокращения числа узлов сетки по рассматриваемому в данной работе алгоритму на 15−20%. Точность численного решения уравнений мелкой воды и уравнения переноса (адвекции) на сфере при этом изменяется не существенно (менее 10%).

Приведем основные результаты работы, являющиеся одновременно и положениями, выносимыми на защиту:

1. Создана вычислительно эффективная двумерная в вертикальной плоскости негидростатическая модель квазисжимаемой атмосферы. Предложен усовершенствованный метод численного решения негидростатических уравнений гидротермодинамики квазисжимаемой атмосферы, основанный на работе [101]. Обсуждаемый в диссертации алгоритм обладает большим запасом устойчивости по сравнению с [101], за счет использования двуслойной по времени полулагранжевой схемы SETTLS [50]. Метод расчета геопотенциала, предложенный в диссертации, имеет меньшую численную ошибку.

2. Результаты численного воспроизведения орографических волн демонстрируют качественное и количественное согласие с работами других авторов, выполненных с использованием моделей сжимаемой атмосферы [53, 85, 130]. Шаг по времени рассматриваемой модели при этом выше используемого в указанных работах.

3. Разработан алгоритм построения оптимальной редуцированной широтно-долготной сетки на сфере для глобальной полулагранжевой конечно-разностной модели атмосферы. На основе численных экспериментов в диссертации показано, что точность решения уравнения переноса (адвекции) на сфере, рассчитанного с помощью полулагранжевого метода, является монотонной функцией относительного уменьшения числа узлов оптимальной редуцированной сетки по сравнению с регулярной сеткой. На основе общепринятых тестов сделан вывод о несущественном ухудшении точности (менее 10%) численного решения уравнений мелкой воды на сфере при сокращении числа узлов сетки на 15−20%.