诚心求教插值问题

opeveu

请问各位大佬一个问题，是有关插值的问题，比如我想通过wrfout文件提取50m高度的风速，通过查看发现50m高度并不是在某一个eta_level上，第一个level和第二个level分别是30m和60m，然后我通过家园里的脚本提取了50m高度的风速，请问这个脚本的插值方法是哪种呢？
这是从家园里找的脚本：
begin
a = addfile("/mnt/disk3/lxg_data/xlwyongmoni/20201112RL=3/wrfout_d04_2020-12-18_00:00:00" ,"r")
times = wrf_user_list_times(a)
ntimes =dimsizes(times)-1
time=-1
u  = wrf_user_getvar(a,"ua",time)
v  = wrf_user_getvar(a,"va",time)
rh=wrf_user_getvar(a, "rh",time)
height  = wrf_user_getvar(a, "z",time)
ter = wrf_user_getvar(a, "ter",time)
Latitude = 45.40
Longitude = 127.6563
hh=100
res = True
res@returnInt = True
point = wrf_user_ll_to_ij(a,Longitude,Latitude,res)  
point =point -1
x = point(0)
y = point(1)
print(x)
print(y)
nheight = conform(height,ter,(/0,2,3/))   
height = height - nheight
u_ph  = wrf_user_intrp3d( u,height,"h", hh,0.,False)  
v_ph  = wrf_user_intrp3d( v,height,"h", hh,0.,False)
rhh=wrf_user_intrp3d(rh,height,"h", hh,0.,False)
spdh= sqrt(u_ph^2 + v_ph^2)
windspdh=spdh(:,x,y)
dirh=wind_direction(u_ph,v_ph,0)
winddirh=dirh(:,x,y)
rhc=rhh(:,x,y)
npts=ntimes+1
fName = "2020-12-18_100m_heightdata.txt"
data  = new( npts, "string")
print("  Time     Wind_dir_hm      Wind_speed_hm ")
print ( times+sprintf("%12.1f",rhc)+sprintf("%12.0f",winddirh)+sprintf("%12.1f", windspdh) )

data = times + sprintf("%12.1f",rhc)+sprintf("%12.0f",winddirh)+sprintf("%12.1f", windspdh)
asciiwrite (fName,data)

end

opeveu

本人用的初始及边界条件数据为FNL数据

opeveu

也就是ncl在处理wrfout文件时，垂直方向的插值法是样条插值法还是分段插值法或者是别的？

Liujiaxin

顶一波，坐等大佬解释！

飘散的梦

隐约记得是样条

chaoschaos

dddddd字数补丁

dashuku

默认是线性插值吧？

opeveu

dashuku 发表于 2022-5-18 14:14
默认是线性插值吧？

垂直方向是样条插值，水平方向是线性插值