E-P通量的Fortran程序

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

想请问一下楼主，在算EP flux的第二项ep2时，有一项是seta的纬向平均对p的偏导。我看到程序里p是从1000hPa到10hPa的排列顺序，而u，v,t都是从10hPa到1000hPa的排列顺序。那么求导数的这项就添加了一个负号，变成正的了，为什么要这么处理呢？
算散度时，第二项也是有一个偏p的求导，程序里似乎又没添加负号。
所以对求通量里的偏p，和求散度里的偏p这项的计算有些无所适从了。。。

cbsdly792 · 发表于 2014-8-20 10:41:52

学习了，感谢分享~谢谢

budapestw · 发表于 2014-8-21 09:32:47

学习一下谢谢楼主

一水天下 · 发表于 2014-9-4 20:55:15

楼主的程序验证过嘛。。对照公式似乎有很多问题。。。

cold_wq · 发表于 2014-10-4 15:12:00

一水天下发表于 2014-9-4 20:55
楼主的程序验证过嘛。。对照公式似乎有很多问题。。。

看你验证的是哪个公式了，我记得这个程序的开头貌似是给出了相应的文献，有哪些问题你可以说出来大家探讨，这个程序是我计算之后，并且和很多别人的计算结果相比较了的

cold_wq · 发表于 2014-10-4 15:14:18

影伊发表于 2014-8-16 16:41
想请问一下楼主，在算EP flux的第二项ep2时，有一项是seta的纬向平均对p的偏导。我看到程序里p是从1000hPa ...

关于这个程序，我当时也是进行了推导的，嗯，因为书上的仅仅给出的是公式，编程的时候很多都是进行了简化的，你再推导下就i知道了，程序没问题这个我敢保证，因为我也用这个程序进行过计算，并且与别人已发表的记过进行过对比，无误

cold_wq · 发表于 2014-10-4 15:15:08

忘却发表于 2014-5-3 20:36
楼主我这几天在算三维波作用通量那个垂直分量总是算的有问题不知道可否借鉴楼主的程序啊？

不好意思，电脑出问题，程序已经找不着了，我当时有传给其它人

cold_wq · 发表于 2014-10-4 15:15:35

柠檬汁儿发表于 2014-3-15 10:23
lz还来么……同求三维for程序

我手上的程序也没有了

longlong · 发表于 2014-11-27 20:18:38

{:eb502:}谢谢分享

Sunnystone24 · 发表于 2015-4-26 10:26:11

怎么点击扣了金币和贡献，没下下载下来呢。。。也没用其他下载工具，就是网页下载啊？？

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

[源代码] E-P通量的Fortran程序