【已解决】位涡插值到等熵面值缺测

是冉冉升起的冉

位涡计算使用pot_vort_isobaric这个函数，插值使用的int2p_n_Wrap，位涡的值也输出来看了都正常，开始以为是0值的影响，0值替换成缺测值后，缺测值使用poisson_grid_fill插值了，位涡值用ismissing函数看过没有缺测，但是插值到等熵面就全部缺测，论坛上以前有人说是位温低层缺测，我选择的层也取高了，也检查了位温的结果。
不知道哪里出了问题！有没有人遇见过这种情况！！
----------------------------------------------------------------------------------------------------------------------------
计算位涡的变量减去了daily climatology想要看异常，但之后插值就会显示缺测。
不看异常就不会缺测！！！所以选择插值完再减去气候态计算的插值，虽然不知道方法对不对。

检查数据，常用的两个代码。

;-----将0设置为缺失值

    u_ave_10_pos=where(u_ave_10_pos.ne.0,u_ave_10_pos,u_ave_10_pos@_FillValue)
;-----------------------缺测值 插值
    nscan = 1500             ; usually *much* fewer
    eps   = 1.e-2            ; variable depended
    gtype = True            ; "lon" does  generate cyclic fields
    guess = 1                ; use zonal means
    relc  = 0.6              ; standard relaxation coef
    opt   := 0
    poisson_grid_fill(u_ave_10_pos, gtype, guess, nscan, eps, relc, opt)