用ncl将wrfout输出的风场插值到等熵面上

清晨迪声

各位大神，我最近用ncl的int2p_n函数将wrfout的风场从sigma面插值到等熵面后，出现了很多缺测值，求支招，谢谢！
以下是我的程序中用来插值的部分：

  u = wrf_user_getvar(f,"ua",0)
  the = wrf_user_getvar(f, "theta",0)

  do k = 0,33
    lvl(k) = avg(the(k,:,:))
   end do         
       
   lvl@units = the@units
   uu = int2p_n(the,u,lvl,1,0)

其中：
u和the的Dimensions and sizes:        [bottom_top | 34] x [south_north | 600] x [west_east | 600]

coliane

同等熵面插值出现缺测 不知道怎么办

coliane

我找到原因了，是因为int2p这个函数要求pin这个变量必须是单调，这就要求输入的位温变量不能有缺测，不然会被认为不单调，而受地形影响，很多地区在低层都有缺测，如青藏高原600hPa以下就是缺测了。所以int2p这个函数做不了插值，才会输出缺测。解决办法就是，我们看位涡一般都只看高层的，在输入位温变量时可以把受地形影响而缺测的部分给截掉，再放入int2p函数进行计算就可以了。

清晨迪声

coliane 发表于 2022-2-25 14:27
我找到原因了，是因为int2p这个函数要求pin这个变量必须是单调，这就要求输入的位温变量不能有缺测，不然会 ...

原来如此，谢谢你~