请教FORTRAN从高斯输出中格式化读入问题

ChemiAndy

今天编译以前组里人写的一个小code能编译但运行出错。以前用低版本的f77编译运行肯定是没问题的，[size=13.63636302948px]现在用了高版本的gfortran就不行。功能很简单，从高斯输出的SCF Done这一行读入能量值和所用的循环数。

code相关部分如下：

1.    10 READ(10,'(A)',END=20)LINE
   
2.    25 FORMAT(' SCF Done:  E(RHF) = ',G19.11,' A.U. after ',I4,
   
3.      &  ' cycles')
   
4.       IF (INDEX(LINE,'SCF Done').NE.0) THEN
   
5.          READ(LINE,25) V,ITER
   
6.       ENDIF

复制代码

出错提示，

1. Fortran runtime error: Constant string in input format
   
2. (' SCF Done:  E(RHF) = ',G19.11,' A.U. after ',I4,            ' cycles')

复制代码

高斯输出文件这一行的内容为：

1. SCF Done:  E(RHF) =  -1.06610864862     A.U. after    1 cycles
   

复制代码

谁能帮着看一下应该怎么改。

ChemiAndy

sobereva 发表于 2014-11-14 11:36
这么写是违背语法的。write的时候能用25号格式，而在读取的时候不能把本来要读的内容直接写在格式里。
如 ...

Sob老大说的很对。我已经按照你的说法改过code了。
原来的code是老板写的。老。板。写。的。。。

ChemiAndy

sobereva 发表于 2014-11-14 11:36
这么写是违背语法的。write的时候能用25号格式，而在读取的时候不能把本来要读的内容直接写在格式里。
如 ...

Sob老大说的很对。我已经按照你的说法改过code了。
原来的code是老板写的。老。板。写。的。。。

sobereva

这么写是违背语法的。write的时候能用25号格式，而在读取的时候不能把本来要读的内容直接写在格式里。
如果你想空过字符部分，应该用比如24x,f19.11,10x,i4 ，x前面的数代表读的时候跳过的字符数

ChemiAndy

搞定了。格式化行的字符都用X编辑符代替就可以了

1.    25 FORMAT(21X,G19.12,12X,I4,7X)

复制代码