有没有人帮我看下这个程序里面的计算是不是对的呀。。怎么和降水相关性不显著呢

kzpyezi

cris 发表于 2013-5-27 16:02

                               
登录/注册后可看大图

你的goto 122，但是122在哪里呢？还有就是你在比较水汽压和饱和水汽压之后，将温度以0.005的减少量逐渐循 ...

122        IF(E(I,J,K,it).lt.ES(I,J,K,it)) then
         T(i,j,k,it)=T(i,j,k,it)-0.005
          Es(I,J,K,it)=6.11*exp(BA*t(I,J,K,it)/(273.16-BB+t(I,J,K,it)))
      goto 122

kzpyezi

cris 发表于 2013-5-27 16:02

                               
登录/注册后可看大图

你的goto 122，但是122在哪里呢？还有就是你在比较水汽压和饱和水汽压之后，将温度以0.005的减少量逐渐循 ...

而且这个位温公式是1000hpa的，我查了的，是对的啊。

cris

kzpyezi 发表于 2013-5-27 16:12

                               
登录/注册后可看大图

我木有查到位温的计算公式，请问你知道么？  冰面的和睡眠的系数是不同的呀

你的位温公式是不是手误写错了呀。应该是 TP(i,j,K,IT)=TT(I,J,K,IT)*(1000./P(i,j,K,IT))**(RD/CPD) ,RD/CPD=0.286.

cris

kzpyezi 发表于 2013-5-27 16:12

                               
登录/注册后可看大图

我木有查到位温的计算公式，请问你知道么？  冰面的和睡眠的系数是不同的呀

IF(T(I,J,K),IT.GE.0)ES(I,J,K,IT)=6.112*10**(7.5*T(I,J,K,IT)/(237.3+T(I,J,K,IT))) !对于水面而言
                IF(T(I,J,K,IT).LT.0)ES(I,J,K,IT)=6.112*10**(9.5*T(I,J,K,IT)/(265.5+T(I,J,K,IT))) !冰面,这是大气物理学那本书上的计算水面饱和水汽压和冰面饱和水汽压的公式。

liutaoz21

kzpyezi 发表于 2013-5-27 11:29

                               
登录/注册后可看大图

贴程序么？还是资料呀？

是贴程序...

蔚蓝的天

顶一个哈哦

wuyula

。。。。。。。。。。。。。学习