K指数脚本gs

gzdsg

ldwloveyyz

'reinit'
'open f:\suoyouctl\suoyoushuju'
'set mpdset hires country1'
'set grads off'
'set lev 850'
'set t 605 972'
'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)'
;

ldwloveyyz

'reinit'
'open f:\suoyouctl\suoyoushuju'
'set mpdset hires country1'
'set grads off'
'set lev 850'
'set t 605 972'
'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)'
;

ldwloveyyz

处理的时FNL资料，如果时其他资料只需把相应温度、湿度变量替换下就可以了