X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=files%2Fpbs8_new.csh;h=0c87f80aa63a47597032223a68dc267ec8674279;hb=8dce3b49dabcc85baae2b796ae055da8866424bb;hp=6d9d029e2c385a1a276a1acbe9c08a163b7fd0bc;hpb=81afcbe06c2b06f623c30db8f3e83c7fca106b97;p=django_unres.git

diff --git a/files/pbs8_new.csh b/files/pbs8_new.csh
index 6d9d029..0c87f80 100755
--- a/files/pbs8_new.csh
+++ b/files/pbs8_new.csh
@@ -67,6 +67,21 @@ 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
@@ -75,11 +90,26 @@ rm ave
 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
+ 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
+  /users2/local/mmtsb/perl/convpdb.pl -setchain A -renumber 1 -out generic_noh plik.pdb |grep ATOM > plik1chain.pdb
+  /users2/local/mmtsb/perl/convpdb.pl -setchain B -renumber 1 -out generic_noh plik2.pdb |grep ATOM >> plik1chain.pdb  
+  foreach i (`seq 1 5`)
+    sed /TER/q MODEL$i.pdb_ > tmp1.pdb
+    /users2/local/mmtsb/perl/convpdb.pl -setchain A -renumber 1 -out generic_noh tmp1.pdb |grep ATOM > MODEL${i}chain.pdb    
+    sed '0,/TER/d' MODEL$i.pdb_ > tmp2.pdb
+    /users2/local/mmtsb/perl/convpdb.pl -setchain B -renumber 1 -out generic_noh tmp2.pdb |grep ATOM >> MODEL${i}chain.pdb    
+    rm tmp1.pdb tmp2.pdb MODEL$i.pdb_
+    ../files/DockQ.py MODEL${i}chain.pdb plik1chain.pdb > dockq_$i.out
+  end
+
+ endif
  /users2/local/mmtsb/perl/convpdb.pl -renumber 1 -out generic tmp.pdb > plik1.pdb
  rm tmp.pdb