tail -q -n +$line file_GB*.stat | awk '{if ( NF == 10 ) print}' >remd_all.stat
cat file_GB*.stat |awk '{if ( NF == 10 ) print}' > remd_all0.stat
endif
-../../files/matplotlib_hist.py $temperatures
+../files/matplotlib_hist.py $temperatures
rm remd_all.stat remd_all0.stat
#CLUSTER WHAM
ln -s file_wham_T*K_ave.pdb ave
../files/cluster2allatom.sh ave
rm ave
+
+#amber min refinement
+foreach i (`seq 1 5`)
+../files/refine.sh MODEL$i.pdb $ssbond
+end
+
/users2/local/pymol_1.6/pymol -c ../files/model.pml
+awk '{printf "%s%s%s\n",substr($0,0,21)," ",substr($0,23)}' plik.pdb > tmp.pdb
+/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
+/users2/local/bin/tmscore MODEL4.pdb plik1.pdb > tmscore4.out
+/users2/local/bin/tmscore MODEL5.pdb plik1.pdb > tmscore5.out
+rm plik1.pdb
+
+if (-e saxs.data) then
+awk -f ../files/saxs_dist.awk MODEL1.pdb > MODEL1_saxs.data
+awk -f ../files/saxs_dist.awk MODEL2.pdb > MODEL2_saxs.data
+awk -f ../files/saxs_dist.awk MODEL3.pdb > MODEL3_saxs.data
+awk -f ../files/saxs_dist.awk MODEL4.pdb > MODEL4_saxs.data
+awk -f ../files/saxs_dist.awk MODEL5.pdb > MODEL5_saxs.data
+../files/plot_saxs.py
+endif
+
#END
touch finished