初学者分享下rgrid2rcm的学习心得

小虹

rgrid2rcm函数是插值函数，所用的插值方法为反距离权重法，可用于规则网格数据往曲线网格数据（如WRF数据）插值。
附上函数说明网址http://www.ncl.ucar.edu/Document/Functions/Built-in/rgrid2rcm.shtml
下面是我将1*1的资料插值到wrfinput的网格中的脚本：

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
begin
;********************************************************
; open file
;********************************************************
  f = addfile ("wrfinput_d01.nc", "r")

  xlon   = f->XLONG(0,:,:)              ;wrfinput的经纬度，由于是兰伯特投影，因此不是均匀网格，经纬度是2维的
  xlat   = f->XLAT(0,:,:)
  
  ff = addfile ("20080127.nc", "r")  ;1*1的数据

lat_0=ff->lat_0(::-1)    ;1-D lat(181)  lon(360) 由于ff文件中纬度是起于90，止于-90，而rgrid2rcm要求经纬度必须是递增，因此反向读入
lon_0=ff->lon_0-180  ;;; ff文件中lon(0)=180,因此这里减去180才与实际情况相对应

en_tem  = ff->t_P1_L100_GLL0(0,1,7,::-1,:)     ; (ensemble0, initial_time0_hours, lv_ISBL0, lat_0, lon_0)数据lat反向读入，与上对应

  en_tem!0="lat_0"      ;赋予数据属性
  en_tem!1="lon_0"
  en_tem&lat_0=lat_0(:)
  en_tem&lon_0=lon_0(:)
  wrf_tem=rgrid2rcm(en_tem&lat_0,en_tem&lon_0,en_tem,xlat,xlon,0)

  print(wrf_tem)

end

kongfeng0824

谢谢分享                     

天边的那边

谢谢分享

卿暮霏

谢谢分享

青山迢迢

感谢楼主分享~

金玉姻缘

楼主你好，我有一点看不懂，请教一下，你写的“  f = addfile ("wrfinput_d01.nc", "r")，xlon   = f->XLONG(0,:,:) ，xlat   = f->XLAT(0,:,:);wrfinput的经纬度，由于是兰伯特投影，因此不是均匀网格，经纬度是2维的 ”，不是很理解这个经纬度为什么是二维的，看了很多nc数据，它经纬度一般都是一维变化的啊，为什么wrf的数据是二维的？还望解释一下，谢谢啦，还有一点就是变量名为什么是XLAT，XLON，给人感觉是它是径向或者纬向网格点的数目，不知道理解的对不对

小虹

金玉姻缘 发表于 2014-11-30 18:05
楼主你好，我有一点看不懂，请教一下，你写的“  f = addfile ("wrfinput_d01.nc", "r")，xlon   = f->XLON ...

你说的经纬度是一维的情况，指的是经纬度随格点呈线性变化的，所以可以用一维数组来表示。wrf的格点是如果采用兰伯特投影，格点与格点之间的距离不是均匀的，你可以找个wrfinput或者wrfout数据读一下看看，画个图就明白了

金玉姻缘

小虹 发表于 2014-11-30 22:17
你说的经纬度是一维的情况，指的是经纬度随格点呈线性变化的，所以可以用一维数组来表示。wrf的格点是如 ...

@小虹 我用rcm2rgrid函数把regcm模式运行的数据插值到自己设定的网格上了，这个和你说的那个二维经纬度坐标一样，现在加深了对这个的理解，谢啦啊

2511624023

谢谢分享，还是很用用的

追风的阳光

小虹 发表于 2014-11-30 22:17
你说的经纬度是一维的情况，指的是经纬度随格点呈线性变化的，所以可以用一维数组来表示。wrf的格点是如 ...

lz你好，我最近也在使用rcm2grid等函数，我想将wrf网格投影到经纬网格上，以便利用再分析数据，可是脚本运行后出现了问题（插值后的值全部都是缺测值），下附脚本：
......................
in=addfile("/home/zssapr/wrfout_d01_2014-07-24-00_f02.nc","r")
;************get variables***********************************
   UA = wrf_user_getvar(in, "ua",-1)
  VA = wrf_user_getvar(in, "va",-1)
  p = wrf_user_getvar(in, "pressure",-1)
  lat2d = in->XLAT(0,:,:)     ;经纬网格为262*448
  lon2d = in->XLONG(0,:,:)  ;最低纬经度范围为112.87至122.7289
  ;print(lon2d(0,:))
  lat2d@units = "degrees_north"
  lon2d@units = "degrees_east"
  ua2 = wrf_user_intrp3d(UA,p,"h",900,0.,False) ;插值到900hPa
  va2 = wrf_user_intrp3d(VA,p,"h",900,0.,False)
  T =in->Times
  scale=1.e05
  wks = gsn_open_wks("ps","flow")
;*********************wrf grid to  rectilinear grid**********
   nlat=448
   lat  = latGau(nlat, "lat", "latitude", "degrees_north")
   lon = fspan(115,120,262) ; create longitudes or read from file
   UA2  = rcm2rgrid_Wrap(lat2d,lon2d,ua2,lat,lon,1)  ;插值到经纬网格上
   VA2  = rcm2rgrid_Wrap(lat2d,lon2d,va2,lat,lon,1)
   print(UA2)
............................
我print（UA2）后，发现输出的全部都是缺测值（9.96921e+36），是怎么一回事，你能提供些建议吗？