爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3221|回复: 1

[求助] 计算ri 出现问题 ,这是为啥呢??结果不对

[复制链接]

新浪微博达人勋

发表于 2013-1-29 15:59:07 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 爱的侍者 于 2013-1-29 21:06 编辑

用matlab计算结果正确
用fortran计算结果不对
贴出一段程序看一下为啥呢
              -51.76    6.7    6.6    0.6   -0.3
这是fortran算出来的
第一个为结果其他四个 分别为sumu(i,1)sumu(i,4)stmt(i,1)stmt(i,2)

if(stmt(i,1).ne.99999.and.stmt(i,2).ne.99999.and.sumu(i,1).ne.99999.and.sumu(i,4).ne.99999)then
ri(i)=(9.8/(stmt(i,1)+273.15))*((stmt(i,2)-stmt(i,1))/60+0.0065)/((sumu(i,4)-sumu(i,1))/60)**2
else
ri(i)=99999
endif
write(25,'(f20.1,1x,f7.1,1x,f7.1,1x,f7.1,1x,f7.1)')ri(i),sumu(i,2),sumu(i,4),stmt(i,1),stmt(i,2)

这样出了问题

正确结果
>> (9.8/(0.6+273.15))*((-0.3-0.6)/60+0.0065)/(((6.6-6.7)/60).^2)
ans =
-109.5452
应该是这个
为啥呢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-1-29 21:31:29 | 显示全部楼层
12345.png
左边matlab算的
结果正确
右边fortran算的结果存在问题 ,这是为啥
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表