24小时热门版块排行榜

>论坛更新日志 (2230)
>导师招生 (802)
>虫友互识 (410)
>休闲灌水 (395)
>招聘信息布告栏 (126)
>论文投稿 (90)
>博后之家 (78)
>文献求助 (49)
>硕博家园 (38)
>教师之家 (35)
>考研 (31)
>学术会议 (28)
>公派出国 (22)
>考博 (19)
>基金申请 (16)
>找工作 (11)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

dongdaozhu

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 62
帖子: 15
在线: 11.9小时
虫号: 3115518
注册: 2014-04-05
专业: 大地测量学

[求助] 读取数据文件精度问题！真心求大神赐教已有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  请问这种情况该怎么解决？

回复此楼