我有一个gromac文件从中提取具体细节。文件格式如下。我只需要在每一列中输入值。在
Generated by trjconv : a bunch of waters t= 0.00000
3000
1SOL OW 1 1.5040 2.7580 0.6820
1SOL HW1 2 1.4788 2.7853 0.7702
1SOL HW2 3 1.4640 2.8230 0.6243
2SOL OW 4 1.5210 0.9510 2.2050
2SOL HW1 5 1.5960 0.9780 2.1520
2SOL HW2 6 1.4460 0.9940 2.1640
1000SOL OW 2998 1.5310 1.7952 2.1981
1000SOL HW1 2999 1.4560 1.7375 2.1836
1000SOL HW2 3000 1.6006 1.7369 2.2286
3.12736 3.12736 3.12736
Generated by trjconv : a bunch of waters t= 9000.00000
3000
1SOL OW 1 1.1579 0.4255 2.1329
1SOL HW1 2 1.0743 0.3793 2.1385
1SOL HW2 3 1.1514 0.4772 2.0526
2SOL OW 4 2.2976 2.5403 0.1574
2SOL HW1 5 2.2161 2.4981 0.1847
2SOL HW2 6 2.3118 2.6092 0.2223
1000SOL OW 2998 2.5927 1.0233 0.4695
1000SOL HW1 2999 2.5927 0.9869 0.3810
1000SOL HW2 3000 2.5365 0.9642 0.5197
3.10941 3.10941 3.10941
我试着按照代码来读它。在
^{pr2}$我想提取下面的数据来浮动变量。在
我的想法是把数据分成p[3],p[4],p[6] 但我失败了。 当我读文件时,我想跳过3000行之后的信息。每3000行后有三行给出信息。我也想跳过前两行。gro文件具有特定的格式。在
http://manual.gromacs.org/current/online/gro.html
输出应如下所示。我要用这些值做一些计算。在
第一套
X = [1.5040, 1.4788, 14640,... , 1.6006]
Y = [2.7580, 2.7853, , 1.7369]
Z = [0.6820, 0.7702, , 2.2286]
第二套
X = [1.1579, 1.0743, ,... , 2.5365]
Y = [0.4255, 0.3793, ,... , 0.9642]
Z = [2.1329, 2.1329, ,... , 0.5197]
请仔细阅读我的问题。我需要一些指导来写上面的程序。我是python的新手。我做过研究,但找不到任何帮助来提取这些数字。在
谢谢你
最后,我能够编写完整的程序来读取我的gromac文件。如果想从grmomac文件中读取数据并从中提取特定值,这可能会有所帮助。此代码调整为读取10行后指定的位置查找。我用这个来读3000行。这是100%成功的工作。在
有人帮我解决具体问题。Thanx也去找他们。这是一个非常好的网站。我喜欢这个。这有助于学习。在
相关问题 更多 >
编程相关推荐