爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 6491|回复: 3

如何将风场图叠加在地图上

[复制链接]
发表于 2017-4-16 19:29:23 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

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
QQ图片20170416193612.png
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2017-4-16 19:39:30 | 显示全部楼层
求助有大神前辈们么
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2017-4-16 19:47:48 | 显示全部楼层
还想请教四维的u,v风格式为short转化float这个步骤是否为多余
密码修改失败请联系微信:mofangbao
发表于 2018-1-2 18:12:15 | 显示全部楼层
请问你知道short转化float这个步骤是否为多余了吗
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表