- 积分
- 1742
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-6-26
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 BDRUSH 于 2014-2-22 00:24 编辑
大家看看我这个图,我有2个问题,第一个问题是里面的网格能去掉么,第二个问题是横纵坐标为啥经纬度不显示出来,需要设置什么吗?
下面是脚本。
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
a1 = addfile("/cygdrive/D/Study/NCL/home/data/uwnd.1961.nc","r")
a2 = addfile("/cygdrive/D/Study/NCL/home/data/vwnd.1961.nc","r")
u = short2flt(a1->uwnd(989,5,:,:))
v = short2flt(a2->vwnd(989,5,:,:))
lon=a1->lon
lat=a1->lat
wks=gsn_open_wks("pdf","1961290-1961300")
res=True
res@vcMinFracLengthF=0.25
res@vcRefMagnitudeF=18.0
res@vcRefLengthF=0.02
res@vcMonoLineArrowColor=False
res@vcMinDistanceF =0.015
res@vfXCStartV=lon(0)
res@vfXCEndV=lon(dimsizes(lon)-1)
res@vfYCStartV=lat(0)
res@vfYCEndV=lat(dimsizes(lat)-1)
res@mpLimitMode="LatLon"
res@mpCenterLonF=-180
res@mpMinLonF=-240
res@mpMaxLonF=-120
res@mpMinLatF=0
res@mpMaxLatF=60
res@mpPerimOn=True
res@tiMainString=":F26:Wind Feild---1961"
res@tiXAxisString=":F25:longitude"
res@tiYAxisString=":F25:latitude"
res@lbLabelFont="times-roman"
res@pmLabelBarDisplayMode="Always"
res@lbPerimOn=False
res@pmLabelBarSide="Bottom"
res@lbOrientation="Horizontal"
res@gsnScale=True
res@vpXF=0.2
res@vpYF=0.8
res@vpWidthF=0.6
res@vpHeightF=0.4
cmap=(/(/1.00,1.00,1.00/),(/0.00,0.00,0.00/),(/.560,.500,.700/),(/.300,.300,.700/),(/.100,.100,.700/),(/.000,.100,.700/),(/.000,0.300,.700/),(/.000,0.500,.500/),(/.000,.700,.100/),(/.060,.680,.000/),(/.550,0.550,.000/),(/.570,0.420,.000/),(/.700,0.285,.000/),(/.700,0.180,.000/),(/.870,0.050,.000/),(/1.00,0.000,.000/),(/.700,0.700,.700/)/)
gsn_define_colormap(wks,cmap)
map=gsn_vector_map(wks,u,v,res)
end
|
|