数值模式 | 如何对WRF变量场的地形缺测进行填补（二）

 "Smarter, albeit slower version of wrf.interplevel()"

 * pressure（等压面）
 * ght_msl / ght_agl（等高度面）
 * theta / theta-e（等位温面）

 from wrf import vinterp
 interp_field = vinterp(
     wrfin,                    # WRF NetCDF 文件对象
     field=tk,                 # 待插值的 3D 场
     vert_coord="pres",        # 目标垂直坐标：等压面
     interp_levels=[850, 500], # 目标层（hPa）
     extrapolate=False,        # 是否在地形以下外推
     field_type="tk",          # 变量类型标识
     log_p=True                # 对数气压插值
 )

 1. 它只在垂直方向外推，不是水平方向的空间插值
 2. 有条件限制：只有当 vert_coord 和 field_type 匹配特定组合时才生效
 3. 如果条件不满足或 field_type=None，则直接使用最低模式层值，不做任何外推

 ! 使用美国标准大气递减率 (US Standard Atmosphere lapse rate)
 tlev = tlhsl + (zlhsl - zlev)*USSALR   ! 温度外推

 ! 对压力坐标：
 zlev = zbot + vt/USSALR*(1. - (vlev/pbot)**EXPON)

 ! 对高度坐标：
 plev = pbot*(1. + USSALR/vt*(zbot - zlev))**EXPONI

 * 利用最低模式层温度、气压、湿度和标准大气递减率
 * 在垂直方向上估算地形以下的物理量
 * 本质上是1D的物理外推，不是 2D 的空间插值

from wrf import getvar, vinterp
t850_1 = vinterp(wrf_file, t, 'pressure', [850],extrapolate=True )