Автор |
Сообщение |
lalex
Участник
Москва+Крым
|
# Дата: 2 Jun 2011 21:02 - Поправил: lalex
Да нет, вроде бы оба Grib2. Раскодировал, понял - который меньше это одноградусный ;)
Вы, я так понимаю, Grib-1 и одноградусную модель предпочитаете?
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 2 Jun 2011 21:24
Вы, я так понимаю, Grib-1 и одноградусную модель предпочитаете?
lalex
На Метеоцентре используется 0.5-градусная ГФС в формате ГРИБ1.
|
lalex
Участник
Москва+Крым
|
# Дата: 3 Jun 2011 14:23 - Поправил: lalex
Ясно, и с inventory тоже вроде начало проясняться. А насчет выбора конкретных серверов: каталог rotating обновляется быстрее, но в нем заблаговременность только 48 часов, вроде бы?
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 3 Jun 2011 16:02
каталог rotating обновляется быстрее, но в нем заблаговременность только 48 часов, вроде бы?
lalex
Нет, почему же.
180 часов.
http://nomad3.ncep.noaa.gov/pub/gfs/rotating/
|
lalex
Участник
Москва+Крым
|
# Дата: 3 Jun 2011 19:45
Хм. В этом каталоге вроде бы одноградусные грибы, матрица 360х181 (на пробу распаковал). Размер 32 Мб.
Полуградусные, я так понял, в каталоге http://nomad3.ncep.noaa.gov/pub/gfs/rotating-0.5/ (если на Nomad3). И вот там уже - только по 48 часов (файлы притом 56 Мб и сразу Grib-2).
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 3 Jun 2011 20:01
Полуградусные, я так понял, в каталоге http://nomad3.ncep.noaa.gov/pub/gfs/rotating-0.5/ (если на Nomad3). И вот там уже - только по 48 часов (файлы притом 56 Мб и сразу Grib-2).
lalex
Да, точно, я был неправ.
Можно брать 0.5-градусную на 180 ч вот тут
http://nomad1.ncep.noaa.gov/cgi-bin/ftp2u_gfs0.5.s h
Я там беру для Метеоцентра данные ГФС. Но система не шибко надёжная, от 1 до 5 дней в месяц данные не обновляются.
В общем, лучше берите СМС - там всё просто и стабильно :)
|
lalex
Участник
Москва+Крым
|
# Дата: 3 Jun 2011 21:52 - Поправил: lalex
"от 1 до 5 дней в месяц данные не обновляются"
Я как раз хотел об этом спросить... Наверно, по этой причине на Метеоцентре висит прогноз за 01.06?
С другой стороны, на основном сервере Nomads вроде бы все давно уже обновилось, по третье июня включительно, прогон за 12 засветился недавно. И заблаговременность полноценная - 192 часа... Ну в общем к этому всему надо пристреляться самостоятельно, как я понял.
http://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/ prod/
насчет СМС внимательно читаю все отзывы, слышал, думаю. Как и насчет других моделей. Но это отдельная тема, пока что.
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 3 Jun 2011 21:57
насчет СМС внимательно читаю все отзывы, слышал, думаю
lalex
Основные копья ломаются относительно кол-ва конвективных осадков по СМС (оно часто завышено). Если брать из СМС обычные осадки, то будет примерно то же самое, что и по ГФС.
Правда, при этом 30% гроз и ливней будут просто-напросто неспрогнозированы (это присуще ГФС), но для многих метеоклубовцев это меньшее зло, чем перестраховочные прогнозы СМС вкупе с 90% предупреждённостью опасных конв. явлений (у ГФС эта предупреждённость 70%).
|
gasha
Участник
Россия, Карелия, Петрозаводск
|
# Дата: 14 Jun 2011 21:37 - Поправил: gasha
Отсюда скачиваю 65 grib2 файлов (0-192 часа). Затем в батнике командой copy /B склеиваю их в один большой гриб(gfs.t00z.mastergrb2f). По размеру он как исходные 65 грибов (gfs.t00z.mastergrb2f00 и т.д.).
После идёт следующий батник:
perl g2ctl.pl gfs.t00z.mastergrb2f >gfs.t00z.mastergrb2f.ctl
gribmap -i gfs.t00z.mastergrb2f.ctl
Так вот в получившемся ctl файле 61 момент времени, а не 65.
dtype grib2
ydef 41 linear 50.000000 0.5
xdef 101 linear 0.000000 0.500000
tdef 61 linear 00Z14jun2011 3hr
Версия g2ctl.pl у меня 0.0.8.2
wgrib2 версии 0.1.9.1
gribmap (This is gribmap for GrADS Version 2.0.a9.oga.1)
Куда пропадают 4 момента времени (183, 186, 189, 192 часа)?
Вот пример скрипта (tmp2m.gs):
V=tmp2m
V2=pressfc
hour=0
'set mpdset hires'
'set map auto'
'set mproj scaled'
'set lat 50 70'
'set lon 0 50'
'open gfs.t00z.mastergrb2f.ctl'
timex=1
while( timex <= 65)
'set t ' timex
'value='V'-273.16'
'value2='V2''
'set gxout shaded'
'd value'
'run cbarn.gs'
'set gxout contour'
'set ccolor 0'
'd value2'
'q time'
times = sublin(result,1)
timeword = subwrd(times,3)
'draw title 'V' / 'V2', valid: 'timeword''
'printim 'V'_'hour'.png white x800 y800'
timex = timex + 1
hour=hour+3
'clear'
endwhile
Генерируемые файлы от tmp2m_0.png до tmp2m_180.png содержат информация, с ними всё хорошо. А вот оставшиеся с tmp2m_183.png по tmp2m_192.png пустые.
grads пишет:
Data Request Warning: Request is completely outside file limits
Data Request Warning: Request is completely outside file limits
Cannot contour grid - all undefined values
Cannot plot color bar: No shading information
Cannot contour grid - all undefined values
и таких 4 пары сообщений.
Что доказывает, что g2ctl.pl проигнорировал 4 момента времени (183-192). Где в нём чего надо исправить, чтобы заработало. Сижу в никаю в этот скрипт, пока не нашёл.
P.S. Вот изменённый скрипт meteogram_gfs.gs, данные для построения беруться из большого гриба) И понятное дело, что строиться метеограмма по 61 моменту, а не по 65 как хотелось бы. Так как после g2ctl.pl получается не совсем правильный ctl
P.P.S. ОТБОЙ ТРЕВОГИ!
Всё хорошо)) Ошибка была при скачивании грибов: как раз последние 4 из 65 качались с параметрами 61-ого. Теперь все проблемы исчезли))
|
Mesocyclon
Участник
Троицк, Ключевка, Челябинская область
|
# Дата: 5 Май 2014 13:23
Как декодировать GRIB1 или GRIB2 без использования сторонних программ? Реализовать, например на С++ , в интернете всё облазил, не нашёл схему декодирования...
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 5 Май 2014 17:39 - Поправил: Corvus
Mesocyclon
Зачем реализовывать это всё с нуля, если есть прекрасные консольные программы (под Винду и Юникс) wgrib и wgrib2?
Наверняка и их исходные коды есть в инете (на C), их можно встроить в любые свои приложения (хотя, ИМХО, гораздо проще запускать wgrib и wgrib2 в виде отдельных задач по планировщику, а потом своими программами "кушать" уже расшифрованные ГРИБы).
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 5 Май 2014 17:40
Кроме того, в инете есть подробное описание кода ГРИБ. Я 12 лет назад писал на Дельфи самопальную программку для декодирования ГРИБ.
|
qwerty
Участник
Владимирская область, г.Киржач
|
# Дата: 5 Май 2014 18:36 - Поправил: qwerty
Скачал гриб ukmet
Отсюда
http://nomads.ncep.noaa.gov/pub/data/nccf/com/ukme t/prod/ukmet.201405/
Создал файл ukm.ctl
Пытался открыть в Градсе в итоге он мне выдал ошибку Can't open description file
В чём дело?
|
Mesocyclon
Участник
Троицк, Ключевка, Челябинская область
|
# Дата: 5 Май 2014 19:18
Corvus
Прочитал что в формате GRIB2 использован алгоритм сжатия JPEG2000, спрашивается зачем было использовать видеоформат.
http://www.gdal.org/frmt_grib.html
либо я не так понял?
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 6 Май 2014 11:41
Mesocyclon
Не забивайте себе голову лишней информацией, лучше изучите программы ВГРИБ и ВГРИБ2, там и исходные коды есть
http://www.cpc.ncep.noaa.gov/products/wesley/wgrib .html
http://www.cpc.ncep.noaa.gov/products/wesley/wgrib 2/
|
TornadoF5
Участник
Харьков, Украина. (Игорь)
|
# Дата: 6 Май 2014 13:36
Corvus
Проверьте пожалуйста почту, я Вам письмо отправлял дня 2 назад :)
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 6 Май 2014 16:24
TornadoF5
Я видел, при наличии времени отвечу на него.
|
met2
Участник
Минск
|
# Дата: 8 Май 2014 15:15
Corvus
TornadoF5
ну или degrib. она еще проще
|
qwerty
Участник
Владимирская область, г.Киржач
|
# Дата: 17 Jul 2014 14:57 - Поправил: qwerty
Почему при открытие ctl файла выскакивает ошибка Can't open description file ?
|
qwerty
Участник
Владимирская область, г.Киржач
|
# Дата: 16 Авг 2014 15:13
Как открыть grib файл?
|
Mesocyclon
Участник
Троицк, Ключевка, Челябинская область
|
# Дата: 16 Авг 2014 16:55
qwerty
попробуйте этой программой
www.giss.nasa.gov/tools/panoply/download_win.html
|
Meteoman16
Участник
Белогорск, Крым
|
# Дата: 17 Янв 2015 14:10
чем можно открыть GRIB2? уже пытался разными программами, не получается,может как-то можно конвертировать в GRIB1, его я могу открыть
|
matrix145
Участник
Курган
|
# Дата: 17 Янв 2015 14:30
Meteoman16
Файл grb2 легко открывается программой panoplay, выше уже кидали ссылку.
http://www.giss.nasa.gov/tools/panoply/download_wi n.html
|
Meteoman16
Участник
Белогорск, Крым
|
# Дата: 17 Янв 2015 14:58
matrix145 скачал, однако открыть не удалось, тормозит
|
matrix145
Участник
Курган
|
# Дата: 17 Янв 2015 15:19
matrix145 скачал, однако открыть не удалось, тормозит
Я хорошо открываю через него файлы, даже большого размера.
|
Meteoman16
Участник
Белогорск, Крым
|
# Дата: 17 Янв 2015 19:12
matrix145
я так понимаю, что работает он с JAva. У меня открывается окно после загрузки пустое, так и стоит
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 18 Янв 2015 01:28
Meteoman16
Нужно установить или обновить Java на Вашем компьютере.
В Файрфоксе это делается через закладку "Плагины".
|