/users2/local/mmtsb/perl/convpdb.pl -out generic_noh z_tleap.pdb > noh.pdb
/users2/local/pymol_1.6/pymol -cq /users2/CASP12/EXAMPLES/all-atom/amber_files/dss.py | awk \
-'{ires[NR]=$2;ss[NR]=$3}END{for (i=2;i<=NR;i++) {if (ss[i-1]!=ss[i]&&ss[i-1]=="L") printf "%s;",ires[i];if (ss[i-1]!=ss[i]&&ss[i]=="L") printf "%s;%s\n",ires[i-1],ss[i-1]}}' \
+'{ires[NR]=$2;ss[NR]=$3}END{for (i=2;i<=NR;i++) {if (ss[i-1]!=ss[i]&&ss[i-1]!="L") printf "%s;%s\n",ires[i-1],ss[i-1];if (ss[i-1]!=ss[i]&&ss[i]!="L") printf "%s;",ires[i];}}' \
> zakresy.csv
$AMBERHOME/bin/ambpdb -p t.top <restrt> stop.pdb
rm $pdb
+cp stop.pdb ${pdb}_
/users2/local/mmtsb/perl/convpdb.pl -amber -out generic_noh stop.pdb > $pdb
rm tmp1.pdb z_tleap.pdb t.top t.crd leap.log noh.pdb zakresy.csv RST.dat restrt mdinfo stop.pdb