fatal:Subscript out of range, error in subscript #0

追风的阳光

附上相关脚本：;*******p 垂直插值*********************  
  minlat   =36.2
  maxlat   =35
  minlon   =119.2
  maxlon   =119.2
  loc =wrf_user_ll_to_ij(f, (/minlon,maxlon/),(/minlat,maxlat/), True)

  plane=(/loc(0,0),loc(1,0),loc(0,1),loc(1,1)/)
  opts          =True
  angle         =0.
  p_plane   = wrf_user_intrp3d(p, p,"v",plane, angle,opts)

；******************水平插值得到计算假相当位温的变量（p_plane从2开始是因为0,1层为缺测层）******************
pressure_levels = (/p_plane(2:, 0) /)   ; pressure levels to plot
  nlevels         = dimsizes(pressure_levels)     ; number of pressure levels
  tc_plane = wrf_user_intrp3d(tc,p,"h",pressure_levels,0.,False)
  rh_plane = wrf_user_intrp3d(rh,p,"h",pressure_levels,0.,False)

；*******************计算假相当位温******************************
tk=tc_plane+273.166
  P = conform(tc_plane,pressure_levels ,0)
  a1 = where(tk .gt. 263.0, 0.622*6.11*exp(17.26*(tk-273.16)/(tk-35.86)), \
                     0.622*6.11*exp(21.87*(tk-273.16)/(tk-7.66)))   ;
  b1= where(tk .gt. 263.0, P-0.278*exp(17.26*(tk-273.16)/(tk-35.86)),\
                 P-0.278*exp(21.87*(tk-273.16)/(tk-7.66)))   ;
  qs1=a1/b1
  q1=qs1*rh_plane
  e1=P*q1/100./(0.62197+q1/100.0)

  tk1=55.0+2840.0/(3.5*log(tk)-log(e1)-4.805)
  pot1=tk*(1000/P)^(0.2854*(1.0-0.28*q1/100.0))
  ept1=pot1*exp(((3376./tk1)-2.54)*q1/100.0*(1.0+0.81*q1/100.0)) ;
  ept1@description = "0se"
  ept1@units = "K"

；*********************赋坐标属性****************************
  copy_VarMeta(avo,ept1)
；***********************垂直俩点插值**************************
ept_plane  = wrf_user_intrp3d(ept1,p,"v",plane,angle,opts)

在垂直插值语句上出现了上述错误，请问是为什么啊。。。。谢谢！


Variable: ept1
Type: float
Total Size: 43194368 bytes
            10798592 values
Number of Dimensions: 3
Dimensions and sizes:        [bottom_top | 92] x [south_north | 448] x [west_east | 262]
Coordinates:
Number Of Attributes: 7
  FieldType :        104
  MemoryOrder :        XY
  stagger :       
  coordinates :        XLONG XLAT
  units :        10-5 s-1
  description :        Absolute Vorticity
  _FillValue :        9.96921e+36