24小时热门版块排行榜    

查看: 989  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

dongdaozhu

新虫 (初入文坛)

[求助] 读取数据文件精度问题!真心求大神赐教 已有1人参与

程序代码
module typedef
  type DATE
    character(len=3) :: Key
        integer :: L,M
    real(8) :: C,S
  end type
  end module
  program resault
   use typedef
   implicit none
   type(DATE) :: s(1913)
   integer i,k
   character( len = 2 ) :: cTemp
!定义变量
do k = 1,2
  write( cTemp,'(i2)' ) k
  open ( 1, file = '' // trim(adjustl( cTemp )) //'.gfc', status = 'old' )
  open ( 2, file = '' // trim(adjustl( cTemp )) //'.txt' )
  do i=1,22
   read(1,*)
  ENDDO
!前22行不需要读取
  do i=23,1913
    read(1,*) s(i)%Key, s(i)%L, s(i)%M, s(i)%C, s(i)%S
end do
do i=23,1913
write (2,*)s(i)
end do
close (2)
close (1)
end do
end program
读取的源文件数据格式如下:
gfc     0    0  0.100000000000E+01  0.000000000000E+00  0.0000E+00  0.0000E+00
gfc     1    0  0.000000000000E+00  0.000000000000E+00  0.0000E+00  0.0000E+00
gfc     2    0  -.484169206174E-03  0.000000000000E+00  0.0000E+00  0.0000E+00
gfc     3    0  0.957185436561E-06  0.000000000000E+00  0.0000E+00  0.0000E+00
gfc     4    0  0.539987318528E-06  0.000000000000E+00  0.0000E+00  0.0000E+00
gfc     5    0  0.686421621455E-07  0.000000000000E+00  0.0000E+00  0.0000E+00
gfc     6    0  -.149983726101E-06  0.000000000000E+00  0.0000E+00  0.0000E+00
gfc     7    0  0.904960683013E-07  0.000000000000E+00  0.0000E+00  0.0000E+00
gfc     8    0  0.494749421093E-07  0.000000000000E+00  0.0000E+00  0.0000E+00
gfc     9    0  0.280128784478E-07  0.000000000000E+00  0.0000E+00  0.0000E+00
gfc    10    0  0.533523562152E-07  0.000000000000E+00  0.0000E+00  0.0000E+00
程序运行后 ;
gfc           0           0   1.00000000000000       0.000000000000000E+000
gfc           1           0  0.000000000000000E+000  0.000000000000000E+000
gfc           2           0 -4.841692061740000E-004  0.000000000000000E+000
gfc           3           0  9.571854365609999E-007  0.000000000000000E+000
gfc           4           0  5.399873185280000E-007  0.000000000000000E+000
gfc           5           0  6.864216214550000E-008  0.000000000000000E+000
gfc           6           0 -1.499837261010000E-007  0.000000000000000E+000
gfc           7           0  9.049606830130000E-008  0.000000000000000E+000
gfc           8           0  4.947494210930000E-008  0.000000000000000E+000
gfc           9           0  2.801287844780000E-008  0.000000000000000E+000
gfc          10           0  5.335235621520000E-008  0.000000000000000E+000
可以看到第四行的与原来的数据不同了  0.957185436561E-06 变为了9.571854365609999E-007  后面多出了9999  请问这种情况该怎么解决?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiangxingge

金虫 (正式写手)

引用回帖:
4楼: Originally posted by jiangxingge at 2014-05-26 16:37:28
REAL(8)  保持了 12位的精度 很不错了!
你控制一下 输出 的 格式 ,不要留那么长的位数,它会自动四舍五入 ,这样 就 跟原来 一模一样了。

比如
write(2,'(a3,i5,i5,2E20.12E2)') s(i)
5楼2014-05-26 16:46:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

dongdaozhu

新虫 (初入文坛)

真心求大神啊!!!
2楼2014-05-25 15:51:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiangxingge

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
REAL(8)  保持了 12位的精度 很不错了!
你控制一下 输出 的 格式 ,不要留那么长的位数,它会自动四舍五入 ,这样 就 跟原来 一模一样了。
4楼2014-05-26 16:37:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见