重磅！副高脊线程序MATLAB版

斥鷃

附件为编写锝比较粗糙的副热带高压动力学脊线的程序。算法是从45°N开始往南搜索找到第一个出现西风转换为东风的位置。函数识别纬度位置的精度完全取决于输入数据精度，且所有输入数据精度如果需要提高识别精度需要进一步改写程序，这里只作为抛砖引玉。注意讨论副高脊线位置，由于其位置随时间和区域变化非常大，所以程序还设置了用于识别脊线经纬度范围的参数，但即使如此仍需小心检验脊线位置识别的准确与否，附上的程序仅为参考。
函数中投影矩阵R可以使用R=makerefmat(起始经度,起始纬度,经度分辨率,纬度分辨率);来生成，例如起始经度为60°E，起始纬度维20°S，分辨率为1.25°的输入资料，投影矩阵R=makerefmat(60,-20,1.25,1.25);更多有关投影矩阵的细节可以参考我的另一个帖子《经纬度和矩阵行列标的转换（基于MAP工具箱）及一个副高西伸脊点程序》。
以下为函数的help：

1. 提取脊线纬度位置以及其上某物理量的程序，函数形式：
   
2. [position output]=westernpts(u,R,lonrange,latrange,inputvar)
   
3. 输入变量分别为： u风速、投影矩阵R、经度范围、纬度范围、脊线上所需物理量，
   
4. 输入变量均需要是经纬度存储的二维变量，行标（第一维）为纬度，列标（第二维）为经度，可以允许缺测值
   
5. 输出变量position为每个经度上的脊线纬度值，output为每个经度上所需物理量的值
   
6. 举例：[pts ridhgt]=ridgeline(uwnd,R,[60 160],[10 45],vwnd);

复制代码

Copy_of_ridgeline.m (1.51 KB, 下载次数: 64)

2017-12-18 20:11 上传

点击文件名下载附件
副高脊线程序 再上一张高大上的三维图，图为7月份平均的副热带高压及南亚高压脊面及脊面上经向风风速的分布，其中红黄色系为南风，绿蓝色系为北风，黑灰背景为地形。

8月份平均的副热带高压及南亚高压脊面及脊面上经向风风速的分布，其中红黄色系为南风，绿蓝色系为北风，黑 ...

扫眉才子

感谢分享，非常6666

看云淡风清

楼主厉害、膜拜下

风之精灵

厉害喽，膜拜大神

mapleleaves99

对此颇有兴趣，学习一下

micolian

楼主厉害，膜拜一下

Annieee

求问大神，判断脊线时用的风速是500hPa层上的吗

kc121

这图画的还不错！学习一下。