- 积分
- 31030
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2019-5-13
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
a = addfile("/mnt/d/ncl/20090605uv.nc","r")
u8 = a->u(8,{25:40},{110:130})/1000
v8 = a->v(8,{25:40},{110:130})/1000
wks = gsn_open_wks("X11","20090605wind_vector")
gsn_define_colormap(wks,"rainbow")
;--------------Set some basic resources-----------------
res = True
res@tiMainString = "2009_0605_07:00wind_vector"
res@gsnLeftString = ""
res@gsnRightString = ""
res@mpShapeMode ="FreeAspect"
res@vpXF = 0.1
res@vpYF = 0.8
res@vpHeightF = 0.7
res@vpWidthF= 0.7
res@mpGeophysicalLineColor = "Black"
res@mpNationalLineColor = "Black"
res@mpGridLineColor = "Black"
res@mpDataBaseVersion = "MediumRes"
res@mpDataSetName = "Earth..4"
res@mpGeophysicalLineThicknessF = 2
res@mpNationalLineThicknessF = 1
res@mpProvincialLineColor = "black"
res@mpProvincialLineThicknessF = 1
res@mpMinLatF=25
res@mpMaxLatF=40
res@mpMinLonF=110
res@mpMaxLonF=130
res@mpFillOn = True
res@mpOceanFillColor ="lightblue"
res@gsnAddCyclic = False
res@mpOutlineOn = True
;----设置填充箭头样式--------
res@vcGlyphStyle = "WindBarb"
res@vcRefAnnoOn = True ;风速单位标注
res@vcWindBarbScaleFactorF = 2.5 ; 风杆符合国内习惯(4m/s一个长杆)
res@vcRefMagnitudeF = 4. ; make vectors larger
res@vcRefLengthF = 0.040 ; ref vec length
plot = gsn_csm_vector_map(wks,u8,v8,res)
数据的分辨率是0.25X0.25的照理说在选定的25:40里纬度里应该有几十个值(确实对u8读取的pirnt的时候也是有的,但是为什么用画出来就只有几个箭头呀
|
|