读取多个文件之后提示的错误

小诸葛

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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
  files = systemfunc ("ls ../../wrfout_d01*")
  f     = addfiles (files + ".nc" , "r")   ；多个文件读取
  f     = addfile("../../*.nc","r")  ；单个文件读取
  wks = gsn_open_wks(“x11”,"wrfout_wind")
  res = True
  res@MainTitle       = " WRF RESULT--WIND "
  res@tiXAxisString   = " Latitude "
  res@tiYAxisString   = " Longitude "
   pltres = True
  mpres  = True
   times  = wrf_user_list_times(f)  ; get times in the file
   ntimes = dimsizes(times)         ; number of times in the file
   do it = 0,ntimes-1
      u  = wrf_user_getvar(f,"ua",it)      ; u averaged to mass points
      v  = wrf_user_getvar(f,"va",it)      ; v averaged to mass points
   end do
     uu  = dim_avg_n(u,0)
     vv  = dim_avg_n(v,0)  
     spd = sqrt(uu*uu + vv*vv)           ; speed in m/sec
     wrf_smooth_2d( spd, 3)              ; smooth speed   
     spd@description = "Wind Speed"
     spd@units = "m/s"
       opts = res
       opts@ContourParameters = (/ 1., 12., 1. /)
       opts@cnFillOn          = True
       opts@cnLevelSelectionMode   ="ManualLevels"         
       gsn_define_colormap(wks,"BlAqGrYeOrReVi200")
       opts@cnLabelMasking    = True
       contour = wrf_contour(f,wks,spd(:,:),opts)
       delete(opts)   
end

    当读取单个文件的时候没有问题，但是当读取多个文件的时候会提示：fatal:Argument type mismatch on argument (0) of (wrf_contour) can not coerce，，错误在contour = wrf_contour(f,wks,spd(:,:),opts)这一行，查询了下解决办法：http://www.ncl.ucar.edu/Support/talk_archives/2011/2451.html，建议将contour = wrf_contour(f,wks,spd(:,:),opts)改为contour = wrf_contour(f(0),wks,spd(:,:),opts)，但是还是提示同样的错误，大家有谁遇到过同样的问题或者哪位达人帮忙解决下，不胜感激！

小诸葛

mofangbao 发表于 2012-3-6 09:40

                               
登录/注册后可看大图

虽然没有人能够解答你的问题，但是我仍然期待楼主能够把解决方法描述一下，这样就不仅仅是你解决了问题， ...

就是将contour = wrf_contour(f,wks,spd(:,:),opts)改为contour = wrf_contour(f[0],wks,spd(:,:),opts)，不是小括号而是中括号

小诸葛

问题已解决

mofangbao

小诸葛 发表于 2012-3-6 09:18

                               
登录/注册后可看大图

问题已解决

虽然没有人能够解答你的问题，但是我仍然期待楼主能够把解决方法描述一下，这样就不仅仅是你解决了问题，还能帮助以后的很多朋友，先谢过了！

gaoyc

今天学习ncl，也遇到了同样的问题。看到你的解决方案，很快解决了我的问题，非常感谢！

liangsheng3460

不错  解决了问题

晓猪王木木

解决了问题，但是不知道原理是啥？

咩大人

遇到了同样的问题  感谢楼主贴出解决方法
希望以后大家解决问题后能把解决方法也贴上来，方便大家

rjys1219

好，学习中...

苏绮

谢谢楼主！问题解决。我的理解是批量读文件以后，f就多了一维，所以f[0]才代表文件，不知道对的吧