2,684
edits
SSX-PepT Se: Difference between revisions
Jump to navigation
Jump to search
→color the 4 clusters
| Line 38: | Line 38: | ||
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} \ | 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 | else{print $0}}' selection_1 iso.pdb >sel_1.pdb | ||
--> do that for all 4 clusters | --> do that for all 4 clusters: | ||
for k in 2 3 4 | |||
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 | |||
do | |||
done | |||
== look at the 4 clusters == | == look at the 4 clusters == | ||