好难过啊 rcm2grid 插值出错

丫土豆wpy

我想把用rcm2grid把 wrf 输出结果插值到格点上，但是插值结果只有一个格点有值,其他都是缺省。不知道为什么啊
求助
t=rcm2rgrid(lat_wrf,lon_wrf,wrf_out,lat,lon,0);

weiguang1233

脚本贴多一点看看，这样什么都看不到

丫土豆wpy

weiguang1233 发表于 2016-10-26 16:59
脚本贴多一点看看，这样什么都看不到

begin
;---Set some needed arrays                           
     
  lines1= asciiread("H500_summer_era.txt",(/480,241/),"float");
  
  lat= asciiread("lat_era.txt",(/241/),"float");
  lon= asciiread("lon_era.txt",(/480/),"float");
  lat_wrf=asciiread("lat.txt",(/232,196/),"float");
  lon_wrf=asciiread("lon.txt",(/232,196/),"float");
  
  diri= "./"
  fils1=systemfunc("ls "+diri+"*H500*exp*.txt")
  
  do j=0,0
  
  lines4= asciiread(fils1(j),(/232,196/),"float"); spatial data

  lines4@_FillValue=-999
  lines4=lines4*1000;
  
  
  t=rcm2rgrid(lat_wrf,lon_wrf,lines4,lat,lon,0);

  
  
  end do
  


  
end

丫土豆wpy

weiguang1233 发表于 2016-10-26 16:59
脚本贴多一点看看，这样什么都看不到

求助求助 找了两天了

丫土豆wpy

weiguang1233 发表于 2016-10-26 16:59
脚本贴多一点看看，这样什么都看不到

求助求助 找了两天了

weiguang1233

丫土豆wpy 发表于 2016-10-26 17:24
求助求助 找了两天了

没看出来，如果排除了你保存数据有问题的情况，可以一步步print，试一下哪一步出了问题。有一点不懂，为什么你要把wrfout保存在*.txt里。而不用ncl官网的例子呢

1. Example 1
   
2. 
   
3. Interpolate to a gaussian T42 (64,128) grid.
   
4.    lat = latGau(64)       ; create gaussian latitudes or read from file
   
5.    lon = fspan(0,359,128) ; create longitudes or read from file
   
6. 
   
7.    f     = addfile ("some_RCM_file.nc", "r")
   
8.    lat2d = f->xlat                         ; size = (nlat,nlon)
   
9.    lon2d = f->xlon                         ; size = (nlat,nlon)
   
10.    x     = f->X
    
11. 
    
12.    xgrd  = rcm2rgrid_Wrap(lat2d,lon2d,x,lat,lon,0)
    
13. 
    
14. 
    
15. If x was of size:
    
16.             x(nlat,nlon) ==>  xgrd(64,128)
    
17.        x(ntim,nlat,nlon) ==>  xgrd(ntim,64,128)
    
18.   x(ntim,klev,nlat,nlon) ==>  xgrd(ntim,klev,64,128)
    
19. 
    

复制代码

我对wrfout，通过rcm2rgrid处理的方法跟ncl官网差不多。坐等大神们回答。

丫土豆wpy

weiguang1233 发表于 2016-10-26 18:38
没看出来，如果排除了你保存数据有问题的情况，可以一步步print，试一下哪一步出了问题。有一点不懂，为 ...

那个 我是先在matlab 里面 先计算了一遍 我在检查检查

Will

我也是被这个问题困扰了很久，建议楼主使用ESMF_regrid_with_weights这个命令，可能有点复杂，多看看官网吧。

Royone

{:5_214:}

chenchaohui2001

这个函数比较慢