- 积分
- 3167
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-10-5
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 oyy. 于 2020-3-5 16:04 编辑
楼主最近在研究T-N波作用通量,费了很久才理清这个玩意儿咋算, 今天有空来整理一下程序。感谢气象家园@伽蓝鸟 大佬 以及同学cj xc zxy的帮助。
先列一下在论坛里找到的资料,供大家参考。
Python计算T-N波作用通量
三维波作用通量ncl脚本分享
NCL 画波作用通量
Programs
下图是我参考的公式
公式参考文献:
1. (球坐标系) 施春华,金鑫,刘仁强,2017.大气动力学中三种 Rossby 波作用通量的特征差异和适用性比较[J] .大气科学学报,40( 6) : 850-855.
2. (直角坐标系)Yang SY, Li T (2016) Intraseasonal variability of air temperature over the mid-high latitude Eurasia in boreal winter. Clim Dyn 47:2155–2175
3. (直角坐标系,主要参考) Gao M-N, Yang J, Wang B, Zhou S-Y, Gong D-Y, Kim SJ(2017) How are heat waves over Yangtze River valley associated with atmosphericquasi-biweekly oscillation? Climate Dynamics:1-17
不同的时间尺度下 对数据的处理也不同. 楼主主要研究 夏季(5-8月)季节内时间尺度(10-30天) 中高纬波列的波作用通量, 程序求的是三维波作用通量的水平分量(垂直方向不考虑)。
流函数是用u v求得,有的程序中是用位势高度hgt来表示流函数。大家可以都试试, 楼主偷懒只算了前一种。
另外发现的一点是,ncl函数uv2sfvpF求得流函数存在问题,跟环流场并不能很好的对应。可以参考基友xc发现的这个网址辐散风和旋转风,这个网页还有很多干货值得学习。
公式中 主要是流函数pthi的扰动场 以及u v的背景态。扰动场(u v先10-30天滤波,然后回归、合成等求得u’, v’。最后求得pthi’)。背景态是相对于季节内时间尺度而言(u v求30天以上低通滤波,然后对时间维求平均)。
如有错漏 欢迎指正。
|
-
-
-
TN.ncl
4.21 KB, 下载次数: 580, 下载积分: 金钱 -5
|