[...]
In Utrecht inderdaad, indien het bevalt dan is de master een serieuze optie. Vind het ook een stuk interessanter dan de mechanica waar ik bij de Werktuigbouw mee te maken heb.
[...]
Ja, dat klinkt eigenlijk ook wel logisch. Ik heb nog gekeken naar de kosten van een HiRLAM licentie bij het KNMI, maar daar was ik snel van genezen :S.
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.
Nu de Noorse versie gebruiken, en eens kijken in hoeverre dat nauwkeurig is. Uiteindelijk is het natuurlijk de bedoeling om dit aan de neerslagsoortradar te koppelen waaraan ik werk.
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.
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
Bericht laatst bijgewerkt: 17-10-2009 20:04
