- 积分
- 6467
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-11-6
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2016-12-12 23:31:00
|
显示全部楼层
您好,这是我的脚本,虽然图出来了,可是对过程还是不太理解,纯粹是按照ncl提示的错误进行修改的,所以图也不知道对不对,我想画纬向风速的纬向垂直剖面图,原来垂直上有27层,插值后就变93了,我也不知道对应的垂直坐标是什么,wrf_user_intrp3d函数里的p就是从wrfout里找的气压值,因为看官网上说p要和数据的维数相同。
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/wrf/WRFUserARW.ncl"
begin
f=addfile("c:/wrfout_d01_2000-01-24_12.nc","r")
u=wrf_user_getvar(f,"U",0)
p=wrf_user_getvar(f,"PB",0)
lon=wrf_user_getvar(f,"XLONG",0)
ps=fspan(1000,50,93)
ps@units="hpa"
wks=gsn_open_wks("png","u")
gsn_define_colormap(wks,"rainbow")
plane = (/30,10,30,60 /)
u_plane=wrf_user_intrp3d(u,p,"v",plane,90.0,False)
u_plane!0="lev"
u_plane!1="lon"
u_plane&lev=ps
u_plane&lon=lon(0,0:71)
res=True
plot=gsn_csm_pres_hgt(wks,u_plane,res)
end
|
|