-setenv UNRES_BIN /users2/czarek/UNRES/run/ADAM/unres-mult-symetr_KCC_ifort_MPICH-NEWCORR-SAXS-NMRAMB-Bfac.exe
+#ORG: setenv UNRES_BIN /users2/czarek/UNRES/run/ADAM/unres-mult-symetr_KCC_ifort_MPICH-NEWCORR-SAXS-NMRAMB-Bfac.exe
+#setenv UNRES_BIN /users2/adam/bin/unres-ms_ifort_MPICH_SC-HCD5-DFA-nmr.exe
+#DZIALA: setenv UNRES_BIN /users2/adam/bin/unres-ms_ifort_MPICH_SC-HCD5_nmr.exe
+setenv UNRES_BIN /users2/adam/bin/unres-ms_ifort_MPICH_SC-HCD5_nmr-SAXS.exe
setenv THETPARPDB $DD/thetaml_ext.5parm
setenv ROTPARPDB $DD/scgauss_ext.parm
setenv ROTPAR $DD/rotamers_AM1_aura_ext.10022007.parm
setenv THETPARPDB $DD/thetaml_ext.5parm
setenv ROTPARPDB $DD/scgauss_ext.parm
setenv ROTPAR $DD/rotamers_AM1_aura_ext.10022007.parm
setenv TORDPAR $DD/pot_tord_G631_DIL_ext.parm
setenv ELEPAR $DD/electr_631Gdp_ext.parm
setenv SIDEPAR $DD/scinter_GB_ext_lip.parm
setenv TORDPAR $DD/pot_tord_G631_DIL_ext.parm
setenv ELEPAR $DD/electr_631Gdp_ext.parm
setenv SIDEPAR $DD/scinter_GB_ext_lip.parm
setenv SCPPAR $DD/scp_ext.parm
setenv PATTERN $DD/patterns.cart
setenv LIPTRANPAR $DD/Lip_tran_initial_ext.parm
setenv SCPPAR $DD/scp_ext.parm
setenv PATTERN $DD/patterns.cart
setenv LIPTRANPAR $DD/Lip_tran_initial_ext.parm
-setenv WHAM_BIN /users2/czarek/UNRES/run/ADAM/wham_ifort_KCC_MPICH-okeanos_NEWCORR-SAXS-NMRAMB-Bfac.exe
+#ORG: setenv WHAM_BIN /users2/czarek/UNRES/run/ADAM/wham_ifort_KCC_MPICH-okeanos_NEWCORR-SAXS-NMRAMB-Bfac.exe
+#FAIL: setenv WHAM_BIN /users2/adam/bin/wham_ifort_MPICH-SC-HCD5-nmr-SAXS.exe
+#DZIALA: setenv WHAM_BIN /users2/adam/unres/bin/wham/wham_ifort_MPICH-SC-HCD5-NMR.exe
+setenv WHAM_BIN /users2/adam/bin/wham_ifort_MPICH-SC-HCD5-nmr-SAXS.exe
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*//'`
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*//'`
setenv headval `echo $i "*" $onestruc|bc -l`
setenv halfi `echo $i/2|bc`
head -n $headval tmp.pdb | tail -n $onestruc >MODEL${halfi}.pdb
setenv headval `echo $i "*" $onestruc|bc -l`
setenv halfi `echo $i/2|bc`
head -n $headval tmp.pdb | tail -n $onestruc >MODEL${halfi}.pdb
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
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
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
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/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
+ foreach i (`seq 1 $numstruc`)
+ /users2/local/bin/tmscore MODEL${i}.pdb plik1.pdb > tmscore${i}.out
+ end
-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
+foreach i (`seq 1 $numstruc`)
+awk -f ../files/saxs_dist.awk MODEL${i}.pdb > MODEL${i}_saxs.data
+end
+#NMR:
+if (-e file_wham_T*K_0001.nmr) then
+ setenv LD_LIBRARY_PATH `pwd`/../files/gnuplot:$LD_LIBRARY_PATH
+ foreach i (`seq -f "%04g" 1 $numstruc`)
+ set pl_cur=`ls -1 file_wham_T*K_${i}.nmr`
+ ln -s $pl_cur _a.nmr
+ ../files/gnuplot/gnuplot ../files/nmr_cont.gnu
+ rm _a.nmr
+ mv _a.png `basename $pl_cur nmr`png
+ end
+endif
+
+