U/V格点数据计算涡度、散度

MeteoInfo

MeteoInfo中的GridData类有很强的运算能力，这样就可以对格点数据进行计算分析。这个例子是用U/V格点数据计算风速、涡度、散度的程序。

#---- Import class libraries
import clr
clr.AddReferenceByPartialName("System.Windows.Forms")
clr.AddReferenceByPartialName("System.Drawing")
from System.Windows.Forms import *
from System.Drawing import *
clr.AddReference("MeteoInfoC.dll")
from MeteoInfoC import *
from MeteoInfoC.Data import *

#---- Define variables
BaseDir = "C:\\Program Files (x86)\\MeteoInfo\\"
MapDir = BaseDir + "Map\\"
DataDir = BaseDir + "Sample\\"

#---- Create MIApp object
myApp = MIApp()
#---- Open layers
myApp.OpenLayer(MapDir + "country1.shp")
myApp.SetLegendBreak("country1.shp",0,Color.Yellow,Color.Gray,1,True,False,True)
#---- Open data
myApp.OpenGrADSData(DataDir + "GrADS\\model.ctl")
myApp.SetDrawType("shaded")
#---- Get grid data
U = myApp.GetGridData("U")
V = myApp.GetGridData("V")
#---- Data calculation
#Speed = DataMath.Magnitude(U,V)
#Speed = DataMath.Sqrt(U * U + V * V)
#vort = DataMath.Hcurl(U,V)
digv = DataMath.Hdivg(U,V)
myApp.Display(digv)
#myApp.SetDrawType("vector")
#myApp.Display("U","V")
myApp.MoveLayerToTop("country1.shp")
myApp.ZoomLonLatEx(0,360,-90.1,90.1)
myApp.SetTitle("MeteoInfo script demo - Grid data calculation")
myApp.MapLayout.PaintGraphics()
#myApp.SaveFigure(".\\try.png")
Application.Run(myApp)





GrADS用同样的数据做出来的涡度、散度图：

咖啡好苦

meteinfo还是用不习惯，毕竟没有普及，学校还都是用的GRADS

mofangbao

我是专门来坐沙发的，MeteoInfo必将吸引更多用户！

artmunich

前几天自己摸索MeteoInfo的时候，还想着问楼主这样的问题呢，马上就去试验一下

MeteoInfo

咖啡好苦 发表于 2011-7-21 11:53

                               
登录/注册后可看大图

meteinfo还是用不习惯，毕竟没有普及，学校还都是用的GRADS

生活需要丰富多彩，就像咖啡有多种口味。MeteoInfo并非要取代GrADS，它是为了让这个领域更多彩而生的。

programbug

好好学习中。。。

MeteoInfo

sawnan 发表于 2012-2-2 16:31

                               
登录/注册后可看大图

请问版主，请问一下涡度散度是用哪个公式算出来的？ 我自己也用原始数据算涡度散度，但是算出来量级与grads ...

MeteoInfo中计算涡度和散度的方法是和GrADS一样的，当初参考了GrADS的代码。GrADS是开源软件，你可以去看看它的相关代码，关键是先弄明白cdiff函数。

sawnan

好的，谢谢。已经找到linux下的grads源码，正在查看中

sawnan

请问版主，请问一下涡度散度是用哪个公式算出来的？ 我自己也用原始数据算涡度散度，但是算出来量级与grads相差不多，等值线的分布相差挺大的。很郁闷。 能提供一下公式或者相关的计算代码么？

王小二

非常不错~~~~