- 积分
- 57
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2021-10-16
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
初学者求助,主要为以下三个问题:
一、请问如何使左图(meteoinfolab绘制)箭头像右图(NCL绘制)那样顺眼一点?
二、如何在左图上添加如右图在右上角的参考图注?
三、图注箭头的值是怎么求的(有设置间隔取值,以及箭头大小和宽度,但不知对应参考箭头的值应该是多少)
四、colorbar上注明的单位怎么移动到如图二位置?
(ps:也不知道这样用quiverm计算水汽通量对不对,grib数据超出上传文件限制大小,无法上传,第一次发帖求助,不知道怎么搞,见谅)
- f1 = addfile('C:/CYGWIN/home/Administrator/year_tcwv.grb')
- pw = f1['Total_column_water_vapour_surface_Mixed_intervals_Average'][:,:,:]
- pwave = mean(pw, axis=0)
- axesm()
- xlim(60,110)
- ylim(10,40)
- levs = [ 1, 2, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50]
- cols = [(204,0,0),(255,51,51),(255,102,51),(255,153,102),(255,204,102),(255,255,153),(255,255,204),(204,255,255),(153,255,255),(102,204,255),(51,153,255),(51,51,255)]
- layer = contourfm(pwave, levs, colors=cols)
- colorbar(layer,label='$\r{kg}*\r{s}^{-2}
- )
- f2=addfile("C:/CYGWIN/home/Administrator/year_flux.grb","r")
- ew=f2['Vertical_integral_of_eastward_water_vapour_flux_entire_atmosphere_Mixed_intervals_Average'][:,:,:]
- sn=f2['Vertical_integral_of_northward_water_vapour_flux_entire_atmosphere_Mixed_intervals_Average'][:,:,:]
- u = mean(ew, axis=0)
- v = mean(sn, axis=0)
- uu = u[::20,::20]
- vv = v[::20,::20]
- layer = quiverm(uu, vv, colors='k', size=3, width=0.8)
- #legend(label='?kg/(m*s)')
- lat=array([27,27,33,33])
- lon=array([92,98,98,92])
- geoshow(lat,lon,displaytype='polygon',facecolor=None,edgecolor='k',size=1)
复制代码
|
-
-
|