用NCL中的rcm2rgrid函数插值不成功

lx1073760154

以下是我的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/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"
; load "/home/lixiang/files/NCLfiles/lxdef/lxdef.ncl"

begin
    ;数据读取与预处理
    ;读取再分析资料获得分析场
    f0 = addfile("/home/lixiang/files/data_download/"+\
    "meiyu0708/analysis/add_200hPa/pl07Asia_meiyu.nc", "r")
    t1 = (7-1)*24 + (8-8)  ;７日,08时（BJT）
    t2 = (8-1)*24 + (8-8)  ;8日,08时（BJT）
    u0 = short2flt(f0->u(t1:t2,:,:,:))
    v0 = short2flt(f0->v(t1:t2,:,:,:))
    u0!1 = "level"
    u0!2 = "latitude"
    u0!3 = "longitude"
    u0&level = f0&level
    u0&latitude = f0&latitude
    u0&longitude = f0&longitude
    copy_VarMeta(u0, v0)
    f = addfile("/home/lixiang/files/WRFfiles/meiyu0708/"+\
    "SKEB/test01/wrfout_d01_2020-07-06_18:00:00"+".nc", "r")
    p = wrf_user_getvar(f,"pressure",-1) ;获取总的气压数据
    u = wrf_user_getvar(f, "ua", -1)
    v = wrf_user_getvar(f, "va", -1)
    u := wrf_user_intrp3d(u, p, "h", 500, 0., False)
    v := wrf_user_intrp3d(v, p, "h", 500, 0., False)
    lat2d = wrf_user_getvar(f, "XLAT", 0)
    lon2d = wrf_user_getvar(f, "XLONG", 0)  
    ugrd = rcm2rgrid(lat2d, lon2d, u(0,:,:), u0&latitude, u0&longitude, 0)
    ; vgrd = rcm2rgrid(lat2d, lon2d, v, u0&latitude, u0&longitude, 1)
    print(ugrd(:,9))
    exit()
    fbindirwrite("u500.grd", u0)
    fbindirwrite("v500.grd", v0)
    ; printVarSummary(ugrd)
end

TangYh

请问解决了吗，我的也是这样

lx1073760154

TangYh 发表于 2020-11-30 09:57
请问解决了吗，我的也是这样

解决了，是范围的问题，要保证插值的经纬度范围包含在被插值的经纬度范围之内。