我有一个输入文件（文件名：L0analysis），如下所示
...... !headers
3 2.1 2.4 2.5 !ntfmr rcut(1) rcut(2) rcut(3)
...... !other data
rcut（i）中的索引i由ntfmr
决定在我的程序中，
......
ncrit=55
open(ncrit,file='Lanalysis')
......!(read the headers. skipped here)
read(ncrit,*) ntfmr, (rcut(i),i=1,ntfmr)
do i=1,ntfmr ! # of network formers and their cutoffs
rcut2(i)=rcut(i)*rcut(i)
write(*,*)'ntfmr',ntfmr,'i=',i,rcut(i)
enddo
.......
.......
我的输出是
ntfmr 3 i = 1 2.10000000000000
ntfmr 3 i = 2 2.40000000000000
ntfmr 3 i = 3 4.41000000000000
我试图初始化rcut（i）= 0。这是行不通的。我找不到任何改变这个值的地方。 我无法弄清楚出了什么问题。任何建议或评论都将受到高度赞赏。