【已解决】Fortran程序出错,好像是精度问题

学而不厌

非常简单的程序，用来求一个数组的和，但是出了问题。程序如下：
program main
   implicit none
   integer,parameter:: nx=28,nt=36,nn=11
   real b(nt),Rd(nx,nt),sum
   integer i,k,year(nt)

   open(21,file='Dec_1981-2016_Rd.txt',form='formatted')  
   do k=1,nt
       read(21,*) year(k),(Rd(i,k),i=1,nx)
   enddo
   close(21)

   open(22,file='sum.txt',form='formatted')  
   sum=0.0
   do k=1,nt
       b(k)=Rd(24,k)
       sum=sum+b(k)
       write(22,*) b(k),sum
   enddo
end

sum.f90 (430 Bytes, 下载次数: 0)

2017-12-25 21:43 上传

点击文件名下载附件     Dec_1981-2016_Rd.txt (16.37 KB, 下载次数: 2)

2017-12-25 21:44 上传

点击文件名下载附件
sum.txt (1.13 KB, 下载次数: 0)

2017-12-25 21:45 上传

点击文件名下载附件


最终答案本来应该是85.3，但是计算结果缺显示为85.30001。查看一下应该实在黄色标识的地方出现错误（2010年），这是为什么？有对其他几列进行计算，发现有的计算对了，有的有问题，尽管误差不大，但是这到底是什么原因造成的。



查了一下书，表控输出一个实型量占11列，其中一位小数点，七位小数。单精度实型变量（real）的有效位数是6~7。
可能是由于表控输出超过了有效位数的原因，造成输出错误的？
如果是这样的话，每次运算结果不是应该都不一样的，超出有效位数的部分应该会发生变化才对吧！但是我的计算结果一直没有变化。
求哪位大神帮忙解惑一下吧？

学而不厌

哪位大神帮忙解惑一下吧

香草拿铁

输入sum的时候定义一下只输出小数点后一位

学而不厌

香草拿铁 发表于 2017-12-26 09:19
输入sum的时候定义一下只输出小数点后一位

主要是想搞明白问题出在哪里了

weinihou

大体是你说的，但这不叫“错误”，是数值计算 一定会碰到的  数值精度的问题

如果真想尽量减少，可以考虑提升下变量精度，比如 real*8 之类的声明法
或者板凳建议的，格式输出一下，眼不见心不烦