NCL双线性插值wrfout数据

最爱milk的味道

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"
;************************************************
begin
  diri = "/cygdrive/H/try/qw/"   

  a=addfile("wrfout_d04_2008-04-19_12_00_00.nc","r")              
  
  u=a->U10(:,:,:)  
  lat=a->XLAT
  lon=a->XLONG
  ;printVarSummary(u)  
  lat@units = "degrees_north"
  lon@units = "degrees_east"      
  ;u@missing_value=  missing
  ;u@_FillValue=  missing
     
LAT=fspan(31.625,34.875,14)      
LON=fspan(119.,122.,13)
;datarw2=new((/14,13/),"float")
datarw2=linint2_Wrap(lon,lat,u,True,LON,LAT,0)  
;datarw2=linint2_Wrap(west_east,south_north,u,True,LON,LAT,0)
  datarw2@missing_value=  missing
  datarw2@_FillValue=  missing
  datarw2!0   ="time"  
  datarw2!1   ="lat"
  datarw2!2   = "lon"
  datarw2&lat =LAT
  datarw2&lon =LON
  LAT@units = "degrees_north"
  LON@units = "degrees_east"      
  print(datarw2@missing_value)

outname="qs_u-0419"  
outfile ="/cygdrive/h/try/qw/"+outname+".nc"     
  system("/bin/rm -rf "+outfile)
  out= addfile(outfile,"c")  
  out->U10=datarw2
end
出错：fatal:Number of dimensions in parameter (0) of (linint2_Wrap) is (3), (1) dimensions were expected
fatal:["Execute.c":8578]:Execute: Error occurred at or near line 25 in file /cygdrive/i/try/qw/chazhi.ncl
请大神赐教！

weiguang1233

看提示是25行什么东西写的有问题呀

最爱milk的味道

weiguang1233 发表于 2015-9-24 08:35
看提示是25行什么东西写的有问题呀

改来改去还是错的，不知道怎么改

zhouxiaosky

多谢分享~~~

回_回

说你的lon的维数有问题啊，目前看不出什么 printVarSummary(lon) 看看维数对不对

最爱milk的味道

回_回 发表于 2015-10-1 22:08
说你的lon的维数有问题啊，目前看不出什么 printVarSummary(lon) 看看维数对不对

printVarSummary(lon)之后不太能看懂，能否加个QQ请教一下

罗夏

多谢，收下了

abc258886@

那个插值方法更好？？？

heohd

感谢你的分享，对我有很大帮助

小莹子

最爱milk的味道 发表于 2015-10-15 08:38
printVarSummary(lon)之后不太能看懂，能否加个QQ请教一下

你好请问你学下了吗？想请教下LAT=fspan(31.625,34.875,14)      
LON=fspan(119.,122.,13)，这两行怎么理解那？