24小时热门版块排行榜    

查看: 1051  |  回复: 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的回帖

dongdaozhu

新虫 (初入文坛)

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

dongdaozhu

新虫 (初入文坛)

真伤心!!
3楼2014-05-26 16:15:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiangxingge

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
REAL(8)  保持了 12位的精度 很不错了!
你控制一下 输出 的 格式 ,不要留那么长的位数,它会自动四舍五入 ,这样 就 跟原来 一模一样了。
4楼2014-05-26 16:37:28
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 dongdaozhu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 东南大学364求调剂 +4 JasonYuiui 2026-03-15 4/200 2026-03-16 08:36 by Linda Hu
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考研] 材料工程专硕274一志愿211求调剂 +5 薛云鹏 2026-03-15 5/250 2026-03-15 20:38 by Logic2024
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +3 Losir 2026-03-12 3/150 2026-03-14 12:11 by 热情沙漠
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 一志愿浙江大学0856材料与化工求调剂 +4 yansheng@211 2026-03-09 5/250 2026-03-14 02:10 by JourneyLucky
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 材料工程专硕,一志愿中国矿业大学,总分314,求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 材料与化工304求B区调剂 +5 邱gl 2026-03-11 6/300 2026-03-13 22:37 by JourneyLucky
[考研] 333求调剂 +3 球球古力 2026-03-11 3/150 2026-03-13 21:27 by JourneyLucky
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 285化工学硕求调剂(081700) +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[考研] 一志愿河海大学085900土木水利专硕279求调剂不挑专业 +4 SunWwWwWw 2026-03-10 8/400 2026-03-13 02:23 by SunWwWwWw
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
[考研] 294 英二数二物化 求调剂 +6 米饭团不好吃 2026-03-09 6/300 2026-03-09 23:55 by barlinike
信息提示
请填处理意见