关于EP通量请教一下

影伊 · 发表于 2014-8-1 21:09:58

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由影伊于 2015-6-4 20:20 编辑

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

说明，求的是行星波1~3波，对u，v，t作了谐波分析。

1991-2014EPD.gmf (183.65 KB, 下载次数: 64)

影伊 · 发表于 2014-8-6 19:23:42

这个问题我已经解决。在搜索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
(出处: 气象家园)’

非常感谢资源共享。

沙颖凯 · 发表于 2014-8-1 21:35:04

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

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

phinlot · 发表于 2014-8-2 10:24:30

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

吴闫521 · 发表于 2014-8-2 10:38:24

不懂啊好高深啊

影伊 · 发表于 2014-8-2 13:52:48

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

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

影伊 · 发表于 2014-8-2 13:56:46

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

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

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

phinlot · 发表于 2014-8-2 15:19:24

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

影伊 · 发表于 2014-8-2 16:33:33

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

谢谢你的解惑。之前我有尝试添加负号，但出图只是箭头朝下，如图

我觉得还是不太对。我理解的应该是整个图像反转，越往高层，箭头朝上。grads里有这样的命令么？（图像反转）或者说这只是个投机取巧的方法？
并没有真正解决问题？

影伊 · 发表于 2014-8-2 16:34:40

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

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

phinlot · 发表于 2014-8-2 18:13:01

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

		自动登录	找回密码
密码			立即注册

[求助] 关于EP通量请教一下

登录后查看更多精彩内容~