MATLAB能算逐日的副高脊线位置指数吗？

Jeremy

之前用GrADs能做出逐日的副高脊线图。

如下图：

但是，我想要逐日的副高脊线位置数据。不知道能不能用MATLAB算出来，不知道各位大神有什么方法算出逐日的副高脊线位置数据。

第一次在MATLAB版发帖，还望各位大神多多指教，小弟不胜感激。

PS:发现气象家园没有求助贴这一项，不知道能不能送贡献和金钱。如果可以的话，回答能解决在下的问题的回帖，小弟会奖励888金钱以表谢意。当然，也欢迎给积极回帖讨论。

门外小爷

大神给解答

Jeremy

门外小爷 发表于 2017-1-2 15:25
大神给解答

兄台会吗？

aoyun

解决了没兄弟，我也急求

Jeremy

aoyun 发表于 2017-3-7 20:58
解决了没兄弟，我也急求

没有啊！发现越想越复杂：
① 算出的脊线需不需要考虑586高压的？

② 双脊线的情况怎么办？

liuqingyi

我最近在做相关的研究，参考了（气象家园里另外一个fortran的程序，见http://bbs.06climate.com/forum.php?mod=viewthread&tid=6555），（国家气候中心西太平洋副热带高压监测业务规定），使用JRA55数据计算了四个副热带高压指数（面积，强度，脊线和西伸脊点）。附件是我的matlab代码，可能会有缺陷，你可以参考一下，有问题我们可以再讨论。

Jeremy

liuqingyi 发表于 2017-11-13 16:19
我最近在做相关的研究，参考了（气象家园里另外一个fortran的程序，见http://bbs.06climate.com/forum.php? ...

非常感谢！

2015sysu

楼主您好~~副高脊线是季节跳动比较明显，直接算汛期一段时间内脊线的平均位置，这样是可以的吗？？还是没有意义？？
其实我是想把每年脊线都画出来，看脊线位置的变动情况。

Jeremy

2015sysu 发表于 2017-12-26 22:16
楼主您好~~副高脊线是季节跳动比较明显，直接算汛期一段时间内脊线的平均位置，这样是可以的吗？？还是没有 ...

可以，用GrADs就能画

2015sysu

Jeremy 发表于 2017-12-27 10:31
可以，用GrADs就能画

非常感谢楼主的回复！！！！
我是参考了论坛里大家共享的程序，根据需求调整成汛期这一时段的。但是输出结果有点问题，588线的西点和国家发布的西伸脊点的经度差得有点多（1993年为例，程序图显示的是140-150E，不过国家公布的74项指标是101E）
所以想请教一下，是我程序调整错了吗？不能直接用ave函数吗？？
还有一个问题是我之前看其他人程序，不太懂的地方，uwnd的作用是什么呀？似乎最终输出的结果并不需要用到uwnd呢？？

程序：
'set t 1'
'define H=ave(hgt.1,t=544,t=549)'
'define U=ave(uwnd.2,t=544,t=549)'
'd H/10.0'
'set gxout contour'
'set clevs 0'
'define ms=cdiff(U,y)'
'd maskout(U,ms)'

图：