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

yestreenstars

猜的：

1. awk '{if($4>10)a[$1]++;b[$1]++}END{for(i in a)if(a[i]==b[i])print i}'

复制代码

Herowinter

回复 11# yestreenstars
目测就这个意思，楼主没说怎么区分不同的原子，
可能是$1，也可能是
795
Generated by VMD这种上面的那个数字表示795表示
不同的原子id。
   

ly5066113

回复 12# Herowinter


795 是总数，数据9楼有例子。

楼主之所以得不到结果，推测最大的可能是需要：

把 $NF 改成 $4
dos2unix file

Herowinter

回复 13# ly5066113
哦，那我们理解错了，我和星辰大大的理解好像
都是1楼那图整段表示第795个原子，原来是每行
表示一个原子，然后这些行会变化1000次，要
找出这些变化中$4一直大于10的行。


   

ly5066113

回复 14# Herowinter


星辰没理解错，只不过他猜测楼主的数据里面没有重复1000次，所以我的代码没结果。