File tree Expand file tree Collapse file tree 4 files changed +19
-10
lines changed Expand file tree Collapse file tree 4 files changed +19
-10
lines changed Original file line number Diff line number Diff line change 24
24
do
25
25
printf " >>> Multiplying vectors by $factor in $f ..."
26
26
if [ -e $f ]; then
27
- awk -v fac=$factor ' {print $1 " " $2 " " $3 " " fac*$4 " " fac*$5 " " fac*$6}' $f >> /tmp/temp.dat
28
- mv -f /tmp/temp.dat $f
27
+ # Generate random file name
28
+ of=" $( mktemp ./XXXXXXXX.tmp) "
29
+ awk -v fac=$factor ' {print $1 " " $2 " " $3 " " fac*$4 " " fac*$5 " " fac*$6}' $f >> ${of}
30
+ mv -f ${of} $f
29
31
printf " Done.\n"
30
32
else
31
33
printf " File doesn't exist, skipping.\n"
Original file line number Diff line number Diff line change 17
17
do
18
18
printf " >>> Removing lines with zero vectors in $f ..."
19
19
if [ -e $f ]; then
20
- awk ' {if($4!=0.0 && $5!=0.0 && $6!=0.0) print $0 }' $f >> /tmp/temp.dat
21
- mv -f /tmp/temp.dat $f
20
+ # Generate random file name
21
+ of=" $( mktemp ./XXXXXXXX.tmp) "
22
+ awk ' {if($4!=0.0 && $5!=0.0 && $6!=0.0) print $0 }' $f >> ${of}
23
+ mv -f ${of} $f
22
24
printf " Done.\n"
23
25
else
24
26
printf " File doesn't exist, skipping.\n"
Original file line number Diff line number Diff line change 23
23
# Check that it is a LAMMPS file
24
24
islmp=$( grep " atom types" $f | wc -l)
25
25
if [ $islmp = 1 ]; then
26
+ # Generate random file name
27
+ of=" $( mktemp ./XXXXXXXX.tmp) "
26
28
# If the keyword "Atoms" is not followed by "# charge", add it
27
29
atomic=$( grep " Atoms" $f | grep " charge" | wc -l)
28
30
if [ $atomic = 0 ]; then
29
31
sed -i ' /Atoms/ c\Atoms # charge' $f
30
32
fi
31
- # Add a column of zero after column 2
32
- awk ' {if(NR>9 && NF>4 && $1>0 && $2>0)
33
+ nrAtom=$( awk ' /^Atom/{print NR}' $f )
34
+ # Add a column of zero after keyword "Atom"
35
+ awk -v n=" ${nrAtom} " ' {if(NR>n && NF>4 && $1>0 && $2>0)
33
36
{print $1 "\t" $2 "\t0.0\t" $(NF-2) "\t" $(NF-1) "\t" $NF}
34
37
else
35
38
{print}
36
- }' $f > /tmp/temp.lmp
37
- mv -f /tmp/temp.lmp $f
39
+ }' $f > ${of}
40
+ mv -f ${of} $f
38
41
printf " Done.\n"
39
42
else
40
43
printf " Not a LAMMPS data file, skipping.\n"
Original file line number Diff line number Diff line change 22
22
# Check that it is a LAMMPS file
23
23
islmp=$( grep " atom types" $f | wc -l)
24
24
if [ $islmp = 1 ]; then
25
+ # Generate random file name
26
+ of=" $( mktemp ./XXXXXXXX.tmp) "
25
27
# If the keyword "Atoms" is not followed by "# atomic" add it
26
28
charge=$( grep " Atoms" $f | grep " atomic" | wc -l)
27
29
if [ $charge = 0 ]; then
32
34
{print $1 "\t" $2 "\t" $(NF-2) "\t" $(NF-1) "\t" $NF}
33
35
else
34
36
{print}
35
- }' $f > /tmp/temp.lmp
36
- mv -f /tmp/temp.lmp $f
37
+ }' $f > ${of}
38
+ mv -f ${of} $f
37
39
printf " Done.\n"
38
40
else
39
41
printf " Not a LAMMPS data file, skipping.\n"
You can’t perform that action at this time.
0 commit comments