计算K指数

'reinit'
 
'open e:/mpv/data/1.ctl'
 
'set lon 90 125'
 
'set lat 10 35'
 
'set lev 850'
 
'define T850=tmpprs-273.16'
 
'define RH=rhprs'
 
'define Td850=T850-((14.55+0.114*T850)*(1-0.01*rh) + pow((2.5+0.007*T850)*(1-0.01*rh),3) + (15.9+0.117*T850)*pow((1-0.01*rh),14))'
 
'set lev 700'
 
'define T700=tmpprs-273.16'
 
'define RH=rhprs'
 
'define Td700=T700-((14.55+0.114*T700)*(1-0.01*rh) + pow((2.5+0.007*T700)*(1-0.01*rh),3) + (15.9+0.117*T700)*pow((1-0.01*rh),14))'
 
'set lev 500'
 
'define T500=tmpprs-273.16'
 
'define K=T850-T500+Td850-(T700-Td700)'
 
'set cint 2.5'
 
'set cthick 6'
 
'd k'

没有验证过，公式比较简单，应该没有问题。K指数计算的关键在于Td的计算，这里采用的是Michael Maxwell 的算法。
如果有误，敬请指正！