在fortran中解析输入文件

Question

这是我的旧文章 thread 的延续.

我有一个来自不同代码的文件，我应该解析以用作我的输入。其中的一个片段如下所示:

GLOBAL SYSTEM PARAMETER 
NQ                 2 
NT                 2 
NM                 2 
IREL               3 
************************************* 
BEXT        0.00000000000000E+00 
SEMICORE  F 
LLOYD     F 
NE                32         0 
IBZINT             2 
NKTAB            936 
XC-POT    VWN       
SCF-ALG   BROYDEN2   
SCF-ITER          29 
SCF-MIX     2.00000000000000E-01 
SCF-TOL     1.00000000000000E-05 
RMSAVV      2.11362995016878E-06 
RMSAVB      1.25411205586140E-06 
EF          7.27534671479201E-01 
VMTZ       -7.72451391270293E-01 
************************************* 

等等。

目前我正在逐行阅读它，如下:

Program  readpot 
use iso_fortran_env 
Implicit None 
integer ::i,filestat,nq 
character(len=120):: rdline 
character(10)::key!,dimension(:),allocatable ::key 
real,dimension(:),allocatable ::val 
i=0 

open(12,file="FeRh.pot_new",status="old") 
readline:do 
  i=i+1 
  read(12,'(A)',iostat=filestat) rdline!(i) 

  if (filestat /= 0) then 
    if (filestat == iostat_end ) then 
      exit readline 
    else 
      write ( *, '( / "Error reading file: ", I0 )' )  filestat 
      stop 
    endif 
  end if 

  if (rdline(1:2)=="NQ") then 
    read(rdline(19:20),'(i)'),nq 
    write(*,*)nq 
  end if 
end do readline 

End Program  readpot 

因此，我必须读取每一行，手动找到与键对应的值列，然后写入(为简洁起见，我只显示了一个值)。我的问题是，这是执行此操作的正确方法吗？或者还有其他更简单的方法？请告诉我。

Answer 1

如果文件没有可变性，您几乎不需要解析它。假设您已为文件中所有感兴趣的数据项声明了变量，并且这些变量的名称显示在文件的行上。例如

  INTEGER :: nq , nt, nm, irel
  REAL:: scf_mix, scf_tol  ! '-' not allowed in Fortran names
  CHARACTER(len=48) :: label, text
  LOGICAL :: semicore, lloyd
  ! Complete this as you wish

然后写一段这样的代码

  OPEN(12,file="FeRh.pot_new",status="old") 
  READ(12,*) ! Not interested in the 1st line
  READ(12,*) label, nq
  READ(12,*) label, nt
  READ(12,*) label, nm
  READ(12,*) label, irel
  READ(12,*) ! Not interested in this line
  READ(12,*) label, bext
  READ(12,*) label, semicore
  ! Other lines to write
  CLOSE(12)

Fortran 的列表定向输入可以理解行中的空格来分隔值。它不会将这些空格读取为字符变量的一部分。这种行为可以改变，但在你的情况下你不需要这样做。请注意，当读入逻辑变量时，它还会将字符 F 理解为 .false.。

我的代码片段只是忽略了标签和解释行。如果你性格紧张，你可以处理它们，也许

IF (label/='NE') STOP

或者任何你想要的。