求大侠解惑：如何利用WRF输出速度变量U、V求和速度以及风向？

kingkong

xyan88 发表于 2016-9-25 10:38
学习了 请问最后是怎么样插值的啊

我是用了matlab的函数scatteredInterpolant，你可以查查它的用法，另外，ncl也能处理

hegc100

使用NCL的函数：
十米高度直接有
u10 = a->U10
v10 = a->V10
其余层高度的风速：
    u  = wrf_user_getvar(a,"ua",-1)        ; u to mass points
    v  = wrf_user_getvar(a,"va",-1)        ; v to mass points
插值函数，可以插值获得非eta层的：
  u_plane  = wrf_user_intrp3d( u,z,"h",height_levels,0.,False)
  v_plane  = wrf_user_intrp3d( v,z,"h",height_levels,0.,False)
由测站经纬度，取得最近计算格点函数：
loc = wrf_user_ll_to_ij(pointa, 105.1874, 23.478, True);
ny = loc(1) - 1
mx = loc(0) - 1
由风速矢量U,V求风向函数：
wrfdirt=wind_direction(U,V)