- 积分
- 15
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-3-15
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
做的是700hPa7个年份雨季结束前两候的平均风场图,纬度0~45,经度40~130做出的图片只有地图地图没有出现风场图,希望前辈们指导下。
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"
begin
u=new((/7,10,73,144/),float)
v=new((/7,10,73,144/),float)
uf_1989_4748=addfile("~/data/NCEP-daily/uwnd/uwnd.1989.nc","r")
uf_1993_5253=addfile("~/data/NCEP-daily/uwnd/uwnd.1993.nc","r")
uf_1999_5051=addfile("~/data/NCEP-daily/uwnd/uwnd.1999.nc","r")
uf_2003_5051=addfile("~/data/NCEP-daily/uwnd/uwnd.2003.nc","r")
uf_2007_4950=addfile("~/data/NCEP-daily/uwnd/uwnd.2007.nc","r")
uf_2009_5152=addfile("~/data/NCEP-daily/uwnd/uwnd.2009.nc","r")
uf_2011_5253=addfile("~/data/NCEP-daily/uwnd/uwnd.2011.nc","r")
vf_1989_4748=addfile("~/data/NCEP-daily/vwnd/vwnd.1989.nc","r")
vf_1993_5253=addfile("~/data/NCEP-daily/vwnd/vwnd.1993.nc","r")
vf_1999_5051=addfile("~/data/NCEP-daily/vwnd/vwnd.1999.nc","r")
vf_2003_5051=addfile("~/data/NCEP-daily/vwnd/vwnd.2003.nc","r")
vf_2007_4950=addfile("~/data/NCEP-daily/vwnd/vwnd.2007.nc","r")
vf_2009_5152=addfile("~/data/NCEP-daily/vwnd/vwnd.2009.nc","r")
vf_2011_5253=addfile("~/data/NCEP-daily/vwnd/vwnd.2011.nc","r")
u0 = short2flt(uf_1989_4748->uwnd(230:239,3,:,:))
u1 = short2flt(uf_1993_5253->uwnd(255:264,3,:,:))
u2 = short2flt(uf_1999_5051->uwnd(245:254,3,:,:))
u3 = short2flt(uf_2003_5051->uwnd(245:254,3,:,:))
u4 = short2flt(uf_2007_4950->uwnd(240:249,3,:,:))
u5 = short2flt(uf_2009_5152->uwnd(250:259,3,:,:))
u6 = short2flt(uf_2011_5253->uwnd(255:264,3,:,:))
v0 = short2flt(vf_1989_4748->vwnd(230:239,3,:,:))
v1 = short2flt(vf_1993_5253->vwnd(255:264,3,:,:))
v2 = short2flt(vf_1999_5051->vwnd(245:254,3,:,:))
v3 = short2flt(vf_2003_5051->vwnd(245:254,3,:,:))
v4 = short2flt(vf_2007_4950->vwnd(240:249,3,:,:))
v5 = short2flt(vf_2009_5152->vwnd(250:259,3,:,:))
v6 = short2flt(vf_2011_5253->vwnd(255:264,3,:,:))
lat = uf_1989_4748->u0&lat
lon = uf_1989_4748->u0&lon
u(0,:,:,:)=u0
u(1,:,:,:)=u1
u(2,:,:,:)=u2
u(3,:,:,:)=u3
u(4,:,:,:)=u4
u(5,:,:,:)=u5
u(6,:,:,:)=u6
v(0,:,:,:)=v0
v(1,:,:,:)=v1
v(2,:,:,:)=v2
v(3,:,:,:)=v3
v(4,:,:,:)=v4
v(5,:,:,:)=v5
v(6,:,:,:)=v6
u_houavg=dim_avg_n(u,1)
v_houavg=dim_avg_n(v,1)
u_y=dim_avg_n(u_houavg,0)
v_y=dim_avg_n(v_houavg,0)
u_y!0="lat"
u_y!1="lon"
u_y&lat=lat
u_y&lon=lon
v_y!0="lat"
v_y!1="lon"
v_y&lat=lat
v_y&lon=lon
wks = gsn_open_wks("x11","wind")
resources = True
resources@gsnDraw = True
resources@gsnFrame = True
nlon = dimsizes(lon)
nlat = dimsizes(lat)
resources@vcRefMagnitudeF=5.
resources@vcRefLengthF=1.
resources@vcMinDistanceF=1
resources@vcRefAnnoOn=True
resources@vcGlyphStyle="WindBarb"
resources@mpLimitMode="LatLon"
resources@mpMinLatF=0.
resources@mpMaxLatF=45.
resources@mpMinLonF=40.
resources@mpMaxLonF=130.
map = gsn_csm_vector_map (wks,u_y,v_y,resources)
end
|
|