- 积分
- 10703
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-11-22
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 我不's玩具 于 2017-4-30 20:21 编辑
做极地投影想画出弧形的区域 如图
经过很多日的尝试,我找到了解决方法
贴在这里吧
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
;
; Open a netCDF file and read some data.
;
a = addfile("$NCARG_ROOT/lib/ncarg/data/cdf/meccatemp.cdf","r")
t = a->t(0,:,:) ; Read first time step
wks = gsn_open_wks("png","draw_order")
res = True
res@cnFillOn =True
res@gsnDraw = False
res@gsnFrame = False
res@gsnMaximize = True
res@gsnAddCyclic =False ; Don't add longitude cyclic pt.
res@mpOutlineOn = True
;res@cnLineDrawOrder = "Draw" ; These resources can be one of
;res@cnFillDrawOrder = "Draw" ; "PreDraw", "Draw" or "PostDraw".
; res@mpOutlineDrawOrder = "Draw"
;res@mpFillDrawOrder = "PreDraw"
plot = gsn_csm_contour_map_polar(wks,t({0:90},:),res)
gres = True
gres@gsFillColor =-1 ;填充透明的颜色
gres@gsEdgesOn =True ;打开边界
gres@gsEdgeThicknessF=10
gres@gsFillIndex =0 ;填充类型
gres@gsLineThicknessF =7 ;
lat = (/72, 35, 35, 72, 72/)
lon = (/-45, -45, 5, 5, -45/)
dum = gsn_add_polygon(wks,plot,lon,lat,gres)
draw(plot)
frame(wks)
end
成功后的图
|
-
箭头指向处
|