| Автор |
Сообщение |
Spasatel
Участник
Талдом
|
# Дата: 12 Дек 2012 16:25
Alter
Спасибо
|
met2
Участник
Минск
|
# Дата: 12 Дек 2012 19:52
Spasatel
А вот такой вопрос. Как часто обновляются и когда прогностические карты различных моделей? Те к примеру GFS и UKMET один раз в сутки (те только срок 0) или чаще?
2 раза. Прогоны 00 и 12 ч. Примерно в 10.30 и 22.30 местного карты на сервере.
TornadoF5
Черкассы добавлены. Вечером посчитаются.
|
TornadoF5
Участник
Харьков, Украина. (Игорь)
|
# Дата: 12 Дек 2012 19:54
met2
Спасибо, а то один обл. центр был обижен :)
|
Alter
Участник
Москва, ЮЗАО
|
# Дата: 27 Дек 2012 01:54
met2
3) Учли пожелания ряда товарищей по фазе осадков. Для уточнения фазы подключили T850 для трех моделей UKMET, GFS, CMC.
Спасибо Вам большое.
Сегодня пожинали первые плоды обновленного таким образом прогноза. Прогноз GFS показывал ледяной дождь при отрицательной температуре - а она в последнее время явно точнее при определении высотных температур -, и действительно, гололедные осадки в количестве около 2 мм. выпали сегодня утром в столичном регионе.
|
met2
Участник
Минск
|
# Дата: 18 Янв 2013 15:56
Внимательнее! Теперь в мультимодельных метеограммах зелёненьким - GEM-модель. Канада, CMC.
T и Td
Давление и средний ветер.
|
met2
Участник
Минск
|
# Дата: 25 Янв 2013 13:19 - Поправил: met2
Взялись за WRF. Сразу напишу: не знаю, осилим ли...
Думаем считать пока раз в три часа с заблаговременностью 48 ч на сетку 15-20 км.
Домен, примерно: западнее Варшавы - Ярославль / севернее Питера - Одесса.
На входе - нулевой реанализ любимого GFS.
Хотим посчитать вот такие поля.
Если есть замечания, предложения. Пишите.
Конвекция и турбулентность
CIN – энергия задерживающего слоя [Дж/кг].
4LFTX – Best LI индекс плавучести (разность поднимающейся воздушной частицы и окр. воздуха от земли до 500 гПа) [K].
SI – Showalter Index – индекс неустойчивости.
SB CAPE – энергия неустойчивости (0-500 гПа) [Дж/кг].
ML CAPE – энергия неустойчивости (в слое 0-180 гПа) [Дж/кг].
KX – K-index – индекс неустойчивости.
SX – Sweat-индекс – индекс опасных, связанных с конвекцией явлений.
TKE – кинетическая энергия турбулентности [Дж/кг].
RI – число Ридчарсона.
Облачность, %
TCDC – общее кол-во облаков.
LCDC – облака нижнего яруса.
MCDC – облака среднего яруса.
HCDC – облака верхнего яруса.
CDCON – облака конвективные.
FG – приземная обл./туман.
Температура, K
TMP2m – температура воздуха на 2 м.
TMAX, TMIN – температура максим. и мин. воздуха на 2 м.
TMPSFC – температура поверхности почвы.
TMP10sm – температура почвы на глубине 10 см.
DPT – точка росы.
DEPR – дефицит точка росы.
VTMP – виртуальная температура
Влажность
SPFH удельная влажность [кг/кг].
RH относительная влажность [%].
POT – потенциальная температура [K].
EPOT – потенциальная псевдоэквивалентная температура [K].
Осадки
CSNOW – фаза осадков, снег [0-нет/1- да].
CRAIN – фаза осадков, дождь [0-нет/1- да].
CFRZR – фаза осадков, замерзающий дождь [(0-нет/1- да].
CPRAT – мгновенная интенсивность конвективных осадков [кг/м2/с].
PRATE – мгновенная интенсивность осадков [кг/м2/с].
APCP – общее количество осадков [кг/м2].
AСPCP – количество конвективных осадков осадков [кг/м2].
Вегетация
VEG – индекс вегетации, %.
Влагозапас
SRWEQ – водность в снегопаде [кг/м2/с].
WEASD – водозапас в снежном покрове [кг/м2].
PWAT – влагозапас в атмосфере [кг/м2]
CWAT – водность облаков [кг/кг].
Давление
PRES – давление на поверхности [Па]
PRMSL – давление, приведенное к уровню моря [Па]
HGT – высота геопотенциала [гп. м].
Ветер
UGRD – u-компонента ветра.
VGRD – v-компонента ветра.
GUST – скорость порывов ветра на высоте 10 м [м/с].
ABSV – вихрь абсолютный [1/c].
ABSD – дивергенция абсолютная [1/c].
DZDT – вертикальные токи [м/с].
Поверхности
-10 см
0 м
2 м
10 м
1000 гПа
975 гПа
950 гПа
925 гПа
900 гПа
875 гПа
850 гПа
800 гПа
750 гПа
700 гПа
600 гПа
500 гПа
400 гПа
300 гПа
200 гПа
100 гПа
Уровень тропопаузы ?
Уровень конвекции ?
Уровень конденсации ?
Уровени изотермы 0 и -22 ?
|
TornadoF5
Участник
Харьков, Украина. (Игорь)
|
# Дата: 25 Янв 2013 13:35
met2
Скажите, планируется проноз WRF в табличном виде, как по остальным моделям?
|
Cumulonimbus incus
Участник
Кишинёв, Молдова
|
# Дата: 25 Янв 2013 13:36
met2
Если Вы сделаете такие карты и прогнозы, то все мы будем у Вас в долгу :-) Успехов в начинаниях, идея отличная!
|
met2
Участник
Минск
|
# Дата: 25 Янв 2013 13:54
Вооще-то, идея выложить в раздачу в grib2 ;) по типу, как это делают GFS/CMC.
|
LESS
Участник
Приокский лесс на юге МО
|
# Дата: 25 Янв 2013 13:58
met2
Уровень тропопаузы ?
Уровень конвекции ?
Уровень конденсации ?
Уровени изотермы 0 и -22 ?
А что знаки вопроса означают?~
Пишите.
Если не секрет, что планируете "подкладывать" в качестве рельефа поверхности?
|
spralex
Участник
г. Конотоп, Сумская обл., Украина
|
# Дата: 25 Янв 2013 14:54 - Поправил: spralex
met2
Отличная идея! Удачи в реализации!
Хотелось бы еще, по-возможности, добавить индекс МКС (MCS):
Формула
где каждый член в уравнении (индекс Li, сдвиг в слое 0-3 км и адвекция температуры на уровне 700 гПа) нормирован. Следует обратить внимание, что данный параметр имеет смысл, если имеются условия для развития конвекции (например, при Li < 0)
и еще, если возможно, псевдо-отражаемость (sim. Reflectivity). Пример: http://www.modellzentrale.de/WRF/00Z/12h/maxDBZ_eu .png
|
TornadoF5
Участник
Харьков, Украина. (Игорь)
|
# Дата: 27 Янв 2013 10:30
met2
Скажите пожалуйста, можно ли как-то избавиться от технической ошибки в высоте СП? Просто по факту в Харькове сейчас 5 см снега, а в таблице целых 32 см!
|
VoDya
Участник
Липецк
|
# Дата: 27 Янв 2013 12:59
Уровени изотермы 0 и -22 ?
met2
а почему -22, а не -10, по теории критичная температура при которой начинается существенное образовании осадков
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 27 Янв 2013 19:29 - Поправил: Corvus
а почему -22
VoDya
Изотермы 0 и -22 используются на метеорадарах для распознавания осадков и гроз. Но неплохо было бы и -10 добавить.
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 27 Янв 2013 19:49
met2
Поверхности
Я бы добавил, как в ГФС, слои (по Т, влажности, ветру)
30-0 mb above ground
60-30 mb above ground
90-60 mb above ground
120-90 mb above ground
150-120 mb above ground
180-150 mb above ground
Они нужны для детального прогноза низкой облачности, тумана и т.п. целей. Ценность этих уровней в том, что они именно относительно рельефа, а не уровня моря.
|
met2
Участник
Минск
|
# Дата: 28 Янв 2013 13:47
TornadoF5
Скажите пожалуйста, можно ли как-то избавиться от технической ошибки в высоте СП? Просто по факту в Харькове сейчас 5 см снега, а в таблице целых 32 см!
К сожалению, никак.
Какой у GFS на входе реанализ, то такой и расчет.
Реанализ - фактические данные, "натянутые" (синтерполированные) на регулярную сетку. Остаётся только догадываться данные с каких метеостанций будут включены в реанализ, какие параметры и виды интерполяции используются для построения регулярной сетки.
Счет модели осуществляется по узлам регулярной сетки.
При формировании выходных grib-файлов результаты счета нужно опять синтерполировать на 0.5-градусную сетку.
Теперь, мы берем гриб-файл со снегом и при записи в базу еще раз интерполируем на Харьков.
Смотрите, по карте рядом с Харьковом, в Белгороде 49 см, и рядом 0.5 см.
В результате - вот они, 30 см. Ровно и красиво размазаны. :)
http://meteoinfo.by/maps/?type=fact&map=SNOW_UKR&d ate=2013012800&time=06
|
TornadoF5
Участник
Харьков, Украина. (Игорь)
|
# Дата: 28 Янв 2013 15:14
met2
Благодарю за подробное объяснение! Теперь всё понятно :)
|
met2
Участник
Минск
|
# Дата: 2 Фев 2013 14:54
Докладываю.
Уважаемый kato http://rain.hyarc.nagoya-u.ac.jp/~kato/work/GMT/in dex-e.html пошел нам на встречу и любезно сделал патч grdbarb для GMT 4.5.9 под linux.
Grdbarb - инструмент, позволяющий отображать ветер так, как принято в метеорологии.
Буквально вчера удалось собрать GMT 4.5.9 под linux CentOS 5 (тот еще гемор). Но в результате - та-дам...
http://meteoinfo.by/maps/?type=cmc&map=TMP850&date =2013020200&time=12
Конечно, мы и раньше наносили ветер на карты с помощью GMT, но то была версия под Windows, что не очень серьезно.
Пока прикрутили ветер к основным поверхностям 925, 850, 700, 500, 300.
Чтобы карта была читабельной, наносим ветер через 1.2°.
|
TornadoF5
Участник
Харьков, Украина. (Игорь)
|
# Дата: 8 Фев 2013 11:03
met2
Уже второй день не поступают снимки на ваш сайт:
http://www.pogoda.by/eumetsat/
|
met2
Участник
Минск
|
# Дата: 8 Фев 2013 14:27
TornadoF5
Поправили.
|
met2
Участник
Минск
|
# Дата: 19 Фев 2013 17:54
Кое чё поменяли в разделе "Карты". Надеюсь, улучшили :)
|
TornadoF5
Участник
Харьков, Украина. (Игорь)
|
# Дата: 19 Фев 2013 18:20
met2
Кое чё поменяли
А что именно. Можно конкретнее пожалуйста?
|
met2
Участник
Минск
|
# Дата: 19 Фев 2013 21:04 - Поправил: met2
TornadoF5
Да по интерфейсу (дурацко звучит). :)
Добавили GEM в быстрое переключениями между картами.
|
met2
Участник
Минск
|
# Дата: 20 Фев 2013 18:42
Часто на почту приходит "Как вы делаете прогнозы на POGODA.BY?" Овечал товарищу. Продублирую на форуме. М.б. кому-то станет интересно.
1. =====================================
За 2-3 месяца сделать сайт и запрограммировать невозможно.
Опишу основные этапы. Препроцессинг.
1. Загрузка данных в формате grib2 с сайта NCEP GFS High Resolution (0.5
degree, 1 week) http://nomad3.ncep.noaa.gov/ncep_data/index.html
Загрузка файлов осуществляется утилитой WGET, которая входит в дистрибутив Linux.
2.Декодировать данные с помощью утилиты DEGRIB
http://www.nws.noaa.gov/mdl/degrib/ . Эта же утилита осуществляет
интерполяцию данных на точку.
3. Записать проитерполированные данные в базу MySQL (входит в комплект Linux). Здесь же можно выполнить некоторые предварительные вычисления. Например, мы вычисляем точку росы из удельной влажности - SPFH.
Постпроцессинг.
При обращении пользователя к странице сайта осуществляется выборка данных из базы MySQL.
И прежде, чем вывести результаты на экран запускаются процедуры,
преобразующие сырые данные в понятный пользователю прогноз погоды.
Процедуры запрограммированы на PHP (входит в комплект Linux). В процедурах реализованы методики и критерии по определению явлений погоды.
Все методики можно найти в небезысвестных учебниках: Синоптической
Метеорологии, Физике Атмосферы.
Алгоритмы, естественно, раскрыть не могу.
2. ===================
Практическая реализация "на коленке". Основной инструмент wgrib2. Работает в Windows. Архив rar.
В архиве находится утилита wgrib2. Она умеет раскодировать файлы в формате grib2 NOAA и CMC с подстановкой географических координат.
Утилита состоит из исполняемого .exe-файла и трёх dll. Тут же находится grb2.bat - это командный файл, который запускается в командной строке Windows (Меню Пуск > Выполнить> cmd.exe).
В папке data лежат исходные grib-файлы.
Запускайте grb2.bat и на выходе получатся текстовые .csv файлы со значениями в конкретной географической точке.
Открыв в текстовом редакторе grb2.bat можно координаты изменить, также можно параметрами поменять формат вывода данных.
Параметры программы и помощь по работе wgrib2:
http://www.cpc.ncep.noaa.gov/products/wesley/wgrib 2/
Если не ошибаюсь, там есть примеры, в которых реализуется подстановка множества координат. Утилита работает в цикле, где ей "скармливаеся файлик" с географическими координатами.
Удачи! Эксперементируйте.
PS
Здесь можно взять грибы модели CMC.
Они удобно разрезаны: 1 параметр, 1 срок, 1 поверхность
http://dd.weatheroffice.gc.ca/model_gem_global/hig h_resolution/grib2/lat_lon/
|
misha
Участник
Украина, восточный Донбасс
|
# Дата: 20 Фев 2013 23:50
met2
Спасибо, почитал.
Ваша команда насчитывает 1-2 или 10-15 замечательных специалистов?
|
met2
Участник
Минск
|
# Дата: 21 Фев 2013 12:00
Ваша команда насчитывает 1-2 или 10-15 замечательных специалистов?
Не припомню, чтоб одновременно работало более 4 человек :)
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 22 Фев 2013 10:31 - Поправил: Corvus
1. Загрузка данных в формате grib2 с сайта NCEP GFS High Resolution (0.5
degree, 1 week) http://nomad3.ncep.noaa.gov/ncep_data/index.html
Загрузка файлов осуществляется утилитой WGET, которая входит в дистрибутив Linux.
2.Декодировать данные с помощью утилиты DEGRIB
http://www.nws.noaa.gov/mdl/degrib/ . Эта же утилита осуществляет
интерполяцию данных на точку.
3. Записать проитерполированные данные в базу MySQL (входит в комплект Linux). Здесь же можно выполнить некоторые предварительные вычисления. Например, мы вычисляем точку росы из удельной влажности - SPFH.
met2
Я с базой данных связываться не стал, всё храню в файлах.
1. Скачиваю ГРИБЫ2 wget'ом.
2. С помощью wgrib2 -d NUM FILENAME -no_header -order we:sn -bin NEWFILENAME
конвертирую ГРИБ2 в простой бинарный сеточный формат (одна переменная за один момент времени в одном файле), откуда данные легко может читать РНР-скрипт.
NUM - это номер переменной в большом ГРИБЕ2, содержащем сразу много переменных за один момент времени.
При этом данные по одной точке считываются из нескольких сотен файликов, но происходит это очень быстро, на всё про всё менее 0.1 секунды.
3. Прогноз создаётся на лету РНР скриптом (по координатам пункта).
Чтение wgrib2'ом из ГРИБА2 точки по координатам пробовал, но работает это медленно - несколько секунд, то есть в десятки и сотни раз медленнее, чем моя технология.
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 22 Фев 2013 10:33
З.ы. Интерполяции по пространству нет. Просто берётся ближайший узел модельной сетки.
|
met2
Участник
Минск
|
# Дата: 23 Фев 2013 01:26
Corvus
Да, у wgrib2 интерполяции нет, а у degrib билинейная.
|
Gennady
Участник
|
# Дата: 23 Фев 2013 17:52
met2
Corvus
Очень полезная для меня информация. Большое спасибо.
|
