Fortran：如何从文本文件的每一行中读取第一个字符？

Question

这是我第一次尝试在Fortran中编程。 我正在尝试编写一个程序，该程序打印斐波那契数列的前1476个项，然后检查每个项的第一个数字，并存储数组中出现的1s，2s，3s，...，9s的数量。

我似乎无法弄清楚的问题是如何阅读每个术语的第一位数字。 我已经尝试了几件事，但是由于对Fortran技术的了解有限而遇到困难。 我将这些术语写到一个文本文件中，其想法是读取每一行的第一位并在数组中累加相应的数字。 有人对此有任何建议吗？

到目前为止，这是我的代码：

（ 编辑 ：我包括了读取文件的代码。现在它仅打印出3.60772951994415996E-313，它看起来像某种地址，因为它不是斐波那契数字之一。而且，这是唯一的事情打印，我希望它将打印出文件的每一行...）

（ 编辑编辑 ：考虑到这一点之后，也许有一种方法可以将写入文本文件的格式设置为仅第一位。是否可以将实数的有效位数设置为一个？：P）

subroutine writeFib(n)
  integer ::  i
  real*8 :: prev, current, newFib
  prev = 0
  current = 1
  do i = 1, n
     newFib = prev + current
     prev = current
     current = newFib
     write(7,*) newFib
  end do
  return
end subroutine

subroutine recordFirstDigits(a)
  integer :: openStat, inputStat
  real*8 :: fibNum
  open(7, file = "fort.7", iostat = openStat)
  if (openStat > 0) stop "*** Cannot open the file ***"
  do
     read(7, *, iostat = inputStat) fibNum
     print *,fibNum
     if (inputStat > 0) stop "*** input error ***"
     if (inputStat < 0) exit ! end of file
  end do
  close(7)
end subroutine

program test
  integer :: k, a(9)
  k = 1476
  call writeFib(k)
  call recordFirstDigits(a)
end program

Answer 1

可以阅读FORMAT（A1）吗？ 已经20年了，所以我不记得确切的语法了。

Answer 2

我想知道为什么未关闭文件7时open语句成功。 我认为您在读写之间需要一个endfile语句和/或一个rewind语句。

Paul Tomblin发布了您解决问题后首先要做的事情。

Answer 3

尽管提出了建议，但也忘记了几件事。 REAL类型的范围，以及一些格式问题。

无论如何，这是一个经过修补的解决方案，可以编译并正常工作，因此请尝试看看是否适合您。 我自由选择了自己的斐波纳契数计算方法。

  program SO1658805
  implicit none

  integer, parameter :: iwp = selected_real_kind(15,310)
  real(iwp) :: fi, fib
  integer :: i
  character(60) :: line
  character(1) :: digit
  integer :: n0=0, n1=0, n2=0, n3=0, n4=0, n5=0, n6=0, n7=0, n8=0, n9=0

  open(unit=1, file='temp.txt', status='replace')
  rewind(1)
  !-------- calculating fibonacci numbers -------
  fi = (1+5**0.5)/2.
  do i=0,1477
    fib = (fi**i - (1-fi)**i)/5**0.5
    write(1,*)fib,i  
  end do
  !----------------------------------------------
  rewind(1)

  do i=0,1477
    read(1,'(a)')line
    line = adjustl(line)
    write(*,'(a)')line


    read(line,'(a1)')digit

     if(digit.eq.' ') n0=n0+1
     if(digit.eq.'1') n1=n1+1
     if(digit.eq.'2') n2=n2+1
     if(digit.eq.'3') n3=n3+1
     if(digit.eq.'4') n4=n4+1
     if(digit.eq.'5') n5=n5+1
     if(digit.eq.'6') n6=n6+1
     if(digit.eq.'7') n7=n7+1
     if(digit.eq.'8') n8=n8+1
     if(digit.eq.'9') n9=n9+1
  end do
  close(1)

  write(*,'("Total number of different digits")')
  write(*,'("Number of digits 0: ",i5)')n0
  write(*,'("Number of digits 1: ",i5)')n1
  write(*,'("Number of digits 2: ",i5)')n2
  write(*,'("Number of digits 3: ",i5)')n3
  write(*,'("Number of digits 4: ",i5)')n4
  write(*,'("Number of digits 5: ",i5)')n5
  write(*,'("Number of digits 6: ",i5)')n6
  write(*,'("Number of digits 7: ",i5)')n7
  write(*,'("Number of digits 8: ",i5)')n8
  write(*,'("Number of digits 9: ",i5)')n9

  read(*,*)

  end program SO1658805

w，...我刚刚读到您需要存储在数组中的位数。 虽然我只是数了数。

哦，好吧，……“留给读者练习……” :-)

Answer 4

我收到“行尾”运行时错误

您不显示! code to read here... ! code to read here...这使得很难猜测您在做什么错了:-)

也许您需要一个循环来读取每一行，然后在没有更多行时跳出循环以继续执行语句。

像这样：

   do
      read(7,*,end=10) fibNumber
   end do
   10 continue

最好还是看一下此revcomp程序中使用的更现代的样式 。

Answer 5

这里有一些提示：

1. 您不需要使用character ， 只需很少的文件i / o即可解决此问题（除非您忘记声明必须创建文件）。
2. 因此，请使用数学来查找统计信息 。 斐波纳契数有很多资源，它们可能提供简化的见解，或者至少是一种独立检查您答案的方法。
   
   • 这是非Fortran语言中的一个复杂提示：
     floor(10^(frac(log_10(7214989861293412))))
     （将其放在Wolfram Alpha中以查看其功能。）
     
   • 一个更简单的提示（针对另一种方法）是，您可以在Fortran中使用循环结构内部的简单算术来做得很好-至少是在解决方案的第一遍。
3. 随时累积统计数据 。 该建议甚至适用于您的角色驱动方法。 （此问题非常适合为您的统计数据提出一个可爱的索引方案，但是有些人讨厌编程中的可爱方案。如果您不担心可爱性……那么只要您的键在Fortran中就可以有关联数组是整数;-)
4. 此问题最重要的方面是将用于计算答案的数据类型 。 例如， 这是您必须打印的最后一个数字 。

干杯，--Jared