怎样读取这个索要列的文件的数据？

popcorn

文件格式如下：
Tracer particle file
   26   16
       x       y       z      vx      vy      vz
      sx      sy      sz     sxy     syz     szx
     efp     rho    rvol  active
  0.51806E-01
  0.00000E+00  0.00000E+00  0.50000E+01  0.00000E+00  0.00000E+00  0.00000E+00
-0.11880E+00 -0.11880E+00 -0.11880E+00  0.00000E+00  0.00000E+00  0.00000E+00
  0.10000E+01  0.17170E+01  0.00000E+00  0.39781E+05
  0.00000E+00  0.00000E+00  0.60000E+01  0.00000E+00  0.00000E+00  0.00000E+00
-0.11880E+00 -0.11880E+00 -0.11880E+00  0.00000E+00  0.00000E+00  0.00000E+00
  0.10000E+01  0.17170E+01  0.00000E+00  0.41901E+05
  0.00000E+00  0.00000E+00  0.70000E+01  0.00000E+00  0.00000E+00 -0.80578E-08
-0.11880E+00 -0.11880E+00 -0.11880E+00  0.00000E+00  0.00000E+00  0.00000E+00
  0.10000E+01  0.17170E+01  0.00000E+00  0.36869E+05
0.10014E+01
  0.00000E+00  0.00000E+00  0.50000E+01  0.00000E+00  0.00000E+00  0.29901E-02
-0.40473E-01 -0.40473E-01 -0.40473E-01  0.00000E+00  0.00000E+00  0.00000E+00
  0.10000E+01  0.17170E+01  0.00000E+00  0.39781E+05
  0.00000E+00  0.00000E+00  0.60000E+01  0.00000E+00  0.00000E+00  0.15986E-02
-0.43063E-01 -0.43063E-01 -0.43063E-01  0.00000E+00  0.00000E+00  0.00000E+00
  0.10000E+01  0.17170E+01  0.00000E+00  0.41901E+05
  0.00000E+00  0.00000E+00  0.70000E+01  0.00000E+00  0.00000E+00 -0.15985E-02
-0.43063E-01 -0.43063E-01 -0.43063E-01  0.00000E+00  0.00000E+00  0.00000E+00
  0.10000E+01  0.17170E+01  0.00000E+00  0.36869E+05
·····················
我要读取的是SZ列的数据存到矩阵里，因为文件数据不是很规则，有些行只有一列数据！望大家帮忙，等着作毕业论文！谢谢各位高手呢！

苍山负雪

一段python代码

2. #! /usr/bin/env python
3. import os,sys
5. ifile = sys.argv[1]; ofile = sys.argv[2]
6. finp = open(ifile,'r'); fout = open(ofile,'w')
7. lines = finp.readlines()
9. for line in lines:
10.     if len(line.split()) > 2:
11.     foo = line.split()[2]
12.     fout.write(foo + '\n')
13.   
14. finp.close(); fout.close()

复制代码

2. d:\test.py inputfile outputfile
4. d:\type outputfile
5. 0.50000E+01
6. -0.11880E+00
7. 0.00000E+00
8. 0.60000E+01
9. -0.11880E+00
10. 0.00000E+00
11. 0.70000E+01
12. -0.11880E+00
13. 0.00000E+00
14. 0.50000E+01
15. -0.40473E-01
16. 0.00000E+00
17. 0.60000E+01
18. -0.43063E-01
19. 0.00000E+00
20. 0.70000E+01
21. -0.43063E-01
22. 0.00000E+00

复制代码

popcorn

谢谢呢！
我试一下。

violetwind

原来苍山负雪是舞蛇高手啊,
真是失敬

苍山负雪

/blush 我只是初学者，献丑献丑。violetwind 多多指教

quartzbj

能不能解释一下？
直接把这段程序放到fortran命令中吗？
楼主试验后怎么样？

苍山负雪 wrote:
一段python代码

2. #! /usr/bin/env python
3. import os,sys
5. ifile = sys.argv[1]; ofile = sys.argv[2]
6. finp = open(ifile,'r'); fout = open(ofile,'w')
7. lines = finp.readlines()
9. for line in lines:
10.     if len(line.split()) > 2:
11.     foo = line.split()[2]
12.     fout.write(foo + '\n')
13.   
14. finp.close(); fout.close()

复制代码

2. d:\test.py inputfile outputfile
4. d:\type outputfile
5. 0.50000E+01
6. -0.11880E+00
7. 0.00000E+00
8. 0.60000E+01
9. -0.11880E+00
10. 0.00000E+00
11. 0.70000E+01
12. -0.11880E+00
13. 0.00000E+00
14. 0.50000E+01
15. -0.40473E-01
16. 0.00000E+00
17. 0.60000E+01
18. -0.43063E-01
19. 0.00000E+00
20. 0.70000E+01
21. -0.43063E-01
22. 0.00000E+00

复制代码

苍山负雪

关于python请搜索一下本版，以前贴过一篇介绍。

不过不能直接在fortran里面用啦。需要安装python的解释器。

让python读入文件后，每一行数据表示成一个字符串。把字符串用空格劈开，形成多个字符串，如果多于两个字符串，那就保留第三个字符串，很简单的。

quartzbj

有时间看来要学一下,好像比用fortran 提取数据简单些吧.