На главную | База 1 | База 2 | База 3

РД 52.04.651-2003

РУКОВОДЯЩИЙ ДОКУМЕНТ

МЕТОДИЧЕСКИЕ УКАЗАНИЯ

Первичная обработка
результатов судовых измерений
метеорологических, актинометрических
и оптических величин

Предисловие

1 РАЗРАБОТАН государственным учреждением «Главная геофизическая обсерватория им. А.И. Воейкова» (ГУ ГГО) Росгидромета, Гидрометеорологическим научно-исследовательским центром Российской Федерации (Гидрометцентр России)

2 РАЗРАБОТЧИКИ В.Ю. Окоренков (руководитель разработки), Р.Г. Тимановская (ответственный исполнитель разработки), Г.П. Резников, В.В. Рудометкина (ГУ ГГО), Р.С. Фахрутдинов (Гидрометцентр России)

3 УТВЕРЖДЕН Руководителем Федеральной службы России по гидрометеорологии и мониторингу окружающей среды (Росгидромет)

4 ЗАРЕГИСТРИРОВАН ЦКБ ГМП за номером 52.04.651-2003

5 ВВЕДЕН ВПЕРВЫЕ

6 Дополнения и возможные уточнения приводимых в настоящих методических указаниях методов обработки и расчетов возможны только с разрешения разработчика настоящих методических указаний

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

СОДЕРЖАНИЕ

1 Область применения. 4

2 Нормативные ссылки. 4

3 Термины и определения. 4

4 Сокращения. 4

5 Общие положения. 5

6 Первичная обработка результатов измерений атмосферного давления и барической тенденции. 6

6.1 Атмосферное давление. 6

6.2 Барическая тенденция. 7

7 Первичная обработка результатов измерений скорости и направления ветра. 8

7.1 Метод расчета скорости и направления истинного ветра по данным измерений скорости кажущегося ветра и его направления. 8

7.2 Метод расчета скорости ветра (кажущегося или истинного), измеряемой анемометрами, снабженными шкалой. 9

8 Первичная обработка результатов измерений температуры воздуха, воды.. 9

8.1 Обработка результатов измерений по жидкостным термометрам температуры воздуха, воды и температуры смоченного термометра. 9

8.2 Обработка результатов измерений температуры воздуха, воды по электрическим термометрам.. 10

9 Методы расчета величин, характеризующих влажность воздуха. 10

9.1 Перечень величин, характеризующих влажность воздуха. 10

9.2 Методы расчета парциального давления водяного пара. 10

9.3 Метод расчета давления насыщенного водяного пара. 12

9.4 Методы расчета точки росы и точки льда (инея) 12

9.5 Методы расчета относительной влажности воздуха. 13

9.6 Метод расчета дефицита насыщения. 13

10 Первичная обработка результатов измерений параметров обледенения судна. 13

11 Первичная обработка результатов визуальных наблюдений за гидрометеорологическими величинами. 14

12 Первичная обработка результатов измерений составляющих радиационного баланса. 14

12.1 Основное выражение для расчета составляющих радиационного баланса. 14

12.2 Метод расчета прямой солнечной радиации на горизонтальную поверхность. 15

12.3 Метод расчета коротковолнового альбедо. 15

12.4 Методы расчета часовых сумм составляющих радиационного баланса. 15

12.5 Методы расчета радиационного баланса подстилающей поверхности по измерениям его коротковолновой и длинноволновой составляющих. 16

12.6 Метод пересчета составляющих радиационного баланса, представленных в разных единицах измеряемой величины.. 17

13 Методы расчета астрономических величин. 17

13.1 Методы расчета местных среднего, истинного солнечного времени и уравнения времени. 17

13.2 Метод расчета высоты Солнца. 17

13.3 Метод расчета склонения Солнца. 18

14 Методы расчета оптических величин вертикального столба атмосферы.. 18

14.1 Перечень оптических величин вертикального столба атмосферы, определяемых по данным измерений прямой солнечной радиации и данным аэрологического зондирования атмосферы.. 18

14.2 Метод расчета коэффициента прозрачности атмосферы.. 18

14.3 Метод расчета натурального коэффициента ослабления прямой солнечной радиации аэрозолем над океаном.. 19

14.4 Метод расчета общего и общего приведенного влагосодержания вертикального столба атмосферы по данным аэрологического зондирования. 20

Приложение А (обязательное) Общая формула перевода данных измерений из единиц промежуточной величины в единицы измеряемой величины.. 21

Приложение Б (справочное) Дополнительные рекомендации к методу расчета направления истинного ветра. 22

Приложение В (обязательное) Методы расчета характеристик влажности воздуха при отрицательных значениях температуры воздуха и отсутствии информации об агрегатном состоянии воды на резервуаре смоченного термометра. 22

Приложение Г (справочное) Оценки общего влагосодержания вертикального столба атмосферы высотой 10 км по данным измерений характеристик влажности воздуха в приводном слое атмосферы.. 24

Приложение Д (обязательное) Тестовый пример расчета общего и общего приведенного влагосодержания вертикального столба атмосферы.. 25

Библиография. 26

Введение

В 1986 г. на НИС и НИСП Росгидромета был внедрен разработанный в ГУ ГГО комплекс программ (КСМАП) для ЭВМ типа ЕС и СМ по обработке результатов судовых метеорологических и актинометрических наблюдений. Он позволял на всех судах, где производились гидрометеорологические наблюдения штатными наблюдателями, унифицировать и автоматизировать процесс первичной обработки, контроля и архивации результатов этих наблюдений [1].

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

В последние годы парк ЕС и СМ ЭВМ заменен на ПЭВМ. Это обстоятельство обусловило необходимость разработки КСМАП применительно к ПЭВМ. С другой стороны, появилась возможность разработки КСМАП силами УГМС для ПЭВМ разной конфигурации.

Накопленный опыт по разработке и эксплуатации КСМАП для ЕС и СМ ЭВМ показал, что и в случае использования ПЭВМ программное обеспечение процедур первичной обработки результатов судовых наблюдений и расчета ряда характеристик приводного слоя атмосферы, составляющих один или несколько автономных блоков в системе КСМАП, должно также базироваться на единых методах и алгоритмах обработки. Это является одним из требований наставления [2]. Кроме того, в этих методах должны быть приведены в соответствие с существующими государственными и отраслевыми стандартами термины и обозначения геофизических величин с целью ликвидации разночтения и понимания их по разным литературным источникам.

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

Поскольку разработкой программного обеспечения, как правило, занимаются программисты, а не гидрометеорологи, было признано целесообразным обобщить все существующие методы первичной обработки, оформив их в виде руководящего документа, что, с одной стороны, значительно упростит работу с ними, а с другой - позволит на местах (в УГМС) самостоятельно разрабатывать программное обеспечение первичной обработки, контроля и архивации судовой гидрометеорологической, актинометрической информации для ПЭВМ.

Регламентирование использования единых методов первичной обработки результатов метеорологических и актинометрических измерений позволит на судах разных ведомств получать сопоставимые во времени и пространстве данные наблюдений, выполняемых в том числе по программе комплексного мониторинга состояния Мирового океана. В результате будет обеспечена сопоставимость получаемых данных с данными судовых наблюдений, архивированными во ВНИИГМИ-МЦД ранее.

РД 52.04.651-2003

РУКОВОДЯЩИЙ ДОКУМЕНТ

МЕТОДИЧЕСКИЕ УКАЗАНИЯ

Первичная обработка результатов судовых измерений
метеорологических, актинометрических
и оптических величин

Дата введения - 2004-08-01

1 Область применения

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

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

Настоящие методические указания обязательны при обработке на судах или в НИУ и УГМС Росгидромета результатов неавтоматизированных судовых метеорологических, актинометрических, оптических измерений и наблюдений.

2 Нормативные ссылки

В настоящих методических указаниях использованы ссылки на следующие стандарты:

ГОСТ 112-78 Термометры метеорологические, стеклянные. Технические условия

ГОСТ 4401-81 Атмосфера стандартная. Параметры

ОСТ 52.04.10-83 Актинометрия. Термины, буквенные обозначения и определения основных величин

ГОСТ 8.524-85 Таблицы психрометрические. Построение, содержание, расчетные соотношения

ГОСТ 8.567-99 Измерение времени и частоты. Термины и определения

3 Термины и определения

В настоящих методических указаниях применяемые термины, обозначения гидрометеорологических, актинометрических и оптических величин соответствуют представленным в действующих РД [2 - 7], журналах КГМ-15 и УКГМ-15А [3, 6].

4 Сокращения

В настоящих методических указаниях приняты следующие сокращения:

ВНИИГМИ-МЦД - Всероссийский научно-исследовательский институт гидрометеорологической информации - Мировой центр данных.

ГУ ГГО - государственное учреждение «Главная геофизическая обсерватория им. А.И. Воейкова».

КСМАП - комплекс судовых метеорологических и актинометрических программ.

НИС - научно-исследовательское судно.

НИСП - научно-исследовательское судно погоды.

ПАП - первичный актинометрический преобразователь.

ПИП - первичный измерительный преобразователь.

ПЭВМ - персональная электронная вычислительная машина.

Росгидромет - Федеральная служба России по гидрометеорологии и мониторингу окружающей среды.

СРВ - составляющие радиационного баланса.

СГМС - судовая гидрометеорологическая станция.

УГМС - межрегиональное территориальное управление по гидрометеорологии и мониторингу окружающей среды.

5 Общие положения

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

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

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

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

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

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

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

5.2 Первичная обработка результатов судовых гидрометеорологических, актинометрических и оптических измерений и наблюдений в конкретный срок состоит из ряда последовательных операций:

- занесение результатов измерений и наблюдений в журналы (книжки) УКГМ-15А, КГМ-15 или непосредственно в ПЭВМ;

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

- технический контроль результатов измерений и наблюдений;

- критический контроль результатов измерений и наблюдений;

- подготовка оперативных сообщений;

- формирование отчетной документации.

5.3 Результаты измерений и наблюдений заносят в журналы КГМ-15, УКГМ-15А в соответствии с требованиями, представленными в [3, 4, 6] или в самих журналах.

5.4 Данные измерений и наблюдений переводят в единицы измеряемых, наблюдаемых величин в соответствии с разделами 6 - 14 настоящих методических указаний, а в цифры кода - в соответствии с [4].

5.5 Технический и критический контроль результатов метеорологических измерений и наблюдений осуществляют в соответствии с требованиями [8, 9].

5.6 Формирование отчетной документации отчетных форм осуществляют в соответствии с требованиями [3, 6].

5.7 Форматы величин, используемых в разных формулах (форматы значений величин), соответствуют представленным в указанных в подразделе 5.3 журналах; форматы выходных данных после расчетов по методам, рассматриваемым в настоящих методических указаниях, округление рассчитываемых величин соответствуют форматам, представленным в наставлениях [2, 3].

5.8 Контроль за выполнением требований настоящих методических указаний возлагается на методистов-метеорологов, судовых инспекторов УГМС.

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

6 Первичная обработка результатов измерений атмосферного давления и барической тенденции

6.1 Атмосферное давление

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

6.1.2 В соответствии с требованиями [2], данные измерений атмосферного давления Ризм (отсчеты по барометру) должны быть приведены к уровню моря и температуре воздуха 0 °С. Если отсчеты по барометру Ризм выражены в гектопаскалях, то такое приведение осуществляют по формуле

Р0 = Ризм + DРш + DРт + DРу,                                               (1)

а если они выражены в миллиметрах ртутного столба, то по формуле

Р0 = 1,3332(Ризм + DРш + DРт + DРу),                                       (2)

где Р0 - атмосферное давление, приведенное к уровню моря и температуре воздуха 0 °С, гПа;

Ризм - отсчет по прибору (барометру, барометру-анероиду), гПа или мм рт.ст.;

DРш - поправка шкалы к отсчету по прибору, гПа или мм рт.ст. Информация о поправках шкалы (если таковые имеются) приводится в свидетельстве о поверке к прибору в виде таблицы, необходимую поправку рассчитывают методом интерполяции;

DРт - температурная поправка для приведения атмосферного давления к температуре 0 °С, которую рассчитывают по формуле, указанной в свидетельстве о поверке, по температуре воздуха ta, измеренной вблизи прибора, гПа или мм рт.ст.;

DРу - поправка на приведение атмосферного давления к уровню моря, которую рассчитывают по формуле

DРу = DРН(Н + DН)                                                      (3)

(DРН = 0,133 гПа/м или 0,1 мм рт.ст./м (соответствует изменению атмосферного давления на 1 м высоты),

Н - высота установки прибора над уровнем моря для измерения атмосферного давления (отсчитывается от положения максимальной ватерлинии [3]), м. Информация о высоте установки прибора имеется в журналах УКГМ-15А, КГМ-15 [3, 6],

DН - разность уровней «закрытого» моря (типа Каспийского) и Мирового океана, м. Эту разность берут со знаком «плюс», если уровень «закрытого» моря выше уровня Мирового океана, и со знаком «минус», если этот уровень ниже уровня Мирового океана (информация о высоте уровня «закрытого» моря имеется в журнале КГМ-15 или в УГМС. Для открытых морей и акватории океанов DН = 0);

1,3332 - численный коэффициент перевода миллиметров ртутного столба в гектопаскали: 1 мм рт.ст. = 1,3332 гПа.

6.1.3 При расчете значений Р0 по формулам (1) и (2) следует помнить, что:

- все поправки рассчитывают с округлением до 0,1 гПа или до 0,1 мм рт.ст. и берут для расчетов со своим знаком;

- все слагаемые в правой части должны быть выражены в одних единицах измерения (гПа или мм рт.ст.).

6.1.4 Приведем примеры обработки результатов измерений атмосферного давления при плавании судна в океане и в Каспийском море.

Примеры

1 Плавание судна проходит в океане. Атмосферное давление измеряют по анероиду № 392890 (извлечение из свидетельства о поверке представлено в таблице 5 наставления [3]), установленному на высоте Н = 10,1 м над максимальной ватерлинией. Отсчет по анероиду Ризм = 741,9 мм рт.ст., отсчет по термометру в рубке ta = 12,4 °С; значение DРу = 10,1 × 0,1 - 1,01 мм рт.ст. (DН = 0). Из свидетельства о поверке следует, что DРш = -0,6 мм рт.ст., а DРт = 0,3 мм рт.ст.

По формуле (2) рассчитывают значение атмосферного давления, приведенное к температуре 0 °С и уровню моря: Р0 = 990,1 гПа (округление до десятых долей гектопаскаля осуществляют в соответствии с требованиями [2]).

2 Исходные данные те же, что и в примере 1, но плавание проходит в Каспийском море, уровень которого на начало 1994 г. был на 26,8 м ниже уровня Мирового океана. По формуле (3) определяют, что DРу = 0,1×(10,1 - 26,8) - 1,67 мм рт.ст.

Значение P0, рассчитанное по формуле (2), оказалось равным 986,5 гПа.

6.2 Барическая тенденция

6.2.1 Барическая тенденция описывается двумя параметрами - ее величиной DР, отражающей количественное изменение атмосферного давления за 3 ч, предшествующие сроку наблюдения t, ч, и характеристикой а, описывающей качественное изменение атмосферного давления за эти 3 ч [3, 4]. По данным судовых наблюдений барическую тенденцию DР рассчитывают по формуле

                                                        (4)

где  - атмосферное давление, измеренное в срок t и приведенное по формуле (1) или (2) к уровню моря и температуре воздуха 0 °С, гПа;

 - атмосферное давление, измеренное в срок t минус 3 ч (t - 3 ч) и приведенное по формуле (1) или (2) к уровню моря и температуре 0 °С, гПа.

6.2.2 Характеристику барической тенденции а не рассчитывают, а представляют в цифрах кода [4]: кодирование осуществляют по значению DР и визуальной оценке характера изменения давления за 3 ч, фиксируемого на диаграммном бланке барографа [3, 4, 6], идентифицируя его с одним из типичных видов хода давления, приведенных в кодовой таблице [4].

7 Первичная обработка результатов измерений скорости и направления ветра

7.1 Метод расчета скорости и направления истинного ветра по данным измерений скорости кажущегося ветра и его направления

7.1.1 В соответствии с [2, 3, 6], на судне скорость V и направление d истинного ветра не измеряют, а рассчитывают по данным измерений на судне (по анеморумбометрам или другим приборам и способам) скорости Vк и направления dк кажущегося ветра и данным о скорости Vc и компасном курсе dc судна.

7.1.2 Расчет скорости истинного ветра V и его направления d осуществляют по формулам

                                     (5)

d = dс + dк ± dки,                                                          (6)

где Vc - скорость хода судна, уз. Для расчетов берут с округлением до 0,1 уз;

Vк - скорость кажущегося ветра, м/с. Для расчетов берут с округлением до 0,01 м/с;

dс - курс судна, целые градусы;

dк - направление кажущегося ветра, отсчитываемое от компасного курса судна по часовой стрелке, целые градусы;

dки - угол между кажущимся и истинным ветром, отсчитываемый по часовой стрелке, целые градусы. Рассчитывают по формуле

dки = arc cos[(Vк - 0,5144Vc cosdк)/V]*.                                         (7)

______________

* В [3] значение dки предлагается рассчитывать через arcsin[(0,5144V0sindк)/V], который при определенных значениях Vc, V и Vк не всегда имеет однозначное решение (что проиллюстрировано в приложении В), а это усложняет расчеты. По этой причине в настоящих методических указаниях предлагается уточненный метод расчета dки.

Если dк < 180°, то значение d рассчитывают по формуле (6), в которой dки берут со знаком «плюс»; если dк ³ 180°, то значение d рассчитывают по формуле (6), в которой dки берут со знаком «минус».

При расчете скорости и направления истинного ветра по формулам (5) и (6) необходимо учитывать следующие ситуации:

а) V = 0 или Vк = Vс, тогда dк = 0;

б) Vк = 0 или V = Vс, тогда d = dc;

в) Vc = 0 или V = Vк, тогда d = dк.

Примечания

1 Если направление кажущегося ветра определялось по отношению к географическому меридиану dкс (например, по ветровому конусу или судовому компасу [3, 6]), то d в формулах (5) и (6) следует рассчитывать из соотношения: dк = dкс - dс.

2 Если угол dк < 0, то для расчета по формуле (6) значение dк = dкс - dс + 360°.

3 Если рассчитанное значение направления истинного ветра оказалось отрицательным (dрасч < 0), то окончательное значение d = 360° + dрасч; если dрасч > 360°, то окончательное значение d = dpacч - 360°.

4 Если на судне скорость ветра измеряют по анемометру интегрирующего типа**, используя показания по его шкале, отнесенные к интервалу времени от включения до выключения анемометра, расчет скорости кажущегося ветра для определения скорости истинного ветра V (м/с) осуществляют в соответствии с подразделом 7.2.

5 В формулах (5) - (7) значения скорости ветра и хода судна округляют до 0,1 м/с и 0,1 уз соответственно, а направления ветра и курса судна - до 1°.

______________

** Анемометр такого типа, по существу, является счетчиком количества оборотов ветроприемника за конкретный интервал времени.

7.2 Метод расчета скорости ветра (кажущегося или истинного), измеряемой анемометрами, снабженными шкалой

При измерении на судне скорости ветра по анемометрам, снабженным шкалой с делениями, V или Vк определяют по изменению показаний анемометра (дел.) за интервал времени (с), измеренный секундомером, т.е. по промежуточной величине Dn (дел./с). В свидетельстве о поверке анемометра приводится таблица перевода Dn в метры за секунду (таблица 1).

Таблица 1 - Извлечение из свидетельства о поверке к ручному анемометру № 424875

Dn, дел./с

Скорость V, м/с

Dn, дел./с

Скорость V, м/с

1

1,3

6

5,8

2

2,2

7

6,7

3

3,1

8

7,6

4

4,0

9

8,6

5

4,9

10

9,4

...

...

20

18,8

7.2.1 Если связь между Dn анемометра и скоростью ветра по данным из свидетельства о поверке окажется линейной, то скорость ветра V или Vк (м/с), определяют по формуле

V = 0,1(V20 - V10)(nt+Dt - nt)Dt-1,                                          (8)

где V20 - скорость ветра из свидетельства о поверке к анемометру, соответствующая значению Dn = (nt + Dt - nt)Dt-1 = 20 дел./с, м/с. Извлечение из одного такого свидетельства приведено в таблице 1;

V10 - скорость ветра из свидетельства о поверке, соответствующая значению Dn = 10 дел./с, м/с;

nt - отсчет по шкале анемометра перед измерением скорости кажущегося (истинного) ветра до момента включения анемометра, дел.;

nt+Dt - отсчет по шкале анемометра через промежуток времени Dt (с), оканчивающийся в момент выключения анемометра, дел.;

0,1 - коэффициент учитывает то, что разность скоростей ветра в первой скобке делится на 10, с/дел.

7.2.2 Если связь между Dn анемометра и скоростью ветра нелинейная, то значение скорости ветра V определяют методом интерполяции по формуле

V = [(Vj - Vt)/(nj - nt)](nt + Dt - nt)Dt-1,                                         (9)

при j > i Dni < (nt + Dt - nt)Dt-1 < Dnj.

В формуле (9) единицы используемых величин те же, что и в формуле (8).

8 Первичная обработка результатов измерений температуры воздуха, воды

8.1 Обработка результатов измерений по жидкостным термометрам температуры воздуха, воды и температуры смоченного термометра

Жидкостные термометры проградуированы, как правило, в единицах измеряемой величины, т.е. в градусах Цельсия (°С) [2, 3], поэтому обработка результатов измерений сводится к введению поправок из свидетельств о поверке в отсчеты по термометрам. В итоге обработку осуществляют по формуле

ta = Xk + Dtk                                                               (10)

или по формуле

t¢ = X¢k + Dt¢k,                                                          (10а)

где ta - температура воздуха, воды, °С;

t¢ - температура смоченного термометра, °С;

Xk = tk (или X¢k = t¢k) - k-й отсчет по термометру, °С.

Для расчетов берут с округлением до 0,1 °С;

Dtk (или Dt¢k) - поправка для k-го отсчета из свидетельства о поверке, °С. Для расчетов берут с округлением до 0,1 °С (со своим знаком).

8.2 Обработка результатов измерений температуры воздуха, воды по электрическим термометрам

Как правило, если для измерений используют неавтоматические измерительные комплексы, температуру по электрическим термометрам отсчитывают в единицах промежуточных величин по показаниям электроизмерительных приборов xt. В этих случаях обработку результатов измерений температуры осуществляют по формуле А.1 (приложение А) с округлением до 0,1 °С.

9 Методы расчета величин, характеризующих влажность воздуха

9.1 Перечень величин, характеризующих влажность воздуха

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

9.2 Методы расчета парциального давления водяного пара

9.2.1 Парциальное давление водяного пара е рассчитывают по температуре воздуха tа, температуре смоченного термометра t¢ и атмосферному давлению Р, используя основную психрометрическую формулу, представленную в ГОСТ 8.524.

Таблица 2 - Перечень измеряемых или рассчитываемых величин характеризующих влажность воздуха

Наименование, условное обозначение и единицы величин

Определение величины по ГОСТ 8.524

по ГОСТ 8.524

по Психрометрическим таблицам [10]

Точка росы td, °C

Точка росы td, °C

Температура воздуха, при которой водяной пар, содержащийся в воздухе, достигает состояния насыщения при данном давлении водяного пара и неизменном атмосферном давлении

Парциальное давление водяного пара е, гПа

Упругость водяного пара е, гПа

Давление водяного пара, содержащегося в воздухе

Давление насыщенного водяного пара Ew над поверхностью воды, гПа

Максимальная упругость водяного пара (упругость насыщения) над водой Ew, гПа

Давление водяного пара, находящегося в состоянии термодинамического равновесия с жидкой фазой химически чистой воды, при плоской поверхности границ раздела фаз в однокомпонентной системе*

Давление насыщенного водяного пара Et над поверхностью льда, гПа

Максимальная упругость водяного пара (упругость насыщения) Et надо льдом, гПа

Давление водяного пара, находящегося в состоянии термодинамического равновесия с твердой фазой химически чистой воды, при плоской поверхности границ раздела фаз в однокомпонентной системе*

Дефицит насыщения d, гПа

Дефицит влажности d, гПа

Разность между давлением насыщенного пара Ew и парциальным давлением водяного пара е при данной температуре

Относительная влажность воздуха f, %

Относительная влажность воздуха f, %

Отношение парциального давления водяного пара е, содержащегося в воздухе, к давлению насыщенного водяного пара Ew при одинаковых температуре и атмосферном давлении

* Иными словами, вне смеси водяного пара с другими парами и газами. Значение парциального давления насыщенного водяного пара Ес, смешанного с газами, например воздухом, несколько отличается при тех же условиях от значения его в однокомпонентной системе, что учитывается в прецизионных измерениях, но игнорируется в метеорологических наблюдениях.

9.2.1.1 Если на резервуаре смоченного термометра вода* (как правило, при температуре смоченного термометра t¢w > 0 °С), то расчет осуществляют по формуле

e = Ec(t¢w) - kwAP(ta - t¢w)(1 + awt¢w),                                      (11)

______________

* Информация об агрегатном состоянии воды на резервуаре смоченного термометра (вода или лед) должна быть в журнале УКГМ-15А в графе «Индекс смоченного термометра» [6].

где е - парциальное давление водяного пара, гПа. Рассчитывают с округлением до 0,01 гПа;

Ec(t¢w) - давление насыщенного водяного пара, содержащегося в воздухе, гПа. Рассчитывают по формуле (14) с округлением до 0,01 гПа;

t¢w - температура смоченного термометра, когда на резервуаре термометра вода, °С. При расчете берут с округлением до 0,1 °С;

kw - коэффициент, учитывающий агрегатное состояние конденсированной фазы воды на резервуаре смоченного термометра: если вода в жидком состоянии, то kw = 1;

А - психрометрический коэффициент для атмосферного воздуха стандартного состава, (°С)-1 (по ГОСТ 4401). Значение коэффициента для конкретного экземпляра аспирационного психрометра приводится в его паспорте; при отсутствии в паспорте такой информации при расчетах принимают А = Аа = 662×10-6 (°С)-1 для психрометров типа М-36 и МВ-4 с номинальным значением скорости вентиляции 2 м/с (по ГОСТ 8.524) и с термометрами типа ТМ-6 (по ГОСТ 112);

Р - атмосферное давление, определяемое из соотношения: Р = Ризм + DРш, где Ризм - давление, измеренное на уровне установки прибора (барометра) с учетом внесения поправки шкалы DРш, гПа. При расчете берут с округлением до 0,1 гПа;

tа - температура воздуха (по сухому термометру), °С. При расчете берут с округлением до 0,1 °С;

aw - коэффициент, учитывающий зависимость от температуры удельной теплоты фазового перехода конденсированной фазы воды в пар и других величин, входящих в выражение для психрометрического коэффициента, (°С)-1. При расчете принимают aw = 0,00115 (°С)-1 (по ГОСТ 4401).

9.2.1.2 Если на резервуаре смоченного термометра лед, то расчет осуществляют по формуле

e = Ec(t¢i) - ki(tа - t¢i),                                                    (12)

где Ес(t¢i) - давление насыщенного водяного пара, содержащегося в воздухе, гПа. Рассчитывают по формуле (14а) с округлением до 0,01 гПа;

t¢i - температура смоченного термометра, когда на резервуаре смоченного термометра лед, °С. При расчете берут с округлением до 0,1 °С;

ki - коэффициент, учитывающий агрегатное состояние конденсированной фазы воды на резервуаре смоченного термометра: вода в твердом состоянии (лед) - ki = 0,8822.

9.2.1.3 Если информация об агрегатном состоянии воды на резервуаре смоченного термометра в данных наблюдений отсутствует при минус 10,0 °С < ta £ 0 °С, то проверяют выполнение условия t¢ < ta. При этом:

- невыполнение условия означает, что вода на резервуаре смоченного термометра находится в твердом состоянии (лед), и расчет значения е осуществляют по формуле (12). В этом случае, хотя е > Ес(t¢i), но d ³ 0, так как (в соответствии с подразделом 9.6) d = E0(t¢w) - e;

- если условие выполняется, то расчет значения е осуществляют по формуле

е = 0,5(е11 + е12),                                                       (13)

где e11 и e12 - значения е, рассчитанные соответственно по формулам (11) и (12).

9.2.2 Расчет е по температуре воздуха ta и относительной влажности воздуха f осуществляют по следующей схеме:

- по формуле (14) рассчитывают с округлением до 0,01 гПа значение Ew(tа), подставив вместо t¢ значение tа;

- по значениям f и Ew(ta) рассчитывают с округлением до 0,01 гПа значение е по формуле е = 0,01fEw(ta).

9.2.3 Расчет е по точке росы td осуществляют по следующей схеме: по формуле (14) рассчитывают значение Ew(t¢w) с округлением до 0,01 гПа, подставив вместо t¢w значение td, которое будет тождественно значению е: Ew(t¢w) = e.

9.3 Метод расчета давления насыщенного водяного пара

Давление насыщенного водяного пара Ес(t¢) рассчитывают с округлением до 0,01 гПа по данным измерений температуры смоченного термометра t¢:

а) если на резервуаре смоченного термометра вода (как правило, при температуре смоченного термометра t¢w > 0 °С), то расчет осуществляют по эмпирической формуле по ГОСТ 8.524

Ес(t¢)* = Ew(t¢w) = Е0exp[awt¢w/(bw + t¢w)],                                      (14)

_____________

* См. сноску к таблице 2.

где E0 = 6,1121 гПа - давление насыщенного водяного пара при ta = 0 °С;

aw = 17,5043 и bw = 241,2 °С - постоянные для воды;

б) если на резервуаре смоченного термометра лед, то расчет осуществляют по эмпирической формуле по ГОСТ 8.524

Ес(t¢)* = Ei(t¢i) = Е0exp[ait¢i/(bi + t¢i)],                                    (14a)

_____________

* См. сноску к таблице 2.

где ai = 22,4893 и bi = 272,881 °С - постоянные для льда.

9.4 Методы расчета точки росы и точки льда (инея)

9.4.1 Расчет точки росы td и точки льда (инея) ti по данным измерений температуры воздуха ta и температуры смоченного термометра t¢w или t¢i осуществляют с округлением до 0,1 °С.

Схемы расчета точки росы и точки льда следующие.

а) схема расчета td:

- по формуле (11) рассчитывают е;

- это значение е подставляют в формулу

                                                 (15)

где e = ln Е0 = 1,8103;

б) схема расчета ti:

- по формуле (12) или (13) рассчитывают е;

- это значение е подставляют в формулу

                                                (15a)

___________

** Здесь е - парциальное давление водяного пара (согласно таблице 2).

9.4.2 Расчет точки росы td по данным измерений относительной влажности воздуха f и температуры воздуха ta осуществляют с округлением до 0,1 °С по формуле по ГОСТ 8.524

                                   (16)

в которой значение ta для расчета берут с округлением до 0,1 °С, а f - до 1 %.

9.5 Методы расчета относительной влажности воздуха

9.5.1 Расчет относительной влажности воздуха f по температуре воздуха tа и температуре смоченного термометра t¢ осуществляют с округлением до 1 % по формуле

f = 100е/Ew(ta),                                                          (17)

в которой значения е и Ew(tа) для расчета берут с округлением до 0,01 гПа. Значение е рассчитывают по формулам (11), (12) или (13), Ew(tа) - по формуле (14), в которой t¢w заменяют на ta.

9.5.2 Расчет относительной влажности воздуха f по температуре воздуха ta и точке росы td осуществляют с округлением до 1 % по формуле по ГОСТ 8.524

                                        (17а)

где Т = 273,15 + ta, К*;

Td = 273,15 + td, К;

G = 6888,2 К;

D = -5,3627.

9.6 Метод расчета дефицита насыщения

Расчет осуществляют с округлением до 0,01 гПа по формуле

d = Ew(ta) - e,                                                         (18)

где Ew(ta) рассчитывают по формуле (14) с заменой t¢ на ta, а е - в соответствии с подразделом 9.2 - в зависимости от состава измеряемых величин.

10 Первичная обработка результатов измерений параметров обледенения судна

Обледенение судна описывается двумя параметрами:

- толщиной отложения льда Ht в срок наблюдения t, целые сантиметры;

- характеристикой обледенения Rs, отражающей скорость нарастания (таяния) отложения льда в сантиметрах за единицу времени (за час).

Расчет значений Ht и Rs осуществляют по формулам

                                                   (19)

______________

* К - обозначение единицы измерения температуры по абсолютной термометрической шкале Кельвина.

Rs = (Ht - Ht-Dt)/Dt,                                                      (20)

где Ht рассчитывают с округлением до 1 см;

Hi - толщина отложения льда в i-й точке измерения (в соответствии с [3, 6] п ³ 3), мм;

Ht-Dt - толщина отложения льда за предыдущий срок наблюдения t - Dt, см;

Dt - промежуток времени (предшествующий сроку наблюдения t), за который определяют характеристику обледенения, ч.

В соответствии с [3,6], Dt = 6 ч, если Rs < 0,7 см/ч, и Dt = 1 ч, если Rs ³ 0,7 см/ч.

11 Первичная обработка результатов визуальных наблюдений за гидрометеорологическими величинами

В соответствии с [2, 3, 6], производят визуальные наблюдения и оценивают:

- ветровые волны и зыбь;

- гидрометеорологические явления, погоду в срок наблюдения и прошедшую погоду, характеристики облачности;

- метеорологическую дальность видимости.

Первичная обработка результатов наблюдений и оценок в этом случае сводится к их кодированию в соответствии с требованиями действующего на момент наблюдений кода КН-01 [4].

12 Первичная обработка результатов измерений составляющих радиационного баланса

12.1 Основное выражение для расчета составляющих радиационного баланса

Как правило, все СРБ, кроме коротковолнового альбедо Ак (альбедо не измеряют, а рассчитывают в соответствии с подразделом 12.3), измеряют в единицах промежуточных величин, выражая их в делениях шкал измерительных приборов. Методы обработки предполагают перевод результатов измерений из таких единиц в общепринятые единицы измеряемой величины.

Общее выражение для преобразования результата измерений (мгновенного значения) любой СРБ и представления его в общепринятых единицах величины имеет вид

                                                   (21)

где Xi - значение СРБ в i-й момент времени, при этом Xi = Si, Qi, Di, , Bi, ,  (соответственно прямая, суммарная, рассеянная и отраженная коротковолновая солнечная радиация, радиационный баланс подстилающей поверхности, коротковолновая и длинноволновая составляющие радиационного баланса (по ОСТ 52.04.10)), кВт/м2;

a - цена деления шкалы измерительного прибора (диаграммного бланка), Y/дел. (Y - наименование единицы промежуточной величины, измеряемой соответствующим прибором, например милливольты). Значение a приводится в свидетельстве о поверке измерительного прибора;

b - коэффициент перевода цены деления шкалы а измерительного прибора в милливольты на деление (мВ/дел.), если в поверочном свидетельстве размерность а приводится не в мВ/дел. Если a выражено в мВ/дел., тогда b = 1;

Сj - корректирующий параметр, зависящий от типа ПАП (j - тип ПАП) и измерительных приборов (например, азимутальная поправка к показаниям пиранометра): при расчетах обычно типу прибора присваивают конкретный номер, например j = 1, 2, ..., J.

Ni - показание измерительного прибора, деления его шкалы (показание вольтметра, потенциометра, ордината диаграммного бланка и т.д.) в i-й момент времени;

 - значение места нуля измерительного прибора в делениях его шкалы в i-й момент времени;

 - коэффициент преобразования (чувствительность) ПАП (актинометрический, радиационный, радиометрический датчик). Индекс ta означает, что коэффициент, как правило, зависит от температуры воздуха. Эта зависимость приводится в свидетельстве о поверке ПАП. В системе Росгидромета  принято выражать в мВ×м2/кВт (по ОСТ 52.04.10) при температуре воздуха t = 20 °С (K20).

С учетом этого значение  в свидетельстве о поверке ПАП может быть приведено в виде постоянного числа для температуры 20 °С (K20) или в виде функции от температуры воздуха tа. Если в свидетельстве о поверке такая функция отсутствует, то приведение коэффициента K20 к температуре воздуха tа, при которой производится измерение СРБ, осуществляют по формуле

                                               (22)

где  - температурный коэффициент ПАП, который для актинометра типа М-3 равен 0,0008, для пиранометра типа М-80М равен 0,0011, а для балансомера типа М-10 равен 0.

Как правило, выражение (22) используют для приведения коэффициента K20 к температуре воздуха tа при расчетах мгновенных значений и часовых сумм прямой солнечной радиации S.

Расчет мгновенных значений СРБ осуществляют с округлением до 0,01 кВт/м2.

12.2 Метод расчета прямой солнечной радиации на горизонтальную поверхность

Предлагаемый метод расчета используют при определении значений В. Расчет прямой солнечной радиации на горизонтальную поверхность S¢ (кВт/м2) осуществляют с округлением до 0,01 кВт/м2 по формуле

S¢ = S sin hO,                                                            (23)

где hO - высота Солнца, градусы угловые (...°) или радианы (рад). Рассчитывают по формуле (34); значение sinhO определяют с округлением до 0,001.

12.3 Метод расчета коротковолнового альбедо

Альбедо выражают в процентах и его расчет осуществляют с округлением до 1 % по формуле

Aк = 100Rк/Q,                                                           (24)

где Q и Rк - соответственно суммарная и отраженная коротковолновая солнечная радиация, измеренная с точностью до 0,01 кВт/м2.

12.4 Методы расчета часовых сумм составляющих радиационного баланса

12.4.1 Часовые суммы СРБ (МДж/м2) (по ОСТ 52.04.10) измеряют или рассчитывают с округлением до 0,01 МДж/м2; в формулах (25) и (26)  выражено в мВ×м2/Вт.

12.4.2 При снятии показаний приборов с дискретностью Dt = 1 с расчет часовых сумм СРБ  (МДж/м2) осуществляют по формуле

                                    (25)

где X - одна из СРБ (X = S, Q, D, Rк, В, Вд, Вк), для которой рассчитывают часовые суммы;

п - номер часа, начиная с которого рассчитывают часовую сумму (n = 0, 1, 2, ..., 23 ч по местному среднему солнечному времени tм, рассчитываемому по формуле (31));

1 и 3600 - крайние номера секунд в часе.

12.4.3 Если показания приборов фиксируют с шагом Dt ¹ 1 с, расчет  осуществляют по формуле

                                (26)

где J - количество отсчетов за 1 ч: J = 3600/Dt;

Dt - дискретность снятия показаний, с.

12.4.4 Метод расчета часовых сумм СРБ (МДж/м2) по среднечасовым значениям, отсчетов по измерительному прибору Nср, места нуля  и  (чувствительности ПАП), как правило, применяют в ситуациях, когда вручную обрабатывают диаграммные бланки с непрерывной регистрацией СРБ. Расчет  в таких ситуациях осуществляют:

______________

* Поскольку чувствительность отдельных ПАП зависит от температуры воздуха, предусмотрен вариант расчета часовых сумм по среднечасовым значениям  ПАП. Это же замечание касается и места нуля измерительного прибора.

а) если Ncp и  выражены в делениях измерительного прибора, a - в мВ/дел., и значения  - в мВ×м2/Вт, по формуле

                                        (27)

б) если Ncp и  выражены в делениях измерительного прибора, a - в мВ/дел.,  - в мВ×м2/кВт, по формуле

                                           (28)

12.5 Методы расчета радиационного баланса подстилающей поверхности по измерениям его коротковолновой и длинноволновой составляющих

В зависимости от наличия данных измерений указанных величин для расчета В используют разные методы:

а) при наличии данных измерений Q, Rк и Вд расчет осуществляют по формуле

B = Q - Rк - Bд;                                                         (29)

б) при наличии данных измерений S, D, Rк и Вд расчет осуществляют по формуле

B = SsinhO + D - Rк - Вд;                                               (30)

При расчетах значений В по формулам (29), (30) все составляющие должны быть выражены в кВт/м2 с округлением до 0,01 кВт/м2.

12.6 Метод пересчета составляющих радиационного баланса, представленных в разных единицах измеряемой величины

В ранее изданных справочниках, таблицах, монографиях и научных статьях СРБ и их суммы могут быть выражены в единицах, отличающихся от использованных в настоящих методических указаниях. Перевод единиц СРБ в единицы величин, использованные в настоящих методических указаниях, и наоборот, следует осуществлять по соотношениям наставления [11]:

1 Вт = 1 Дж/с;

1 кал/(см2×мин) = 0,698 кВт/м2;

1 кВт/м2 = 1000 Вт/м2;

1 кал/см2 = 0,0419 МДж/м2;

1 кВт/м2 = 0,1 Вт/см2;

1 ккал/см2 = 41,9 МДж/м2;

1 МДж/м2 = 100 Дж/см2;

1 кал/см2 = 4,19 Дж/см2.

13 Методы расчета астрономических величин

13.1 Методы расчета местных среднего, истинного солнечного времени и уравнения времени

Местное среднее солнечное время tм (ч и мин) связано с временем tUTC* (ч и мин) tUTC* º tгр - гринвичское время [5]) соотношением

_____________

* UTC (universal time coordinated) - Международная шкала координированного времени, ГОСТ 8.567 (эквивалентно среднему гринвичскому времени (СГВ)).

tм = tUTC + l/15,                                                        (31)

где l - долгота пункта наблюдения, градусы и минуты угловые (соответственно ...° и ...¢).

Истинное солнечное время tи.о.в. связано с местным средним солнечным соотношением

tи.с.в. = tм + t                    при tм £ 24 ч или                                   (32)

tи.c.в. = tм - 24 + t             при tм > 24 ч**,

______________

** При расчете tм на ЭВМ бывают ситуации, когда рассчитанное значение tм > 24 ч.

где t - уравнение времени, мин. Значение t рассчитывают по формуле

t = 0,0172 + 0,4281cos q0 - 7,3515sin q0 - 3,3495cos 2q0 - 9,3619sin 2q0.     (33)

Здесь q0 = 2pdn/365 рад (p = 3,1415) или q0 = 360dn/365°; dn - номер дня от начала года, изменяющийся от 1 (1 января) до 365 (31 декабря) для обычного года и до 366 для високосного года.

13.2 Метод расчета высоты Солнца

Расчет высоты Солнца hO в конкретной точке Мирового океана для любого момента истинного солнечного времени tи.с.в. конкретного дня N месяца М {hO [tи.с.в.(NM)]} осуществляют с округлением до 0,1° или 0,001 рад по формуле

              (34)

где  - широта пункта наблюдения в момент времени tи.с.в. N-го дня М-го месяца, ...°. Для расчетов используют с округлением до 0,1°;

dNM - склонение Солнца в местный полдень N-го дня М-го месяца, ...°. Для расчетов используют с округлением до 0,1°;

Wk - часовой угол Солнца в радианах или градусах, отсчитываемый от времени наступления истинного полдня (tи.с.в. = 12 ч) в пункте наблюдения, т.е. при tи.с.в. = 12 ч значение Wk = 0.

При расчете Wk следует помнить, что промежутку времени Dt = 1 ч соответствует значение DW¢k = 0,262 рад º 15° (DWk = 15°).

С учетом сказанного:

- при выражении часового угла Солнца в радианах: Wk = 0,262 (tи.с.в. - 12),

- при выражении часового угла Солнца в градусах: Wk = 15,0 (tи.с.в. - 12).

Примечание - Выбор единиц для Wk при расчетах высоты Солнца hO по формуле (34) зависит от выбора единиц измерения j и d (градусы или радианы).

13.3 Метод расчета склонения Солнца

Склонение Солнца d (в радианах) может быть определено из следующего приближения, согласно руководству [5]:

d @ 0,006918 - 0,399912cos q0 + 0,070257sin q0 - 0,006758cos q0 + 0,000908sin 2q0.    (35)

Максимальная ошибка расчета d по формуле (35) составляет 35².

14 Методы расчета оптических величин вертикального столба атмосферы

14.1 Перечень оптических величин вертикального столба атмосферы, определяемых по данным измерений прямой солнечной радиации и данным аэрологического зондирования атмосферы

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

- коэффициент прозрачности атмосферы Р2 (метод расчета представлен в подразделе 14.2);

- натуральный коэффициент ослабления прямой солнечной радиации аэрозолем m¢а (метод расчета представлен в подразделе 14.3).

Исходными данными для расчета перечисленных оптических величин являются:

- результаты измерений прямой солнечной радиации Sизм в момент времени tи.c.в. с округлением до 0,01 кВт/м2;

- координаты пункта измерений S: широта j(tи.с.в.) и долгота l(tи.с.в.), используемые в расчетах с округлением до 0,1° или 0,001 рад;

- высота Солнца hʘ в момент времени tи.с.в. (метод расчета представлен в подразделе 13.2);

- склонение Солнца d (метод расчета представлен в подразделе 13.3),

- данные аэрологического зондирования атмосферы для расчета общего влагосодержания вертикального столба атмосферы W (расчет осуществляют по методу, представленному в подразделе 14.4); если информация о данных аэрологического зондирования атмосферы отсутствует, значение W может быть получено из других источников, например по эмпирическим соотношениям, представленным в приложении Г [12].

14.2 Метод расчета коэффициента прозрачности атмосферы

Коэффициент прозрачности атмосферы Р2 - безразмерная величина. Его рассчитывают с округлением до 0,001 по данным измерений прямой солнечной радиации S на перпендикулярную поверхность, измеренной в условиях открытого диска Солнца [2]).

Расчет Р2 осуществляют по формуле из наставления [11]

Р2 = (Sr/l,367)C,                                                         (36)

где S - измеренная прямая солнечная радиация, кВт/м2;

r - поправочный множитель для приведения значений S к среднему расстоянию между Солнцем и Землей, который рассчитывают по формуле из руководства [5]

r = 1/(1,00011 + 0,034222cos q0 + 0,00128sin q0 + 0,000719cos 2q0 +0,000077sin 2q0),            (37)

где q0 = 2pdn/365 рад или q0 = 360 dn/365°;

dn - номер дня от начала года, изменяющийся от 1 (1 января) до 365 (31 декабря) для обычного года и до 366 для високосного года;

С - поправочный множитель для приведения значений Sr к высоте Солнца hO = 30°, который рассчитывают по формуле

С = (sin hO + 0,205)/1,41.                                                 (38)

14.3 Метод расчета натурального коэффициента ослабления прямой солнечной радиации аэрозолем над океаном

Натуральный коэффициент ослабления прямой солнечной радиации аэрозолем m¢а рассчитывают с округлением до 0,01 по формуле из методических указаний [7]

                  (39)

где  = 0,696 кВт/м2 и  = 0,574 кВт/м2 - подозонные солнечные постоянные для спектральных интервалов Dl1 = 0,3 ... 0,8 мкм и Dl2 = 0,8 ... 4,5 мкм соответственно;

 - интегральная функция пропускания солнечной радиации для рэлеевского рассеяния в зависимости от высоты Солнца hO в i-й момент времени tм. Ее определяют по таблице 3 с округлением до 0,001;

 - интегральная функция пропускания для водяного пара, которую определяют по таблице 4 с округлением до 0,01:  - оптическая масса атмосферы при  (), Wпр - общее приведенное влагосодержание вертикального столба атмосферы (г/см2), которое рассчитывают по методу согласно подразделу 14.4 с округлением до 0,01 г/см2;

 - измеренная прямая солнечная радиация на перпендикулярную поверхность в i-й момент времени tм, кВт/м2. Для расчетов берут с округлением до 0,01 кВт/м2.

Таблица 3 - Значения функции Т(hʘ) для разных высот Солнца hʘ

Высота Солнца hʘ, ...°

Десятки

Единицы

0

1

2

3

4

6

6

7

8

9

30

0,728

0,733

0,737

0,742

0,747

0,752

0,756

0,761

0,766

0,770

40

0,775

0,778

0,780

0,782

0,786

0,788

0,790

0,732

0,795

0,798

60

0,800

0,802

0,803

0,805

0,806

0,808

0,810

0,811

0,813

0,814

60

0,816

0,817

0,818

0,819

0,820

0,822

0,823

0,824

0,825

0,826

70

0,827

0,828

0,828

0,828

0,829

0,830

0,830

0,830

0,831

0,832

80

0,832

0,832

0,832

0,833

0,833

0,833

0,833

0,833

0,834

0,834

90

0,834

Таблица 4 - Значения функции Р(m, Wпр) для разных значений произведения Wnp

Целые Wпр, см

Десятые доли Wпр, см*

0,0

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

0

1,00

0,91

0,88

0,92

0,85

0,84

0,82

0,82

0,82

0,81

1

0,80

0,80

0,80

0,80

0,79

0,79

0,78

0,78

0,78

0,77

2

0,77

0,77

0,77

0,76

0,76

0,75

0,75

0,75

0,75

0,75

3

0,74

0,74

0,74

0,74

0,74

0,74

0,73

0,73

0,73

0,73

4

0,72

0,72

0,72

0,72

0,72

0,72

0,71

0,71

0,71

0,71

5

0,71

0,70

0,70

0,70

0,70

0,69

0,69

0,69

0,69

0,69

6

0,68

0,68

0,68

0,68

0,68

0,68

0,68

0,68

0,68

0,68

7

0,67

0,67

0,67

0,67

0,67

0,67

0,67

0,67

0,67

0,67

8

0,66

0,66

0,66

0,66

0,66

0,66

0,66

0,66

0,66

0,66

9

0,66

0,66

0,66

0,66

0,65

0,65

0,65

0,65

0,65

0,65

* Значения W и Wпр рассчитывают в г/см3, при этом следует помнить, что слой осажденной воды высотой 1 см эквивалентен 1 г/см2, т.е. численно W (см) @ W (г/см2).

14.4 Метод расчета общего и общего приведенного влагосодержания вертикального столба атмосферы по данным аэрологического зондирования

14.4.1 Общее влагосодержание W (г/см2) вертикального столба атмосферы высотой Н = 10 км рассчитывают с округлением до 0,1 г/см2 по формуле

                                              (40)

а общее приведенное влагосодержание Wnp (г/см2) - по формуле

                                     (41)

где а(hk) - профиль абсолютной влажности атмосферы, г/м3. Метод расчета профиля представлен в 14.4.2;

P(hk) - профиль атмосферного давления (данные аэрологического зондирования атмосферы), гПа;

Р(0) - атмосферное давление вблизи водной поверхности (данные аэрологического зондирования атмосферы), гПа;

DhJ - мощность (толщина) J-го слоя атмосферы между двумя уровнями измерений атмосферного давления, температуры и влажности воздуха, м;

k - порядковый номер уровня измерений аэрологическим зондом;

0,0001 - коэффициент перевода значений W и Wпр из граммов на метр квадратный (г/м2) в граммы на сантиметр квадратный (г/см2).

14.4.2 В результате аэрологического зондирования атмосферы получают распределение с высотой (профиль) температуры воздуха t(hk) (°C), относительной влажности воздуха f(hk) (%) и атмосферного давления Р(hk) (гПа). Профиль абсолютной влажности воздуха в атмосфере а(hk) рассчитывают по формуле из методических указаний [7]

а(hk) = 2,167f(hk)E(hk)/[273,2 + t(hk)],                                       (42)

где E(hk) - давление насыщенного водяного пара на уровне измерения k, гПа. Его рассчитывают по выражению из ГОСТ 8.524 (см. также формулу (14) в разделе 9 настоящих методических указаний)

Е(hk) = Е0ехр{at(hk)/[b + t(hk)]},                                          (43)

где Е0 = 6,1121 гПа;

a и b - постоянные, которые принимают значения соответственно aw = 17,5043 и bw = 241,2 °С для воды во всем диапазоне значений температуры t(hk), в том числе и отрицательных*.

______________

* Это обусловлено тем, что все ПИП влажности (в том числе и радиозондовые) калибруются по отношению к жидкой фазе воды, включая область отрицательных значений температуры.

14.4.3 Метод расчета W и Wпр состоит в следующем:

а) рассчитывают общее влагосодержание атмосферы DWJ (г/м2) с округлением до 0,01 г/м2 в слое мощностью DhJ по формуле

DWJ = 0,5[a(hk) + a(hk+1)](hk+1 - hk)

при k, J = 1, 2, ..., N;                                                     (44)

б) рассчитывают общее приведенное влагосодержание атмосферы (г/м2), с округлением до 0,01 г/м2 в слое мощностью DhJ (м) по формуле

                                               (45)

в) рассчитывают общее W и общее приведенное Wпр влагосодержание с округлением до 0,01 г/м2 по формулам

                                                            (46)

                                                         (47)

Как следует из [12], расчет DWJ от уровня моря до высоты 1 км можно осуществлять с шагом DhJ = 100 м, а выше 1 км - с шагом DhJ = 1 км, что практически не отражается на точности расчета W при меньших значениях DhJ.

14.4.4 В приложении Д приведен тестовый пример расчета W и Wnp в соответствии с 14.4.3.

Приложение А
(обязательное)

Общая формула перевода данных измерений из единиц промежуточной величины в единицы измеряемой величины

Общая формула перевода результатов измерений из единиц промежуточной величины в единицы измеряемой величины имеет следующий вид:

Xi = bz[(xt - xJ)f(m)]F(n),                                                   (A.1)

где Xi - результат измерения любой гидрометеорологической величины с помощью ПИП или ПАП (в i-й момент времени, i = i, ..., I) в принятых единицах этой величины;

bz - знак результата измерения, сохраняющий постоянство с момента времени z;

xi - результат измерения в единицах промежуточной величины по ПИП или ПАП в i-й момент времени;

xJ - значение места нуля ПИП или ПАП, сохраняющее постоянство в течение J-го интервала времени;

f(т) - функция перевода из промежуточных величин в значения измеряемой величины для конкретного ПИП или ПАП, которая зависит от аргумента т;

F(п) - функция поправок к результатам измерений величины для конкретного ПИП или ПАП, зависящая от значения измеряемой величины.

Как правило, значения функций f(m) и F(п) задаются в виде таблиц или в аналитическом виде и приводятся в свидетельствах о поверке к ПИП, ПАП и измерительным приборам либо к измерительным системам, состоящим из ПИП (ПАП) и измерительного прибора, либо к измерительным каналам автоматических гидрометеорологических станций, состоящим из электрической цепи, на входе которой используется ПИП (ПАП), а на выходе - измерительное устройство.

В случаях, когда измерения производят по приборам непосредственно в единицах измеряемой величины (например, атмосферное давление измеряют в гектопаскалях), формула (А.1) принимает вид

Xi = bz(xi - xJ)F(n),                                                       (A.2)

где значения F(п), хi и xJ берут в единицах измеряемой величины.

Приложение Б
(справочное)

Дополнительные рекомендации к методу расчета направления истинного ветра*

_____________

* Метод расчета представлен в наставлении [3].

Б.1 В наставлении [3] и в 7.1.2 направление истинного ветра d предлагается рассчитывать по одной и той же формуле (6)

d = dc + dк ± dки,                                                        (Б.1)

где dc - компасный курс судна;

dк - направление кажущегося ветра (в соответствии с 7.1.2);

dки - угол между кажущимся (наблюденным на судне) и истинным ветром (в соответствии с 7.1.2).

Однако в 7.1.2 dки предлагается рассчитывать по формуле (7), которая имеет однозначное решение, а в [3] - через arcsin[(0,5144Vcsindк)/V], который при определенных значениях Vс и V (скорость перемещения судна и скорость истинного ветра соответственно) не всегда имеет однозначное решение.

Так, если Vcsindк < V < Vc, то существуют два решения:

а)  < 90°,

б)  = 180° - .                                                                                                   (Б.2)

С учетом сказанного при условии, когда Vcsin dк < V < Vс, направление истинного ветра рассчитывают по формуле

d = dc + dк + 180° - .                                              (Б.3)

Без учета этого условия расчеты направления истинного ветра d не всегда будут правильными.

Б.2 Иллюстрацией к сказанному является следующий пример.

Пример - При dс = 260°, Vс = 12,5 уз, dк = 40°, Vк = 2,5 м/с, скорость истинного ветра V, рассчитанная по формуле (5), оказалась равной 4,8 м/с, а направление d, рассчитанное по формуле (Б.1), оказалось равным 0°. Если dки рассчитывать через arcsin[(0,5144Vсsindк)/V], то

dки = arc sin [(0,5144Vcsin dк)/V] = 60°.                             (Б.4)

Расчеты dки, выполненные по формуле (7) и графическим способом с помощью ветрочета КСМО-1 [3], показали, что dки = 120°.

Расчеты dки, выполненные по формуле (Б.4) с учетом проверки выполнения соотношения Vcsindк < V < Vc показали, что dки также равно 120°.

Приложение В
(обязательное)

Методы расчета характеристик влажности воздуха при отрицательных значениях температуры воздуха и отсутствии информации об агрегатном состоянии воды на резервуаре смоченного термометра

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

Погрешность психрометрического метода увеличивается с уменьшением температуры tа. Так, в области повышенных значений относительной влажности воздуха f при ta = 10 °С погрешность измерения температуры t¢ смоченного термометра Dt¢ = 0,1 °С приводит к погрешности измерения f, равной 1 % (Df = 1 %), тогда как при ta = -10 °С та же погрешность Dt¢ = 0,1 °С эквивалентна Df = 3 %. По этой причине указанное значение tа = -10 °С принято за рекомендуемый нижний предел температурного диапазона применения психрометрического метода измерения характеристик влажности воздуха, хотя никаких ограничений физического характера для его использования при более низкой температуре не существует.

К еще большей погрешности приводит отсутствие информации об агрегатном состоянии воды на резервуаре смоченного термометра при отрицательных значениях температуры. Продолжая рассмотрение выше указанного примера, допустим, что tа = t¢ = -10 °С. Если известно, что вода находится в жидком (переохлажденном) состоянии, то f = 100 %; если она в твердом состоянии (лед), то

f = e/Ew(ta) = Et(-10)/Ew(-10) = 2,60/2,86 @ 0,91,                            (B.1)

т.е. 91 %. В случае отсутствия информации об агрегатном состоянии воды на резервуаре смоченного термометра приходится принять гипотезу об одном из альтернативных состояний, и если гипотеза не оправдывается, то возникает погрешность, в приводимом примере составляющая 9 %: Df = 100 - 91. Она, как видно, превышает погрешность, обусловленную погрешностью измерения t¢, в 3 раза.

Способ снижения этой погрешности при отсутствии информации об агрегатном состоянии воды на резервуаре смоченного термометра основан на использовании среднего арифметического из двух возможных крайних значений величины. В рассматриваемом примере парциальное давление водяного пара е, рассчитанное для воды по формуле (11), составляет 2,86 гПа [е = Ew(t¢) = Ew(-10) = 2,86 гПа], а для льда, рассчитанное по формуле (12), составляет 2,60 гПа [е = Ei(t¢) = Ei(-10) = 2,60 гПа]. При отсутствии информации об агрегатном состоянии воды наилучшей оценкой е будет среднее (виртуальное) из полученных значений: е = 0,5(2,86 + 2,60) = 2,73 гПа. Обобщая изложенное, получим формулу (13) раздела 9 настоящих методических указаний. Далее с полученным «виртуальным» значением е следует обращаться, как с реальным. При этом необходимо помнить, что при вычислении относительной влажности f, дефицита насыщения d при отрицательном значении температуры воздуха ta давление насыщенного водяного пара Ew(tа) берут по отношению к воде в жидком состоянии. С учетом этого «виртуальное» значение f в рассматриваемом примере будет равно 95,5 % (f = 2,73/2,86), а дефицит насыщения d = Ew(ta) - е = 2,86 - 2,73 = 0,13 гПа.

Очевидно, что полученная таким образом оценка е никогда не совпадет с фактическим значением, несмотря на то, что она наилучшая в рассматриваемых условиях. Но при этом и погрешность никогда не превзойдет половины той погрешности, которая появляется вследствие произвольного предположения об агрегатном состоянии воды на резервуаре смоченного термометра при отсутствии такой информации. Оцениваемая «половинная» погрешность может иметь как положительный, так и отрицательный знак, и носит случайный характер. В рассматриваемом примере реализуется самая большая по абсолютному значению погрешность: Df = ±4,5 % (если фактически на резервуаре жидкая вода, то Df = 95,5 - 100 = -4,5 %, а если лед, то Df = 95,5 - 91,0 = +4,5 %), В середине рабочего поддиапазона отрицательных значений температуры воздуха (от 0 °С до минус 10 °С), т.е. при ta = -5 °С, наибольшее значение Df = ±2,5 %.

Возвращаясь к началу рассмотрения примера, можно видеть, что при температуре воздуха ta = -10 °С и его относительной влажности f = 91 % (при этом на резервуаре смоченного термометра лед) психрометрическая разность tа - t¢ = 0, т.е. ta = t¢ = -10 °С. Физически это означает, что в данных условиях имеет место термодинамическое равновесие между льдом на резервуаре термометра и водяным паром в воздухе, т.е. не происходит ни испарения, ни конденсации.

Очевидно, что если относительная влажность окажется выше указанного равновесного значения, т.е. f > 91 %, то будет происходить конденсация водяного пара на ледяной поверхности, что вызовет повышение температуры смоченного термометра вследствие выделения теплоты фазового перехода (пара в лед). Тогда окажется, что t¢ > tа, т.е. психрометрическая разность ta - t¢ приобретет отрицательное значение. Такое превышение t¢ над ta достигает наибольшего значения при f = 100 % и в данном примере составляет 0,3 °С, т.е. при f > 91 % увеличению f на каждые 3 % будет соответствовать в среднем приращение t¢ на 0,1 °С.

Таким образом, сам факт превышения t¢ над tа (или, иначе говоря, отрицательный знак психрометрической разности ta - t¢) свидетельствует о твердом агрегатном состоянии воды на резервуаре смоченного термометра. И поскольку конденсация на лед является процессом, обратным испарению, то, естественно, расчет значений е следует выполнять в этом случае по формуле (12). Изложенное выше приводит к выводу о том, что в подобных ситуациях, встречающихся в практике измерений, превышение t¢ над ta не следует воспринимать как ошибку, особенно в тех случаях, когда оно не достигает больших значений.

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

Формат (устройство) таблиц такой же, как и входящих в состав [10] таблиц, рассчитанных отдельно для каждого из двух возможных агрегатных состояний воды на резервуаре смоченного термометра (см. с. 32 - 36 в [10]), и отличается от них только сигнальной надписью в верхних углах страниц (вместо «вода!» или «лед!» надписано: «агрегатное состояние неизвестно!»). Отсюда видно, что и правила пользования рассматриваемыми таблицами остаются неизменными.

Такие таблицы являются, следовательно, расширением существующих психрометрических таблиц [10], при очередном издании которых это может быть учтено.

Приложение Г
(справочное)

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

Как следует из раздела 14, для оценки натурального коэффициента ослабления прямой солнечной радиации аэрозолем а необходима, наряду с другими данными, информация об общем влагосодержании W вертикального столба атмосферы высотой 10 км. Наилучший способ получения такой информации - это расчет W по данным аэрологического зондирования атмосферы. Однако его выполняют с борта единичных судов, в то время как характеристики влажности воздуха в приводном слое могут определяться с борта любого судна, на котором организована и функционирует СГМС.

Учитывая повышенное по сравнению с сушей влагосодержание атмосферы над поверхностью моря и тот факт, что наибольший вклад в общее влагосодержание вертикального столба атмосферы вносит влагосодержание ее приводного слоя, были найдены эмпирические соотношения между W и отдельными характеристиками влажности воздуха приводного слоя атмосферы. С целью выявления наиболее корректного эмпирического соотношения для определения W в трудах ГГО [12] проведено сопоставление значений W, рассчитанных по данным аэрологического зондирования атмосферы, и значений Wi, рассчитанных по эмпирическим соотношениям, представленным в таблице Г.1.

Таблица Г.1 - Сведения об эмпирических соотношениях для расчета значений Wi

Вид эмпирического соотношения

Исследуемая акватория Мирового океана

Единица величины, используемая в приводимых соотношениях

W

е

Количество облаков п

Коэффициент g0

1 W1 = 0,25е - 2,16

Тропическая зона Атлантики

г/см2

гПа

-

-

2 W2 = 0,15e

Тропическая зона Атлантики

г/см2

гПа

-

-

3 W3 = (1,55 + 0,046п)е1,075

Разные акватории океанов

мм

мм рт.ст.

балл

-

4 W4 = 1,79пе1,08 (при п ¹ 0)

Разные акватории океанов

мм

мм рт.ст.

балл

-

5 W5 = 0,7159е + 0,3144е2 + 0,0002426e3

Разные акватории океанов

мм

гПа

-

-

Результаты сопоставления показали, что по эмпирическим соотношениям 1 - 5 можно получать надежные значения W: средние расхождения между общим влагосодержанием вертикального столба атмосферы W, определенным по данным аэрологического зондирования, и рассчитанными по указанным соотношениям значениями Wi составили от 4 % до 12 % от значения W для конкретных условий.

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

Wпр = 0,69W + 0,27,

в котором значение W (г/см2) рассчитывают по любому соотношению из таблицы ГЛ.

Следует помнить, что 1 см осажденной воды соответствует 1 г/см2, поэтому значения W, рассчитанные в миллиметрах, необходимо перевести в граммы на сантиметр квадратный по соотношению: W [г/см2] = 0,1W [мм].

Приложение Д
(обязательное)

Тестовый пример расчета общего и общего приведенного влагосодержания вертикального столба атмосферы

Д.1 Исходные данные для расчета W и Wпр, являющиеся данными аэрологического зондирования атмосферы с борта НИС «Муссон» 2 августа 1985 г. в одной из точек Северной Атлантики, с результатами расчета представлены в таблице Д.1.

Д.2 Переводим åDW º W и åDWпр º Wпр из граммов на метр квадратный в граммы на сантиметр квадратный по формулам

                           (Д.1)

                         (Д.2)

Таблица Д.1

i

J

Исходные данные

Результаты расчета

h км

t(hi)C

f(hi) %

Р(hi) гПа

Е(hi) гПа

е(hi) гПа

а(hi) г/м3

DW г/м2

 г/м2

1

-

0

20,2

85

1000,0

23,25

19,76

14,59

-

-

2

1

0,1

20,1

85

999,4

23,11

19,64

14,51

1455

1454

3

2

0,2

20,0

84

999,0

22,96

19,29

14,27

1439

1438

4

3

0,3

19,8

82

998,0

22,68

18,60

13,75

1401

1398

5

4

0,4

19,0

80

990,0

21,60

17,28

12,82

1328

1318

6

5

0,5

18,5

75

984,0

20,94

15,71

11,67

1224

1204

7

6

0,7

17,9

75

978,0

20,18

15,13

11,26

2293

2256

8

7

0,9

15,0

78

965,0

16,81

13,11

9,86

2112

2066

9

8

1,0

14,8

73

900,0

16,60

12,12

9,12

949

916

10

9

2,0

10,0

65

810,0

12,15

7,90

6,04

7585

6827

11

10

3,0

2,0

50

750,0

7,04

3,52

2,77

4405

3668

12

11

4,0

-4,0

42

680,0

4,56

1,92

1,54

2160

1620

13

12

6,0

-16,0

42

430,0

1,79

0,75

0,63

2170

1476

14

13

9,0

-30,0

25

50,0

0,52

0,13

0,12

1125

484

15

14

10,0

-35,0

15

20,0

0,32

0,05

0,04

80

4

Сумма å

29726

26129

Библиография

[1]

Методические указания по производству метеорологических и актинометрических наблюдений на научно-исследовательских судах. Ч. III. Требования для подготовки данных для ввода в ЭВМ. - Л.: Ртп. ГГО, 1986

[2]

РД 52.04.316-92

Наставление гидрометеорологическим станциям и постам. Гидрометеорологические наблюдения на морских станциях. Вып. 9, ч. II. Гидрометеорологические наблюдения на судовых станциях, производимые штатными наблюдателями. Кн. 1. Общие методические требования к организации и обеспечению гидрометеорологических и актинометрических наблюдений на судах

[3]

РД 52.04.585-97

Наставление гидрометеорологическим станциям и постам. Гидрометеорологические наблюдения на морских станциях. Вып. 9, ч. III. Гидрометеорологические наблюдения, производимые штурманским составом на морских судах

[4]

РД 52.27.133-89

Международный код FM 12-1X SYNOP и 13-1Х SHIP для передачи данных приземных гидрометеорологических наблюдений с наземных и морских станций

[5]

Руководство по метеорологическим приборам и методам наблюдений. - Женева: ВМО, 1983. - № 8

[6]

РД 52.04,316-92

Наставление гидрометеорологическим станциям и постам. Гидрометеорологические наблюдения на морских станциях. Вып. 9, ч. II. Гидрометеорологические наблюдения на судовых станциях, производимые штатными наблюдателями. Кн. 2. Методики выполнения гидрометеорологических наблюдений на судах

[7]

РД 52.04.129-87

Методические указания. Измерение аэрозольной мутности атмосферы по данным интегрального коротковолнового излучения в условиях НИС

[8]

Р 52.04.605-99

Рекомендации. Оценка работ судовых гидрометеорологических станций

[9]

РД 52.04.614-2000

Наставление гидрометеорологическим станциям и постам. Вып. 3, ч. II. Обработка материалов метеорологических наблюдений

[10]

Психрометрические таблицы. - 2-е изд. - Л.: Гидрометеоиздат, 1981

[11]

РД 52.04.562-98

Наставление гидрометеорологическим станциям и постам. Вып. 5. Актинометрические наблюдения. Ч. 1. Актинометрические наблюдения на станциях

[12]

Кессель В.Н.

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

 

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

Лист регистрации изменений

Номер изменения

Номер страницы

Номер документа

Подпись

Дата

измененной

замененной

новой

аннулированной

внесения изменения

введения изменения