| Автор |
Сообщение |
isamegrelo
Участник
Zugdidi
|
# Дата: 19 Апр 2020 12:32
CorvusCorax
'qu850 = spfhprs(lev=850)*ugrdprs(lev=850)'
'qv850 = spfhprs(lev=850)*vgrdprs(lev=850)'
'vimf850 = -1*hdivg(qu850,qv850)*100000'
Удельная влажность умножается на U- и V- ветер, а ведь функция hdivg возвращает значение только двух компонент u,v.
В классическом понимании это уже не дивергенция ветра :)
|
isamegrelo
Участник
Zugdidi
|
# Дата: 19 Апр 2020 12:45
CorvusCorax
В общем индекс довольно продвинутый, но его можно усовершенствовать путем подбора новых коэффициентов из множественной регрессионной модели.
|
CorvusCorax
Автор сайта
|
# Дата: 19 Апр 2020 13:13
функция hdivg возвращает значение только двух компонент u,v.
В классическом понимании это уже не дивергенция ветра :)
isamegrelo
В метеорологии дивергенцию обычно считают по двум горизонтальным осям без учёта вертикальной оси, пренебрегая ею. Я сомневаюсь, что глобальные модели могут див. по верт. оси посчитать. Скорее только мезомасштабные типа ВРФ, при разрешении 2 км и менее.
|
CorvusCorax
Автор сайта
|
# Дата: 19 Апр 2020 13:16
нужны данные о грозовых очагах в радиусе 16 км от Зугдиди за последние 5 лет.
http://meteocenter.asia/ts.php?p=37279
Вы можете экспортировать архив в Excel в таком формате?
Дата - Очаги
01.07.2015 - 2
isamegrelo
Попробую завтра сделать.
|
isamegrelo
Участник
Zugdidi
|
# Дата: 19 Апр 2020 13:22
Попробую завтра сделать.
CorvusCorax
Спасибо.
|
isamegrelo
Участник
Zugdidi
|
# Дата: 19 Апр 2020 13:31 - Поправил: isamegrelo
CorvusCorax
'tc = TMPPRS(LEV=850)-273'
'rh = RHPRS(LEV=850)'
'TD850 = tc-((14.55+0.114*tc)*(1-0.01*rh) + pow((2.5+0.007*tc)*(1-0.01*rh),3) + (15.9+0.117*tc)*pow((1-0.01*rh),14))'
'D850 = tc-TD850'
'tc = TMPPRS(LEV=700)-273'
'rh = RHPRS(LEV=700)'
'TD700 = tc-((14.55+0.114*tc)*(1-0.01*rh) + pow((2.5+0.007*tc)*(1-0.01*rh),3) + (15.9+0.117*tc)*pow((1-0.01*rh),14))'
'D700 = tc-TD700'
'T850 = TMPPRS(LEV=850)-273'
'T500 = TMPPRS(LEV=500)-273'
'tc = TMPPRS(LEV=500)-273'
'rh = RHPRS(LEV=500)'
'TD500 = tc-((14.55+0.114*tc)*(1-0.01*rh) + pow((2.5+0.007*tc)*(1-0.01*rh),3) + (15.9+0.117*tc)*pow((1-0.01*rh),14))'
'D500 = tc-TD500'
Интересно, а зачем в списке 2 переменные T850 и T500? Они далее не используются в расчете, только tc, при TD850 > T850, TD700 > T700, TD500 > T500
|
isamegrelo
Участник
Zugdidi
|
# Дата: 19 Апр 2020 15:35 - Поправил: isamegrelo
Одна из главных причин, из-за которой индексы не всегда точно прогнозируют грозы кроется в пренебрежительном отношении к CIN.
Ни в одном известном индексе он не учитывается, а зря. При CIN <-100 конвекция тормозится, а при <-300 вообще блокируется.
|
CorvusCorax
Автор сайта
|
# Дата: 19 Апр 2020 19:01
зачем в списке 2 переменные T850 и T500?
isamegrelo
Я уже не помню за давностью времени. Это же отрывок из большого скрипта. Возможно, далее Т850 и Т500 использовались для расчёта других параметров.
|
CorvusCorax
Автор сайта
|
# Дата: 19 Апр 2020 19:07 - Поправил: CorvusCorax
в пренебрежительном отношении к CIN.
Ни в одном известном индексе он не учитывается, а зря
isamegrelo
CIN важен для дневных внутримассовых ливней и гроз, когда конвекция термическая, то есть начинается от земли. Для тропической зоны это, наверное, актуально.
Но в умеренной зоне примерно половина гроз (а может и больше) связаны с фронтами, причём возникают они часто вечером, ночью или даже утром. Во всех этих случаях конвекция свободная, начинается с некоторого уровня в нижней или средней тропосфере, при этом CIN нет.
|
isamegrelo
Участник
Zugdidi
|
# Дата: 20 Апр 2020 12:20
Это же отрывок из большого скрипта
CorvusCorax
Верно, я не обратил внимание на последнюю строку.
'p = 6.5*T850 - 3.25*T500 - 3.36*D850 - 3.36*D700 - 0.105*D500 - 2.6*no4lftxsfc + 0.52*vimf + 0.46*shear73 - 43.5'
|
isamegrelo
Участник
Zugdidi
|
# Дата: 20 Апр 2020 12:26 - Поправил: isamegrelo
Решил заменить удельную влажность на влагосодержание и умножить на 0.00001, чтобы преобразовать большие числа.
'define t850 = tmpprs(lev=850) - 273.15'
'define rh = rhprs(lev=850)'
'define td850 = t850 - ((14.55 + 0.114 * t850) * (1 - 0.01 * rh) + pow((2.5 + 0.007 * t850) * (1 - 0.01 * rh),3) + (15.9 + 0.117 * t850) * pow((1 - 0.01 * rh),14))'
'define d850 = t850 - td850'
'define t700 = tmpprs(lev=700) - 273.15'
'define rh = rhprs(lev=700)'
'define td700 = t700 - ((14.55 + 0.114 * t700) * (1 - 0.01 * rh) + pow((2.5 + 0.007 * t700) * (1 - 0.01 * rh),3) + (15.9 + 0.117 * t700) *pow((1 - 0.01 * rh),14))'
'define d700 = t700 - td700'
'define t500 = tmpprs(lev=500) - 273.15'
'define rh = rhprs(lev=500)'
'define td500 = t500 - ((14.55 + 0.114 * t500) * (1 - 0.01 * rh) + pow((2.5 + 0.007 * t500) * (1 - 0.01 * rh),3) + (15.9 + 0.117 * t500) * pow((1 - 0.01 * rh),14))'
'define d500 = t500 - td500'
'define qu2m = pwatclm * ugrd10m'
'define qv2m = pwatclm * vgrd10m'
'define vimf2m = -1 * hdivg(qu2m,qv2m) * 100000'
'define qu850 = pwatclm * ugrdprs(lev=850)'
'define qv850 = pwatclm * vgrdprs(lev=850)'
'define vimf850 = -1 * hdivg(qu850,qv850) * 100000'
'define qu700 = pwatclm * ugrdprs(lev=700)'
'define qv700 = pwatclm * vgrdprs(lev=700)'
'define vimf700 = -1 * hdivg(qu700,qv700) * 100000'
'define vimf = 1200 * (vimf2m + vimf850 + vimf700)'
'define du = ugrdprs(lev=300) - ugrdprs(lev=700)'
'define dv = vgrdprs(lev=300) - vgrdprs(lev=700)'
'define shear73 = sqrt(du * du + dv * dv)'
'define p = 6.5 * t850 - 3.25 * t500 - 3.36 * d850 - 3.36 * d700 - 0.105 * d500 - 2.6 * no4lftxsfc + 0.52 * vimf + 0.46 * shear73 - 43.5'
'd p*0.00001'
|
CorvusCorax
Автор сайта
|
# Дата: 20 Апр 2020 13:06
Вы можете экспортировать архив в Excel в таком формате?
Дата - Очаги
01.07.2015 - 2
isamegrelo
Готово:
30-05-2011 - 81
31-05-2011 - 1
02-06-2011 - 4
18-06-2011 - 12
19-06-2011 - 11
20-06-2011 - 6
23-06-2011 - 8
02-07-2011 - 1
12-07-2011 - 24
15-07-2011 - 3
16-07-2011 - 5
23-07-2011 - 21
12-08-2011 - 45
13-08-2011 - 8
16-08-2011 - 28
01-09-2011 - 2
02-09-2011 - 91
03-09-2011 - 21
04-09-2011 - 79
09-09-2011 - 6
23-09-2011 - 2
24-09-2011 - 2
29-09-2011 - 28
30-09-2011 - 26
02-10-2011 - 7
12-10-2011 - 5
13-10-2011 - 47
19-10-2011 - 10
25-10-2011 - 16
07-11-2011 - 4
29-11-2011 - 4
10-01-2012 - 25
11-04-2012 - 1
13-05-2012 - 8
26-05-2012 - 2
01-06-2012 - 10
07-06-2012 - 22
08-06-2012 - 4
15-06-2012 - 3
23-06-2012 - 2
25-06-2012 - 16
26-06-2012 - 2
27-06-2012 - 29
28-06-2012 - 41
29-06-2012 - 111
30-06-2012 - 60
01-07-2012 - 20
09-07-2012 - 1
10-07-2012 - 19
14-07-2012 - 3
15-07-2012 - 31
17-07-2012 - 1
18-07-2012 - 5
29-07-2012 - 5
30-07-2012 - 1
01-08-2012 - 47
09-08-2012 - 48
10-08-2012 - 70
13-08-2012 - 128
14-08-2012 - 14
15-08-2012 - 3
16-08-2012 - 7
19-08-2012 - 134
29-08-2012 - 4
21-09-2012 - 17
22-09-2012 - 19
23-09-2012 - 24
04-10-2012 - 2
05-10-2012 - 20
09-10-2012 - 11
11-10-2012 - 11
12-10-2012 - 58
26-10-2012 - 1
01-11-2012 - 8
03-11-2012 - 3
05-12-2012 - 69
30-12-2012 - 47
25-01-2013 - 15
13-03-2013 - 3
23-04-2013 - 2
28-05-2013 - 52
01-06-2013 - 16
03-06-2013 - 4
06-06-2013 - 3
10-06-2013 - 53
17-06-2013 - 12
18-06-2013 - 91
19-06-2013 - 242
20-06-2013 - 32
21-06-2013 - 3
01-07-2013 - 4
02-07-2013 - 2
12-07-2013 - 20
14-07-2013 - 14
21-07-2013 - 27
28-07-2013 - 6
29-07-2013 - 70
30-07-2013 - 7
08-08-2013 - 25
12-08-2013 - 1
13-08-2013 - 20
19-08-2013 - 32
21-08-2013 - 3
23-08-2013 - 14
24-08-2013 - 38
25-08-2013 - 5
27-08-2013 - 4
30-08-2013 - 19
31-08-2013 - 110
01-09-2013 - 18
05-09-2013 - 1
15-09-2013 - 95
21-09-2013 - 8
24-09-2013 - 13
28-09-2013 - 179
29-09-2013 - 146
03-10-2013 - 4
27-11-2013 - 16
08-12-2013 - 2
09-12-2013 - 2
13-01-2014 - 40
05-04-2014 - 6
08-04-2014 - 1
28-04-2014 - 3
31-05-2014 - 1
02-06-2014 - 4
03-06-2014 - 14
10-06-2014 - 14
11-06-2014 - 18
12-06-2014 - 11
13-06-2014 - 3
29-06-2014 - 1
06-07-2014 - 4
07-07-2014 - 120
08-07-2014 - 163
14-07-2014 - 24
20-07-2014 - 1
24-07-2014 - 6
25-07-2014 - 9
27-07-2014 - 40
04-08-2014 - 2
05-08-2014 - 2
06-08-2014 - 7
10-08-2014 - 24
18-08-2014 - 1
19-08-2014 - 2
20-08-2014 - 21
23-08-2014 - 6
26-08-2014 - 2
05-09-2014 - 1
07-09-2014 - 1
10-09-2014 - 52
11-09-2014 - 199
15-09-2014 - 1
16-09-2014 - 34
28-09-2014 - 73
29-09-2014 - 68
04-10-2014 - 11
18-10-2014 - 1
07-02-2015 - 4
03-04-2015 - 4
17-05-2015 - 12
31-05-2015 - 12
06-06-2015 - 1
12-06-2015 - 28
15-06-2015 - 11
19-06-2015 - 1
24-06-2015 - 3
25-06-2015 - 48
28-06-2015 - 7
08-07-2015 - 12
11-07-2015 - 16
26-07-2015 - 3
05-08-2015 - 4
06-08-2015 - 102
16-08-2015 - 16
19-08-2015 - 12
20-08-2015 - 10
21-08-2015 - 33
24-08-2015 - 33
25-08-2015 - 11
31-08-2015 - 2
14-09-2015 - 60
03-10-2015 - 3
04-10-2015 - 3
07-10-2015 - 1
01-11-2015 - 8
08-11-2015 - 2
11-11-2015 - 31
04-12-2015 - 3
03-01-2016 - 14
20-04-2016 - 28
15-05-2016 - 2
07-06-2016 - 10
08-06-2016 - 24
24-06-2016 - 16
28-06-2016 - 3
30-06-2016 - 15
01-07-2016 - 2
03-07-2016 - 5
04-07-2016 - 8
05-07-2016 - 1
06-07-2016 - 4
07-07-2016 - 10
10-07-2016 - 20
11-07-2016 - 4
18-07-2016 - 4
19-07-2016 - 2
28-07-2016 - 19
30-07-2016 - 43
31-07-2016 - 2
10-08-2016 - 6
11-08-2016 - 42
14-08-2016 - 25
16-08-2016 - 47
17-08-2016 - 1
18-08-2016 - 21
26-08-2016 - 1
31-08-2016 - 140
01-09-2016 - 40
07-09-2016 - 7
14-09-2016 - 26
15-09-2016 - 60
21-09-2016 - 1
24-09-2016 - 110
28-09-2016 - 52
03-12-2016 - 11
04-12-2016 - 3
07-12-2016 - 8
08-05-2017 - 1
29-05-2017 - 2
30-05-2017 - 15
31-05-2017 - 15
10-06-2017 - 12
12-06-2017 - 3
15-06-2017 - 8
16-06-2017 - 20
19-06-2017 - 4
04-07-2017 - 2
07-07-2017 - 359
19-07-2017 - 55
22-07-2017 - 4
23-07-2017 - 119
24-07-2017 - 5
25-07-2017 - 1
14-08-2017 - 44
15-08-2017 - 2
17-08-2017 - 4
22-08-2017 - 15
23-08-2017 - 8
28-08-2017 - 34
29-08-2017 - 1
30-08-2017 - 3
31-08-2017 - 8
01-09-2017 - 11
06-09-2017 - 6
07-09-2017 - 32
23-09-2017 - 34
24-09-2017 - 1
26-09-2017 - 2
27-09-2017 - 12
14-10-2017 - 22
15-10-2017 - 20
22-10-2017 - 79
23-10-2017 - 45
23-11-2017 - 20
05-12-2017 - 57
09-05-2018 - 5
10-05-2018 - 2
21-05-2018 - 1
29-05-2018 - 1
01-06-2018 - 34
15-06-2018 - 1
18-06-2018 - 2
20-06-2018 - 4
22-06-2018 - 57
29-06-2018 - 3
04-07-2018 - 10
05-07-2018 - 4
07-07-2018 - 1
15-07-2018 - 4
16-07-2018 - 8
21-07-2018 - 1
24-07-2018 - 1
01-08-2018 - 5
02-08-2018 - 31
03-08-2018 - 30
04-08-2018 - 4
07-08-2018 - 143
08-08-2018 - 40
09-08-2018 - 9
10-08-2018 - 33
12-08-2018 - 12
19-08-2018 - 10
31-08-2018 - 21
06-09-2018 - 29
10-09-2018 - 3
14-09-2018 - 53
16-09-2018 - 1
17-09-2018 - 4
23-09-2018 - 15
25-09-2018 - 1
26-09-2018 - 2
28-09-2018 - 1
29-09-2018 - 1
02-10-2018 - 2
03-10-2018 - 63
24-10-2018 - 4
25-10-2018 - 1
12-11-2018 - 12
11-12-2018 - 2
11-04-2019 - 1
09-05-2019 - 1
07-06-2019 - 2
12-06-2019 - 15
17-06-2019 - 22
18-06-2019 - 70
20-06-2019 - 7
24-06-2019 - 4
27-06-2019 - 65
28-06-2019 - 183
29-06-2019 - 32
08-07-2019 - 1
12-07-2019 - 27
14-07-2019 - 1
15-07-2019 - 15
19-07-2019 - 1
21-07-2019 - 3
24-07-2019 - 9
25-07-2019 - 7
30-07-2019 - 2
04-08-2019 - 29
17-08-2019 - 9
20-08-2019 - 2
30-08-2019 - 29
02-09-2019 - 11
06-09-2019 - 58
07-09-2019 - 18
17-09-2019 - 37
28-09-2019 - 11
09-10-2019 - 1
11-10-2019 - 1
12-10-2019 - 5
30-10-2019 - 1
30-11-2019 - 9
27-12-2019 - 1
13-02-2020 - 3
|
isamegrelo
Участник
Zugdidi
|
# Дата: 20 Апр 2020 13:21
CorvusCorax
Огромнейшее спасибо!
|
CorvusCorax
Автор сайта
|
# Дата: 20 Апр 2020 13:43
Решил заменить удельную влажность на влагосодержание и умножить на 0.00001, чтобы преобразовать большие числа.
'define qu2m = pwatclm * ugrd10m'
'define qv2m = pwatclm * vgrd10m'
isamegrelo
Удельную влажность нельзя заменять на влагосодержание!
Фишка в том, чтобы считать конвергенцию влаги именно для каждого уровня.
Уд. влажность разная для каждого конкретного уровня.
А влагосодержание - это одно число, сумма всех влагосодержаний по уровням над данной точкой атмосферы.
|
isamegrelo
Участник
Zugdidi
|
# Дата: 20 Апр 2020 20:46
CorvusCorax
Кто сформировал индекс GFP? Как подбирались коэффициенты?
Можно улучшить его путем добавления новых уровней 925, 600 и 400. А также сделать табличный прогноз на 48 часов.
|
CorvusCorax
Автор сайта
|
# Дата: 21 Апр 2020 10:50
Кто сформировал индекс GFP? Как подбирались коэффициенты?
isamegrelo
Я.
Начало:
http://meteoclub.ru/index.php?action=vthread&forum =16&topic=3212&page=4#7
Продолжение:
http://meteoclub.ru/index.php?action=vthread&forum =16&topic=3212&page=4#19
Окончательный вид:
http://meteoclub.ru/index.php?action=vthread&forum =16&topic=3212&page=4#24
|
CorvusCorax
Автор сайта
|
# Дата: 21 Апр 2020 10:53
сделать табличный прогноз на 48 часов.
isamegrelo
Это тогда же и сделано, в 2012 году.
В прогнозах на сайте Метеоцентр.Азия есть колонка PTS (вероятность грозы), она как раз рассчитывается по методу GFP
http://meteocenter.asia/?m=g&p=35952
http://meteocenter.asia/?m=g&p=37279
|
isamegrelo
Участник
Zugdidi
|
# Дата: 21 Апр 2020 14:22
есть колонка PTS (вероятность грозы)
CorvusCorax
Это хорошо, летом буду оценивать качество такого прогноза.
|
isamegrelo
Участник
Zugdidi
|
# Дата: 21 Апр 2020 14:36
Что касается моего поиска эталона для предсказания гроз, то он оказался провальным. При равных метеоусловиях в Зугдиди вероятность грозы 50 на 50.
|
CorvusCorax
Автор сайта
|
# Дата: 22 Апр 2020 10:38
При равных метеоусловиях в Зугдиди вероятность грозы 50 на 50.
isamegrelo
Гроза - локальное явление. Даже при идеальных условиях грозы в большинстве случаев отмечаются не сплошь, а на какой-то части территории.
Удаление фактической грозы 16 км от точки - слишком строгий критерий, он для всех методов будет давать оправдываемость не выше 50-60% (за исключением супергрозовых регионов типа экваториальной Африки или Индонезии).
Для оценки методов лучше брать более мягкий критерий удаления фактической грозы от точки, например 50 км.
|
isamegrelo
Участник
Zugdidi
|
# Дата: 22 Апр 2020 13:50 - Поправил: isamegrelo
CorvusCorax
На этой карте вероятность грозы от индекса GFP? Во многих местах при >90% нет грозоотметок)
|
CorvusCorax
Автор сайта
|
# Дата: 22 Апр 2020 15:16
На этой карте вероятность грозы от индекса GFP? Во многих местах при >90% нет грозоотметок)
isamegrelo
1. Да.
2. Это естественно. Ни один метод не является панацеей. Он даёт фоновые условия, благоприятные для гроз. А дальше действует столько факторов, что и перечислить сложно, а не то что учесть.
Метеорология тем и отличается от математики, что в ней 2+2 не всегда 4, а часто 3 или 5, иногда даже 0 или 10 :)
|
CorvusCorax
Автор сайта
|
# Дата: 22 Апр 2020 15:21
К тому же в умеренной зоне, в отличие от тропиков, грозы кратковременны (обычно длятся полчаса - час), тогда как методы показывают благоприятные для гроз условия на протяжении многих часов (в среднем 6-12 ч, а порой и сутки напролёт).
Поэтому "мгновенная" карта лишь частично отражает ситуацию.
Вот сегодня спорадические грозы были на Южном Урале (там, где вероятность по карте высокая), а сейчас есть очаг на севере Узбекистана.
|
isamegrelo
Участник
Zugdidi
|
# Дата: 23 Апр 2020 14:24 - Поправил: isamegrelo
Ни один метод не является панацеей.
CorvusCorax
Может быть лучше исходить из абсолютных значений индекса?
А вероятность грозы только в виде контурных линий.
Вот как на этой карте, CAPE (заливка), а LI (контур)
|
isamegrelo
Участник
Zugdidi
|
# Дата: 23 Апр 2020 16:12
В GrADS можно создавать такие графики, прогноз индекса SWEAT для Зугдиди.
|
CorvusCorax
Автор сайта
|
# Дата: 24 Апр 2020 10:02
вероятность грозы только в виде контурных линий.
isamegrelo
Моя карта рассчитана не только на метеорологов, но и на широкий круг пользователей (пилоты и т.д.), в этом случае заливка вероятности гроз нагляднее и проще для быстрого восприятия, не надо вглядываться в изолинии.
|
isamegrelo
Участник
Zugdidi
|
# Дата: 24 Апр 2020 13:07
CorvusCorax
А можно сделать карту отдельно?
|
isamegrelo
Участник
Zugdidi
|
# Дата: 25 Апр 2020 14:14
Оказывается индекс GFP не ограничивается сотней, он почти как SWEAT, может быть и 500 :)
|
isamegrelo
Участник
Zugdidi
|
# Дата: 25 Апр 2020 14:15
Из этого следует что вероятность грозы 100% на карте Метеоцентра возможно равняется всего 20%
|
isamegrelo
Участник
Zugdidi
|
# Дата: 25 Апр 2020 14:17 - Поправил: isamegrelo
Вот точная копия этого скрипта:
http://meteoclub.ru/index.php?action=vthread&forum =16&topic=7521#8
'define tc850=tmpprs(lev=850)-273'
'define rh850=rhprs(lev=850)'
'define td850=tc850-((14.55+0.114*tc850)*(1-0.01*rh850)+po w((2.5+0.007*tc850)*(1-0.01*rh850),3)+(15.9+0.117* tc850)*pow((1-0.01*rh850),14))'
'define d850=tc850-td850'
'define tc700=tmpprs(lev=700)-273'
'define rh700=rhprs(lev=700)'
'define td700=tc700-((14.55+0.114*tc700)*(1-0.01*rh700)+po w((2.5+0.007*tc700)*(1-0.01*rh700),3)+(15.9+0.117* tc700)*pow((1-0.01*rh700),14))'
'define d700=tc700-td700'
'define tc500=tmpprs(lev=500)-273'
'define rh500=rhprs(lev=500)'
'define td500=tc500-((14.55+0.114*tc500)*(1-0.01*rh500)+po w((2.5+0.007*tc500)*(1-0.01*rh500),3)+(15.9+0.117* tc500)*pow((1-0.01*rh500),14))'
'define d500=tc500-td500'
'define qu2m=spfh2m*ugrd10m'
'define qv2m=spfh2m*vgrd10m'
'define vimf2m=-1*hdivg(qu2m,qv2m)*100000'
'define e850=6.11*exp((17.2694*(tmpprs(lev=850)-273))/(tmp prs(lev=850)-35.86))'
'define spfh850=0.622*e850*(rhprs(lev=850)/100)/850'
'define qu850=spfh850*ugrdprs(lev=850)'
'define qv850=spfh850*vgrdprs(lev=850)'
'define vimf850=-1*hdivg(qu850,qv850)*100000'
'define e700=6.11*exp((17.2694*(tmpprs(lev=700)-273))/(tmp prs(lev=700)-35.86))'
'define spfh700=0.622*e700*(rhprs(lev=700)/100)/700'
'define qu700=spfh700*ugrdprs(lev=700)'
'define qv700=spfh700*vgrdprs(lev=700)'
'define vimf700=-1*hdivg(qu700,qv700)*100000'
'define vimf=1200*(vimf2m+vimf850+vimf700)'
'define du=ugrdprs(lev=300)-ugrdprs(lev=700)'
'define dv=vgrdprs(lev=300)-vgrdprs(lev=700)'
'define shear=sqrt(du*du+dv*dv)'
'define t850=tmpprs(lev=850)-273'
'define t500=tmpprs(lev=500)-273'
'define gfp=6.5*t850-3.25*t500-3.36*d850-3.36*d700-0.105*d 500-2.6*no4lftxsfc+0.52*vimf+0.46*shear-43.5'
'd gfp'
|
