Работа с GrADS

Метеоцентр.Азия - наш новый сайт с высокодетализированными прогнозами погоды по пунктам России и мира
Облегчённая версия Метеоклуба (для КПК)
Карта активных участников Метеоклуба (105 кБ)    Таблица дней рождения активных участников Метеоклуба
Клуб любителей метеорологии (группа ВКонтакте)

Работа с GrADS

Сейчас в Метеоклубе:
Участников - 0
Максимальное одновременное количество посетителей: 0 [1 Янв 1970 03:00]
Гостей - / Участников -

 - Начало - Ответить - Статистика - Регистрация - Поиск -
МЕТЕОКЛУБ : независимое сообщество любителей метеорологии (Европа и Азия) : ФОРУМ О ПОГОДЕ И ПРИРОДЕ / Компьютерная техника и интернет в метеорологии / Работа с GrADS
. 1 . 2 . 3 . 4 . 5 . 6 . 7 . 8 . 9 . 10 . >>
Автор Сообщение
Kib
Участник
Письмо
# Дата: 10 Окт 2009 18:11


Скажите пожалуйста, как открыть файл с расширением .nc. Сразу использовать GrADS не получается, поскольку сначала нужно составить описание

zzann
Участник
Письмо
Севастополь
# Дата: 10 Окт 2009 19:39


вместо open ... нужно писать sdfopen ...
если новая версия ГРАДСа, то откроет .nc без ctl.
у меня версия 1.9.0 открывает .nc

Kib
Участник
Письмо
# Дата: 11 Окт 2009 08:17 - Поправил: Kib


У меня при команде sdfopen ... на экран выводятся только аргументы: широта, долгота, уровень, время. Хотелось бы знать имя содержащегося в файле метеопараметра, чтобы построить карту

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 11 Окт 2009 09:06


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

Для отображения списка переменных после команды sdfopen...
нужно дать команду:
q file

dmitpav
Участник
Письмо
Ukraine
# Дата: 11 Янв 2010 22:41


Подскажите, пожалуйта, есть ли какие-то пособия на русском языке по GrADS? Необходимо сформировать файлы для него из текстовых данных... Учебник на официальном сайте (150 ст. на английском языке) пока что не очень поддаётся пониманию

Kib
Участник
Письмо
# Дата: 19 Янв 2010 19:11 - Поправил: Kib


У меня есть. Дал один знакомый. Вы можете указать ваш e-mail на форуме или прислать письмо на imagtriada@mail.ru. Тогдя я Вам пришлю инструкцию по электронной почте

Kib
Участник
Письмо
# Дата: 20 Янв 2010 09:59


Я уже не раз скачивал из Интернета .nc и особых проблем не возникало. Сейчас при попытке открыть файл sdfopen... выводится сообщение:

Scanning self-describing file: ...
gadsdf: SDF file has no discernable X coordinate

Kib
Участник
Письмо
# Дата: 23 Янв 2010 11:48


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

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 23 Янв 2010 14:42


NCEP придумал какую-то засаду для защиты информации
Kib

Маловероятно. Непохоже на них.

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

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 23 Янв 2010 14:52


Похожая ситуация описана тут:

https://lists.ictp.it/pipermail/regcnet/2006/00024 5.html

Mesocyclon
Участник
Письмо
Троицк, Ключевка, Челябинская область
# Дата: 25 Фев 2012 16:26 - Поправил: Mesocyclon


Скрипты для GRADS по моделям GEFS FNMOC ансамбль
ссылка

Образцы
GEFS


FNMOC


Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 17 Июл 2012 15:15


Полезные советы по работе с ГРАДС (в том числе ссылка на перевод мануала по ГРАДС 2008 года)

http://forum.meteoinfo.ru/viewtopic.php?f=13&t=48

Совет по внедрению русских надписей на карты ГРАДС

http://forum.meteoinfo.ru/viewtopic.php?f=13&t=51

SergejKo
Участник
Письмо
Russian Federation
# Дата: 27 Июл 2012 23:54


Corvus
Спасибо! Давно интересовался, как можно делать надписи на русском в Градсе.

Gennady
Участник
Письмо
# Дата: 28 Июл 2012 08:49


Corvus
Совет по внедрению русских надписей на карты ГРАДС

Спасибо!

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 28 Июл 2012 15:24


Gennady
SergejKo

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

Например, чтобы написать на карте "Фактические грозовые очаги", в скрипте надо написать "Fakti]eskie grozov#e o]agi".

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 31 Июл 2012 12:47


Оказывается, у НАСА есть модель Земли

GEOS-5 Earth System Modeling and Data Assimilation

http://geos5.org/wiki/index.php?title=GEOS-5_Earth _System_Modeling_and_Data_Assimilation

В том числе она включает и модель атмосферы. Самое интересное, что эти данные доступны для использования в ГРАДС через ДОДС-сервер:

http://opendap.nccs.nasa.gov:9090/dods/GEOS-5

В частности, глобальные прогностические метеоданные на 5 дней с разрешением около 0.5 градуса:

http://opendap.nccs.nasa.gov:9090/dods/GEOS-5/lega cy_fp/0.5_deg/fcast/tavg2d_met_x.latest.info

http://opendap.nccs.nasa.gov:9090/dods/GEOS-5/lega cy_fp/0.5_deg/fcast/inst3d_met_p.latest.info

Надо разбираться, клон ли это ГФС, или принципиально новая глобальная модель, ещё неведомая нам :)

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 31 Июл 2012 13:07 - Поправил: Corvus


Некоторые прогностические характеристики облачности и осадков (включая КАПЕ) есть даже с разрешением 0.25 градуса

http://opendap.nccs.nasa.gov:9090/dods/GEOS-5/fp/0 .25_deg/fcast/inst3_2d_met_Nx.latest.info

И ряд других метеовеличин тоже с таким разрешением

http://opendap.nccs.nasa.gov:9090/dods/GEOS-5/fp/0 .25_deg/fcast/tavg1_2d_slv_Nx.latest.info

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 31 Июл 2012 13:56


Глобальная модель FNMOC, одноградусная. Как я понимаю, это ансамбль на основе модели НОГАПС.

http://nomads.ncep.noaa.gov:9090/dods/fens/fens201 20731/fens_all_00z.info

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 31 Июл 2012 14:01


Глобальная модель GENS, одноградусная, - ансамбль на основе ГФС.

http://nomads.ncep.noaa.gov:9090/dods/gens/gens201 20731/gep_all_00z.info

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 31 Июл 2012 14:04


Оценка интенсивности осадков со спутников Метеосат-7 и Метеосат-9. Данные через каждые полчаса. Формат ГРИБ2.

http://oiswww.eumetsat.int/~idds/html/grib.html

Mesocyclon
Участник
Письмо
Троицк, Ключевка, Челябинская область
# Дата: 1 Авг 2012 16:26 - Поправил: Mesocyclon


Corvus

Некоторые прогностические характеристики облачности и осадков (включая КАПЕ) есть даже с разрешением 0.25 градуса

Спасибо! Модель абсолютно другая, это видно по картам.

Скрипт:

'reinit'

'sdfopen http://opendap.nccs.nasa.gov:9090/dods/GEOS-5/fp/0 .25_deg/fcast/inst3_2d_met_Nx.latest'

'q ctlinfo'
_ctl = result

maps = 62

tsize = subwrd(_tdef,2)
_t1 = 1 ;* 2nd half of timeseries
_t2 = 45
tsize = _t2 - _t1 + 1
'set t '_t1' '_t2
'q dims'
times = sublin(result,5)
_time1 = subwrd(times,6)
_time2 = subwrd(times,8)
_tdim = _time1' '_time2

tincr = subwrd(_tdef,5)
_tdef = 'tdef 'tsize' linear '_time1' 'tincr

huh = subwrd(_tdef,4)

i = 3
while ( i<maps )

setmap(1)
'set t 'i

'set mpdset hires'
'set mpt * off'
'set mpt 0 1 1 9'
'set mpt 1 58 1 4'


'set gxout shaded'
'set csmooth on'
'set cint 0.2'
'set cmin 0.0'
'set cmax 15.0'
'set clevs 0.5 1 2 4 8 16 32'
'set ccols 0 4 11 13 3 7 8 2 6 0'
'mask=sumg(1,t=1,t-0,2)-sumg(1,t=1,t-1,2)'
_trr = 'tloop(prectot*10000-mask*const(10000*prectot(t-1) ,0,-u))'
'd '_trr
'set gxout contour'
'set cstyle 3'
'set ccolor 2'
'set line 9 3 6'
'set clab off'
'set clevs 0.5 1 2 4 8 16 32'

'set csmooth on'
'd '_trr

'set gxout contour'
'set csmooth on'
'set ccolor 6'
'set clab on'
'set cthick 10'
'set clevs 100 300 600 1000 1500 2000 3000'
'set ccols 4 11 13 3 7 8 2 6 13 7'
'd cape'

'cbarm'

'q dims'
times = sublin(result,5)
hub = subwrd(times,6)

'set string 1 bc 6'
'set strsiz 0.15 0.25'

'set strsiz 0.12'
'set string 1 r 3 90' ; 'draw string 9.9 4.6 mm'
'set string 1 r 3 0' ; 'draw string 9.45 0.6 Valid: 'hub
'set string 1 l 3 0' ; 'draw string 0.15 0.2 Data: GEOS-5 0.25 deg model, run: 'huh
'set string 1 r 11 0'
'set strsiz 0.14'
'set string 1 l 6 0' ; 'draw string 0.15 0.6 3hr Precip (mm) & CAPE (j/kg)'
'set string 1 l 5 0'
'set strsiz 0.1 0.1'
drawone(33.0497,68.9667 ,Murmansk )
.
.
.
.
drawone(55.1000,51.6831 ,Orenburg )
'printim gfs_slp_'i'.png x900 y675'

i = i+1
endwhile


function setmap(args)
'set mproj nps'
'set mpvals 45 82 44 65'
'set lon 32 97'
'set lat 41 70'
'set mpdset hires'
'set display color white'
'c'
'set grads off'
'set parea 0.02 9.54 0.8 8.50'
'set rgb 99 1 1 1'
'set line 99 1 6'


return

function drawone(arg1, arg2, arg3)
'q w2xy 'arg1' 'arg2
x=subwrd(result,3)
y=subwrd(result,6)
'draw string 'x' 'y+0.05' .'arg3
return

Сразу не понятно в каких единицах измеряются осадки. В скрипте исходные единицы я увеличил на 10000.

Образец



Видно явное завышение.

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 1 Авг 2012 19:26


Mesocyclon

Осадки там суммарные, как в японской модели. То есть для пересчёта на 3 ч надо из текущего срока вычитать предыдущий. Вы учли это?

В чём смысл этих строк?

'mask=sumg(1,t=1,t-0,2)-sumg(1,t=1,t-1,2)'
_trr = 'tloop(prectot*10000-mask*const(10000*prectot(t-1) ,0,-u))'


Mesocyclon
Участник
Письмо
Троицк, Ключевка, Челябинская область
# Дата: 1 Авг 2012 19:48


Corvus

Это как раз и есть пересчёт на 3 часа, у ГФС эта фишка работает.

Всёравно кажется что кол-во осадков завышается.

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 2 Авг 2012 09:26


Mesocyclon

Всё непонятно, зачем нужно это:

'mask=sumg(1,t=1,t-0,2)-sumg(1,t=1,t-1,2)'

ИМХО, вполне достаточно было бы написать так:

_trr = '10000*(prectot-prectot(t-1))'

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 3 Авг 2012 15:49


Сервис для экспорта данных вышеуказанной модели ГЕОС-5 в разных форматах файлов (ГРИБ, ГРАДС, НЕТ-ЦДФ и т.д.)

http://portal.nccs.nasa.gov/cgi-lats4d/webform.cgi ?&i=GEOS-5/fp/0.25_deg/fcast/inst3_2d_met_Nx.lates t

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 3 Авг 2012 15:57


http://acmg.seas.harvard.edu/geos/wiki_docs/geos5/ GEOS5_V1_File_Specification.pdf


Описание переменных этой модели

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 3 Авг 2012 18:20 - Поправил: Corvus


Список переменных для файлов tavg2d_met_x (1 file per time, 8 files per day: 01:30, 04:30, 07:30, 10:30, 13:30, 16:30,
19:30, 22:30 GMT)

Variable Name Description Units
ALBEDO Surface albedo fraction
ALBNIRDF Diffuse beam NIR surface albedo fraction
ALBNIRDR Direct beam NIR surface albedo fraction
ALBVISDF Diffuse beam VIS surface albedo fraction
ALBVISDR Direct beam VIS surface albedo fraction
BSTAR Surface buoyancy scale m s-2
CLDHGH High-level (above 400 hPa) cloud fraction fraction
CLDLOW Low-level (1000-700 hPa) cloud fraction fraction
CLDMID Mid-level (700-400 hPa) cloud fraction fraction
CLDTOT Total cloud fraction fraction
DISPH Displacement Height m
DTG Total change in skin temperature K
EFLUX Latent heat flux (positive upward) W m-2
EMIS Surface emissivity dimensionless
EVAP Surface evaporation kg m-2 s-1
FRLAKE Fraction of lake type in grid box fraction
FRLAND Fraction of land type in grid box fraction
FRLANDICE Fraction of land ice type in grid box fraction
FROCEAN Fraction of ocean in grid box fraction
GRN Vegetation greenness fraction fraction
GWETROOT Root zone soil wetness fraction
GWETTOP Top soil layer wetness fraction
HFLUX Sensible heat flux (positive upward) W m-2
LAI Leaf area index m2 m-2
LWGDWN Surface downward longwave flux W m-2
LWGNETCLR Net surface downward longwave flux assuming clear skyW m-2
LWGNET Net surface downward longwave flux at the ground W m-2
LWGUP Longwave flux emitted from surface (upward) W m-2
LWI Surface types 0=water, 1=land, 2=ice

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 3 Авг 2012 18:20


LWTUP Upward longwave flux at top of atmosphere W m-2
LWTUPCLR Upward longwave flux at top of atmosphere assuming
clear sky
W m-2
PARDF Surface downward photosynthetically active radiation
diffuse flux
W m-2
PARDR Surface downward photosynthetically active radiation
beam flux
W m-2
PBLH Planetary boundary layer height m
PRECANV Surface precipitation flux from anvils kg m-2 s-1
PRECCON Surface precipitation flux from convection kg m-2 s-1
PRECLSC Surface precipitation flux from large-scale kg m-2 s-1
PRECSNO Surface snowfall flux kg m-2 s-1
PRECTOT Total surface precipitation flux kg m-2 s-1
PS Time averaged surface pressure Pa
QV10M Specific humidity interpolated to 10 m above the
displacement height
kg kg-1

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 3 Авг 2012 18:21


QV2M Specific humidity interpolated to 2 m above the
displacement height
kg kg-1
RHOA Surface air density kg m-3
SLP Sea level pressure Pa
SNOMAS Snow mass as liquid water equivalent depth m
SNODP Snow depth m
SWGDWN Surface downward shortwave flux W m-2
SWGDWNCLR Surface downward shortwave flux assuming clear sky W m-2
SWGNET Net surface downward shortwave flux W m-2
SWGNETCLR Net surface downward shortwave flux assuming clear
sky
W m-2
SWTDWN Incident shortwave radiation at top of atmosphere W m-2
SWTUP Top of atmosphere outgoing shortwave flux W m-2
SWTUPCLR Top of atmosphere outgoing shortwave flux assuming
clear sky
W m-2
T10M Temperature interpolated to 10 m above the displacement
height
K
T2M Temperature interpolated to 2 m above the displacement
height
K
TAUGWX Eastward (zonal) gravity wave surface stress N m-2
TAUGWY Northward (meridional) gravity wave surface stress N m-2
TAUHGH Optical thickness of high clouds dimensionless
TAULOW Optical thickness of low clouds dimensionless
TAUMID Optical thickness of mid-level clouds dimensionless
TAUTOT Optical thickness of all clouds dimensionless
TAUX Eastward (zonal) surface wind stress N m-2
TAUY Northward (meridional) surface wind stress N m-2
TO3 Total Column Ozone Dobson
TPW Total precipitable water kg m-2
TROPP Tropopause pressure Pa
TROPQ Tropopause specific humidity kg kg-1

Corvus
Автор сайта
Письмо
Владимир (г. Байконур)
# Дата: 3 Авг 2012 18:22


TROPT Tropopause temperature K
TSKIN Skin temperature K
TTO3 Tropospheric Total Ozone Column Dobson
U10M Eastward (zonal) wind at 10 m above displacement
height
m s-1
U2M Eastward (zonal) wind at 2 m above the displacement
height
m s-1
U50M Eastward (zonal) wind at 50 m above displacement
height
m s-1
USTAR Surface velocity scale m s-1
V10M Northward (meridional) wind at 10 m above the
displacement height
m s-1
V2M Northward (meridional) wind at 2 m above the
displacement height
m s-1
V50M Northward (meridional) wind at 50 m above
displacement height
m s-1
Z0H Roughness length, sensible heat m

. 1 . 2 . 3 . 4 . 5 . 6 . 7 . 8 . 9 . 10 . >>
Ваш ответ

          Отменить *Что это?

 » Логин  » Пароль 
 
 
Полезная информация:



Поддержка: miniBB forum software © 2001-2017