【求助】关于确定极涡边缘的方法

dieam

又看了一下，之前有些步骤考虑的过于简单了。

1. 计算出Epv的之后，插值到等效纬度中，每个Epv等值线包围的面积与每个纬度包围的面积是相同的

2. 在确定极涡边缘的时候，计算Epv梯度最大值，并结合纬向西风最大值区域作为约束

3. 极涡边界区域为等Epv线围绕极涡边缘的局部最大凹凸曲率的位置

                               
登录/注册后可看大图

                               
登录/注册后可看大图

                               
登录/注册后可看大图

在 Nash et al. (1996) An objective determination of the polar vortex using Ertel's potential vorticity 文章中，提到极涡边界为最大纬向西风区域中位涡经向梯度最大处所对应的等效纬度

最近在用NCL编写这个程序，遇到一些问题想和各位大佬讨论一下

1. 首先我在计算位温的时候，考虑到地形的影响只计算600hPa以上的部分，之后将纬向风以及PV插值到240K~740K的等熵面上。这个过程采用int2p_n_Wrap函数

pv     = int2p_n_Wrap(theta(level|:,time|:,lat|:,lon|:), Epv(level|:,time|:,lat|:,lon|:),\
                       pout, linlog, pdim)

不知道这个考虑是否正确呢

2. 关于最大纬向西风区域的确定【这个应该是我比较困惑的部分】

关于纬向西风最大值的话，考虑采用dim_max函数，即

u_maxlat = dim_max(u(lon|:,lat|:))

但这种方法返回的是每个经度上最大西风的纬度，即相当于急流轴，那怎样体现最大纬向西风区域呢？

在文中也有提到，“Since the wind is not single-valued around the vortex edge, we averaged the wind along Epv (Ertel's potential vorticity) isolines to determine the location of the jet maximum relative to the area under the Epv isolines."

那么，沿Epv等值线的平均风，是不是我要确定逐个Epv的值，然后和这个值对应的所有纬向风的平均值？

例如，某一根Epv等值线为5 （10^6 m^2 kg^-1 s^-1），将其他值设为缺测，采用where函数对应到纬向风数据中，留下与5 PVU对应的部分

u_5  = u
u_5  = where(pv.eq.5, u, u@_FillValue)

3. Ertel's位涡经向梯度最大位置与最大西风是接近的，那是否可以只结算最大梯度的位置来确定极涡边界呢？

以上是我个人的一些想法，也许会有很多错误。

欢迎大家一起讨论，也请计算过极涡边界的大佬赐教