应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 读取数据文件精度问题！真心求大神赐教
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1058  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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  请问这种情况该怎么解决？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-05-24 17:54:19
已阅   回复此楼   关注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 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600，320分求调剂 +15 大馋小子 2026-04-04 16/800 2026-04-06 06:34 by houyaoxu
[考研] 材料求调剂 +11 一样YWY 2026-04-05 11/550 2026-04-05 23:36 by 来看流星雨10
[考研] 288求调剂，一志愿华南理工大学071005 +6 ioodiiij 2026-04-04 6/300 2026-04-05 10:09 by guoweigw
[考研] 070300化学学硕311分求调剂 +10 梁富贵险中求 2026-04-04 12/600 2026-04-05 09:37 by guoweigw
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
[考研] 085601，一志愿厦大334复试被刷求调剂 +13 曾仰之 2026-04-03 15/750 2026-04-04 20:13 by dongzh2009
[考研] 306求调剂 +3 hyb上名工 2026-04-02 3/150 2026-04-04 18:12 by 热情沙漠
[考研] 0710生物学336分求调剂 +6 kiyy 2026-04-01 8/400 2026-04-04 10:10 by kiyy
[考研] 兽医调剂 +3 wh119216 2026-04-02 3/150 2026-04-03 19:34 by zrongyan
[考研] 293求调剂 +5 末未mm 2026-04-02 6/300 2026-04-03 15:20 by 王保杰33
[考研] 282求调剂 不挑专业 求收留 +7 Yam. 2026-03-30 8/400 2026-04-03 14:12 by zhangdingwa
[考研] 366求调剂 +7 sbdnd 2026-04-03 7/350 2026-04-03 12:40 by cymywx
[考研] 330求调剂 +3 白神呜呼呼 2026-04-02 3/150 2026-04-03 10:15 by 蓝云思雨
[考研] 318求调剂，计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[考研] 07生物学求调剂 一志愿同济大学359分 +3 LAMC. 2026-03-30 3/150 2026-04-02 10:26 by 18828373951
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 11:35 by wangjy2002
[考研] 合肥区域性重点一本招收调剂 +4 6266jl 2026-03-30 8/400 2026-03-31 18:43 by 6266jl
[考研] 考研调剂求助 +7 13287130938 2026-03-31 7/350 2026-03-31 16:39 by 690616278
[考研] 哈尔滨工业大学材料与化工专硕378求调剂 +3 塔比乌斯 2026-03-30 3/150 2026-03-30 22:55 by 无际的草原
[考研] 调剂 +4 GK72 2026-03-30 4/200 2026-03-30 20:32 by dick_runner
信息提示
关闭
请填处理意见
关闭