| 查看: 1252 | 回复: 4 | ||
[求助]
读取数据文件精度问题！真心求大神赐教 已有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 请问这种情况该怎么解决？ |
回复此楼» 猜你喜欢
投稿咨询
已经有7人回复
-
今天放榜没戏了吧
已经有11人回复
-
我面上完蛋了
已经有5人回复
-
让我中一个面上吧！
已经有12人回复
-
基金啊基金
已经有4人回复
-
放榜前的不淡定
已经有13人回复
-
时间戳今天，20号变了
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有24人回复
-
估计是周四
已经有3人回复
-
科研孤儿太难了
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于磁性的问题，跪求大神赐教
已经有5人回复
- 求大神赐教模拟陶瓷透波率的方法
已经有3人回复
- 有关氧化石墨烯团聚及分散问题，急！求大神帮助！
已经有23人回复
- 这些错误该怎么改啊 ？求大神赐教
已经有3人回复
- 投稿被拒，求大神分析
已经有4人回复
- fluent中面上的受力（force）和力矩（moment）问题，求大神解惑
已经有10人回复
- CIF有B类错误，求帮忙修改
已经有13人回复
- 求大神分析下，好担心
已经有8人回复
- 求大神帮助解决合成中的问题
已经有13人回复
- 小弟研一，求大神推荐些大数据研究方向？
已经有9人回复
- 求大神帮忙用1stOpt 拟合复数数据
已经有13人回复
- 求大神帮忙设计一个工装夹具
已经有6人回复
- 在ansys中做齿轮瞬态温度场分析，遇到问题，求大神指点
已经有9人回复
- Pd XPS数据处理，不会啊，求大神
已经有11人回复
- 求大神告诉下模拟软件
已经有4人回复
- 2,6-二羟甲基-4-叔丁基苯酚的合成，求大神赐教！！！
已经有7人回复
- 计算模拟编程用，求大神推荐笔记本
已经有16人回复
- ansys 13 fluent 并行计算问题
已经有13人回复
2楼2014-05-25 15:51:12
3楼2014-05-26 16:15:32
jiangxingge
金虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2498.6
- 红花: 5
- 帖子: 322
- 在线: 55.4小时
- 虫号: 3205681
- 注册: 2014-05-14
- 性别: GG
- 专业: 石油、天然气地质学
4楼2014-05-26 16:37:28
jiangxingge
金虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2498.6
- 红花: 5
- 帖子: 322
- 在线: 55.4小时
- 虫号: 3205681
- 注册: 2014-05-14
- 性别: GG
- 专业: 石油、天然气地质学
5楼2014-05-26 16:46:53
