发布MeteoInfo 1.0.8

MeteoInfo

发布MeteoInfo 1.0.8

修改了气象数据对话框，能够同时处理多个文件：


增加了一个脚本编辑器，可以方便地编写MeteoInfo脚本程序（Tools -> Script菜单）。这里除了可以运行之前利用MIApp编写的脚本以外，还可以对正在运行的MeteoInfo桌面软件进行控制。MeteoInfo桌面软件主窗体的对象是mipy，控制图形操作的是mipy.MapDocument，控制气象数据操作的是mipy.MeteoDataset。详细的东西以后慢慢探讨，这里给一个利用多个NC文件计算水汽通量散度的例子：



脚本程序也可以通过命令行窗口来运行：


脚本代码：

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. from MeteoInfoC.Layout import *
    
12. 
    
13. #---- Define variables
    
14. BaseDir = "C:\\Program Files (x86)\\MeteoInfo\"
    
15. MapDir = BaseDir + "Map\"
    
16. DataDir = "D:\\Temp\\nc\"
    
17. 
    
18. #---- Open data
    
19. DataAir = MeteoDataInfo()
    
20. DataUwnd = MeteoDataInfo()
    
21. DataVwnd = MeteoDataInfo()
    
22. DataRhum = MeteoDataInfo()
    
23. DataAir.OpenNCData(DataDir + "air.2011.nc")
    
24. DataUwnd.OpenNCData(DataDir + "uwnd.2011.nc")
    
25. DataVwnd.OpenNCData(DataDir + "vwnd.2011.nc")
    
26. DataRhum.OpenNCData(DataDir + "rhum.2011.nc")
    
27. 
    
28. #---- Meteo dataset
    
29. mDataset = mipy.MeteoDataset
    
30. mipy.MapDocument.ActiveMapFrame.RemoveMeteoLayers()
    
31. mDataset.RemoveAllMeteoData()
    
32. 
    
33. #---- Add meteo data
    
34. mDataset.AddMeteoData(DataAir)
    
35. mDataset.AddMeteoData(DataUwnd)
    
36. mDataset.AddMeteoData(DataVwnd)
    
37. mDataset.AddMeteoData(DataRhum)
    
38. 
    
39. #---- Calculation
    
40. #---- Set date to June 23
    
41. tIdx = 171
    
42. DataAir.TimeIndex = tIdx;
    
43. DataUwnd.TimeIndex = tIdx;
    
44. DataVwnd.TimeIndex = tIdx;
    
45. DataRhum.TimeIndex = tIdx;
    
46. #---- Set level to 700hPa
    
47. lIdx = 3
    
48. DataAir.LevelIndex = lIdx;
    
49. DataUwnd.LevelIndex = lIdx;
    
50. DataVwnd.LevelIndex = lIdx;
    
51. DataRhum.LevelIndex = lIdx;
    
52. #---- Get grid data
    
53. air = DataAir.GetGridData("air")
    
54. uwnd = DataUwnd.GetGridData("uwnd")
    
55. vwnd = DataVwnd.GetGridData("vwnd")
    
56. rhum = DataRhum.GetGridData("rhum")
    
57. #---- Calculate
    
58. prs = 700
    
59. g = 9.8
    
60. es = 6.112*DataMath.Exp(17.67*(air-273.16)/(air-29.65))
    
61. qs = 0.62197*es/(prs-0.378*es)
    
62. q = qs*rhum/100
    
63. qhdivg = DataMath.Hdivg(q*uwnd/g,q*vwnd/g)
    
64. qv = rhum*es/100
    
65. uv = DataMath.Magnitude(uwnd, vwnd)
    
66. uvq = uv*qv/(9.8*1000)
    
67. 
    
68. #---- Plot data
    
69. mDataset.SetDrawType("shaded")
    
70. mDataset.Display(qhdivg)
    
71. #---- Set MapLayout
    
72. aTime = DataAir.GetTime()
    
73. mipy.MapDocument.MapLayout.AddText("Water Vapor Flux Divergence - " + aTime.ToString("yyyy-MM-dd"), 320, 20)
    
74. mipy.MapDocument.MapLayout.PaintGraphics()

复制代码

棉花糖

更新很快

mofangbao

这次更新很给力！

善人/jw

顶起顶起~~

royalosyin

更新的可真够快的啊。

MeteoInfo

mofangbao 发表于 2012-9-21 12:23

                               
登录/注册后可看大图

这次更新很给力！

代码编辑控件还是用了ICsharpCode.TextEditor，表现还不错。代码提示、自动完成等功能实在没精力弄了。

♂雨已~○

王老师真是厉害啊

远山

呵呵，自己赶紧先藏一个带脚本编辑的版本

river

以前没用过，不过看论坛上感觉功能很强大，下下来学习一下

vortexroc

MeteoInfo 发表于 2012-9-21 17:18

                               
登录/注册后可看大图

代码编辑控件还是用了ICsharpCode.TextEditor，表现还不错。代码提示、自动完成等功能实在没精力弄了。

还请老师详细讲讲mipy的用法，比如添加地图，用mDataset.Display(qhdivg)出图后修改等值线的间隔和标注，这些操作是否能直接用mipy完成，还是要再定义miapp对象来操作？