关于提取wrfout固定高度站点数据的问题

羽陌轻寒

各位大佬，我在论坛查找用ncl提取wrfout风速数据的问题时，发现了有楼主分享了相关的脚本，具体链接如下：
http://bbs.06climate.com/forum.php?mod=viewthread&tid=59010&page=1&extra=#pid1020606


可是我在运行过程当中报错了，显示第六十一行出错，维度不匹配、第六十一行语句如下：windspdh=spdh(:,x,y）
上下文语句如下：
;-------------------interpolation----------------------
;wrf_user_intrp3d(三维数组，垂直数组（压力\高度），插值信息，TRUE=从point A点到point B的横截面图;否则为False。)
u_ph  = wrf_user_intrp3d( u,height,"h", hh,0.,False)
v_ph  = wrf_user_intrp3d( v,height,"h", hh,0.,False)
tch   = wrf_user_intrp3d( tc,height,"h", hh,0.,False)
tdh   = wrf_user_intrp3d( td,height,"h", hh,0.,False)
ph   = wrf_user_intrp3d(p,height,"h", hh,0.,False)
pbh  =wrf_user_intrp3d(pb,height,"h", hh,0.,False)
;rhh=wrf_user_intrp3d(rh,height,"h", hh,0.,False)
;-------------------------------------------------------------
spdh= sqrt(u_ph^2 + v_ph^2)
windspdh=spdh(:,x,y)
dirh=wind_direction(u_ph,v_ph,0)
winddirh=dirh(:,x,y)
;print(winddirh)
tcc=tch(:,x,y)
tdc=tdh(:,x,y)
phc=ph(:,x,y)+pbh(:,x,y)  ;units:pa
;psfcp=psfc(:,x,y)
; rhc=rhh(:,x,y)

不知道有没有好心的大佬帮忙解释一下，报错是出于什么原因呢？是对wrfout文件多久输出一次有要求吗？spdh确实是只有一维数据，x，y分别代表什么呢？ 摸索了这么久的ncl语法，对这种多维的还是没搞太清楚。

lxm666

请问老师，这个问题解决了吗？

小怪兽...

您看下原贴，xy好像是指定经纬度附近的那个格网点