ncl画wrf时间高度剖面出错

601337653

向大家请教一个问题，我想画四维（time，lev，lat，lon）wrfout数据的时间，高度剖面图，脚本如下，中间做了一些别的运算，可是ncl报错如下，看不懂是什么意思，反复看了脚本也没看出什么问题，快要疯了，大家帮帮忙！
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("e://chanchu_18km.nc","r")
;a=addfile("c://18km/18km_owa.nc","r")
u=wrf_user_getvar(f,"ua",ispan(18,54,1))
v=wrf_user_getvar(f,"va",ispan(18,54,1))
w=wrf_user_getvar(f,"wa",ispan(18,54,1))           ;垂直速度m/s
T=wrf_user_getvar(f,"tk",ispan(18,54,1))
q=wrf_user_getvar(f,"QVAPOR",ispan(18,54,1))
p=wrf_user_getvar(f,"pressure",ispan(18,54,1))     ;hpa
z=wrf_user_getvar(f,"z",ispan(18,54,1))
rh=wrf_user_getvar(f,"rh",ispan(18,54,1))          ;相对湿度
theta=wrf_user_getvar(f,"th",ispan(18,54,1))
   ;定义相关常数
cp=1004.07    ;干空气定压比热容
R=287.05      ;干空气气体常数
Rv=461.51     ;水汽气体常数
L=2501000     ;水的蒸发（凝结）潜热
dt=30*60    ;时间步长
nt=37         ;时间维
nl=39         ;高度维
nx=144        ;水平维
ny=171
dx=18000       ;格距
dy=18000
nyy=61
nxx=61
nr=30
ntt=360
rho=p/(2.87*T)
z1=new((/nt,nl,ny,nx/),float)
do i=0,nt-1
z1(i,:,:,:)=z
end do
delete(i)
tyc=(/49,49,51,51,53,54,54,54,54,54,54,53,55,57,56,55,55,57,59,63,65,67,70,73,76,80,83,87,90,94,98,103,108,114,121,126,125/)
txc=(/82,79,75,71,69,67,65,61,60,57,55,53,53,50,47,45,45,45,44,42,41,39,40,40,40,40,41,41,42,43,45,47,49,51,54,57,62/)
   ;熵
s=new((/nt,nl,ny,nx/),"float")
rh1=where(rh.eq.0,0.01,rh)
s=cp*log(theta)+L*q/T-Rv*q*log(rh1/100)-R*log(1000)
   ;熵平流项
Sx=rho*s*u
Sy=rho*s*v
Sz=rho*s*w
ds=-center_finite_diff_n(Sx,dx,False,0,2)-center_finite_diff_n(Sy,dy,False,0,1)-center_finite_diff_n(Sz,z1,False,0,0)
owa=new((/nt,nl,nyy,nxx/),"float")  
do m=0,nt-1         
do i=0,nl-1
  do j=0,nyy-1
   do k=0,nxx-1
      owa(m,i,j,k)=ds(m,i,tyc(m)-(nr-1)+j,txc(m)-(nr-1)+k)
    end do
   end do
end do
end do
owa1=dim_sum_n(owa,3)
owa2=dim_sum_n(owa1,2)
wks=gsn_open_wks("eps","aaaa")
gsn_define_colormap(wks,"BlueRed")
res=True
res@gsnMaximize=True
res@cnFillOn=True
res@cnLinesOn=False
res@cnLineLabelsOn=False
res@lbBoxLinesOn=False
res@lbOrientation="Vertical"
res@lbLabelBarOn=True
plot=gsn_csm_contour(wks,owa2,res)
end

601337653

问题解决了，是我自己粗心，差分的时候把维度搞错了。。

xibao

家园里有帖子的，自己认真看一看呢。

xq33118

xibao 发表于 2018-3-11 14:39
家园里有帖子的，自己认真看一看呢。

请问帖子在哪？

xq33118

xibao 发表于 2018-3-11 14:39
家园里有帖子的，自己认真看一看呢。

请问帖子在哪？