Свободный метео-софт

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

Свободный метео-софт

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

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

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

Троицк, Ключевка, Челябинская область
# Дата: 15 Jun 2012 17:23


Немного подправил скрипт(выше)

Скачать можно здесь
http://cliware.narod2.ru/scripts/ts.gs

Образцы:





Mesocyclon
Участник

Троицк, Ключевка, Челябинская область
# Дата: 15 Jun 2012 17:26 - Поправил: Mesocyclon


.

Cumulonimbus incus
Участник

Кишинёв, Молдова
# Дата: 18 Jun 2012 14:44 - Поправил: Cumulonimbus incus


Mesocyclon

Спасибо! Буду пробовать.
Только что нашёл функцию для расчёта лапласиана поля скалярных величин.

function laplace(field)
*------------------------------------------------- ---------------------
*------------------------------------------------- ----------------------
* Leo Kroon (leo.kroon@wur.nl) / Wageningen University MAQ
* 26/sep/03
* Function to calculate the Laplacian of a scalar field
* simple version
*
* Arguments:
* field = name of the scalar variable (input)
*
* laplac = name of the field holding the laplacian (output)
*------------------------------------------------- ----------------------
*-------------------- BEGINNING OF FUNCTION ----------------------------
*------------------------------------------------- ----------------------
*
* Define constants
*
"define degtorad = 3.14159265359/180"
"define earthradius = 6.37e6"
"define tcos = cos(lat*degtorad)"
*
* Horizontal increments
*
"define dy = cdiff(lat,y)*degtorad"
"define dx = cdiff(lon,x)*degtorad"
*
* Define first derivatives
*
"define dhy = cdiff("field",y)"
"define dhx = cdiff("field",x)"
"define dhdy = tcos*dhy/dy"
"define dhdx = dhx/dx"
*
* Define second derivatives
*
"define c2hy = cdiff(dhdy,y)"
"define c2hx = cdiff(dhdx,x)"
"define d2hy = tcos*c2hy/dy"
"define d2hx = c2hx/dx"
*
* Laplacian in spherical coordinates
*
"define laplac = (1/(earthradius*earthradius*tcos*tcos))*(d2hx+d2hy )"
*
* Variable laplac now holds 3-D laplacian field and its name is returned
* for use by the user. All other variables are undefined, to free up space.
*
"undefine degtorad"
"undefine earthradius"
"undefine tcos"
"undefine dy"
"undefine dx"
"undefine dhy"
"undefine dhx"
"undefine dhdy"
"undefine dhdx"
"undefine c2hy"
"undefine c2hx"
"undefine d2hy"
"undefine d2hx"
Return(laplac)


На сайте http://www.scs.netai.net/ с завтрашнего дня будут выкладываться обновлённые карты. В новой версии при расчётах UI (т.е. индекса Пескова) дефицит точки росы на уровне 500 гПа будет заменён средним дефицитом на уровнях 850-500 гПа. Об этом писал Corvus: http://meteoclub.ru/index.php?action=vthread&forum =9&topic=547&page=4#22.

bullterrier
Участник

Пермь
# Дата: 27 Jun 2012 20:42


Mesocyclon

Я в Вашем скрипте, который строит карты осадков, заменил maps=62 на maps=24, i=3 на i=0 и скрипт отказался работать. Подскажите, в чем у меня ошибка?

Mesocyclon
Участник

Троицк, Ключевка, Челябинская область
# Дата: 27 Jun 2012 21:56


bullterrier
Параметр i не должен быть равен 0, так как у GFS осадки считаеются с 3-х часов. Этот параметр вообще трогать не нужно.

bullterrier
Участник

Пермь
# Дата: 27 Jun 2012 22:07


Mesocyclon

Что-то я ступил)) Спасибо!

bullterrier
Участник

Пермь
# Дата: 28 Jun 2012 14:11 - Поправил: bullterrier


Mesocyclon
Cumulonimbus incus

Спасибо Вашим постам, благодаря им я освоил Grads а вечер и ночь с нуля. Сделал, например, расчет индексов K и TQ

'define t7 = (tmpprs(lev=700)-273.15)-(dewpt(tmpprs(lev=700),rh prs(lev=700))-273.15)'
'define t8 = (tmpprs(lev=850)-273.15)+(dewpt(tmpprs(lev=850),rh prs(lev=850))-273.15)'
'define tq = t8-1.7*(tmpprs(lev=700)-273.15)'
'define ki = t8-t7-(tmpprs(lev=500)-273.15)'
'set gxout shade2b'
'set cmin 20'
'set cmax 40'
'd ki'
'cbarn'
'set gxout contour'
'set cmin 10'
'set cint 1'
'd tq'
'draw title TQ Index (lines), K index (shaded)'
'set string 1 c 1 0'
'set strsiz 0.12'
'draw string 5.5 0.3 GFS Model Run: 'runtime' Valid: 'timeword''
'printim tqmd_'timex'.png white'
'c'


Cumulonimbus incus
Участник

Кишинёв, Молдова
# Дата: 28 Jun 2012 15:41 - Поправил: Cumulonimbus incus


bullterrier

Похожим образом и я осваивал его :-) Целый месяц имел дело с очень известным в этой ветке скриптом meteogram_gfs.gs. Изменял, пытался понять, что к чему. Короче, с этим скриптом и документацией осилил я ГрАДС. Сейчас для меня написание скриптов не мучение, как было вначале, а радость. И благодаря этому я могу делать карты таких индексов, которых нигде в Интернете нет.

bullterrier
Участник

Пермь
# Дата: 28 Jun 2012 17:56


Cumulonimbus incus

А где можно посмотреть имена переменных модели GFS?

Еще дополнил скрипт вероятностью грозы и прогнозом диаметра града по методу Заводченкова А.Ф.:

'set gxout shade2b'
'set cmin 0'
'd (5*ki-100+(25/3)*(6-no4lftxsfc))/2'
'cbarn'
'define g1 = 0.17*(tmpprs(lev=850)-273.15)'
'define g2 = 0.45*(-273.15+dewpt(tmpprs(lev=850),rhprs(lev=850) ))'
'define g3 = 0.05*(tmpprs(lev=850)-tmpprs(lev=400))'
'define gz = g1+g2+g3-0.58*no4lftxsfc-10.1'
'set gxout contour'
'set cmin 0'
'set cint 0.5'
'd gz'
'draw title pts (%) & hail diameter (sm)'
'set string 1 c 1 0'
'set strsiz 0.12'
'draw string 5.5 0.3 GFS Model Run: 'runtime' Valid: 'timeword''
'printim gzmd_'timex'.png white'
'c'


Cumulonimbus incus
Участник

Кишинёв, Молдова
# Дата: 28 Jun 2012 22:39


bullterrier

Список переменных ГФС-ки на DODS-сервере можно найти тут. Дату и время можно ставить любые за последний месяц. Правда, варианта для времени есть всего 4 - 00, 06, 12, 18, т.е. время прогонов ГФС.

http://nomads.ncep.noaa.gov:9090/dods/gfs_hd/gfs_h d20120628/gfs_hd_00z.info

Есть и другие серверы, вроде http://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/ prod/ или http://nomads.ncdc.noaa.gov/data/gfs4/ (очень полезный архив, рекомендую. Но с ним без пол-литра не разберёшься :-)).

Gennady
Участник

# Дата: 29 Jun 2012 04:04


bullterrier
Mesocyclon
Cumulonimbus incus

Спасибо Вашим постам, благодаря им я освоил Grads

Присоединяюсь к Вашей оценке постов Mesocyclon
Cumulonimbus incus. Спасибо им.

bullterrier
Участник

Пермь
# Дата: 29 Jun 2012 17:06


Cumulonimbus incus

Спасибо! У меня еще вот такой вопрос: как сделать так, чтобы ветер отображался в таком виде, как здесь: URL?

Corvus
Автор сайта

###Corvus###
# Дата: 29 Jun 2012 18:15


как сделать так, чтобы ветер отображался в таком виде, как здесь: URL?
bullterrier


'set csmooth on'
'set gxout shaded'
'v = mag(UGRD10m,VGRD10m)'
'd v'
'set gxout stream'
'set ccolor 0'
'set cstyle 2'
'd UGRD10m;VGRD10m'


bullterrier
Участник

Пермь
# Дата: 29 Jun 2012 19:05


Corvus

Большое спасибо!
В общем, я дополнил скрипт, который сделали Cumulonimbus incus и Mesocyclon. Теперь этот скрипт строит полный набор прогностических карт по модели ГФС на трое суток для Пермского края и прилегающих территорий. Строятся следующие карты:
- конвективных индексов (CAPE, TQ, Вайтинга, Totals Totals, MCS, Фатеева, Томпсона и др.)
- вероятности грозы и прогноза диаметра града по методу Заводченкова А.Ф.
- температуры и геопотенциала для уровней 1000, 850, 700, 500 гПа
- осадков и давления
- тропопаузы
- температуры и ветра у земли.

Спасибо всем, кто мне в этом помогал!
Скачать тут: URL

Corvus
Автор сайта

###Corvus###
# Дата: 10 Авг 2012 09:10


Добавлю по карте тумана. Есть такая книжечка "Прогноз погоды для авиации". Автор Джозеф. Дж. Джорж; стр. 296.
Очень интересное исследование. ;)

met2

Хе-хе, а ведь это перевод (1966 года) легендарной книги, в которой был предложен К-индекс прогноза гроз (его в СССР почему-то называли индексом Вайтинга, и до сих пор в СНГ так называют).

George, J. J., 1960: Weather Forecasting for Aeronautics. Academic Press, 673 pp.

Corvus
Автор сайта

###Corvus###
# Дата: 10 Авг 2012 16:02


"Прогноз погоды для авиации". Автор Джозеф. Дж. Джорж

Эта книга есть на Алиб.Ру за 2000 руб.

http://www.alib.ru/5_dzhozef_dzh_dzhordzhprognoz_p ogody_dlya_aviacii_w1t35707fcb5206ff6ef2d325a66710 2737629.html

George, J. J., 1960: Weather Forecasting for Aeronautics. Academic Press, 673 pp.



А эта на Амазоне за 25...72 доллара

http://www.amazon.com/gp/offer-listing/B000XB70OK/ sr=/qid=/ref=olp_tab_used?ie=UTF8&colid=&coliid=&c ondition=used&me=&qid=&seller=&sr=

Corvus
Автор сайта

###Corvus###
# Дата: 10 Авг 2012 16:13


George, J. J., 1960: Weather Forecasting for Aeronautics. Academic Press, 673 pp.


О, среди редакторов книги есть мифический Вайтинг (R. M. Whiting) :)

http://scinet.dost.gov.ph/union/ShowSearchResult.p hp?s=2&f=&p=&x=&page=&sid=1&id=Weather+forecasting +for+aeronautics&Mtype=BOOKS

Corvus
Автор сайта

###Corvus###
# Дата: 10 Авг 2012 16:20


http://oiswww.eumetsat.org/IPPS/html/GE/EUM.kml

Оценка мгновенной интенсивности осадков с геостационарных спутников (я недавно давал ссылку на эти же данные в формате ГРИБ2).

Впечатляет, если открыть этот файлик в Гугле Земля :)

Или вставить этот адрес в строку поиска в Картах Гугле.

spralex
Участник

г. Конотоп, Сумская обл., Украина
# Дата: 10 Авг 2012 16:46 - Поправил: spralex


Corvus

среди редакторов книги есть мифический Вайтинг

Тайну можно считать разгаданной? :)

Оценка мгновенной интенсивности осадков с геостационарных спутников

Это датчик MPE на спутнике MSG, есть еще вот такие карты (с наложением на ИК снимок):


Не всегда точен, но в случаях конвективных очагов для анализа помогает, особенно где нет радарных данных.

Еще бы вот эти данные (Convective Rainfall Rate проект SAF ) с лучшей привязкой найти, они поточнее:...


Viper
Участник

Украина,
# Дата: 10 Авг 2012 17:33


spralex
Можно ссылочку?

spralex
Участник

г. Конотоп, Сумская обл., Украина
# Дата: 10 Авг 2012 17:47


Viper

Ответил в личку.

Mesocyclon
Участник

Троицк, Ключевка, Челябинская область
# Дата: 17 Авг 2012 18:03 - Поправил: Mesocyclon


Скрипт для построения карты осадков по новой модели GEOS-5 c разрешением 0.25 градуса и шагом 1 час!!!
скачать тут
http://cliware.narod2.ru/scripts/GEOS5-PRECIP-1H_5 DAY.gs

Текстовый вариант

'reinit'
'sdfopen http://opendap.nccs.nasa.gov:9090/dods/GEOS-5/fp/0 .25_deg/fcast/tavg1_2d_flx_Nx.latest'
'q ctlinfo'
_ctl = result
maps = 121
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'
_trr = '3600*(prectot-prectot(t-1))'
'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'
'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)'
'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

пример


Mesocyclon
Участник

Троицк, Ключевка, Челябинская область
# Дата: 17 Авг 2012 18:28


Скрипт для построения карты давления на уровне моря по модели GEOS-5 c разрешением 0.25 градуса и шагом 1 час.
скачать тут
http://cliware.narod2.ru/scripts/GEOS5-SEA-1H_5DAY .gs
Текстовый вариант
'reinit'
'sdfopen http://opendap.nccs.nasa.gov:9090/dods/GEOS-5/fp/0 .25_deg/fcast/inst1_2d_hwl_Nx.latest'
'q ctlinfo'
_ctl = result
maps = 121
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 contour'
'set csmooth on'
'set ccolor 1'
'set clab on'
'set cint 2'
'd slp/100'
'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 Sea level pressure(mb)'
'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

пример


Mesocyclon
Участник

Троицк, Ключевка, Челябинская область
# Дата: 4 Окт 2012 21:00


Corvus
У Вас на сайте лежат исходники delphi. Как я понял это расшифровщик SYNOP?
http://meteocenter.net/currentmaps/bitmap/

Corvus
Автор сайта

###Corvus###
# Дата: 4 Окт 2012 21:24


Mesocyclon

Да, я когда-то начинал его делать на Делфи (в 2001 году), но потом переключился на веб-технологии (язык РНР), с тех пор они так и валяются без дела.

Так же там заготовки для программы рисования кольцовок. Всё это я потом тоже реализовал на РНР.

Mesocyclon
Участник

Троицк, Ключевка, Челябинская область
# Дата: 13 Окт 2012 11:43 - Поправил: Mesocyclon


Небольшая программа для загрузки карт с сайта http://www2.wetter3.de/Archiv/
Скачивает карты одного параметра за один год.
Папка с картами создаётся там же где находится сама программа.
Скачать программу
http://www.cliware.narod2.ru/2/maps.exe



Во время скачивания форма немного подвисает.

Например, если Вы запустили программу с рабочего стола и скачиваете карты за 2005 год, то на рабочем столе создастся папка с именем 2005. Туда будут сохраняться все изображения по мере скачивания.

Замечания и предложения по улучшению программы принимаются :)))

Mesocyclon
Участник

Троицк, Ключевка, Челябинская область
# Дата: 14 Окт 2012 09:08


последняя версия программы
http://cliware.narod2.ru/2/mapsdownload.exe

Mesocyclon
Участник

Троицк, Ключевка, Челябинская область
# Дата: 16 Окт 2012 17:14 - Поправил: Mesocyclon


Написал аналогичную программу по загрузке карт реанализа с сайта wetterzentrale.de.

Скриншот


Скачать программу
http://cliware.narod2.ru/2/wetterzentrale.exe

Просьба отписаться о работе проги...

artur
Участник

# Дата: 16 Окт 2012 17:53


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

Mesocyclon
Участник

Троицк, Ключевка, Челябинская область
# Дата: 16 Окт 2012 17:56


artur
Спасибо! :)

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

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

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


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