NCL中条件语句产生缺失值，NCL无法确定分支。。。

zhaolinanshiwo

如题有大神知道吗？附上运行的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/wrf/WRFUserARW.ncl"

begin

a = addfile("./wrfout_d03_2013-04-05_12:00:00.nc","r")

wks =gsn_open_wks("png","plt_HeightLevel")

res = True

res@MainTitle = "REAL-TIME WRF"

res@Footer = False

pltres = True

mpres = True

times =wrf_user_getvar(a,"times",-1)

ntimes = dimsizes(times)

height_levels = (/250.,2000./)

nlevels = dimsizes(height_levels)

do it = 0,ntimes-1,2

       print("Workingon time:" + times(it))

       res@TimeLabel= times(it)

       tc= wrf_user_getvar(a,"tc",it)

       u  = wrf_user_getvar(a,"ua",it)

       v  = wrf_user_getvar(a,"va",it)

       p  = wrf_user_getvar(a,"pressure",it)

       z  = wrf_user_getvar(a,"z",it)

       rh= wrf_user_getvar(a,"rh",it)

       dolevel = 0,nlevels-1

              height= height_levels(level)

              

              p_plane  = wrf_user_intrp3d(p,z,"h",height,0.,False)

              tc_plane= wrf_user_intrp3d(tc,z,"h",height,0.,False)

              rh_plane= wrf_user_intrp3d(rh,z,"h",height,0.,False)

              u_plane  = wrf_user_intrp3d( u,z,"h",height,0.,False)

              v_plane  = wrf_user_intrp3d(v,z,"h",height,0.,False)

              u_plane= u_plane*1.94386

              v_plane= v_plane*1.94386

              u_plane@units= "kts"

              v_plane@units= "kts"

              opts= res

              opts@cnLineColor= "Blue"

              opts@gsnContourLineThicknessesScale= 3.0

              contour_tc= wrf_contour(a,wks,tc_plane,opts)

              delete(opts)

              opts= res

              opts@cnLineColor= "Blue"

              opts@gsnContourLineThicknessesScale= 3.0

              contour_p= wrf_contour(a,wks,p_plane,opts)

              delete(opts)

              opts= res

              opts@cnFillOn= True

              opts@ContourParameters= (/10.,90.,10./)

              opts@cnFillColors=(/"White","White","White","White","Chartreuse","Green","Green3","Green4"/)

              contour_rh= wrf_contour(a,wks,rh_plane,opts)

              delete(opts)

              opts= res

              opts@FieldTitle= "Wind"

              opts@NumVectors= 47

              vector= wrf_vector(a,wks,u_plane,v_plane,opts)

              delete(opts)

              plot=wrf_map_overlays(a,wks,(/contour_rh,contour_tc,contour_p,vector/),pltres,mpres)

              

       enddo

end do

end

zhaolinanshiwo

自顶

禾栗儿

遇到了同样的问题，楼主解决了吗