K指数计算方法

传说中的谁

很久没发帖了，看到最近家园里原创帖子比较少，所以整个简单的来糊弄一下大家，千万别贱笑哈
资料还是FNL，下载地址：http://dss.ucar.edu/datasets/ds083.2/
脚本如下：

1. 
   
2. 'reinit'
   
3. 'open e:/mpv/data/1.ctl'
   
4. 'set lon 90 125'
   
5. 'set lat 10 35'
   
6. 'set lev 850'
   
7. 'define T850=tmpprs-273.16'
   
8. 'define RH=rhprs'
   
9. '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))'
   
10. 'set lev 700'
    
11. 'define T700=tmpprs-273.16'
    
12. 'define RH=rhprs'
    
13. '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))'
    
14. 'set lev 500'
    
15. 'define T500=tmpprs-273.16'
    
16. 'define K=T850-T500+Td850-(T700-Td700)'
    
17. 'set cint 2.5'
    
18. 'set cthick 6'
    
19. 'd k'

复制代码

结果图

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

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'

笨笨不笨

支持原创，今晚哥睡沙发，呼呼……{:soso_e108:}

zsmnantong

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

zsmnantong

可不可以多给些物理量的脚本啊，如SI指数、CAPE值等。最好还有湿位涡、Q矢量。要求有点过分了
本人刚学习，又急需写文章。
这个网站好像要注册的吧，我试了几次都没注册上，提示邮件地址有什么错误的，记不得了。
1*1还是0.5*0.5啊？

Mouse

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

                               
登录/注册后可看大图

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

好学是一种好的学习态度，对于一些常用物理量的计算，要自己动手多写一写才会理解个中含义

mofangbao

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

                               
登录/注册后可看大图

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

超版很辛苦的...感动了

言深深

好样的再补充两句，要不会被当做图片党的{:soso_e153:}

mofangbao

zhoudog 发表于 2011-8-19 08:22

                               
登录/注册后可看大图

呵呵~ 偶也好久没发帖了~ 近来太忙~有空一定来~

有空来就行，像传说这样一来就是经典的，开个玩笑啦，加油~

runtarzan

笨笨不要笨笨

传说中的谁

runtarzan 发表于 2011-8-19 10:22

                               
登录/注册后可看大图

笨笨不要笨笨

是笨笨嫑笨笨，我的新浪围脖，欢迎围观