- 积分
- 15
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2021-5-15
- 最后登录
- 1970-1-1

|

楼主 |
发表于 2021-5-16 16:20:32
|
显示全部楼层
我的脚本是这样的,我也没看出来具体是那里出了错
fn1 = 'D:\\m\\uv\\20080820.000'
fn2= 'D:\\m\\height\\20080820.000'
fn3='D:\\m\\temp\\20080820.000'
f2 = addfile_micaps(fn2)
f1 = addfile_micaps(fn1)
f3 = addfile_micaps(fn3)
t1 = f1.gettime(0)
u = f1['U'][0,0,0:14,2:31]
v = f1['V'][0,0,0:14,2:31]
H=f2['var'][0,0,0:14,2:31]
T=f3['var'][0,0,0:14,2:31]
speed = sqrt(u*u+v*v)
levs=arange(-30,30,0.5)
axesm()
layer2 = contourfm( T,levs,cmap='MPL_coolwarm')
layer3 = contour(H)
geoshow('D:\MeteoInfo\MeteoInfo\map\cn_province')
geoshow('country')
layer = quiverm(u, v, speed, 10, size=8)
quiverkey(layer, 0.85, 0.80, 15, bbox={'edge':True, 'fill':True})
colorbar( layer2)
title(u'850hPa环流场',fontname=u'黑体',fontsize=20,bold=False,color='k') |
|