24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1059  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 312求调剂 +4 gtw1 2026-04-06 4/200 2026-04-06 07:44 by 无际的草原
[考研] 一志愿南昌大学,085600,344分求调剂 +8 调剂上岸玘 2026-04-05 9/450 2026-04-06 07:09 by houyaoxu
[考研] 085600,321分求调剂 +11 大馋小子 2026-04-04 12/600 2026-04-06 06:29 by houyaoxu
[考研] 085600,320分求调剂 +7 大馋小子 2026-04-01 8/400 2026-04-05 21:19 by 学员8dgXkO
[考研] 0703化学321分求调剂 +17 三dd. 2026-03-30 18/900 2026-04-05 18:07 by 蓝云思雨
[考研] 0860 求调剂 一志愿国科大 348 分 +3 WiiiP 2026-04-03 3/150 2026-04-05 17:43 by Ecowxq666!
[考研] 考研生物学考A区211,初试322,科目生化和生物综合,求调剂 +6 。。。54 2026-04-03 6/300 2026-04-05 14:54 by JOKER0401
[考研] 专硕310求调剂 +5 捞捞我…. 2026-04-04 6/300 2026-04-04 23:33 by barlinike
[考研] 368求调剂 +5 今华习 2026-04-03 7/350 2026-04-04 18:47 by imissbao
[考研] 11408 一志愿西电,277分求调剂 +4 zhouzhen654 2026-04-03 4/200 2026-04-04 18:10 by 猪会飞
[考研] 291求调剂 +4 迷蒙木木 2026-04-01 5/250 2026-04-04 15:59 by sihailian3
[考研] 求调剂,一志愿南京航空航天大学 ,080500材料科学与工程学硕 +10 @taotao 2026-04-03 10/500 2026-04-04 09:01 by T可可西里T
[考研] 322求调剂 +6 FZAC123 2026-04-03 6/300 2026-04-03 22:23 by 科研小专家
[考研] 311求调剂 +11 勇敢的小吴 2026-04-02 11/550 2026-04-03 21:46 by qlm5820
[考研] 学硕288调剂!!! +3 小王xw123 2026-04-03 3/150 2026-04-03 21:20 by 啵啵啵0119
[考研] 材料科学与工程339求调剂 +12 hyz0119 2026-03-31 13/650 2026-04-03 18:33 by ls刘帅
[考研] 考研调剂 +8 不爱喝饮料 2026-04-03 8/400 2026-04-03 16:40 by Mistake-J
[考研] 321求调剂 +17 y-yh 2026-04-01 20/1000 2026-04-03 12:57 by y-yh
[考研] 085600 一志愿9 总分351 求调剂学校 +7 czhcz 2026-03-31 9/450 2026-04-01 19:24 by 无际的草原
[考研] 349求调剂 +6 吃的不少 2026-04-01 6/300 2026-04-01 17:55 by JYD2011
信息提示
请填处理意见