Parsing 用fortran解析输入文件

这是我父亲的延续

我有一个来自不同代码的文件，我应该

parse

将其用作输入。 其中的一个片段如下所示：

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 

因此，我必须读取每一行，手动查找对应于键的value列，然后将其写入（为了简单起见，我只显示了一个值）。

---

我的问题是，这是正确的方法吗？还是有其他更简单的方法？请让我知道

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

  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

---

或者你想要的任何东西。

Hi@HighPerformanceMark，是的，格式和顺序总是一样的。他们总是在场。它的输出是另一个代码，所以位置是不可变的。如果这20行是代码中唯一的行，那么使用Fortran将数据转储到一个额外的文件中并使用该文件将其读入代码可能会更容易。