提取轨迹文件中z值一直大于10的原子

lionsingasong

如图，想对一个轨迹文件进行处理，该轨迹共包含795个原子，共1000步；
现想提取出在1000步中Z坐标一直大于10的原子。
求问如何写脚本。

jimmy14k

输出第四列的值大于10的行??

ly5066113

目测是这个意思：

1. awk '$NF>10{a[$1]++}END{for(i in a)if(a[i]==1000)print i}' file

复制代码

lionsingasong

回复 2# jimmy14k


    嗯，z值就是第四个值。
   

lionsingasong

回复 2# jimmy14k

唔，可能是我没有表示清楚，我想要的是在1000次计算中Z的值都大于10的，并不是只要1次大于10。


   

jimmy14k

awk '$4>10' file

lionsingasong

回复 3# ly5066113


    唔，可能是我没有表示清楚，我想要的是在1000次计算中Z的值都大于10的，并不是只要1次大于10。
    Z值是一行坐标的第4个值。

jimmy14k

回复 3# ly5066113

3楼竟然目测对了???????


   

lionsingasong

回复 8# jimmy14k
不对，他的那个语句运行出来 没有值
他的语句我没看懂
ps：我截图里面的ob虽然是一样的 但是是区分开的 可以理解为
a 1 2 3
b 1 2 3
c 1 2 0
。。。。

a 1 2 4
b 1 2 2
c 1 2 1
。。。。

这样重复1000次，比如说我要找出在1000次中第四个值一直大于2.5的那行

   

lionsingasong

回复 3# ly5066113


    你的语句我试了，没有输出我想要的结果
ps：我截图里面的ob虽然是一样的 但是是区分开的 可以理解为
795
VMD
a 1 2 3
b 1 2 3
c 1 2 0
。。。。
795
VMD
a 1 2 4
b 1 2 2
c 1 2 1
。。。。

这样重复1000次，比如说我要找出在1000次中第四个值一直大于2.5的那行