if ( `grep -c pdbref file.inp` ) then
awk '{printf "%s%s%s\n",substr($0,0,21)," ",substr($0,23)}' plik.pdb > tmp.pdb
+ if ( -f "plik1ter.pdb" ) then
+ grep -v TER plik1ter.pdb | awk '{printf "%s%s%s\n",substr($0,0,21)," ",substr($0,23)}' > tmp.pdb
+ endif
/users2/local/mmtsb/perl/convpdb.pl -renumber 1 -out generic tmp.pdb > plik1.pdb
rm tmp.pdb
+
/users2/local/bin/tmscore MODEL1.pdb plik1.pdb > tmscore1.out
/users2/local/bin/tmscore MODEL2.pdb plik1.pdb > tmscore2.out
/users2/local/bin/tmscore MODEL3.pdb plik1.pdb > tmscore3.out