Originally posted by snoopyzhao at 2011-05-10 06:50:40:
部分的程序之外，这部分能看出来是反复读取和写入，可是怎么判断循环结束还是看不懂。另外
160 160 160没有看见读，怎么就能写入了呢。
这个程序真是奇妙啊。
也就高手能写出这样的程序。
要是我自己，打死 ...

Originally posted by Gina88 at 2011-05-10 10:40:20:
snoopyzhao ：其它的我都看明白了，下面一行我上网查了查，trim是去掉字符串前后的空格之后的字符串。
index是看后面的字符串在前面字符串中的位置，如果这个位置不为0，就离开。如果为0就继续。
可是其它行 ...

Originally posted by snoopyzhao at 2011-05-10 10:45:37:
你知道 subline 中的字符串是啥？ '160 160 160'，我把前后的空格都去掉是为了更好的匹配……至于要不要去掉，我还真没有试，我只是觉得这种特征字符串的前后不应该有空格....

Originally posted by Gina88 at 2011-05-10 11:06:29:
谢谢您！
这句话应该不能去掉的，如果去掉就不能实现密度矩阵相减的功能了。
但是我又不明白 if (index (line, trim(subline)) /= 0) exit这句话为什么判断的就是160 160 160那行呢，而不是其它行呢。

16 ...

Originally posted by snoopyzhao at 2011-05-10 11:12:34:
这样吧，

index('    a b c d        ', 'a b c d') 的结果是啥？

index('    a b c d        ', '    a b c d') 的结果是啥？

index('    b b c d        ', 'a b c d') 的结果又是啥？

这跟空格 ...

Originally posted by snoopyzhao at 2011-05-10 11:12:34:
这样吧，

index('    a b c d        ', 'a b c d') 的结果是啥？

index('    a b c d        ', '    a b c d') 的结果是啥？

index('    b b c d        ', 'a b c d') 的结果又是啥？

这跟空格 ...

16楼: Originally posted by snoopyzhao at 2011-05-05 09:19:38:
再贴个完整的程序吧，呵呵……

[code]
program cell_POSCAR

!announcement begin***********************************
implicit none
real(8)::cellcon(3,3)       !lattice constance
real(8),allocat ...