t = temperatuur °C op 850hPa
td = dauwpunt °C op 850hPa
dwpk = dauwpunt °K in 850hPa
lev = 850 (hPa)
vapr= 6,112*exp((17,67*td)/(td+243,5)) = dampdruk
e= vapr*1,001+(lev-100)/900*0,0034
mixr= 0,62197*(e/(lev-e))*1000 = mengverhouding
tlcl= 1/(1/(dwpk-56)+log(t/dwpk)/800)+56 = temperatuur van LCL
theta=t*(1000/lev)0,286 = potentiële temperatuur
thetae=theta*exp((3,376/Tlcl-0,00254)*mixr*1,0+0,00081*mixr) = equivalente potentiële temperatuur
thetaw=-0,00000002610722610731*thetae4+0,000045737373737237*thetae3-0,0303049883448038*thetae2+9,103774048*thetae-746.83 = potentiële natte-bol temperatuur
Theta-W en Theta-E staan heel sterk in verhouding tot elkaar, zodanig dat je met een vierdegraads polynoom de boel kan uitrekenen. Nauwkeuriger is een iteratietechniek waarbij de uitkomst 'gegokt' wordt en de fout telkens minder wordt. Dat is vaak vertragend in een programma of zelfs onmogelijk te implementeren (zoals in GrADS). Deze formule levert resultaten op die op binnen 0,2°C nauwkeurig zijn bij waarden tussen -40 en +40°C.
Gr. Ben
td = dauwpunt °C op 850hPa
dwpk = dauwpunt °K in 850hPa
lev = 850 (hPa)
vapr= 6,112*exp((17,67*td)/(td+243,5)) = dampdruk
e= vapr*1,001+(lev-100)/900*0,0034
mixr= 0,62197*(e/(lev-e))*1000 = mengverhouding
tlcl= 1/(1/(dwpk-56)+log(t/dwpk)/800)+56 = temperatuur van LCL
theta=t*(1000/lev)0,286 = potentiële temperatuur
thetae=theta*exp((3,376/Tlcl-0,00254)*mixr*1,0+0,00081*mixr) = equivalente potentiële temperatuur
thetaw=-0,00000002610722610731*thetae4+0,000045737373737237*thetae3-0,0303049883448038*thetae2+9,103774048*thetae-746.83 = potentiële natte-bol temperatuur
Theta-W en Theta-E staan heel sterk in verhouding tot elkaar, zodanig dat je met een vierdegraads polynoom de boel kan uitrekenen. Nauwkeuriger is een iteratietechniek waarbij de uitkomst 'gegokt' wordt en de fout telkens minder wordt. Dat is vaak vertragend in een programma of zelfs onmogelijk te implementeren (zoals in GrADS). Deze formule levert resultaten op die op binnen 0,2°C nauwkeurig zijn bij waarden tussen -40 en +40°C.
Gr. Ben
Bericht laatst bijgewerkt: 20-11-2008 09:55
