wrf输出结果怎么用rcm2rgrid插值？求指导！

hhumuer

打算用rcm2rgrid 对wrf输出降雨进行插值，写了如下脚本（it 是时间步长），但是，提示错误，求指导。或者哪位高手可以把之前的例子贴一个出来不～～～～

脚本代码：
    rain_exp = wrf_user_getvar(a,"RAINNC",it)
    rain_con = wrf_user_getvar(a,"RAINC",it)


; rcm2rgrid  **********************************************
        nlat=48
        nlon=64
        lat=latGau(nlat,"lat","latitude", "degrees_north")
        lon=fspan(0,nlon,1)
        lat2d=a->xlat
        lon2d=a->xlon

        rain_exp_xgrd =rcm2rgrid(lat2d,lon2d,rain_exp,lat,lon,0)
        rain_con_xgrd =rcm2rgrid(lat2d,lon2d,rain_con,lat,lon,0)

        rain_exp(it, nlat, nlon)==>rain_exp_xgrd(it,64,128)
        rain_con(it, nlat, nlon)==>rain_con_xgrd(it,64,128)
        rain_tot = rain_exp + rain_con
;*******************************************************************

错误提示：
fatal:syntax error: line 80 in file 2gridPrecipitation.ncl before or near =>rain_exp_xgrd
        rain_exp(it, nlat, nlon)==>rain_exp_xgrd
-----------------------------------------------^

fatal:syntax error: line 81 in file 2gridPrecipitation.ncl before or near =>rain_con_xgrd
        rain_con(it, nlat, nlon)==>rain_con_xgrd
-----------------------------------------------^

fatal:Syntax Error in block, block not executed
fatal:error at line 177 in file 2gridPrecipitation.ncl

zpxie

是你函数参数的问题，xgrid=rgrid2rcm(lat, lon, rain, lat2d, lon2d, 1) 刚好你lat和lat2d位置反了 建议使用韩式之前认真仔细看看函数说明

hhumuer

zpxie 发表于 2014-11-14 22:29
是你函数参数的问题，xgrid=rgrid2rcm(lat, lon, rain, lat2d, lon2d, 1) 刚好你lat和lat2d位置反了 建议使 ...

不是噢，我用的是rcm2rgrid,你说的是rgrid2rcm,其实问题是，不管你用哪一个函数，里面的lat lon都要有合理的属性

hxyj

请问楼主，该问题解决了吗？如果解决了，能否告知一下如何解决的？谢谢

sjx8868

rain_exp_xgrd =rcm2rgrid(lat2d,lon2d,rain_exp,lat,lon,0)

天生我才

学习学习

cibaouc

        lat=latGau(nlat,"lat","latitude", "degrees_north")，不用全球模式，用这个貌似不对
        lon=fspan(0,nlon,1)

planck

cibaouc 发表于 2020-4-28 20:09
lat=latGau(nlat,"lat","latitude", "degrees_north")，不用全球模式，用这个貌似不对
        lo ...

官网个例是这样写的，介绍说是适用wrf

yq奇

planck 发表于 2020-5-20 16:48
官网个例是这样写的，介绍说是适用wrf

官网网址可以贴一下吗？十分感谢！