- 积分
- 563
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2020-10-20
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
求助下,用ncl出WRF结果的图时怎么才能只涂色地图里的,去除掉地图以外海的部分的填色?如图怎么把山东省外面的海的部分填色去除或者填成白色?
这是我的部分脚本(出PBL高度的)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
mpres@mpOutlineOn = False
mpres@mpFillOn = False
mpres@mpDataSetName = "Earth..4"
mpres@mpDataBaseVersion = "MediumRes"
mpres@mpOutlineOn = True
mpres@mpMaskAreaSpecifiers = (/"China"/)
mpres@mpOutlineSpecifiers = (/"China","China:Provinces"/)
mpres@mpGeophysicalLineColor = "Black"
mpres@mpLimbLineColor = "Blue"
mpres@mpNationalLineColor = "Black"
mpres@mpPerimLineColor = "Black"
mpres@mpUSStateLineColor = "Black"
mpres@mpNationalLineThicknessF = 2
mpres@mpProvincialLineThicknessF = 1
mpres@mpGridLonSpacingF = 1.0
mpres@mpGridSpacingF = 1
mpres@mpGridLineThicknessF = 0.01
mpres@mpNationalLineThicknessF = 2
mpres@mpProvincialLineThicknessF = 1
mpres@mpGeophysicalLineThicknessF= 2
mpres@mpGridAndLimbOn=True
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
do it = 0,ntimes-1,1 ; TIME LOOP
print("Working on time: " + times(it) )
res@TimeLabel = times(it) ; Set Valid time to use on plots
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
opts = res
opts@cnFillOn = True
opts@cnFillPalette = "BlueYellowRed"
opts@ContourParameters = (/ -700., 700., 100./)
gsn_define_colormap(wks,"BlueRedGray")
contour_pbl = wrf_contour(a,wks,pbl(0,:,:),opts)
delete(opts)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
plot = wrf_map_overlays(a,wks,(/contour_pbl/),pltres,mpres)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
end do
end
|
-
PBL高度图
|