计算涡度、涡度平流、温度平流、散度、水汽通量、水汽通量散度、比湿等的FORTRAN程序

3870383 · 发表于 2016-1-23 17:19:34

@river 楼主请教下，这程序中有错误，不知道仔细看了没有，比如
!------------------------------涡度平流--------------------------------------------
vorop=1.e+36
DO IT=1,NT
DO IZ=1,NZ

DO i=1+2,NX-2
DO j=1+2,NY-2
vorop(i,j,iz,it)=u(i,j,iz,it)*((voro(i+1,j,iz,it)-voro(i-1,j,iz,it))/(2*R*cos(lat(j))*Delta*pi/180))+v(i,j,iz,it)*((voro(i,j+1,iz,it)-vorg(i,j-1,iz,it))/(2*R*Delta*pi/180))
ENDDO
ENDDO

ENDDO
ENDDO

按照您的定义，voro、vorg分别是相对涡度和地转涡度，这里求平流时似乎不太对奥

river · 发表于 2016-1-23 21:17:38

3870383 发表于 2016-1-23 17:19
@river 楼主请教下，这程序中有错误，不知道仔细看了没有，比如
!------------------------------涡度平 ...

这个老早以前的了，你看看天原上的公式吧，这个都是根据天原上的公式编的

d1268y · 发表于 2016-1-27 17:30:38

{:eb502:}

d1268y · 发表于 2016-1-27 17:38:39

谢谢分享！

大西洋的风 · 发表于 2016-1-29 12:06:25

谢谢，正在琢磨

voro(i,j,iz,it)=1./(R*2)*( (v(i+1,j,iz,it)-v(i-1,j,iz,it))/(Delta*pi/180*cos(lat(j)))-(u(i,j+1,iz,it)-u(i,j-1,iz,it))/(Delta*pi/180)+2*u(i,j,iz,it)*tan(lat(j)) )

复制代码

这个不需要改动吧

river · 发表于 2016-1-30 17:41:46

大西洋的风发表于 2016-1-29 12:06
谢谢，正在琢磨

这个不需要改动吧

不需要，就是书上的公式

river · 发表于 2016-1-30 17:41:58

大西洋的风发表于 2016-1-29 12:06
谢谢，正在琢磨

这个不需要改动吧

不需要，就是书上的公式

joyxysetu · 发表于 2016-1-30 21:10:13

回复下赚金币

qiulin_first · 发表于 2016-3-10 09:37:47

谢谢楼主，最喜欢这样的帖子了

qiulin_first · 发表于 2016-3-10 09:38:47

对了，里面有多少个物理量啊，想下一个比较全的程序。

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

[源代码] 计算涡度、涡度平流、温度平流、散度、水汽通量、水汽通量散度、比湿等的FORTRAN程序