| Автор |
Сообщение |
Cumulonimbus incus
Участник
Кишинёв, Молдова
|
# Дата: 18 Jun 2012 15:17
Помогите с формулой для индекса Фауста! Весь форум излазил, все руководства в Интернете посмотрел, а там всё таблицы да графики. А для ГрАДСа нужна формула.
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 18 Jun 2012 21:09
Cumulonimbus incus
Формулы как таковой нет. При дефиците ноль градусов кривые на графике - влажные адиабаты. При дефиците больше нуля поправка пропорциональная дефициту.
|
Cumulonimbus incus
Участник
Кишинёв, Молдова
|
# Дата: 18 Jun 2012 21:49
Corvus
Я понял, но сразу же возникает вопрос: как эти товарищи его считают?
http://portal.chmi.cz/files/portal/docs/meteo/oa/p tu_grafy.html
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 18 Jun 2012 23:09
Cumulonimbus incus
Считать можно, просто скрипт будет громоздкий.
Например, по слоям толщиной в 50 гПа посчитать влажноад. градиент и так от 850 подняться до 500 гПа. Потом ввести поправку на средний дефицит в слое 850-500 гПа.
|
TornadoF5
Участник
Харьков, Украина. (Игорь)
|
# Дата: 19 Jun 2012 20:28
Интересная ситуация с индексами вышла. По прогнозу завтра в Одесской и Николаевской обл. К-инекс будет зашкаливать (около 37-40), при том, что в районе ХФ на Слобожанщине он какой-то вялый (около 30). Основная неустойчивость и осадки ожидаются на северо-востоке страны (в зоне вронта), зато К-индекс что-то вообще непоказательный.
Вопрос: с чем может быть связаны такие значения К-индекса, ведь гроз в тех местах завтра не будет? Обычно он хорошо коррелирует с другими параметрами и правильно отражает грозовые очаги, а сейчас всё наоборот. Ниже привожу прогност. карты К-индекса, ЛИ и АД на завтра по Вознесенску (Николаевская обл.) на срок 15z:
 
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 19 Jun 2012 20:33
TornadoF5
Я сделал вывод (и реализовал на сайте Азия), что вероятность грозы, определённую по К, нужно обязательно корректировать по КФ (индексу Фатеева).
Как-то так:
if( $v['KF'] >= 20 ) $v['PK'] = $v['PK'] + 20;
if( $v['KF'] >= 10 && $v['KF'] < 20 ) $v['PK'] = $v['PK'] + 10;
else if( $v['KF'] < 0 ) $v['PK'] = $v['PK'] - 20;
|
_Ilya_
Участник
г. Белгород
|
# Дата: 19 Jun 2012 20:34
TornadoF5
Да, смотрел этот индекс. Если ему верить, завтра будет дежа-вю вчерашнего, с позволения сказать, фронта.
|
TornadoF5
Участник
Харьков, Украина. (Игорь)
|
# Дата: 19 Jun 2012 20:36
_Ilya_
Если ему верить, завтра будет дежа-вю вчерашнего
Не надо нам такого опять! :)
|
TornadoF5
Участник
Харьков, Украина. (Игорь)
|
# Дата: 19 Jun 2012 20:38
Corvus
Как-то так:
if( $v['KF'] >= 20 ) $v['PK'] = $v['PK'] + 20;
if( $v['KF'] >= 10 && $v['KF'] < 20 ) $v['PK'] = $v['PK'] + 10;
else if( $v['KF'] < 0 ) $v['PK'] = $v['PK'] - 20;
Что это?
|
Sunspot
Участник
Ульяновск
|
# Дата: 19 Jun 2012 21:16
Что это?
TornadoF5
Если вероятность грозы, вычисленная по индексу Фатеева больше или равна 20, тогда вероятность по К=К+20, если по KF в интервале 10..20, тогда по К=К+10, иначе если KF<0 то К=К-20. Как то так, наверное. Извините, если не так понял, Владимир сам лучше объяснит.
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 19 Jun 2012 21:22
Как то так, наверное.
Sunspot
Да, почти так :)
Если индекс Фатеева больше или равен 20, тогда вероятность грозы по К увеличиваем на 20%, если индекс Фатеева в интервале 10..20, увеличиваем на 10%, иначе если индекс Фатеева<0 то вероятность грозы уменьшаем на 20%.
|
TornadoF5
Участник
Харьков, Украина. (Игорь)
|
# Дата: 19 Jun 2012 21:28
Sunspot
Corvus
Спасибо за разъяснения! :)
|
Cumulonimbus incus
Участник
Кишинёв, Молдова
|
# Дата: 19 Jun 2012 21:32 - Поправил: Cumulonimbus incus
TornadoF5
Есть ещё один показатель - средний дефицит точки росы в слое 850-500 гПа. По Фаусту, если он превышает 14 - гроз не будет. В "Руководстве по краткосрочным прогнозам погоды" 1986 г. вообще указаны значения 8...10 в качестве предельных.
Завтра он в указанных Вами областях будет равен 14...22, и лишь вечером местами упадёт до 12. Так что я бы дал по этому региону слабые ливневые осадки с возможными спорадическими грозовыми разрядами, учитывая, что индекс Пескова будет в районе 0...2.
|
spralex
Участник
г. Обнинск, Калужская обл.
|
# Дата: 20 Jun 2012 23:31
На известном ресурсе (модель WRF, тестируемая УкрНИГМИ) появились индексы MCAPE и MCIN.
Пример на субботу:
|
Cumulonimbus incus
Участник
Кишинёв, Молдова
|
# Дата: 21 Jun 2012 21:51
Corvus
Да, скрипт будет сложный. Представляю, сколько времени ГрАДС будет его считать...
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 21 Jun 2012 22:25
Cumulonimbus incus
Попытка не пытка. Главное, грамотно написать цикл, который будет считать Т по слоям с учётом влажноадиабатического градиента.
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 8 Jul 2012 12:41
Можно рассчитать F (индекс Фауста) приближённо.
1. Считаем ЭПТ850 (эквивалентно-потенциальная температура) по Т850, Тд850 и давлению 850 гПа.
2. Считаем ЭПТ500 по Т500, Тд500 и давлению 500 гПа.
3. Считаем псевдоЛИ = ЭПТ850 - ЭПТ500.
4. Рассчитываем поправку к псевдоЛИ в зависимости от среднего дефицита в слое 850-500 гПа:
if($d<6) $r = 0;
else if($d==6) $r = 2;
else if($d==7 || $d==8) $r = 3;
else if($d==9 || $d==10) $r = 4;
else if($d==11 || $d==12 || $d==13) $r = 5;
else if($d==14 || $d==15) $r = 6;
else if($d>=16) $r = 9;
5. Эту поправку прибавляем к псевдоЛИ.
Так я сделал расчёт F в модели НАЕФС.
|
bullterrier
Участник
Пермь
|
# Дата: 8 Jul 2012 14:27
Corvus
Подскажите пожалуйста, верна ли эта формула для подсчета Lifted Index?
Ew = 6.1078 * exp([(9.5939 * Tw) - 307.004]/[(0.556 * Tw) + 219.522])
E = Ew - 0.35 * (Td - Tw)
Tp = -1 * {[ln(E/6.1078) * 219.522] + 307.004} / {[ln(E/6.1078) * 0.556] - 9.59539}
LI= Tc(500mb) - Tp(500mb)
Tc(500mb) - это температура на уровне 500 гПа.
Tp(500mb) - это температура частицы, поднятой от земли до уровня 500 гПа.
Td - точка росы.
Tw - приземная температура (?)
Просто хотелось бы сделать его расчет по данным модели WRF, которую гоняют у нас в ПГНИУ.
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 8 Jul 2012 20:39
bullterrier
А можно ссылку на первоисточник?
Похоже, Tw - температура смоченного термометра, но не уверен.
|
bullterrier
Участник
Пермь
|
# Дата: 8 Jul 2012 21:35 - Поправил: bullterrier
А можно ссылку на первоисточник?
Corvus
URL
Мне LI нужен для подсчета других индексов и прогноза диаметра града по методу Заводченкова. Кстати, этот метод очень себя хорошо показал по модели GFS, на картах все очень наглядно получается.
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 9 Jul 2012 07:20
Мне LI нужен для подсчета других индексов и прогноза диаметра града по методу Заводченкова
bullterrier
Я понимаю, он много для чего нужен, включая индекс Пескова.
Жаль, что на настоящий первоисточник ссылки нет (а только на наш же форум), это бы многое прояснило. У меня опасение, что эти формулы не совсем те.
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 9 Jul 2012 09:09 - Поправил: Corvus
Ну вот, как я и опасался, это смесь бульдога с носорогом.
Вот это расчёт точки росы по психрометру:
'Saturation Vapor Pressure Wet
Ew = 6.1078 * exp([(9.5939 * Tw) - 307.004]/[(0.556 * Tw) + 219.522])
'Saturation Vapor Pressure Dry
Es = 6.1078 * exp([(9.5939 * Td) - 307.004]/[(0.556 * Td) + 219.522])
E = Ew - 0.35 * (Td - Tw) 'Actual Vapor Pressure
Relative Humidity
RH = (E / Es) * 100
Dew Point
Tp = -1 * {[ln(E/6.1078) * 219.522] + 307.004} / {[ln(E/6.1078) * 0.556] - 9.59539}
http://www.aprweather.com/pages/calc.htm
Где Tw - Т смоченного термометра, Td - Т сухого термометра, Tp - точка росы.
А в формуле для ЛИ Тр - это Т поднимающейся частицы и вычисляется она совсем иначе.
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 9 Jul 2012 09:21
Интересный индекс попался
DCI = T(850 mb) + Td(850 mb) - LI(sfc-500 mb)
http://www.crh.noaa.gov/lmk/soo/docu/indices.php
При его значении более 30 высока вероятность сильных гроз и ливней.
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 9 Jul 2012 09:30
Расчёт температуры смоченного термометра (нужен в некоторых методах прогноза шквалов и других конвективных явлений):
Wet Bulb Temperature
Variables:
Tc = Temperature in Degrees C
RH = Relative Humidity in form 88 (not 0.88)
Optional Variable (for more accuracy):
P = Pressure
or
Constant (with up to 0.2° inaccuracy):
P = 985
Tdc = ((Tc - (14.55 + 0.114 * Tc) * (1 - (0.01 * RH)) - ((2.5 + 0.007 * Tc) *
(1 - (0.01 * RH))) ^ 3 - (15.9 + 0.117 * Tc) * (1 - (0.01 * RH)) ^ 14))
E = (6.11 * 10 ^ (7.5 * Tdc / (237.7 + Tdc)))
WBc = (((0.00066 * P) * Tc) + ((4098 * E) / ((Tdc + 237.7) ^ 2) * Tdc)) / ((0.00066 * P) + (4098 * E) / ((Tdc + 237.7) ^ 2))
http://www.aprweather.com/pages/calc.htm
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 9 Jul 2012 09:48 - Поправил: Corvus
Бангладешский сайт про прогнозирование торнадо, анализируются метеоусловия наблюдавшихся случаев
http://bangladeshtornadoes.org/
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 9 Jul 2012 10:41
Считаем псевдоЛИ = ЭПТ850 - ЭПТ500.
Точнее сказать, это будет не псевдоЛИ, а индекс Шоултера, т.к. исходный уровень 850 гПа.
При расчёте ЛИ в качестве исходного берётся тот уровень (включая и приземный), на котором ЭПТ максимальная.
|
bullterrier
Участник
Пермь
|
# Дата: 9 Jul 2012 13:57
Corvus
Спасибо большое. Придется тогда без Lifted Index обойтись. А в WRF оказывается есть переменная, прогнозирующая град. Но ее надо тестировать, смотреть, как она работает. Хорошая модель, только вот переменных CAPE/LI/CIN не хватает для полного счастья.:)
|
bullterrier
Участник
Пермь
|
# Дата: 9 Jul 2012 14:09
А вероятность грозы буду считать по методу Вайтинга: 20-0% 40-100%, и сверху накладывать индексы Фатеева и TQ. Для первого раза, по моему, будет достаточно.
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 9 Jul 2012 16:46
Спасибо большое. Придется тогда без Lifted Index обойтись.
bullterrier
Ну почему, его можно рассчитать аналогично индексу Шоултера:
1. Считаем ЭПТ850 (эквивалентно-потенциальная температура) по Т850, Тд850 и давлению 850 гПа.
2. Считаем ЭПТ500 по Т500, Тд500 и давлению 500 гПа.
3. SI = ЭПТ850 - ЭПТ500.
Только ЭПТ брать не для 850 гПа, а максимальную из приземной и значений на всех уровнях от земли до 700 гПа.
|
Corvus
Автор сайта
###Corvus###
|
# Дата: 11 Jul 2012 14:03 - Поправил: Corvus
Кстати, тот индекс, который в хСССР принято называть индексом Вайтинга, на Западе называют К-индексом. Предложен он неким Джорджем в 1960 году в книге:
George, J. J., 1960: Weather and Forecasting for Aeronautics. Academic Press, 673 pp.
|
