Leuke analyse Bas!
Dank je
Even voor de anderen: per 100 K hPa (geen mb hè, foei!) smeltlaag heb je zo'n 5.2 mm neerslag nodig om de laag naar het vriespunt te laten afkoelen door het smelten en sublimeren van de sneeuwvlokken.
mb en hPa zijn toch volledig equivalent? Of zie ik iets over het hoofd...
Hier dan de kaartjes op basis van GFS06z. Init is vandaag 6 UTC dus het eerste frame correspondeert met morgen 12 UTC en het laatste met 3 UTC in de nacht van donderdag op vrijdag.
Het scriptje heb ik eronder gezet.
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