| Автор |
Сообщение |
Corvus
Автор сайта
###Corvus###
|
# Дата: 6 Апр 2013 06:43
Cumulonimbus incus
А чужой код рассчитывает 112.139, будто такого условия и нет. Кто прав?
Скорее всего, у них условие не учитывается. Расчёт с условиями (когда должна быть выбрана та или иная формула) в ГРАДС сделать непросто, вот они и забили на это.
|
Cumulonimbus incus
Участник
Кишинёв, Молдова
|
# Дата: 6 Апр 2013 12:14
Corvus
Теперь вопрос со СВИТом можно закрывать. Я сделал карту этого индекса на 12 СГВ сегодняшнего дня и сравнил её с готовой картой на AccuWeather.com. Результат - почти полное совпадение:
 
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 6 Апр 2013 22:30
Cumulonimbus incus
Теперь вопрос со СВИТом можно закрывать. Я сделал карту этого индекса на 12 СГВ сегодняшнего дня и сравнил её с готовой картой на AccuWeather.com. Результат - почти полное совпадение
Хорошо! :)
Читали, кстати, моё сообщение про учёт сдвига ветра при прогнозе КД?
http://meteoclub.ru/index.php?action=vthread&forum =16&topic=3212&page=11#7
|
Cumulonimbus incus
Участник
Кишинёв, Молдова
|
# Дата: 17 Апр 2013 10:38
misha
Товарищ капитан, прошу дать справедливую оценку следующим картам. По Вашему решению могу пойти под трибунал.
(Шучу, конечно :-))
   
|
misha
Участник
Украина, восточный Донбасс
|
# Дата: 17 Апр 2013 22:16
Cumulonimbus incus
1. не видно географии (контур геогр. ориентиров)
2. слишком насыщенная заливка, из-за чего всю карту синхронно читать не удаётся
3. линии тока белого цвета теряются на бледной заливке.
Определённо скажу, что читать такую карту я если и смогу, то с большим дискомфортом.
|
TornadoF5
Участник
Харьков, Украина. (Игорь)
|
# Дата: 18 Апр 2013 18:37
Cumulonimbus incus
Отправил вам письмо.
|
bullterrier
Участник
Пермь
|
# Дата: 3 Май 2013 10:49 - Поправил: bullterrier
Алгоритм расчета SB CAPE по ПТсм в OpenGrads по данным модели GEM
*Сначала вычисляем давление на уровне конденсации(spr-приземное давление)
'pl = plcl(t2m+273.15, rh2m, spr)'
*Затем расчитываем коэффициент m, который показывает,, выше или ниже уровня конденсации находится каждый из уровней
'm1000 = pl-1000'
'm985 = pl-985'
'm970 = pl-970'
'm950 = pl-950'
'm925 = pl-925'
'm900 = pl-900'
'm875 = pl-875'
'm850 = pl-850'
'm800 = pl-800'
'm750 = pl-750'
'm700 = pl-700'
'm650 = pl-650'
'm600 = pl-600'
'm550 = pl-550'
'm500 = pl-500'
'm450 = pl-450'
'm400 = pl-400'
'm350 = pl-350'
'm300 = pl-300'
'm275 = pl-275'
'm250 = pl-250'
'm225 = pl-225'
'm200 = pl-200'
'm175 = pl-175'
'm150 = pl-150'
'm100 = pl-100'
'm50 = pl-50'
*Затем перевращаем положительные значения m в 1 при помощи сигмоидальной функции
'n1000 = (2/(1+exp((-1000)*m1000)))-1'
'n985 = (2/(1+exp((-1000)*m985)))-1'
'n970 = (2/(1+exp((-1000)*m970)))-1'
'n950 = (2/(1+exp((-1000)*m950)))-1'
'n925 = (2/(1+exp((-1000)*m925)))-1'
'n900 = (2/(1+exp((-1000)*m900)))-1'
'n875 = (2/(1+exp((-1000)*m875)))-1'
'n850 = (2/(1+exp((-1000)*m850)))-1'
'n800 = (2/(1+exp((-1000)*m800)))-1'
'n750 = (2/(1+exp((-1000)*m750)))-1'
'n700 = (2/(1+exp((-1000)*m700)))-1'
'n650 = (2/(1+exp((-1000)*m650)))-1'
'n600 = (2/(1+exp((-1000)*m600)))-1'
'n550 = (2/(1+exp((-1000)*m550)))-1'
'n500 = (2/(1+exp((-1000)*m500)))-1'
'n450 = (2/(1+exp((-1000)*m450)))-1'
'n400 = (2/(1+exp((-1000)*m400)))-1'
'n350 = (2/(1+exp((-1000)*m350)))-1'
'n300 = (2/(1+exp((-1000)*m300)))-1'
'n275 = (2/(1+exp((-1000)*m275)))-1'
'n250 = (2/(1+exp((-1000)*m250)))-1'
'n225 = (2/(1+exp((-1000)*m225)))-1'
'n200 = (2/(1+exp((-1000)*m200)))-1'
'n175 = (2/(1+exp((-1000)*m175)))-1'
'n150 = (2/(1+exp((-1000)*m150)))-1'
'n100 = (2/(1+exp((-1000)*m100)))-1'
'n50 = (2/(1+exp((-1000)*m50)))-1'
*А отрицательные значения превращаем в 0
'nn1000 = (abs(n1000)+n1000)/2'
'nn985 = (abs(n985)+n985)/2'
'nn970 = (abs(n970)+n970)/2'
'nn950 = (abs(n950)+n950)/2'
'nn925 = (abs(n950)+n925)/2'
'nn900 = (abs(n900)+n900)/2'
'nn875 = (abs(n875)+n875)/2'
'nn850 = (abs(n850)+n850)/2'
'nn800 = (abs(n800)+n800)/2'
'nn750 = (abs(n750)+n750)/2'
'nn700 = (abs(n700)+n700)/2'
'nn650 = (abs(n650)+n650)/2'
'nn600 = (abs(n600)+n600)/2'
'nn550 = (abs(n550)+n550)/2'
'nn500 = (abs(n500)+n500)/2'
'nn450 = (abs(n450)+n450)/2'
'nn400 = (abs(n400)+n400)/2'
'nn350 = (abs(n350)+n350)/2'
'nn300 = (abs(n300)+n300)/2'
'nn275 = (abs(n275)+n275)/2'
'nn250 = (abs(n250)+n250)/2'
'nn225 = (abs(n225)+n225)/2'
'nn200 = (abs(n200)+n200)/2'
'nn175 = (abs(n175)+n175)/2'
'nn150 = (abs(n150)+n150)/2'
'nn100 = (abs(n100)+n100)/2'
'nn50 = (abs(n50)+n50)/2'
*Вычисляем li по птсм для частицы, поднявшейся от земли до всех слоев
'l1000 = 'li('t2m','dewpoint','spr','q1000')
'l985 = 'li('t2m','dewpoint','spr','q985')
'l970 = 'li('t2m','dewpoint','spr','q970')
'l950 = 'li('t2m','dewpoint','spr','q950')
'l925 = 'li('t2m','dewpoint','spr','q925')
'l900 = 'li('t2m','dewpoint','spr','q900')
'l875 = 'li('t2m','dewpoint','spr','q875')
'l850 = 'li('t2m','dewpoint','spr','q850')
'l800 = 'li('t2m','dewpoint','spr','q800')
'l750 = 'li('t2m','dewpoint','spr','q750')
'l700 = 'li('t2m','dewpoint','spr','q700')
'l650 = 'li('t2m','dewpoint','spr','q650')
'l600 = 'li('t2m','dewpoint','spr','q600')
'l550 = 'li('t2m','dewpoint','spr','q550')
'l500 = 'li('t2m','dewpoint','spr','q500')
'l450 = 'li('t2m','dewpoint','spr','q450')
'l400 = 'li('t2m','dewpoint','spr','q400')
'l350 = 'li('t2m','dewpoint','spr','q350')
'l300 = 'li('t2m','dewpoint','spr','q300')
'l275 = 'li('t2m','dewpoint','spr','q275')
'l250 = 'li('t2m','dewpoint','spr','q250')
'l225 = 'li('t2m','dewpoint','spr','q225')
'l200 = 'li('t2m','dewpoint','spr','q200')
'l175 = 'li('t2m','dewpoint','spr','q175')
'l150 = 'li('t2m','dewpoint','spr','q150')
'l100 = 'li('t2m','dewpoint','spr','q100')
'l50 = 'li('t2m','dewpoint','spr','q50')
*Делим на каждом уровне LI на ПТсм слоя, в который вторглась частица, и умножаем на поправку nn (чтобы обрубить значения, лежащие ниже уровня конденсации)
'c1000 = (l1000/(q1000+273.15))*(-1)*nn1000'
'c985 = (l985/(q985+273.15))*(-1)*nn985'
'c970 = (l970/(q970+273.15))*(-1)*nn970'
'c950 = (l950/(q950+273.15))*(-1)*nn950'
'c925 = (l925/(q925+273.15))*(-1)*nn925'
'c900 = (l900/(q900+273.15))*(-1)*nn900'
'c875 = (l875/(q875+273.15))*(-1)*nn875'
'c850 = (l850/(q850+273.15))*(-1)*nn850'
'c800 = (l800/(q800+273.15))*(-1)*nn800'
'c750 = (l750/(q750+273.15))*(-1)*nn750'
'c700 = (l700/(q700+273.15))*(-1)*nn700'
'c650 = (l650/(q650+273.15))*(-1)*nn650'
'c600 = (l600/(q600+273.15))*(-1)*nn600'
'c550 = (l550/(q550+273.15))*(-1)*nn550'
'c500 = (l500/(q500+273.15))*(-1)*nn500'
'c450 = (l450/(q450+273.15))*(-1)*nn450'
'c400 = (l400/(q400+273.15))*(-1)*nn400'
'c350 = (l350/(q350+273.15))*(-1)*nn350'
'c300 = (l300/(q300+273.15))*(-1)*nn300'
'c275 = (l275/(q275+273.15))*(-1)*nn275'
'c250 = (l250/(q250+273.15))*(-1)*nn250'
'c225 = (l225/(q225+273.15))*(-1)*nn225'
'c200 = (l200/(q200+273.15))*(-1)*nn200'
'c175 = (l175/(q175+273.15))*(-1)*nn175'
'c150 = (l150/(q150+273.15))*(-1)*nn150'
'c100 = (l100/(q100+273.15))*(-1)*nn100'
'c50 = (l50/(q50+273.15))*(-1)*nn50'
*Зануляем отрицательные значения
'cc1000 = (abs(c1000)+c1000)/2'
'cc985 = (abs(c985)+c985)/2'
'cc970 = (abs(c970)+c970)/2'
'cc950 = (abs(c950)+c950)/2'
'cc925 = (abs(c950)+c925)/2'
'cc900 = (abs(c900)+c900)/2'
'cc875 = (abs(c875)+c875)/2'
'cc850 = (abs(c850)+c850)/2'
'cc800 = (abs(c800)+c800)/2'
'cc750 = (abs(c750)+c750)/2'
'cc700 = (abs(c700)+c700)/2'
'cc650 = (abs(c650)+c650)/2'
'cc600 = (abs(c600)+c600)/2'
'cc550 = (abs(c550)+c550)/2'
'cc500 = (abs(c500)+c500)/2'
'cc450 = (abs(c450)+c450)/2'
'cc400 = (abs(c400)+c400)/2'
'cc350 = (abs(c350)+c350)/2'
'cc300 = (abs(c300)+c300)/2'
'cc275 = (abs(c275)+c275)/2'
'cc250 = (abs(c250)+c250)/2'
'cc225 = (abs(c225)+c225)/2'
'cc200 = (abs(c200)+c200)/2'
'cc175 = (abs(c175)+c175)/2'
'cc150 = (abs(c150)+c150)/2'
'cc100 = (abs(c100)+c100)/2'
'cc50 = (abs(c50)+c50)/2'
*считаем энергию неустойчивости для каждого слоя по формулам площади трапеции
'a985 = ((cc1000+cc985)*(10*(gp985-gp1000)))/2'
'a970 = ((cc985+cc970)*(10*(gp970-gp985)))/2'
'a950 = ((cc970+cc950)*(10*(gp950-gp970)))/2'
'a925 = ((cc950+cc925)*(10*(gp925-gp950)))/2'
'a900 = ((cc925+cc900)*(10*(gp900-gp925)))/2'
'a875 = ((cc900+cc875)*(10*(gp875-gp900)))/2'
'a850 = ((cc875+cc850)*(10*(gp850-gp875)))/2'
'a800 = ((cc850+cc800)*(10*(gp800-gp850)))/2'
'a750 = ((cc800+cc750)*(10*(gp750-gp800)))/2'
'a700 = ((cc750+cc700)*(10*(gp700-gp750)))/2'
'a650 = ((cc700+cc650)*(10*(gp650-gp700)))/2'
'a600 = ((cc650+cc600)*(10*(gp600-gp650)))/2'
'a550 = ((cc600+cc550)*(10*(gp550-gp600)))/2'
'a500 = ((cc550+cc500)*(10*(gp500-gp550)))/2'
'a450 = ((cc500+cc450)*(10*(gp450-gp500)))/2'
'a400 = ((cc450+cc400)*(10*(gp400-gp450)))/2'
'a350 = ((cc400+cc350)*(10*(gp350-gp400)))/2'
'a300 = ((cc350+cc300)*(10*(gp300-gp350)))/2'
'a275 = ((cc300+cc275)*(10*(gp275-gp300)))/2'
'a250 = ((cc275+cc250)*(10*(gp250-gp275)))/2'
'a225 = ((cc250+cc225)*(10*(gp225-gp250)))/2'
'a200 = ((cc225+cc200)*(10*(gp200-gp225)))/2'
'a175 = ((cc200+cc175)*(10*(gp175-gp200)))/2'
'a150 = ((cc175+cc150)*(10*(gp150-gp175)))/2'
'a100 = ((cc150+cc100)*(10*(gp100-gp150)))/2'
'a50 = ((cc50+cc100)*(10*(gp50-gp100)))/2'
*суммируем все площади, получаем CAPE
'define capesfc = (9.81*(a985+a970+a950+a925+a900+a875+a850+a800+a750+a700+a650+a600+a550+a500+a450+a400+a350+a300+a275+a250+a225+a200+a175+a150+a100+a50))-0.001'
|
bullterrier
Участник
Пермь
|
# Дата: 3 Май 2013 10:51
Для расчета MUCAPE надо посчитать SB CAPE, и затем посчитать CAPE для частицы, поднявшейся от всех слоев, начиная от 1015 гПа и до 700 гПа, из них выбрать максимальное
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 26 Май 2013 09:57
Коллекция скриптов для ГРАДС, включает расчёт Q-вектора и изэнтропический анализ
http://www.met.wau.nl/education/atd/practical/func tions/library.html
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 26 Май 2013 09:58
Скрипт для расчёта Q-вектора, фронтогенеза, геострофического и агеострофического ветра
http://gradsaddict.blogspot.com/2013/05/script-dyn amic-calculates-temperature.html
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 26 Май 2013 09:59
Вообще полезный сайт для пользователей ГРАДС
http://gradsaddict.blogspot.com/
|
qwerty
Участник
Владимирская область, г.Киржач
|
# Дата: 16 Jul 2013 16:44
А сложный язык у этой программы?
|
qwerty
Участник
Владимирская область, г.Киржач
|
# Дата: 16 Jul 2013 19:38
Там какой то сервер надо скачать, а то без него карты не строятся. Чё это такое? Как устанавливать?
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 16 Jul 2013 19:59
Чё это такое?
qwerty
Не "чё", а что. У нас принято на русском общаться.
Поищите в гугле, полно документации по ГРАДС на английском языке. Автоперевод сайтов в Гугле и Яндексе тоже никто ещё не отменял.
|
qwerty
Участник
Владимирская область, г.Киржач
|
# Дата: 19 Jul 2013 19:44
А как запустить на компьютере dods-сервер? А то без него невозможно брать данные в интернете. Градс пишёт ошибку
|
Cumulonimbus incus
Участник
Кишинёв, Молдова
|
# Дата: 19 Jul 2013 20:21
qwerty
Писать скрипт. Другого способа нет.
'reinit'
prompt 'Enter GFS model run date (yyyymmdd) --> '
pull date
prompt 'Enter GFS model run hour (hh) --> '
pull hour
'sdfopen http://nomads.ncep.noaa.gov:9090/dods/gfs_hd/gfs_h d'date'/gfs_hd_'hour'z'
Его нужно сохранить в папке c:\Program Files\OpenGrADS\Contents\Resources\Scripts\ как текстовый файл с расширением .gs. Запускать из командной строки, просто введя имя файла.
Перед сохранением уберите пробелы в ссылке.
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 19 Jul 2013 20:40
qwerty
А как запустить на компьютере dods-сервер?
ДОДС-сервер надо не запустить, а соединиться с ним через интернет. ДОДС-сервер - это источник, распространяющий модельные данные.
Комп пользователя может требовать запустить Х-сервер, он нужен для работы на виндовом компе юниксовой программы (которой является ГРАДС). Как скачать и установить Х-сервер, написано на сайтах про ГРАДС.
|
qwerty
Участник
Владимирская область, г.Киржач
|
# Дата: 20 Jul 2013 00:17
Cumulonimbus incus
У меня при подключение к адресу пишет ошибку DAP-сервера
При заходе на адрес вот выдаёт
|
Igor
Участник
МО, Домодедово
|
# Дата: 20 Jul 2013 00:32
qwerty
А разве пишется "hour`z", а не "hour`s"?
|
Cumulonimbus incus
Участник
Кишинёв, Молдова
|
# Дата: 20 Jul 2013 10:05
qwerty
Заходить на сайт нужно не через браузер, а через командную строку самого ГрАДСа, запустив скрипт.
|
qwerty
Участник
Владимирская область, г.Киржач
|
# Дата: 20 Jul 2013 17:36
Cumulonimbus incus
Да тоже самое
|
Cumulonimbus incus
Участник
Кишинёв, Молдова
|
# Дата: 20 Jul 2013 18:46
qwerty
18-часового прогона ещё нет. 06-часовой последний на данный момент.
|
qwerty
Участник
Владимирская область, г.Киржач
|
# Дата: 20 Jul 2013 19:56
Cumulonimbus incus
Всё равно какие то ошибки идут
|
qwerty
Участник
Владимирская область, г.Киржач
|
# Дата: 21 Jul 2013 13:54
Разобрался с градсом немножко. Смог построить карту Li по Европе.
У меня есть по карте вопрос. Как наносить название?
|
qwerty
Участник
Владимирская область, г.Киржач
|
# Дата: 21 Jul 2013 15:07
Понял как ставить надписи. Команда draw title
|
qwerty
Участник
Владимирская область, г.Киржач
|
# Дата: 22 Jul 2013 16:24
А как построить карту Т850 в градсе?
|
Cumulonimbus incus
Участник
Кишинёв, Молдова
|
# Дата: 22 Jul 2013 16:31
qwerty
't850 = tmpprs(lev=850)-273.15'
'd t850'
|
qwerty
Участник
Владимирская область, г.Киржач
|
# Дата: 23 Jul 2013 18:18
А как можно сделать такую же заливку, вертикальную шкалу и цвет изолоний, как здесь?
http://www.meteogiornale.it/mappe/gfs/europa/tmp85 0_0.png
|
qwerty
Участник
Владимирская область, г.Киржач
|
# Дата: 23 Jul 2013 21:14
И как толщину у изолиний менять?
|
Cumulonimbus incus
Участник
Кишинёв, Молдова
|
# Дата: 25 Jul 2013 21:39
qwerty
как можно сделать такую же заливку, вертикальную шкалу и цвет изолиний
Сложно, но можно. Поставьте себе программку Instant Eyedropper. Потом наводите курсор на нужный цвет, и его значение в RGB автоматом в буфер помещается.
|
