shell 脚本求助 ，循环处理大量文件

yestreenstars

回复 8# hanhuanhuan

你能给个具体的例子来说明一下么？我看了你的主题，又看了你的回复，感觉很矛盾，比如，你主题说：

先是删除前九行，然后对每个文件的第二列进行赋值，接着提取文件前四列

你说你要提取的是文件的前四列，但你在10楼的脚本中却是这样写的：

print $1 ,$3 ,$4 ,$5 ,$6="1.27"

你在8楼的回复中是这样写的：

比如原文件是0001.lammmpstrj，处理后输出成0001.txt文件

我的理解是不修改源文件，把修改后的结果重定向到前缀相同的txt文件，但从你10楼的脚本来看你却修改了源文件：

sed -i '1,9d' $i

看得我头都晕了~

我建议你按照我下面的格式举个例子：

处理前的.lammpstrj文件的内容：
XXX

处理后的.lammpstrj文件的内容：
XXX

根据.lammpstrj得到的.txt文件的内容：
XXX
   

hanhuanhuan

回复 11# yestreenstars
                                                     
谢谢你的回复。我就是想把处理前的文件里前九行描述性的去掉（我用sed在原文件上删除了 ），然后再根据文件中第二列的值（如果是1，$6赋值为1.27，否则赋值1.58），其实我就是提取了原子id以及原子的坐标，加上$6的赋值那一列。

   

hanhuanhuan

回复 12# hanhuanhuan


   

yestreenstars

回复 12# hanhuanhuan

你直接把文本贴到回复里吧，不要用附件，就算用你也别用PDF啊，用TXT不行吗？
   

yxl_ht

三楼回复的多好啊，循环一加，判断加上就能用了啊

hanhuanhuan

回复 14# yestreenstars


处理前：
   
ITEM: TIMESTEP
1236000
ITEM: NUMBER OF ATOMS
16000
ITEM: BOX BOUNDS pp pp ss
-6.7918 6.7918
-6.94724 6.94724
-732.802 734.56
ITEM: ATOMS id type x y z c_5 c_6
176 1 -5.82881 -6.53303 -731.446 -2.42114 0.136341
1106 1 -5.97314 -3.6675 -698.105 -2.79397 0.0938119
1564 1 -2.26312 -4.02711 -731.233 -2.34784 0.311021
634 1 -3.14313 -2.23052 -729.546 -2.88232 0.0540038
982 1 -0.694826 -6.39106 -727.128 -3.27515 0.356942
599 1 -3.6459 -5.92627 -729.823 -2.86464 0.132928
2917 2 -2.5816 -4.36907 -727.462 -6.02845 0.297569
1783 2 -5.66376 -4.47462 -728.803 -5.59957 0.101155
2572 1 -1.59191 -2.13556 -727.987 -2.81018 0.10121
8704 1 -6.19179 -0.552736 -369.505 -2.92728 0.151246
1973 1 -2.6704 -0.0161747 -707.155 -3.16277 0.177424
2434 1 -0.351478 -3.02733 -369.567 -2.8353 0.0386969
1031 2 -0.435456 -1.80943 -730.288 -5.32585 0.162662
2389 1 -3.12658 -0.0973931 -728.272 -2.93912 0.251039

处理后：
176  -5.82881 -6.53303 -731.446  1.27
1106  -5.97314 -3.6675 -698.105 1.27
1564  -2.26312 -4.02711 -731.233  1.27
634  -3.14313 -2.23052 -729.546  1.27
982  -0.694826 -6.39106 -727.128 1.27
599  -3.6459 -5.92627 -729.823  1.27
2917  -2.5816 -4.36907 -727.462  1.58
1783  -5.66376 -4.47462 -728.803  1.58
2572  -1.59191 -2.13556 -727.987 1.27
8704  -6.19179 -0.552736 -369.505 1.27
1973  -2.6704 -0.0161747 -707.155  1.27
2434  -0.351478 -3.02733 -369.567  1.27
1031  -0.435456 -1.80943 -730.288  1.58
2389  -3.12658 -0.0973931 -728.272  1.27

q1208c

回复 10# hanhuanhuan


    #!/bin/bash
for i in ` find file -type f -name "*.lammpstrj" `
do
a=`sed -i '1,9d' $i`
cat $a | while read line
  do
echo $line | awk '{if ($2=="1" print $1 ,$3 ,$4 ,$5 ,$6="1.27" ; else print $1 ,$3 ,$4 ,$5 ,$6="158"}'
  done
done


sed -i 的结果是 "" 的, 所以, 你的 $a 就是个 "", 那下面的while 在等你的输入, 因为你用了 read.

当然不会有结果了.

yestreenstars

回复 16# hanhuanhuan

这样不就好多了吗？一对比就知道你的需求了，比你说一大堆有用。

try:

1. #!/bin/bash
   
2. while read file
   
3. do
   
4.         awk 'NR>9{print $1,$3,$4,$5,$2==1?"1.27":"1.58"}' "$file" > "${file%.*}.txt"
   
5. done <"$(find test -type f -name '*.lammpstrj')"

复制代码

hanhuanhuan

回复 18# yestreenstars
恩恩，是的，第一次在这里发帖求助，不懂规矩， 。谢谢你啦


   

hanhuanhuan

回复 17# q1208c
恩恩，查了下用法，是有问题。谢谢你，