2,684
edits
SSX-PepT Se: Difference between revisions
Jump to navigation
Jump to search
→color the 4 clusters
| Line 37: | Line 37: | ||
color the first cluster: | color the first cluster: | ||
<pre> | |||
grep xtal XSCALE.1.INP |sed s/'INPUT_FILE=..\/tutorial_data\/PepT_SE_145\/xtal_//' |sed 's/.HKL//' |awk '{print $1+1}' > selection_1 | |||
awk 'FNR==NR {a[$1]=0 ;next}{ if($2 in a){printf "%6s%5i%3s%6s%2s%4s%12.3f%8.3f%8.3f%11s\n",$1,$2,"S",$4,$5,$6,$7,$8,$9,$10} \ | |||
else{print $0}}' selection_1 iso.pdb >sel_1.pdb | |||
</pre> | |||
--> do that for the remaining 3 clusters: | --> do that for the remaining 3 clusters: | ||
<pre> | |||
for k in 2 3 4 | |||
do | |||
grep xtal XSCALE.$k.INP |sed s/'INPUT_FILE=..\/tutorial_data\/PepT_SE_145\/xtal_//' |sed 's/.HKL//' |awk '{print $1+1}' > selection_$k | |||
awk 'FNR==NR {a[$1]=0 ;next}{ if($2 in a){printf "%6s%5i%3s%6s%2s%4s%12.3f%8.3f%8.3f%11s\n",$1,$2,"S",$4,$5,$6,$7,$8,$9,$10} \ | |||
else{print $0}}' selection_$k iso.pdb >sel_$k.pdb | else{print $0}}' selection_$k iso.pdb >sel_$k.pdb | ||
done | |||
</pre> | |||
== look at the 4 clusters == | == look at the 4 clusters == | ||