Работа с GrADS

Полезные ссылки:
Метеоцентр.Азия - наш сайт с высокодетализированными прогнозами погоды по пунктам Pоссии и мира
Облегчённая версия Метеоклуба (для смартфонов)

Работа с GrADS

Сейчас в Метеоклубе:
Участников - 4 [ Alexandross, Doof, Spasatel, Mischel ]
Максимальное одновременное количество посетителей: 308 [2 Ноя 2013 22:42]
Гостей - 296 / Участников - 12

 - Начало - Ответить - Статистика - Pегистрация - Поиск -

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

Владимирская область, г.Киржач
# Дата: 25 Jul 2013 23:47 - Поправил: qwerty


Cumulonimbus incus

Получилось. Только как сделать шкалу вертикально, как там незнаю



qwerty
Участник

Владимирская область, г.Киржач
# Дата: 6 Авг 2013 20:03


Никто не знает, почему в этом скрипет идёт ошибка на 2 строке?

function main(args)
if (args = '')
prompt 'Enter forecast date (example, 20110802) --> '
pull date
prompt 'Enter forecast hour (example, 00 or 06 or 12 or 18) --> '
pull hour
else
date = subwrd(args,1)
hour = subwrd(args,2)
endif
'sdfopen http://nomads.ncep.noaa.gov:9090/dods/gfs/gfs20130 806/gfs_06z'
'set mpdset hires'
'set map auto'
'set mproj scaled'
'set lat 40 70'
'set lon 10 45'
'set t 1'
'q time'
times = sublin(result,1)
runtime = subwrd(times,3)
timex = 1
while (timex <= 65)
'set t ' timex
'define tdd5 = tmpprs(lev=500) - dewpt(tmpprs(lev=500),rhprs(lev=500))'
'define conv = -100000*fish_div(ugrd10m,vgrd10m)'
'define shear73 = shear(ugrdprs(lev=300),vgrdprs(lev=300),ugrdprs(le v=700),vgrdprs(lev=700))'
'define ui = -0.4*no4lftxsfc - 0.05*tdd5 + 0.4*conv - 0.07*shear73'
'set gxout shade2b'
'd ui'
'cbarn'
'q time'
times = sublin(result,1)
timeword = subwrd(times,3)
'draw title Peskov Index\GFS Model Run: 'runtime' Valid: 'timeword''
'printim ui_'timeword'.png white'
'c'
timex = timex + 1
endwhile

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 14 Авг 2013 16:37 - Поправил: qwerty


А вообще, как в градсе автоматизировать процесс создания карт?

Alter
Участник

Москва, ЮЗАО
# Дата: 14 Авг 2013 18:04 - Поправил: Alter


qwerty
if (args = '')
Я не знаю, какой здесь синтаксис, но равенство в сравнении величин в с-подобных языках обычно записываются знаком ==. В случае же вышеприведенного сравнения будет возвращен результат присваивания переменной args значения ``(удалось присвоить - true, не удалось - false), что может быть не всегда корректно и вызывать ошибку.
Впрочем, здесь вроде Basic, значит, мои выводы, похоже, не к месту :)

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 24 Авг 2013 18:25 - Поправил: qwerty


***

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 25 Авг 2013 14:16


Как сделать такую же проекцию карты?



Cumulonimbus incus
Участник

Кишинёв, Молдова
# Дата: 27 Авг 2013 22:25 - Поправил: Cumulonimbus incus


qwerty
Как сделать такую же проекцию карты?

Как-то так:

'set mproj nps'
'set lon -45 75'
'set lat 25 75'
'set mpvals -14 42 30 70'


Первая команда задаёт тип проекции (полярная стереографическая), 2 и 3 - общую площадь построения, а последняя - тот участок, что будет виден на экране.



http://iges.org/grads/gadoc/gradcomdsetmpvals.html

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 27 Авг 2013 22:40


Cumulonimbus incus

А вы как линии ветра сделали?

Cumulonimbus incus
Участник

Кишинёв, Молдова
# Дата: 27 Авг 2013 22:41


qwerty

'set gxout barb' (оперение) или 'set gxout stream' (линии тока).

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 27 Авг 2013 22:43


Cumulonimbus incus

А сам скрипт можете выложить?

Cumulonimbus incus
Участник

Кишинёв, Молдова
# Дата: 27 Авг 2013 23:29


qwerty

Карта не моя, это просто пример.

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 27 Авг 2013 23:59


Cumulonimbus incus

А вот как сделать, чтоб стрелки ветра не кучно находились и не загораживали остальное на карте?

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 15 Сен 2013 15:53


А как в Градсе аномалию к нормам подсчитавать?

Cumulonimbus incus
Участник

Кишинёв, Молдова
# Дата: 15 Сен 2013 16:08


qwerty

Нужен ГРИБ с нормами, причём желательно с тем же шагом сетки. Где такое достать, не знаю.

Corvus
Автор сайта

###Corvus###
# Дата: 15 Сен 2013 18:00


Аномалии некоторых метеовеличин есть в готовом виде среди параметров ГФС, например

gpa1000mb
** 1000 mb geopotential height anomaly [gpm]

gpa500mb
** 500 mb geopotential height anomaly [gpm]

http://nomads.ncep.noaa.gov:9090/dods/gfs_hd/gfs_h d20130915/gfs_hd_06z.info

Corvus
Автор сайта

###Corvus###
# Дата: 15 Сен 2013 18:01 - Поправил: Corvus


Если построить карту параметра (gpa500mb - gpa1000mb), это будет по сути карта аномалий средней Т нижнего 5-километрового слоя атмосферы (отличие будет только в коэффициенте).

Тан = (0.1*gpa500mb - 0.1*gpa1000mb)/2

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 15 Сен 2013 18:04


Cumulonimbus incus

А какая команда управляет цветом цифр на изотермах?

Cumulonimbus incus
Участник

Кишинёв, Молдова
# Дата: 15 Сен 2013 18:08


qwerty
А какая команда управляет цветом цифр на изотермах?

Никогда не видел такой команды.

А вот как сделать, чтоб стрелки ветра не кучно находились и не загораживали остальное на карте?

view-source:http://iges.org/grads/gadoc/gradfuncsk ip.html

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 15 Сен 2013 18:11


Cumulonimbus incus



Чтобы цвет значений давления на АТ-850 была, как тут

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 16 Сен 2013 19:59


А как в Градсе качество изображения улучшить?

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 21 Сен 2013 16:04


http://cfs.ncep.noaa.gov/pub/raid0/cfsv2/climo_cfs r_time/sd/

А там часом не нормы за 1981-2010 годов? А то использовать хочу их для построения карт

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 22 Сен 2013 18:40


Corvus
Если построить карту параметра (gpa500mb - gpa1000mb), это будет по сути карта аномалий средней Т нижнего 5-километрового слоя атмосферы (отличие будет только в коэффициенте).

Тан = (0.1*gpa500mb - 0.1*gpa1000mb)/2


И хорошая корреляция с очагами холода и тепла у земли?

dimon2344
Участник

# Дата: 16 Окт 2013 18:10


Мне необходимо сделать наноску станционных данных на карту в виде символов погоды. Как правильно это сделать?

Corvus
Автор сайта

###Corvus###
# Дата: 19 Окт 2013 11:30


qwerty
Тан = (0.1*gpa500mb - 0.1*gpa1000mb)/2

И хорошая корреляция с очагами холода и тепла у земли?



Более-менее хорошая. Но более надёжно использовать не ОТ500/1000, а Т850.

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 16 Ноя 2013 16:30 - Поправил: qwerty


http://www.wetterzentrale.de/pics/Rtavn062.png

Какая проекция у этой картинки и как сделать такое местоположение надписей, как тут?

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 20 Ноя 2013 19:06


А какая проекция нужна для построение карты Арктики?

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 1 Дек 2013 18:36 - Поправил: qwerty


Cumulonimbus incus
Совет по отображению новых государственных границ в ГрАДСе.

1. Скачиваем файл http://www.diva-gis.org/data/misc/countries_shp.zi p и распаковываем его в рабочую папку ГрАДСа. По умолчанию путь такой: c:\OpenGrADS\Contents\Resources\SampleDatasets\.



Файл удалён оттуда. Можете перезалить?

Cumulonimbus incus
Участник

Кишинёв, Молдова
# Дата: 1 Дек 2013 18:54


qwerty

Вот ссылка:

http://yadi.sk/d/4C_i9MOADcQcg

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 1 Дек 2013 18:57 - Поправил: qwerty


Cumulonimbus incus

А можете написать, какая проекция нужна для построение карты северного полюса? Нужно написать карту ЦПВ

Cumulonimbus incus
Участник

Кишинёв, Молдова
# Дата: 1 Дек 2013 19:01


qwerty

ЦПВ - это, как я понимаю, циркумполярный вихрь? Тогда так:

'set mproj nps'
'set lat 60 90'
'set lon -180 180'

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

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

 » Логин  » Пароль 
 
 


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