Mercurial > repos > marcel > caddsuite_linux_x86_64

annotate CADDSuite/galaxyconfigs/countResidues.sh @ 0:bac3c274238f

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Migrated tool version 0.93 from old tool shed archive to new tool shed repository
author marcel
date Tue, 07 Jun 2011 16:43:30 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
1 #!/bin/bash
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
2
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
3 IFS=$'\n';
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
4 residue_list=""
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
5 hetsys_list=""
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
6 hoh_list=""
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
7
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
8 for line in $(cat $1); do
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
9 if [ "${line:0:4}" = "ATOM" ];then
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
10 residue_list+="${line:23:6}\n"
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
11 else
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
12 if [ "${line:0:6}" = "HETATM" ];then
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
13 s2=`echo $line | grep HOH`
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
14 if [ "$s2" != "" ];then
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
15 hoh_list+="${line:21:8}\n"
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
16 else
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
17 s3=`echo $line | grep WAT`
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
18 if [ "$s3" != "" ];then
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
19 hoh_list+="${line:21:8}\n"
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
20 else
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
21 hetsys_list+="${line:21:8}\n"
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
22 fi
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
23 fi
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
24 fi
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
25 fi
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
26 done
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
27
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
28 IFS=$' ';
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
29 let num_res=`echo -e "\n"$residue_list | uniq | wc -l`
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
30 let num_hoh=`echo -e "\n"$hoh_list | uniq | wc -l`
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
31 let num_hetsys=`echo -e "\n"$hetsys_list | uniq | wc -l`
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
32
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
33 # empty line at begin and end of lists -> substract 2
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
34 let num_res-=2
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
35 let num_hoh-=2
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
36 let num_hetsys-=2
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
37 if [ $num_res -lt 0 ]; then
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
38 let num_res=0
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
39 fi
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
40 if [ $num_hoh -lt 0 ]; then
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
41 let num_hoh=0
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
42 fi
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
43 if [ $num_hetsys -lt 0 ]; then
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
44 let num_hetsys=0
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
45 fi
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
46
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
47 echo -n "$num_res "
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
48 if [ $num_res -eq 1 ]; then
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
49 echo "residue"
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
50 else
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
51 echo "residues"
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
52 fi
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
53 echo -n "$num_hoh water "
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
54 if [ $num_hoh -eq 1 ]; then
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
55 echo "molecule"
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
56 else
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
57 echo "molecules"
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
58 fi
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
59 echo -n "$num_hetsys "
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
60 if [ $num_hetsys -eq 1 ]; then
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
61 echo "heterosystem"
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
62 else
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
63 echo "heterosystems"
bac3c274238f Migrated tool version 0.93 from old tool shed archive to new tool shed repository
marcel
parents:
diff changeset
64 fi