码迷,mamicode.com
首页 > 其他好文 > 详细

MeteoInfoLab脚本示例:计算涡度、散度

时间:2015-07-07 16:15:41      阅读:603      评论:0      收藏:0      [点我收藏+]

标签:

用U/V分量数据计算涡度和散度,计算涡度的函数是hcurl,计算散度的函数是hdivg,参数都是U, V。

脚本程序:

f = addfile(D:/Temp/GrADS/model.ctl)
u = f[U][0,0,:,:]
v = f[V][0,0,:,:]
vort = hcurl(u, v)
divg = hdivg(u, v)
axesm()
mlayer = shaperead(D:/Temp/map/country1.shp)
geoshow(mlayer)
#layer = contourfm(vort, 20)
layer = contourfm(divg, 20)
#title(‘Vorticity‘)
title(Divergence)
colorbar(layer, orientation=horizontal, aspect=50)

技术分享

技术分享

MeteoInfoLab脚本示例:计算涡度、散度

标签:

原文地址:http://www.cnblogs.com/yaqiang/p/4626931.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!