X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=files%2Fpbs8.csh;h=650f7f9d63e1ff01473b50afac80eb17210db993;hb=521c866952798209eb432cc8f39434b04eeb56bf;hp=ae80fd0a45514f43766524d210476adb92211578;hpb=2fca6b6df0ebd2df200cddb1e51049afa13399dd;p=django_unres.git

diff --git a/files/pbs8.csh b/files/pbs8.csh
index ae80fd0..650f7f9 100755
--- a/files/pbs8.csh
+++ b/files/pbs8.csh
@@ -53,7 +53,7 @@ else
  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
@@ -67,10 +67,57 @@ setenv CLUSTER_WHAM_BIN /users2/czarek/UNRES/git3_tmp/build/bin/cluster_wham-mul
 
 $MPIRUN -machinefile $PBS_NODEFILE -np 2 $CLUSTER_WHAM_BIN
 
+if `awk '{cap=toupper($0); if (cap!=$0) {print 1;exit}}' file.seq` then
+#no allatom conversion/refinementkx for proteins with D-aminoacids
+ ln -s file_wham_T*K_ave.pdb ave
+ sed -n '/ENERGY/,/ENDMDL/p' ave >tmp.pdb
+ setenv numstruc `grep ENERGY tmp.pdb|wc -l`
+ setenv allline `cat tmp.pdb|wc -l`
+ setenv onestruc `echo $allline "/" $numstruc| bc -l|sed 's/\.0*//'`
+ foreach i (2 4 6 8 10)
+   setenv headval `echo $i "*" $onestruc|bc -l`
+   setenv halfi `echo $i/2|bc`
+   head -n $headval tmp.pdb | tail -n $onestruc >MODEL${halfi}.pdb
+ end
+ rm ave tmp.pdb
+else
+# proteins with L-aminoacids only
 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
+endif
+
 /users2/local/pymol_1.6/pymol -c ../files/model.pml
 
+if ( `grep -c pdbref file.inp` ) then
+ 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
+ if (-e saxs.data) then
+  awk -f ../files/saxs_dist.awk plik1.pdb > plik_saxs.data
+ endif
+# rm plik1.pdb
+endif
+
+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