数据范围不一致（ncl）

追风的阳光

后面是我添加的语句，谢谢！   
。。。。。。。。
    tk=T
     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
     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"
     levels      = ispan(200,350,10)
     rh_res                                = res
     rh_res@cnLineColor                    = "black"
     rh_res@cnLevelSelectionMode           = "ExplicitLevels"
     rh_res@cnLevels                       = levels
     rh_res@cnLineLabelBackgroundColor     = -1    ; transparent
     rh_res@cnLineThicknessF               = 3.5
    rh_res@cnLineLabelsOn                 =True
     ;rh_res@cnHighLabelsOn                 = True
     ;rh_res@cnLowLabelsOn                  = True
     rh_res@cnHighLabelBackgroundColor     = -1
     rh_res@cnLowLabelBackgroundColor      = -1
     rh_res@cnLineDrawOrder                = "draw";"Postdraw"
     rh_res@cnInfoLabelOn                  =False
     ;print(T)
     printVarSummary(ept1)

     contour_rh = gsn_csm_contour(wks,ept1,rh_res)
。。。。。。。。。。。。