计算K指数

MeteoInfo

参照传说中的谁写的K指数计算方法的帖子（K指数计算方法
http://bbs.06climate.com/forum.php?mod=viewthread&tid=1854&fromuid=106
），写了一个MeteoInfo的脚本程序来实现同样的功能，对比了结果还不错。


脚本程序如下：

1. #---- Import class libraries
   
2. import clr
   
3. clr.AddReferenceByPartialName("System.Windows.Forms")
   
4. clr.AddReferenceByPartialName("System.Drawing")
   
5. from System.Windows.Forms import *
   
6. from System.Drawing import *
   
7. clr.AddReference("MeteoInfoC.dll")
   
8. from MeteoInfoC import *
   
9. from MeteoInfoC.Data import *
   
10. from MeteoInfoC.Data.MeteoData import *
    
11. 
    
12. #---- Define variables
    
13. BaseDir = "C:\\Program Files (x86)\\MeteoInfo\"
    
14. MapDir = BaseDir + "Map\"
    
15. DataDir = "F:\\Temp\\grib\"
    
16. 
    
17. #---- Create MIApp object
    
18. myApp = MIApp()
    
19. #---- Open layers
    
20. myApp.OpenLayer(MapDir + "country1.shp")
    
21. myApp.SetLegendBreak("country1.shp",0,Color.Yellow,Color.Black,1,True,False,True)
    
22. #---- Open data
    
23. dataInfo = MeteoDataInfo()
    
24. dataInfo.OpenGRIBData(DataDir + "fnl_20110416_00_00")
    
25. 
    
26. #---- Get temperature grid data - 850 hPa
    
27. dataInfo.LevelIndex = 20
    
28. T850 = dataInfo.GetGridData("Temperature@pressure") - 273.16
    
29. #---- Get relative humidity grid data - 850 hPa
    
30. dataInfo.LevelIndex = 15
    
31. rh = dataInfo.GetGridData("Relative_humidity@pressure")
    
32. #---- Calculate Td
    
33. Td850 = T850-((14.55+0.114*T850)*(1-0.01*rh)+DataMath.Pow((2.5+0.007*T850)*(1-0.01*rh),3)+ \
    
34.   (15.9+0.117*T850)*DataMath.Pow((1-0.01*rh),14))
    
35. #---- Get temperature grid data - 700 hPa
    
36. dataInfo.LevelIndex = 17
    
37. T700 = dataInfo.GetGridData("Temperature@pressure") - 273.16
    
38. #---- Get relative humidity grid data - 700 hPa
    
39. dataInfo.LevelIndex = 12
    
40. rh = dataInfo.GetGridData("Relative_humidity@pressure")
    
41. #---- Calculate Td
    
42. Td700 = T700-((14.55+0.114*T700)*(1-0.01*rh) + DataMath.Pow((2.5+0.007*T700)*(1-0.01*rh),3) + \
    
43.   (15.9+0.117*T700)*DataMath.Pow((1-0.01*rh),14))
    
44. #---- Get temperature grid data - 500 hPa
    
45. dataInfo.LevelIndex = 13
    
46. T500 = dataInfo.GetGridData("Temperature@pressure") - 273.16
    
47. #---- Calculate K
    
48. K = T850-T500+Td850-(T700-Td700)
    
49. #K.SaveAsSurferASCIIFile(DataDir + "K_test.dat")
    
50. 
    
51. #---- Plot data
    
52. myApp.MeteoDataInfo = dataInfo
    
53. myApp.SetLegendScheme(DataDir + "k1.lgs")
    
54. myApp.SetDrawType("contour")
    
55. myApp.Display(K)
    
56. myApp.MoveLayerToTop("country1.shp")
    
57. myApp.ZoomLonLatEx(90,125,10,35)
    
58. #myApp.ZoomLonLatEx(0,360,-90.1,90.1)
    
59. aTime = dataInfo.GetTime()
    
60. text = "K Index - " + aTime.ToString("yyyy-MM-dd")
    
61. myApp.MapLayout.AddText(text, 320, 20)
    
62. myApp.MapLayout.ActiveLayoutMap.GridXDelt = 5
    
63. myApp.MapLayout.ActiveLayoutMap.GridYDelt = 5
    
64. myApp.MapLayout.ActiveLayoutMap.DrawGridLine = False
    
65. #myApp.ProjectLayers("+proj=moll+lon_0=105")
    
66. myApp.MapLayout.PaintGraphics()
    
67. Application.Run(myApp)

复制代码

mofangbao

MeteoInfo的插件功能再强大些就好啦，这样大家写的一些算法程序就能够容易给软件增加更丰富的功能了，像blog有很多插件可以让大家选择安装的那样，不过我是站着说话不腰疼的，哈哈，实现起来还是需要时间的，而且愿意写插件分享的人的多少也很重要

冬日情愫+

楼主您好，我注意到您的脚本程序中计算露点温度的公式，
Td850 = T850-((14.55+0.114*T850)*(1-0.01*rh)+DataMath.Pow((2.5+0.007*T850)*(1-0.01*rh),3)+
  (15.9+0.117*T850)*DataMath.Pow((1-0.01*rh),14))
经常在程序中看到这个公式但是它是来自哪里的呀？

MeteoInfo

冬日情愫+ 发表于 2013-3-10 17:26

                               
登录/注册后可看大图

楼主您好，我注意到您的脚本程序中计算露点温度的公式，
Td850 = T850-((14.55+0.114*T850)*(1-0.01*rh)+D ...

原始出处我也不知道，我没学过气象。

BDRUSH

好复杂的样子啊

夜深沉

脚本在哪里加？我是编程菜鸟~~

Smilely_怡

我想知道K指数里的蒸发量用潜在蒸发比较好还是用实际蒸发量比较好？急急急~

oydl1985

看着比GRADS的GS文件复杂！

MeteoInfo

oydl1985 发表于 2016-12-1 16:49
看着比GRADS的GS文件复杂！

这个帖子很老了，建议看此贴：MeteoInfoLab脚本示例：计算K指数
http://bbs.06climate.com/forum.p ... 846&fromuid=106
(出处: 气象家园)

还有：MeteoInfoLab脚本汇总贴
http://bbs.06climate.com/forum.p ... 151&fromuid=106
(出处: 气象家园)