计算理查逊数的一个小程序，一直看不出错在哪里。。

syl136 · 发表于 2016-11-29 16:03:10

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

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

x

本帖最后由 syl136 于 2016-11-29 20:15 编辑

登录/注册后可看大图

计算理查逊数结果，结果很奇怪，数值很大。。。拜托大家帮忙看看哪里有问题，谢谢了。
变量单位应该都没有问题：
pt       Potential Temperature (K)
t       Temperature (K)
height    Model height (m)
pressure Model pressure (Pa)
u          x-wind component (m s-1)
v          y-wind component (m s-1)

pt是计算的位温

RI=-9999.0
pt=t*(1000/p)**0.286
do k=2,k1-1
do j=2,j1-1
do i=2,i1-1
N2=((9.8/pt(i,j,k))*((pt(i,j,k-1)-pt(i,j,k+1))/(2*(hgt(i,j,k-1)-hgt(i,j,k+1)))))
uz=(u(i,j,k-1)-u(i,j,k+1))/(2*(hgt(i,j,k-1)-hgt(i,j,k+1)))
vz=(v(i,j,k-1)-v(i,j,k+1))/(2*(hgt(i,j,k-1)-hgt(i,j,k+1)))
VWS=uz*uz+vz*vz
if(VWS==0)then
RI(i,j,k)= -9999
else
RI(i,j,k)= N2/VWS
end if
enddo
enddo
enddo

jjpzhzz · 发表于 2016-11-29 19:28:21

同问。。。。。

wow7788905 · 发表于 2016-11-29 19:29:19

11111111111

Lighting · 发表于 2016-11-29 19:55:52

之前计算的也是出现一些问题。FORTRAN版的程序论坛好像已经有现成的了

syl136 · 发表于 2016-11-29 20:03:19

jjpzhzz 发表于 2016-11-29 19:28
同问。。。。。

...你们！！！还给我！！！

syl136 · 发表于 2016-11-29 20:03:41

wow7788905 发表于 2016-11-29 19:29
11111111111

！！！！！我能骂人嘛？！！

syl136 · 发表于 2016-11-29 20:04:33

Lighting 发表于 2016-11-29 19:55
之前计算的也是出现一些问题。FORTRAN版的程序论坛好像已经有现成的了

哦哦，谢谢您，我只找到了GRADS的程序，fortran的好像没看到。。。我再去看看

freekiller · 发表于 2016-11-29 20:12:14

{:eb301:}{:eb301:}{:eb301:}{:eb301:}

syl136 · 发表于 2016-11-29 20:17:46

freekiller 发表于 2016-11-29 20:12

放一回金钱。。全被办公室的抢完了。。

freekiller · 发表于 2016-11-29 20:19:23

syl136 发表于 2016-11-29 20:17
放一回金钱。。全被办公室的抢完了。。

，限制级别低的不准抢

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

[求助] 计算理查逊数的一个小程序，一直看不出错在哪里。。

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

点评

回帖奖励 +30 金钱

评分

回帖奖励 +30 金钱

评分