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

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)