Originally posted by snoopyzhao at 2011-05-08 17:37:30:
你最好还是把运行程序所有的代码+数据文件打包放在附件里吧，要不然，别人也没有办法测试啊，仅仅靠看的功力还是不够的……

Originally posted by Gina88 at 2011-05-08 17:54:09:
附件是我上传的CHGCAR，您可以将它复制一下，一份命名为CHGCAR1，另一份命名为CHGCAR2

program test
implicit none
integer :: i, k
integer :: line_len ! length of element line without leading and tailing space.
integer :: atomsum
character(len=80) :: line
character(len=2), dimension(:), allocatable :: element
integer, dimension(:), allocatable :: atomnum

open(unit=12, file = 'CHGCAR')

do i=1,5
read(12,*)
end do

read(12,'(a)') line
line = adjustl(line)
line_len = len_trim(line)

atomsum=1
do i=1,line_len
   if(line(i:i) == ' ' .and. line(i:i+1) /= ' ') atomsum=atomsum+1
end do

write(*,*) atomsum

allocate(element(atomsum))
allocate(atomnum(atomsum))

k = 1
do i=1,line_len
   if(line(i:i) /= ' ' .and. line(i:i+1) == ' ') then
       element(k) = line(i:i)
       k = k + 1
    end if
   if(line(i:i) /= ' ' .and. line(i:i+1) /= ' ') then
       element(k) = line(i:i+1)
       line(i+1:i+1) = ' '
       k = k + 1
    end if
end do


read(12, *) (atomnum(i), i=1,atomsum)

write(*, *) (atomnum(i), i=1,atomsum)
write(*, '(a5)') (element(i), i=1,atomsum)

end program test

Originally posted by snoopyzhao at 2011-05-08 19:31:46:
你还是没有给全，呵呵，不过不要紧，你的问题在于element不能那么读, read(12,*) 这种只对数字有用，所以，关于字符串，只能人工解析。

下面是一个实现，我没有在你的那个程序上改，你自己看看能不能理解… ...

Originally posted by snoopyzhao at 2011-05-08 17:37:30:
你最好还是把运行程序所有的代码+数据文件打包放在附件里吧，要不然，别人也没有办法测试啊，仅仅靠看的功力还是不够的……

Originally posted by Gina88 at 2011-05-09 00:06:17:
这个软件的目的就是其它部分保持CHGCAR1的不变，
就是想实现384  160   24以下的那部分数据和CHGCAR2相应部分相减。

在元素那行，元素的个数不定，
元素下一行表示相应元素原子的个数。

我的程序按照您 ...

integer,allocatable::atomnum(:)

do i=1,100
  atomnum(i)=0
end do

Originally posted by Gina88 at 2011-05-09 00:06:17:
这个软件的目的就是其它部分保持CHGCAR1的不变，
就是想实现384  160   24以下的那部分数据和CHGCAR2相应部分相减。

在元素那行，元素的个数不定，
元素下一行表示相应元素原子的个数。

我的程序按照您 ...

Originally posted by snoopyzhao at 2011-05-09 06:54:14:
384 160 24 以下那部分数据指的是啥？包括 augmentation occupancies 部分的数据吗？

给一套完整的测试数据，给一个明确的问题的描述。这不是一个很难的问题，应该就很快能够解决了……

Originally posted by Gina88 at 2011-05-09 09:51:06:
en，相减的数据不包含augmentation occupancies的那部分，
相减的是384 160 24以下，augmentation occupancies以上的那部分数据。
即CHGCAR1=CHGCAR1-CHGCAR2。

谢谢大侠了！

数据就是EIGENVAL，还有C ...