Re : Leuk! Waar haal je de data vandaan? yr.no?
Bericht van: Sjoerd (Leidschendam, ZH) , 17-10-2009 20:10[...]
Dat zullen wel heel forse prijzen zijn ja... Ik vraag me af wie dan wel zo´n licentie hebben. Naar mijn weten gebruikt Weeronline gewoon GFS en draait MC haar (naar eigen wensen aangepaste) MM5 model.
Interessante vraag ja, ik zou het niet weten eigenlijk. Jezelf uit de markt prijzen?
[...]
Mooie plannen! Houd er wel rekening mee dat je om juiste voorspellingen te kunnen doen allereerst onderscheid zal moeten maken tussen situaties met frontale winterse neerslag, en situaties met buien in polaire lucht. Je zult voor beide naar andere parameters moeten kijken.
Polaire buien weet ik niet zoveel vanaf, frontale winterse neerslag heb ik me toevallig een tijdje geleden een klein beetje in verdiept. Wat me erg aansprak is het fenomeen 'melting effect'. Dit is het effect dat voor het smelten van sneeuw warmte benodigd is, welke wordt onttrokken aan de omgeving. Dit heeft tot gevolg dat in situaties van langdurige neerslag, met niet te veel warmteadvectie, het vriesniveau langzaam daalt. Als dit maar lang genoeg doorgaat kan het vriesniveau uiteindelijk de grond bereiken, waardoor natte sneeuw over kan gaan in drogere sneeuw. En er zodoende een pak sneeuw kan vallen, terwijl dit op basis van het verticale temperatuurprofiel in eerste instantie niet te verwachten was. Zie onderstaand artikel.
Hierover kwam ik een erg mooie presentatie van Alwin tegen, het ging over NESO; het aanpassen van Hirlam om de neerslagsoort te voorspellen. Toevallig deze week tegengekomen toen ik op zoek was naar het effect van neerslag op de temperatuur.
Voor de verwachting van neerslagsoort bij frontale neerslag zou ik dan ook gebruik maken van de zogenaamde smeltlaag. De intensiteit van de smeltlaag kun je berekenen door het gebied met positieve natteboltemperaturen verticaal te integreren, waardoor de eenheid voor smeltlaag [K mb] wordt. Hiervoor heb je uiteraard wel data nodig met een voldoende verticale resolutie. Naast de smeltlaag, kun je ook de dikte van de vrieslaag berekenen. En door situaties te bekijken waarbij een voldoende intense vrieslaag onder een smeltlaag ligt, zou je dan zelfs ijzel moeten kunnen voorspellen.
Onderstaand een klein GRADS scriptje dat ik een tijdje geleden heb geschreven. Een groot deel van de berekeningen is het uitrekenen van de natteboltemperatuur uit temperatuur, relatieve luchtvochtigheid en luchtdruk. Verder wordt in een while-loop (sorry voor de opmaak, die verdwijnt met kopieren), de som genomen over de positieve waarden van de natteboltemperatuur in de onderste luchtlaag, met stappen van 30 hPa in hoogte.
Groeten,
Bas
t=10
while (t<16)
'set t 't
i=0
'c'
'define prslvl=15'
'define somtemp=0'
'define prs = 0'
'define twsom = 0'
'define twpos = 0'
while (i<180)
'define T = tmp'i+30'_'i'mb-273.15'
'define avp = (rh'i+30'_'i'mb/100)*(0.611*exp((17.27*T)/(T+237.3)))'
'define Td = (116.9+237.3*log(avp))/(16.78-log(avp))'
'define gamma = 0.00000066*(pressfc-prslvl)'
'define delta = (4098*avp)/((Td+237.3)*(Td+237.3))'
'define Tw =((gamma*T)+(delta*Td))/(gamma+delta)'
'delta = (4098*avp)/((0.5*(Tw+Td)+237.3)*(0.5*(Tw+Td)+237.3))'
'define Twnew =((gamma*T)+(delta*Td))/(gamma+delta)'
'twpos = const(maskout(Twnew,Twnew),0,-u)'
'twsom = twsom + Twpos'
'prs = pressfc-prslvl'
'prslvl = prslvl+30'
i = i+30
endwhile
'set grid off'
'd (twsom*180)'
'draw title Dikte van smeltlaag [K hPa] Init+'3*t' uur'
'printim smeltlaag'3*t'.png x640 y480'
t=t+1
endwhile
Groeten,
Bas
