关于EP通量请教一下

影伊

（已解决）关于二维E-P通量的计算及应用[Fortran程序]这个问题有些拖延了。。。

说明，求的是行星波1~3波，对u，v，t作了谐波分析。
1991-2014EPD.gmf (183.65 KB, 下载次数: 65)

2015-6-4 20:18 上传

点击文件名下载附件


首先，EP通量的公式，文献参考(Edmon et al.,1980;Andrews et al.,1987)

E-P通量

E-P通量散度

其中，需要注意的是

F（p）分母项

在绘图时，
散度项处理成： ，并且数值上乘以24*3600，则单位转换成m/s/d
E-P通量的两项处理成：

如果同时画平流层和对流层波作用通量的传播时，还要再做下面的处理。
（2）
(http://www.esrl.noaa.gov/psd/data/epflux/)


但是在处理强的天气事件时，如平流层爆发性增温（SSW），可以不用（2）处理。


文献参考（球面E-P通量的计算及应用）


现在新的问题想请教，grads画图，这个箭头的设置如何调整，才能美观呢

liding

为什么文献中通量量级达到十的八次方？

影伊

这个问题我已经解决。在搜索EP通量的过程中看到了如下几篇相关的帖。
我主要参考了Fortran ‘E-P通量的Fortran程序
http://bbs.06climate.com/forum.php?mod=viewthread&tid=17342&fromuid=33140
(出处: 气象家园)’

另外还有其他绘图软件的程序参考如（a）ncl ：‘E-P通量ncl计算绘图脚本
http://bbs.06climate.com/forum.php?mod=viewthread&tid=19309&fromuid=33140
(出处: 气象家园)’

                                                     （b）matlab ‘问老外要来的EP通量程序
http://bbs.06climate.com/forum.php?mod=viewthread&tid=17788&fromuid=33140
(出处: 气象家园)’

还有这篇帖子里也提到了EP通量。
‘请问有没有求EP通量散度的程序？
http://bbs.06climate.com/forum.p ... 7&fromuid=33140
(出处: 气象家园)’


非常感谢资源共享。

沙颖凯

借楼问个关于对p偏导做差分计算的问题。

在等压面p坐标值不均匀（比如1000, 925, 850, 700, 500hPa）的情形下对p做差分有什么好方法，如何考虑这种不均匀，目前我想到的方法是把数据线性插值到标准的1000, 995, 990, 985, ..., 505, 500hPa层次，再按传统的中央差分处理。不知道有没有更好的方法。

phinlot

楼主的疑问指的是？ P坐标下箭头向下对应于Z坐标箭头向上，楼主是指这个吗

吴闫521

不懂啊好高深啊

影伊

phinlot 发表于 2014-8-2 10:24
楼主的疑问指的是？ P坐标下箭头向下对应于Z坐标箭头向上，楼主是指这个吗

是的，看到的文献中箭头是向上的。grads出图纵坐标是p坐标，所以我不知道该怎么处理了。如果把p坐标转化成z坐标，已经改变了偏p，这个是不是出错了？

影伊

沙颖凯 发表于 2014-8-1 21:35
借楼问个关于对p偏导做差分计算的问题。

在等压面p坐标值不均匀（比如1000, 925, 850, 700, 500hPa）的 ...

把1000~500平均分成n等份，然后插值到每等份中？我也不是很懂，抱歉啦

phinlot

不是出错！ EP flux 垂直分量 楼主算出来应该是负值吧，在P 坐标系下负值实际对应的就是poleward heat flux 和 upward propagating Rossby wave。 之所以画出来的EP flux 垂直分量指向下 而不是向我们通常理解的指向上，我想原因是： 负值代表的是指向垂直坐标的负方向，软件默认为该负方向为指向下。但是在P坐标系下，垂直坐标负方向实际是指向上的！所以只要在画图前把EP flux 垂直分量加个负号就好了。

影伊

phinlot 发表于 2014-8-2 15:19
不是出错！ EP flux 垂直分量 楼主算出来应该是负值吧，在P 坐标系下负值实际对应的就是poleward heat flux ...

谢谢你的解惑。之前我有尝试添加负号，但出图只是箭头朝下，如图 我觉得还是不太对。我理解的应该是整个图像反转，越往高层，箭头朝上。grads里有这样的命令么？（图像反转）或者说这只是个投机取巧的方法？
并没有真正解决问题？

影伊

吴闫521 发表于 2014-8-2 10:38
不懂啊好高深啊

刚开始我也只知道个概念，现在遇到这个问题了，也在学习中。。

phinlot

不好意思不能理解楼主说的“越往高层，箭头朝上”是什么意思。。。heat flux在低层大于高层，所以向上的箭头肯定是在低层更明显。不知楼主理想中或者见到的文献里是个什么样子的图？？不好意思没帮上楼主