- 积分
- 393
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-7-15
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2018-5-15 00:07:38
|
显示全部楼层
绘图部分代码:
;*********************************************************
wks = gsn_open_wks("pdf","wind")
res = True
gsn_define_colormap(wks,"amwg_blueyellowred")
pltres = True
mpres = True
; Plotting options for Wind Vectors
opts = res
; opts@FieldTitle = "Wind" ; overwrite Field Title
opts@NumVectors = 30 ; density of wind barbs
opts@vcLineArrowColor = "red" ; change vector color
opts@gsLineColor = 1
vector = wrf_vector(a,wks,uu(15,:,:),vv(15,:,:),opts)
;delete(opts)
opts1 = res
; opts1@FieldTitle = "Wind" ; overwrite Field Title
opts1@NumVectors = 30 ; density of wind barbs
opts1@vcMonoLineArrowColor = True
opts1@vcLineArrowColor = "blue"
vector1 = wrf_vector(a,wks,uu(60,:,:),vv(60,:,:),opts1)
;delete(opts1)
; MAKE PLOTS
plot = wrf_map_overlays(a,wks,(/vector,vector1/),pltres,mpres)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
end
|
|