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

最新推荐文章于 2023-02-15 16:22:38 发布

weixin_30588675

最新推荐文章于 2023-02-15 16:22:38 发布

阅读量1.7k

点赞数

文章标签： matlab python java

原文链接：http://www.cnblogs.com/yaqiang/p/4626931.html

版权

本文介绍了一个利用GrADS软件中的addfile、hcurl、hdivg等函数计算二维风场涡度和散度的脚本程序，并通过shapefile展示计算结果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

用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)