新手求教脚本问题

wiliiwin

貌似是这样吧

1. awk  '/^HETATM/{a[$3]++;if(a[$3]>=2){$3=$3(a[$3]-1)}else{$3=$3}}1'  aa-51612432042.TXT
   
2. REMARK  99                                                                     
   
3. REMARK  99 MOE v2013.08 (Chemical Computing Group Inc.) Thu Sep 18 16:38:41 2014
   
4. HETATM 1 C <1> 0 36.477 2.216 13.201 1.00 0.00 C
   
5. HETATM 2 C1 <1> 0 36.077 0.948 12.746 1.00 0.00 C
   
6. HETATM 3 S <1> 0 38.174 2.324 13.429 1.00 0.00 S
   
7. HETATM 4 C2 <1> 0 35.558 3.261 13.445 1.00 0.00 C
   
8. HETATM 5 C3 <1> 0 37.240 0.083 12.617 1.00 0.00 C
   
9. HETATM 6 C4 <1> 0 34.675 0.791 12.541 1.00 0.00 C
   
10. HETATM 7 C5 <1> 0 38.426 0.715 12.946 1.00 0.00 C
    
11. HETATM 8 C6 <1> 0 34.199 3.044 13.238 1.00 0.00 C
    
12. HETATM 9 CL <1> 0 37.159 -1.587 12.190 1.00 0.00 CL
    
13. HETATM 10 C7 <1> 0 33.755 1.816 12.787 1.00 0.00 C
    
14. HETATM 11 CL1 <1> 0 34.082 -0.710 11.965 1.00 0.00 CL
    
15. HETATM 12 C8 <1> 0 39.779 0.209 13.005 1.00 0.00 C
    
16. HETATM 13 S1 <1> 0 40.344 -0.472 11.483 1.00 0.00 S
    
17. HETATM 14 O <1> 0 40.492 0.241 14.001 1.00 0.00 O
    
18. HETATM 15 C9 <1> 0 41.323 -1.846 12.123 1.00 0.00 C
    
19. HETATM 16 C10 <1> 0 41.426 -2.888 11.045 1.00 0.00 C
    
20. HETATM 17 O1 <1> 0 40.536 -3.894 11.277 1.00 0.00 O
    
21. HETATM 18 O2 <1> 0 42.176 -2.811 10.081 1.00 0.00 O
    
22. HETATM 19 C11 <1> 0 40.421 -4.889 10.247 1.00 0.00 C
    
23. HETATM 20 C12 <1> 0 39.413 -4.422 9.218 1.00 0.00 C
    
24. HETATM 21 H <1> 0 35.906 4.230 13.792 1.00 0.00 H
    
25. HETATM 22 H1 <1> 0 33.489 3.845 13.425 1.00 0.00 H
    
26. HETATM 23 H2 <1> 0 32.692 1.659 12.622 1.00 0.00 H
    
27. HETATM 24 H1 <1> 0 42.325 -1.503 12.393 1.00 0.00 H
    
28. HETATM 25 H2 <1> 0 40.865 -2.285 13.016 1.00 0.00 H
    
29. HETATM 26 H11 <1> 0 41.392 -5.104 9.786 1.00 0.00 H
    
30. HETATM 27 H21 <1> 0 40.078 -5.813 10.719 1.00 0.00 H
    
31. HETATM 28 H12 <1> 0 39.756 -3.506 8.727 1.00 0.00 H
    
32. HETATM 29 H22 <1> 0 38.458 -4.179 9.695 1.00 0.00 H
    
33. HETATM 30 H3 <1> 0 39.250 -5.190 8.459 1.00 0.00 H
    
34. CONECT    1    2    2    3    4
    
35. CONECT    2    1    1    5    6
    
36. CONECT    3    1    7
    
37. CONECT    4    1   21    8    8
    
38. CONECT    5    2    7    7    9
    
39. CONECT    6    2   10   10   11
    
40. CONECT    7    3    5    5   12
    
41. CONECT    8    4    4   22   10
    
42. CONECT    9    5
    
43. CONECT   10    6    6    8   23
    
44. CONECT   11    6
    
45. CONECT   12    7   13   14   14
    
46. CONECT   13   12   15
    
47. CONECT   14   12   12
    
48. CONECT   15   13   24   25   16
    
49. CONECT   16   15   17   18   18
    
50. CONECT   17   16   19
    
51. CONECT   18   16   16
    
52. CONECT   19   17   26   27   20
    
53. CONECT   20   19   28   29   30
    
54. CONECT   21    4
    
55. CONECT   22    8
    
56. CONECT   23   10
    
57. CONECT   24   15
    
58. CONECT   25   15
    
59. CONECT   26   19
    
60. CONECT   27   19
    
61. CONECT   28   20
    
62. CONECT   29   20
    
63. CONECT   30   20
    
64. END

复制代码

wiliiwin

husong0811 发表于 2014-10-21 17:04
回复 8# wiliiwin

不是mod那一列看有没有重复的 重复的后面加上数字吗  其他我也没有看到区别啊

wiliiwin

回复 17# husong0811
是在原来的列啊 没有变啊 只是在mod列加了数字而已   <1>的列没有变化吧