K指数计算方法

传说中的谁

zsmnantong 发表于 2011-8-19 06:13

                               
登录/注册后可看大图

超版主昨晚没睡吧，呵呵！今天来的早，正好用的上，谢谢分享啊！

没有啦，折腾这么个东西不至于一个晚上的。

传说中的谁

zsmnantong 发表于 2011-8-19 06:22

                               
登录/注册后可看大图

可不可以多给些物理量的脚本啊，如SI指数、CAPE值等。最好还有湿位涡、Q矢量。要求有点过分了
本 ...

其实我也在尝试编这些算法，编出来会跟大家共享的。湿位涡还好，地转Q矢量也不算复杂，但半地转Q矢量、非地转Q矢量、非地转湿Q矢量比较麻烦。另外，plotskew.gs里面包含了一些算法，比如CAPE，可以参考一下。GRADS安装目录的lib文件夹里有一些现成的gs，可以看看，很有用的。

传说中的谁

另外，即使有现成的gs，建议大家还是去看看公式，这样有助于理解其物理意义。

zsmnantong

是要理解物理意义，好好学习，向各位大侠看齐

天目神眉

笨笨依然是这样的给里啊

传说中的谁

天目神眉 发表于 2011-8-19 15:05

                               
登录/注册后可看大图

笨笨依然是这样的给里啊

俺是传说！

天目神眉

传说中的谁 发表于 2011-8-19 15:06

                               
登录/注册后可看大图

俺是传说！

不是同一个人

传说中的谁

天目神眉 发表于 2011-8-19 15:16

                               
登录/注册后可看大图

不是同一个人

改天跟你学习转模式，我装了之后还没正经转过呢

天目神眉

传说中的谁 发表于 2011-8-19 15:17

                               
登录/注册后可看大图

改天跟你学习转模式，我装了之后还没正经转过呢

哈哈哈哈，好，一起学习

zsmnantong

跟一下，K指数Gs。今天才调好的，学习超版主的公式。2.5*2.5的资料http://www.esrl.noaa.gov/psd/data/gridded/reanalysis/
'reinit'
'sdfopen f:\ncep\air.2011.nc'
'sdfopen f:\ncep\rhum.2011.nc'
'set lon 80 130'
'set lat 15 50'
'set time 12z12july2011'
'set grads off'
'set grid on'
'set lev 850'
'define T850=air.1-273.16'
'define Td850=T850-((14.55+0.114*T850)*(1-0.01*rhum.2) + pow((2.5+0.007*T850)*(1-0.01*rhum.2),3) + (15.9+0.117*T850)*pow((1-0.01*rhum.2),14))'
'set lev 700'
'define T700=air.1-273.16'
'define Td700=T700-((14.55+0.114*T700)*(1-0.01*rhum.2) + pow((2.5+0.007*T700)*(1-0.01*rhum.2),3) + (15.9+0.117*T700)*pow((1-0.01*rhum.2),14))'
'set lev 500'
'define T500=air.1-273.16'
'define K=T850-T500+Td850-(T700-Td700)'
'set cint 2.5'
'set cthick 4'
'd k'
'set mpdset cnriver'
'draw  map '
'printim f:\picture\K.png white'
'c'