annotate mayachemtools/docs/modules/html/code/SDFileUtil.html @ 9:ab29fa5c8c1f draft default tip

Uploaded
author deepakjadmin
date Thu, 15 Dec 2016 14:18:03 -0500
parents 73ae111cf86f
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1 <html>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
2 <head>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
3 <title>MayaChemTools:Code:SDFileUtil.pm</title>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
4 <meta http-equiv="content-type" content="text/html;charset=utf-8">
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
5 <link rel="stylesheet" type="text/css" href="../../../css/MayaChemToolsCode.css">
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
6 </head>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
7 <body leftmargin="20" rightmargin="20" topmargin="10" bottommargin="10">
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
8 <br/>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
9 <center>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
10 <a href="http://www.mayachemtools.org" title="MayaChemTools Home"><img src="../../../images/MayaChemToolsLogo.gif" border="0" alt="MayaChemTools"></a>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
11 </center>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
12 <br/>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
13 <pre>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
14 <a name="package-SDFileUtil-"></a> 1 <span class="k">package </span><span class="i">SDFileUtil</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
15 2 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
16 3 <span class="c"># $RCSfile: SDFileUtil.pm,v $</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
17 4 <span class="c"># $Date: 2015/02/28 20:47:18 $</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
18 5 <span class="c"># $Revision: 1.49 $</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
19 6 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
20 7 <span class="c"># Author: Manish Sud &lt;msud@san.rr.com&gt;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
21 8 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
22 9 <span class="c"># Copyright (C) 2015 Manish Sud. All rights reserved.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
23 10 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
24 11 <span class="c"># This file is part of MayaChemTools.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
25 12 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
26 13 <span class="c"># MayaChemTools is free software; you can redistribute it and/or modify it under</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
27 14 <span class="c"># the terms of the GNU Lesser General Public License as published by the Free</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
28 15 <span class="c"># Software Foundation; either version 3 of the License, or (at your option) any</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
29 16 <span class="c"># later version.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
30 17 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
31 18 <span class="c"># MayaChemTools is distributed in the hope that it will be useful, but without</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
32 19 <span class="c"># any warranty; without even the implied warranty of merchantability of fitness</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
33 20 <span class="c"># for a particular purpose. See the GNU Lesser General Public License for more</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
34 21 <span class="c"># details.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
35 22 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
36 23 <span class="c"># You should have received a copy of the GNU Lesser General Public License</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
37 24 <span class="c"># along with MayaChemTools; if not, see &lt;http://www.gnu.org/licenses/&gt; or</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
38 25 <span class="c"># write to the Free Software Foundation Inc., 59 Temple Place, Suite 330,</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
39 26 <span class="c"># Boston, MA, 02111-1307, USA.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
40 27 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
41 28
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
42 29 <span class="k">use</span> <span class="w">strict</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
43 30 <span class="k">use</span> <span class="w">Exporter</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
44 31 <span class="k">use</span> <span class="w">Carp</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
45 32 <span class="k">use</span> <span class="w">PeriodicTable</span> <span class="q">qw(IsElement)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
46 33 <span class="k">use</span> <span class="w">TimeUtil</span> <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
47 34
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
48 35 <span class="k">use</span> <span class="w">vars</span> <span class="q">qw(@ISA @EXPORT @EXPORT_OK %EXPORT_TAGS)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
49 36
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
50 37 <span class="i">@ISA</span> = <span class="q">qw(Exporter)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
51 38 <span class="i">@EXPORT</span> = <span class="q">qw(GenerateCmpdAtomLine GenerateCmpdBondLine GenerateCmpdChargePropertyLines GenerateCmpdCommentsLine GenerateCmpdCountsLine GenerateCmpdAtomAliasPropertyLines GenerateCmpdIsotopePropertyLines GenerateCmpdDataHeaderLabelsAndValuesLines GenerateCmpdMiscInfoLine GenerateCmpdRadicalPropertyLines GenerateCmpdMolNameLine GenerateEmptyCtabBlockLines GenerateMiscLineDateStamp GetAllAndCommonCmpdDataHeaderLabels GetCmpdDataHeaderLabels GetCmpdDataHeaderLabelsAndValues GetCmpdFragments GetCtabLinesCount GetUnknownAtoms GetInvalidAtomNumbers MDLChargeToInternalCharge InternalChargeToMDLCharge MDLBondTypeToInternalBondOrder InternalBondOrderToMDLBondType MDLBondStereoToInternalBondStereochemistry InternalBondStereochemistryToMDLBondStereo InternalSpinMultiplicityToMDLRadical MDLRadicalToInternalSpinMultiplicity IsCmpd3D IsCmpd2D ParseCmpdAtomLine ParseCmpdBondLine ParseCmpdCommentsLine ParseCmpdCountsLine ParseCmpdMiscInfoLine ParseCmpdMolNameLine ParseCmpdAtomAliasPropertyLine ParseCmpdChargePropertyLine ParseCmpdIsotopePropertyLine ParseCmpdRadicalPropertyLine ReadCmpdString RemoveCmpdDataHeaderLabelAndValue WashCmpd)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
52 39 <span class="i">@EXPORT_OK</span> = <span class="q">qw()</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
53 40 <span class="i">%EXPORT_TAGS</span> = <span class="s">(</span><span class="w">all</span> <span class="cm">=&gt;</span> <span class="s">[</span><span class="i">@EXPORT</span><span class="cm">,</span> <span class="i">@EXPORT_OK</span><span class="s">]</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
54 41
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
55 42 <span class="c"># Format data for compounds count line...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
56 <a name="GenerateCmpdCountsLine-"></a> 43 <span class="k">sub </span><span class="m">GenerateCmpdCountsLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
57 44 <span class="k">my</span><span class="s">(</span><span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="cm">,</span> <span class="i">$ChiralFlag</span><span class="cm">,</span> <span class="i">$PropertyCount</span><span class="cm">,</span> <span class="i">$Version</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
58 45
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
59 46 <span class="k">if</span> <span class="s">(</span><span class="i">@_</span> == <span class="n">5</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
60 47 <span class="s">(</span><span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="cm">,</span> <span class="i">$ChiralFlag</span><span class="cm">,</span> <span class="i">$PropertyCount</span><span class="cm">,</span> <span class="i">$Version</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
61 48 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
62 49 <span class="k">elsif</span> <span class="s">(</span><span class="i">@_</span> == <span class="n">3</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
63 50 <span class="s">(</span><span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="cm">,</span> <span class="i">$ChiralFlag</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
64 51 <span class="i">$PropertyCount</span> = <span class="n">999</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
65 52 <span class="i">$Version</span> = <span class="q">&quot;V2000&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
66 53 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
67 54 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
68 55 <span class="s">(</span><span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
69 56 <span class="i">$ChiralFlag</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
70 57 <span class="i">$PropertyCount</span> = <span class="n">999</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
71 58 <span class="i">$Version</span> = <span class="q">&quot;V2000&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
72 59 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
73 60 <span class="k">if</span> <span class="s">(</span><span class="i">$AtomCount</span> &gt; <span class="n">999</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
74 61 <span class="w">croak</span> <span class="q">&quot;Error: SDFileUtil::GenerateCmpdCountsLine: The atom count, $AtomCount, exceeds maximum of 999 allowed for CTAB version 2000. The Extended Connection Table (V3000) format in MDL MOL and SD files is not supported by the current release of MayaChemTools...&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
75 62 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
76 63 <span class="i">$Line</span> = <span class="k">sprintf</span> <span class="q">&quot;%3i%3i%3i%3i%3i%3i%3i%3i%3i%3i%3i%6s&quot;</span><span class="cm">,</span> <span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="cm">,</span> <span class="n">0</span><span class="cm">,</span> <span class="n">0</span><span class="cm">,</span> <span class="i">$ChiralFlag</span><span class="cm">,</span> <span class="n">0</span><span class="cm">,</span> <span class="n">0</span><span class="cm">,</span> <span class="n">0</span><span class="cm">,</span> <span class="n">0</span><span class="cm">,</span> <span class="n">0</span><span class="cm">,</span> <span class="i">$PropertyCount</span><span class="cm">,</span> <span class="i">$Version</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
77 64
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
78 65 <span class="k">return</span> <span class="s">(</span><span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
79 66 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
80 67
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
81 68 <span class="c"># Generate comments line...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
82 <a name="GenerateCmpdCommentsLine-"></a> 69 <span class="k">sub </span><span class="m">GenerateCmpdCommentsLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
83 70 <span class="k">my</span><span class="s">(</span><span class="i">$Comments</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
84 71 <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
85 72
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
86 73 <span class="i">$Line</span> = <span class="s">(</span><span class="k">length</span><span class="s">(</span><span class="i">$Comments</span><span class="s">)</span> &gt; <span class="n">80</span><span class="s">)</span> ? <span class="k">substr</span><span class="s">(</span><span class="i">$Comments</span><span class="cm">,</span> <span class="n">0</span><span class="cm">,</span> <span class="n">80</span><span class="s">)</span> <span class="co">:</span> <span class="i">$Comments</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
87 74
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
88 75 <span class="k">return</span> <span class="i">$Line</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
89 76 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
90 77
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
91 78 <span class="c"># Generate molname line...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
92 <a name="GenerateCmpdMolNameLine-"></a> 79 <span class="k">sub </span><span class="m">GenerateCmpdMolNameLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
93 80 <span class="k">my</span><span class="s">(</span><span class="i">$MolName</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
94 81 <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
95 82
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
96 83 <span class="i">$Line</span> = <span class="s">(</span><span class="k">length</span><span class="s">(</span><span class="i">$MolName</span><span class="s">)</span> &gt; <span class="n">80</span><span class="s">)</span> ? <span class="k">substr</span><span class="s">(</span><span class="i">$MolName</span><span class="cm">,</span> <span class="n">0</span><span class="cm">,</span> <span class="n">80</span><span class="s">)</span> <span class="co">:</span> <span class="i">$MolName</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
97 84
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
98 85 <span class="k">return</span> <span class="i">$Line</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
99 86 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
100 87
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
101 88 <span class="c"># Generate data for compounds misc info line...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
102 <a name="GenerateCmpdMiscInfoLine-"></a> 89 <span class="k">sub </span><span class="m">GenerateCmpdMiscInfoLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
103 90 <span class="k">my</span><span class="s">(</span><span class="i">$ProgramName</span><span class="cm">,</span> <span class="i">$UserInitial</span><span class="cm">,</span> <span class="i">$Code</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
104 91 <span class="k">my</span><span class="s">(</span><span class="i">$Date</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
105 92
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
106 93 <span class="k">if</span> <span class="s">(</span>!<span class="s">(</span><span class="k">defined</span><span class="s">(</span><span class="i">$ProgramName</span><span class="s">)</span> &amp;&amp; <span class="i">$ProgramName</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
107 94 <span class="i">$ProgramName</span> = <span class="q">&quot;MayaChem&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
108 95 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
109 96 <span class="k">if</span> <span class="s">(</span>!<span class="s">(</span><span class="k">defined</span><span class="s">(</span><span class="i">$UserInitial</span><span class="s">)</span> &amp;&amp; <span class="i">$UserInitial</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
110 97 <span class="i">$UserInitial</span> = <span class="q">&quot; &quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
111 98 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
112 99 <span class="k">if</span> <span class="s">(</span>!<span class="s">(</span><span class="k">defined</span><span class="s">(</span><span class="i">$Code</span><span class="s">)</span> &amp;&amp; <span class="i">$Code</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
113 100 <span class="i">$Code</span> = <span class="q">&quot;2D&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
114 101 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
115 102
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
116 103 <span class="k">if</span> <span class="s">(</span><span class="k">length</span><span class="s">(</span><span class="i">$ProgramName</span><span class="s">)</span> &gt; <span class="n">8</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
117 104 <span class="i">$ProgramName</span> = <span class="k">substr</span><span class="s">(</span><span class="i">$ProgramName</span><span class="cm">,</span> <span class="n">0</span><span class="cm">,</span> <span class="n">8</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
118 105 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
119 106 <span class="k">if</span> <span class="s">(</span><span class="k">length</span><span class="s">(</span><span class="i">$UserInitial</span><span class="s">)</span> &gt; <span class="n">2</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
120 107 <span class="i">$UserInitial</span> = <span class="k">substr</span><span class="s">(</span><span class="i">$UserInitial</span><span class="cm">,</span> <span class="n">0</span><span class="cm">,</span> <span class="n">2</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
121 108 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
122 109 <span class="k">if</span> <span class="s">(</span><span class="k">length</span><span class="s">(</span><span class="i">$Code</span><span class="s">)</span> &gt; <span class="n">2</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
123 110 <span class="i">$Code</span> = <span class="k">substr</span><span class="s">(</span><span class="i">$Code</span><span class="cm">,</span> <span class="n">0</span><span class="cm">,</span> <span class="n">2</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
124 111 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
125 112 <span class="i">$Date</span> = <span class="i">GenerateMiscLineDateStamp</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
126 113
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
127 114 <span class="i">$Line</span> = <span class="q">&quot;${UserInitial}${ProgramName}${Date}${Code}&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
128 115
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
129 116 <span class="k">return</span> <span class="i">$Line</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
130 117 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
131 118
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
132 119 <span class="c"># Generate data for compounds misc info line...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
133 <a name="GenerateEmptyCtabBlockLines-"></a> 120 <span class="k">sub </span><span class="m">GenerateEmptyCtabBlockLines</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
134 121 <span class="k">my</span><span class="s">(</span><span class="i">$Date</span><span class="cm">,</span> <span class="i">$Lines</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
135 122
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
136 123 <span class="k">if</span> <span class="s">(</span><span class="i">@_</span> == <span class="n">1</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
137 124 <span class="s">(</span><span class="i">$Date</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
138 125 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
139 126 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
140 127 <span class="i">$Date</span> = <span class="i">GenerateMiscLineDateStamp</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
141 128 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
142 129 <span class="c"># First line: Blank molname line...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
143 130 <span class="c"># Second line: Misc info...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
144 131 <span class="c"># Third line: Blank comments line...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
145 132 <span class="c"># Fourth line: Counts line reflecting empty structure data block...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
146 133 <span class="i">$Lines</span> = <span class="q">&quot;\n&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
147 134 <span class="i">$Lines</span> .= <span class="q">&quot; MayaChem${Date}2D\n&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
148 135 <span class="i">$Lines</span> .= <span class="q">&quot;\n&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
149 136 <span class="i">$Lines</span> .= <span class="i">GenerateCmpdCountsLine</span><span class="s">(</span><span class="n">0</span><span class="cm">,</span> <span class="n">0</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span> . <span class="q">&quot;\n&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
150 137 <span class="i">$Lines</span> .= <span class="q">&quot;M END&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
151 138
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
152 139 <span class="k">return</span> <span class="i">$Lines</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
153 140 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
154 141
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
155 142 <span class="c"># Generate SD file data stamp...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
156 <a name="GenerateMiscLineDateStamp-"></a> 143 <span class="k">sub </span><span class="m">GenerateMiscLineDateStamp</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
157 144 <span class="k">return</span> <span class="i">TimeUtil::SDFileTimeStamp</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
158 145 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
159 146
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
160 147 <span class="c"># Generate data for compound atom line...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
161 148 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
162 <a name="GenerateCmpdAtomLine-"></a> 149 <span class="k">sub </span><span class="m">GenerateCmpdAtomLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
163 150 <span class="k">my</span><span class="s">(</span><span class="i">$AtomSymbol</span><span class="cm">,</span> <span class="i">$AtomX</span><span class="cm">,</span> <span class="i">$AtomY</span><span class="cm">,</span> <span class="i">$AtomZ</span><span class="cm">,</span> <span class="i">$MassDifference</span><span class="cm">,</span> <span class="i">$Charge</span><span class="cm">,</span> <span class="i">$StereoParity</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
164 151 <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
165 152
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
166 153 <span class="k">if</span> <span class="s">(</span>!<span class="k">defined</span> <span class="i">$MassDifference</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
167 154 <span class="i">$MassDifference</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
168 155 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
169 156 <span class="k">if</span> <span class="s">(</span>!<span class="k">defined</span> <span class="i">$Charge</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
170 157 <span class="i">$Charge</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
171 158 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
172 159 <span class="k">if</span> <span class="s">(</span>!<span class="k">defined</span> <span class="i">$StereoParity</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
173 160 <span class="i">$StereoParity</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
174 161 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
175 162 <span class="i">$Line</span> = <span class="k">sprintf</span> <span class="q">&quot;%10.4f%10.4f%10.4f %-3s%2i%3i%3i 0 0 0 0 0 0 0 0 0&quot;</span><span class="cm">,</span> <span class="i">$AtomX</span><span class="cm">,</span> <span class="i">$AtomY</span><span class="cm">,</span> <span class="i">$AtomZ</span><span class="cm">,</span> <span class="i">$AtomSymbol</span><span class="cm">,</span> <span class="i">$MassDifference</span><span class="cm">,</span> <span class="i">$Charge</span><span class="cm">,</span> <span class="i">$StereoParity</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
176 163
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
177 164 <span class="k">return</span> <span class="i">$Line</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
178 165 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
179 166
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
180 167 <span class="c"># Generate data for compound bond line...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
181 168 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
182 <a name="GenerateCmpdBondLine-"></a> 169 <span class="k">sub </span><span class="m">GenerateCmpdBondLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
183 170 <span class="k">my</span><span class="s">(</span><span class="i">$FirstAtomNum</span><span class="cm">,</span> <span class="i">$SecondAtomNum</span><span class="cm">,</span> <span class="i">$BondType</span><span class="cm">,</span> <span class="i">$BondStereo</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
184 171 <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
185 172
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
186 173 <span class="k">if</span> <span class="s">(</span>!<span class="k">defined</span> <span class="i">$BondStereo</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
187 174 <span class="i">$BondStereo</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
188 175 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
189 176 <span class="i">$Line</span> = <span class="k">sprintf</span> <span class="q">&quot;%3i%3i%3i%3i 0 0 0&quot;</span><span class="cm">,</span> <span class="i">$FirstAtomNum</span><span class="cm">,</span> <span class="i">$SecondAtomNum</span><span class="cm">,</span> <span class="i">$BondType</span><span class="cm">,</span> <span class="i">$BondStereo</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
190 177
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
191 178 <span class="k">return</span> <span class="i">$Line</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
192 179 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
193 180
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
194 181 <span class="c"># Generate charge property lines for CTAB block...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
195 182 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
196 <a name="GenerateCmpdChargePropertyLines-"></a> 183 <span class="k">sub </span><span class="m">GenerateCmpdChargePropertyLines</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
197 184 <span class="k">my</span><span class="s">(</span><span class="i">$ChargeValuePairsRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
198 185
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
199 186 <span class="k">return</span> <span class="i">_GenerateCmpdGenericPropertyLines</span><span class="s">(</span><span class="q">&#39;Charge&#39;</span><span class="cm">,</span> <span class="i">$ChargeValuePairsRef</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
200 187 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
201 188
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
202 189 <span class="c"># Generate isotope property lines for CTAB block...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
203 190 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
204 <a name="GenerateCmpdIsotopePropertyLines-"></a> 191 <span class="k">sub </span><span class="m">GenerateCmpdIsotopePropertyLines</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
205 192 <span class="k">my</span><span class="s">(</span><span class="i">$IsotopeValuePairsRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
206 193
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
207 194 <span class="k">return</span> <span class="i">_GenerateCmpdGenericPropertyLines</span><span class="s">(</span><span class="q">&#39;Isotope&#39;</span><span class="cm">,</span> <span class="i">$IsotopeValuePairsRef</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
208 195 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
209 196
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
210 197 <span class="c"># Generate radical property line property lines for CTAB block...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
211 198 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
212 <a name="GenerateCmpdRadicalPropertyLines-"></a> 199 <span class="k">sub </span><span class="m">GenerateCmpdRadicalPropertyLines</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
213 200 <span class="k">my</span><span class="s">(</span><span class="i">$RadicalValuePairsRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
214 201
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
215 202 <span class="k">return</span> <span class="i">_GenerateCmpdGenericPropertyLines</span><span class="s">(</span><span class="q">&#39;Radical&#39;</span><span class="cm">,</span> <span class="i">$RadicalValuePairsRef</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
216 203 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
217 204
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
218 205 <span class="c"># Generate atom alias property line property lines for CTAB block...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
219 206 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
220 207 <span class="c"># Atom alias property line format:</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
221 208 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
222 209 <span class="c"># A aaa</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
223 210 <span class="c"># x...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
224 211 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
225 212 <span class="c"># aaa: Atom number</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
226 213 <span class="c"># x: Atom alias in next line</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
227 214 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
228 <a name="GenerateCmpdAtomAliasPropertyLines-"></a> 215 <span class="k">sub </span><span class="m">GenerateCmpdAtomAliasPropertyLines</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
229 216 <span class="k">my</span><span class="s">(</span><span class="i">$PropertyValuePairsRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
230 217 <span class="k">my</span><span class="s">(</span><span class="i">$Index</span><span class="cm">,</span> <span class="i">$AtomNum</span><span class="cm">,</span> <span class="i">$AtomAlias</span><span class="cm">,</span> <span class="i">$Line</span><span class="cm">,</span> <span class="i">@PropertyLines</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
231 218
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
232 219 <span class="i">@PropertyLines</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
233 220
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
234 221 <span class="k">for</span> <span class="s">(</span><span class="i">$Index</span> = <span class="n">0</span><span class="sc">;</span> <span class="i">$Index</span> &lt; <span class="i">$#</span>{<span class="i">$PropertyValuePairsRef</span>}<span class="sc">;</span> <span class="i">$Index</span> += <span class="n">2</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
235 222 <span class="i">$AtomNum</span> = <span class="i">$PropertyValuePairsRef</span>-&gt;[<span class="i">$Index</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
236 223 <span class="i">$AtomAlias</span> = <span class="i">$PropertyValuePairsRef</span>-&gt;[<span class="i">$Index</span> + <span class="n">1</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
237 224
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
238 225 <span class="i">$Line</span> = <span class="q">&quot;A &quot;</span> . <span class="k">sprintf</span> <span class="q">&quot;%3i&quot;</span><span class="cm">,</span> <span class="i">$AtomNum</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
239 226
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
240 227 <span class="k">push</span> <span class="i">@PropertyLines</span><span class="cm">,</span> <span class="i">$Line</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
241 228 <span class="k">push</span> <span class="i">@PropertyLines</span><span class="cm">,</span> <span class="i">$AtomAlias</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
242 229 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
243 230
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
244 231 <span class="k">return</span> <span class="i">@PropertyLines</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
245 232 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
246 233
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
247 234 <span class="c"># Generate data header labels and values lines...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
248 235 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
249 <a name="GenerateCmpdDataHeaderLabelsAndValuesLines-"></a> 236 <span class="k">sub </span><span class="m">GenerateCmpdDataHeaderLabelsAndValuesLines</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
250 237 <span class="k">my</span><span class="s">(</span><span class="i">$DataHeaderLabelsRef</span><span class="cm">,</span> <span class="i">$DataHeaderLabelsAndValuesRef</span><span class="cm">,</span> <span class="i">$SortDataLabels</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
251 238 <span class="k">my</span><span class="s">(</span><span class="i">$DataLabel</span><span class="cm">,</span> <span class="i">$DataValue</span><span class="cm">,</span> <span class="i">@DataLabels</span><span class="cm">,</span> <span class="i">@DataLines</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
252 239
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
253 240 <span class="k">if</span> <span class="s">(</span>!<span class="k">defined</span> <span class="i">$SortDataLabels</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
254 241 <span class="i">$SortDataLabels</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
255 242 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
256 243
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
257 244 <span class="i">@DataLines</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
258 245 <span class="i">@DataLabels</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
259 246 <span class="k">if</span> <span class="s">(</span><span class="i">$SortDataLabels</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
260 247 <span class="k">push</span> <span class="i">@DataLabels</span><span class="cm">,</span> <span class="k">sort</span> <span class="i">@</span>{<span class="i">$DataHeaderLabelsRef</span>}<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
261 248 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
262 249 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
263 250 <span class="k">push</span> <span class="i">@DataLabels</span><span class="cm">,</span> <span class="i">@</span>{<span class="i">$DataHeaderLabelsRef</span>}<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
264 251 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
265 252 <span class="k">for</span> <span class="i">$DataLabel</span> <span class="s">(</span><span class="i">@DataLabels</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
266 253 <span class="i">$DataValue</span> = <span class="q">&#39;&#39;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
267 254 <span class="k">if</span> <span class="s">(</span><span class="k">exists</span> <span class="i">$DataHeaderLabelsAndValuesRef</span>-&gt;{<span class="i">$DataLabel</span>}<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
268 255 <span class="i">$DataValue</span> = <span class="i">$DataHeaderLabelsAndValuesRef</span>-&gt;{<span class="i">$DataLabel</span>}<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
269 256 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
270 257 <span class="k">push</span> <span class="i">@DataLines</span><span class="cm">,</span> <span class="s">(</span><span class="q">&quot;&gt; &lt;${DataLabel}&gt;&quot;</span><span class="cm">,</span> <span class="q">&quot;$DataValue&quot;</span><span class="cm">,</span> <span class="q">&quot;&quot;</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
271 258 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
272 259 <span class="k">return</span> <span class="i">@DataLines</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
273 260 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
274 261
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
275 262 <span class="c"># Parse data field header in SD file and return lists of all and common data field</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
276 263 <span class="c"># labels.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
277 <a name="GetAllAndCommonCmpdDataHeaderLabels-"></a> 264 <span class="k">sub </span><span class="m">GetAllAndCommonCmpdDataHeaderLabels</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
278 265 <span class="k">my</span><span class="s">(</span><span class="i">$SDFileRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
279 266 <span class="k">my</span><span class="s">(</span><span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$CmpdString</span><span class="cm">,</span> <span class="i">$Label</span><span class="cm">,</span> <span class="i">@CmpdLines</span><span class="cm">,</span> <span class="i">@DataFieldLabels</span><span class="cm">,</span> <span class="i">@CommonDataFieldLabels</span><span class="cm">,</span> <span class="i">%DataFieldLabelsMap</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
280 267
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
281 268 <span class="i">$CmpdCount</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
282 269 <span class="i">@DataFieldLabels</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
283 270 <span class="i">@CommonDataFieldLabels</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
284 271 <span class="i">%DataFieldLabelsMap</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
285 272
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
286 273 <span class="k">while</span> <span class="s">(</span><span class="i">$CmpdString</span> = <span class="i">ReadCmpdString</span><span class="s">(</span><span class="i">$SDFileRef</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
287 274 <span class="i">$CmpdCount</span>++<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
288 275 <span class="i">@CmpdLines</span> = <span class="k">split</span> <span class="q">&quot;\n&quot;</span><span class="cm">,</span> <span class="i">$CmpdString</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
289 276 <span class="c"># Process compound data header labels and figure out which ones are present for</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
290 277 <span class="c"># all the compounds...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
291 278 <span class="k">if</span> <span class="s">(</span><span class="i">@DataFieldLabels</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
292 279 <span class="k">my</span> <span class="s">(</span><span class="i">@CmpdDataFieldLabels</span><span class="s">)</span> = <span class="i">GetCmpdDataHeaderLabels</span><span class="s">(</span>\<span class="i">@CmpdLines</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
293 280 <span class="k">my</span><span class="s">(</span><span class="i">%CmpdDataFieldLabelsMap</span><span class="s">)</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
294 281 <span class="c"># Setup a map for the current labels...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
295 282 <span class="k">for</span> <span class="i">$Label</span> <span class="s">(</span><span class="i">@CmpdDataFieldLabels</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
296 283 <span class="i">$CmpdDataFieldLabelsMap</span>{<span class="i">$Label</span>} = <span class="q">&quot;PresentInSome&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
297 284 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
298 285 <span class="c"># Check the presence old labels for this compound; otherwise, mark &#39;em new...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
299 286 <span class="k">for</span> <span class="i">$Label</span> <span class="s">(</span><span class="i">@DataFieldLabels</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
300 287 <span class="k">if</span> <span class="s">(</span>!<span class="i">$CmpdDataFieldLabelsMap</span>{<span class="i">$Label</span>}<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
301 288 <span class="i">$DataFieldLabelsMap</span>{<span class="i">$Label</span>} = <span class="q">&quot;PresentInSome&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
302 289 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
303 290 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
304 291 <span class="c"># Check the presence this compound in the old labels; otherwise, add &#39;em...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
305 292 <span class="k">for</span> <span class="i">$Label</span> <span class="s">(</span><span class="i">@CmpdDataFieldLabels</span> <span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
306 293 <span class="k">if</span> <span class="s">(</span>!<span class="i">$DataFieldLabelsMap</span>{<span class="i">$Label</span>}<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
307 294 <span class="c"># It&#39;s a new label...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
308 295 <span class="k">push</span> <span class="i">@DataFieldLabels</span><span class="cm">,</span> <span class="i">$Label</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
309 296 <span class="i">$DataFieldLabelsMap</span>{<span class="i">$Label</span>} = <span class="q">&quot;PresentInSome&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
310 297 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
311 298 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
312 299 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
313 300 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
314 301 <span class="c"># Get the initial label set and set up a map...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
315 302 <span class="i">@DataFieldLabels</span> = <span class="i">GetCmpdDataHeaderLabels</span><span class="s">(</span>\<span class="i">@CmpdLines</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
316 303 <span class="k">for</span> <span class="i">$Label</span> <span class="s">(</span><span class="i">@DataFieldLabels</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
317 304 <span class="i">$DataFieldLabelsMap</span>{<span class="i">$Label</span>} = <span class="q">&quot;PresentInAll&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
318 305 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
319 306 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
320 307 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
321 308 <span class="c"># Identify the common data field labels...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
322 309 <span class="i">@CommonDataFieldLabels</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
323 310 <span class="k">for</span> <span class="i">$Label</span> <span class="s">(</span><span class="i">@DataFieldLabels</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
324 311 <span class="k">if</span> <span class="s">(</span><span class="i">$DataFieldLabelsMap</span>{<span class="i">$Label</span>} <span class="k">eq</span> <span class="q">&quot;PresentInAll&quot;</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
325 312 <span class="k">push</span> <span class="i">@CommonDataFieldLabels</span><span class="cm">,</span> <span class="i">$Label</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
326 313 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
327 314 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
328 315 <span class="k">return</span> <span class="s">(</span><span class="i">$CmpdCount</span><span class="cm">,</span> \<span class="i">@DataFieldLabels</span><span class="cm">,</span> \<span class="i">@CommonDataFieldLabels</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
329 316 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
330 317
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
331 318 <span class="c"># Parse all the data header labels and return &#39;em as an list...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
332 319 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
333 320 <span class="c"># Format:</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
334 321 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
335 322 <span class="c">#&gt; Data header line</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
336 323 <span class="c">#Data line(s)</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
337 324 <span class="c">#Blank line</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
338 325 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
339 326 <span class="c"># [Data Header] (one line) precedes each item of data, starts with a greater than (&gt;) sign, and</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
340 327 <span class="c"># contains at least one of the following:</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
341 328 <span class="c"># The field name enclosed in angle brackets. For example: &lt;melting.point&gt;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
342 329 <span class="c"># The field number, DTn , where n represents the number assigned to the field in a MACCS-II database</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
343 330 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
344 331 <span class="c">#Optional information for the data header includes:</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
345 332 <span class="c"># The compound&#146;s external and internal registry numbers. External registry numbers must be enclosed in parentheses.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
346 333 <span class="c"># Any combination of information</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
347 334 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
348 335 <span class="c">#The following are examples of valid data headers:</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
349 336 <span class="c">#&gt; &lt;MELTING.POINT&gt;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
350 337 <span class="c">#&gt; 55 (MD-08974) &lt;BOILING.POINT&gt; DT12</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
351 338 <span class="c">#&gt; DT12 55</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
352 339 <span class="c">#&gt; (MD-0894) &lt;BOILING.POINT&gt; FROM ARCHIVES</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
353 340 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
354 341 <span class="c">#Notes: Sometimes last blank line is missing and can be just followed by $$$$</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
355 342 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
356 <a name="GetCmpdDataHeaderLabels-"></a> 343 <span class="k">sub </span><span class="m">GetCmpdDataHeaderLabels</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
357 344 <span class="k">my</span><span class="s">(</span><span class="i">$CmpdLines</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
358 345 <span class="k">my</span><span class="s">(</span><span class="i">$CmpdLine</span><span class="cm">,</span> <span class="i">$Label</span><span class="cm">,</span> <span class="i">@Labels</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
359 346
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
360 347 <span class="i">@Labels</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
361 348 <span class="j">CMPDLINE:</span> <span class="k">for</span> <span class="i">$CmpdLine</span> <span class="s">(</span><span class="i">@$CmpdLines</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
362 349 <span class="k">if</span> <span class="s">(</span><span class="i">$CmpdLine</span> !~ <span class="q">/^&gt;/</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
363 350 <span class="k">next</span> <span class="j">CMPDLINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
364 351 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
365 352 <span class="c"># Does the line contains field name enclosed in angular brackets?</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
366 353 <span class="s">(</span><span class="i">$Label</span><span class="s">)</span> = <span class="i">$CmpdLine</span> =~ <span class="q">/&lt;.*?&gt;/g</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
367 354 <span class="k">if</span> <span class="s">(</span>!<span class="k">defined</span><span class="s">(</span><span class="i">$Label</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
368 355 <span class="k">next</span> <span class="j">CMPDLINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
369 356 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
370 357 <span class="i">$Label</span> =~ <span class="q">s/(&lt;|&gt;)//g</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
371 358 <span class="k">push</span> <span class="i">@Labels</span><span class="cm">,</span> <span class="i">$Label</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
372 359 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
373 360 <span class="k">return</span> <span class="s">(</span><span class="i">@Labels</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
374 361 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
375 362
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
376 363 <span class="c"># Parse all the data header labels and values</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
377 <a name="GetCmpdDataHeaderLabelsAndValues-"></a> 364 <span class="k">sub </span><span class="m">GetCmpdDataHeaderLabelsAndValues</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
378 365 <span class="k">my</span><span class="s">(</span><span class="i">$CmpdLines</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
379 366 <span class="k">my</span><span class="s">(</span><span class="i">$CmpdLine</span><span class="cm">,</span> <span class="i">$CurrentLabel</span><span class="cm">,</span> <span class="i">$Label</span><span class="cm">,</span> <span class="i">$Value</span><span class="cm">,</span> <span class="i">$ValueCount</span><span class="cm">,</span> <span class="i">$ProcessingLabelData</span><span class="cm">,</span> <span class="i">@Values</span><span class="cm">,</span> <span class="i">%DataFields</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
380 367
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
381 368 <span class="i">%DataFields</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
382 369 <span class="i">$ProcessingLabelData</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
383 370 <span class="i">$ValueCount</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
384 371 <span class="j">CMPDLINE:</span> <span class="k">for</span> <span class="i">$CmpdLine</span> <span class="s">(</span><span class="i">@$CmpdLines</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
385 372 <span class="k">if</span> <span class="s">(</span><span class="i">$CmpdLine</span> =~ <span class="q">/^\$\$\$\$/</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
386 373 <span class="k">last</span> <span class="j">CMPDLINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
387 374 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
388 375 <span class="k">if</span> <span class="s">(</span><span class="i">$CmpdLine</span> =~ <span class="q">/^&gt;/</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
389 376 <span class="c"># Does the line contains field name enclosed in angular brackets?</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
390 377 <span class="s">(</span><span class="i">$Label</span><span class="s">)</span> = <span class="i">$CmpdLine</span> =~ <span class="q">/&lt;.*?&gt;/g</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
391 378 <span class="k">if</span> <span class="s">(</span><span class="k">defined</span> <span class="i">$Label</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
392 379 <span class="i">$CurrentLabel</span> = <span class="i">$Label</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
393 380 <span class="i">$CurrentLabel</span> =~ <span class="q">s/(&lt;|&gt;)//g</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
394 381 <span class="i">$ProcessingLabelData</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
395 382 <span class="i">$ValueCount</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
396 383
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
397 384 <span class="k">if</span> <span class="s">(</span><span class="i">$CurrentLabel</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
398 385 <span class="i">$ProcessingLabelData</span> = <span class="n">1</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
399 386 <span class="i">$DataFields</span>{<span class="i">$CurrentLabel</span>} = <span class="q">&#39;&#39;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
400 387 <span class="k">next</span> <span class="j">CMPDLINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
401 388 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
402 389 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
403 390 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
404 391 <span class="k">if</span> <span class="s">(</span>!<span class="i">$ProcessingLabelData</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
405 392 <span class="c"># Data line containing no &lt;label&gt; as allowed by SDF format. Just ignore it...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
406 393 <span class="k">next</span> <span class="j">CMPDLINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
407 394 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
408 395 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
409 396 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
410 397 <span class="k">if</span> <span class="s">(</span>!<span class="i">$ProcessingLabelData</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
411 398 <span class="k">next</span> <span class="j">CMPDLINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
412 399 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
413 400 <span class="k">if</span> <span class="s">(</span>!<span class="s">(</span><span class="k">defined</span><span class="s">(</span><span class="i">$CmpdLine</span><span class="s">)</span> &amp;&amp; <span class="k">length</span><span class="s">(</span><span class="i">$CmpdLine</span><span class="s">)</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
414 401 <span class="c"># Blank line terminates value for a label...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
415 402 <span class="i">$CurrentLabel</span> = <span class="q">&#39;&#39;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
416 403 <span class="i">$ValueCount</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
417 404 <span class="i">$ProcessingLabelData</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
418 405 <span class="k">next</span> <span class="j">CMPDLINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
419 406 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
420 407 <span class="i">$ValueCount</span>++<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
421 408 <span class="i">$Value</span> = <span class="i">$CmpdLine</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
422 409
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
423 410 <span class="k">if</span> <span class="s">(</span><span class="i">$ValueCount</span> &gt; <span class="n">1</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
424 411 <span class="i">$DataFields</span>{<span class="i">$CurrentLabel</span>} .= <span class="q">&quot;\n&quot;</span> . <span class="i">$Value</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
425 412 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
426 413 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
427 414 <span class="i">$DataFields</span>{<span class="i">$CurrentLabel</span>} = <span class="i">$Value</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
428 415 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
429 416 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
430 417 <span class="k">return</span> <span class="s">(</span><span class="i">%DataFields</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
431 418 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
432 419
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
433 420 <span class="c"># Return an updated compoud string after removing data header label along with its</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
434 421 <span class="c"># value from the specified compound string...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
435 422 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
436 <a name="RemoveCmpdDataHeaderLabelAndValue-"></a> 423 <span class="k">sub </span><span class="m">RemoveCmpdDataHeaderLabelAndValue</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
437 424 <span class="k">my</span><span class="s">(</span><span class="i">$CmpdString</span><span class="cm">,</span> <span class="i">$DataHeaderLabel</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
438 425 <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="cm">,</span> <span class="i">$PorcessingDataHeaderLabel</span><span class="cm">,</span> <span class="i">@CmpdLines</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
439 426
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
440 427 <span class="i">@CmpdLines</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
441 428 <span class="i">$PorcessingDataHeaderLabel</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
442 429
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
443 430 <span class="j">CMPDLINE:</span> <span class="k">for</span> <span class="i">$Line</span> <span class="s">(</span><span class="k">split</span> <span class="q">&quot;\n&quot;</span><span class="cm">,</span> <span class="i">$CmpdString</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
444 431 <span class="k">if</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^&gt;/</span> &amp;&amp; <span class="i">$Line</span> =~ <span class="q">/&lt;$DataHeaderLabel&gt;/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
445 432 <span class="i">$PorcessingDataHeaderLabel</span> = <span class="n">1</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
446 433 <span class="k">next</span> <span class="j">CMPDLINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
447 434 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
448 435
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
449 436 <span class="k">if</span> <span class="s">(</span><span class="i">$PorcessingDataHeaderLabel</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
450 437 <span class="c"># Blank line indicates end of fingerprints data value...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
451 438 <span class="k">if</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^\$\$\$\$/</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
452 439 <span class="k">push</span> <span class="i">@CmpdLines</span><span class="cm">,</span> <span class="i">$Line</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
453 440 <span class="i">$PorcessingDataHeaderLabel</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
454 441 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
455 442 <span class="k">elsif</span> <span class="s">(</span>!<span class="k">length</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
456 443 <span class="i">$PorcessingDataHeaderLabel</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
457 444 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
458 445 <span class="k">next</span> <span class="j">CMPDLINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
459 446 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
460 447
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
461 448 <span class="c"># Track compound lines without fingerprints data...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
462 449 <span class="k">push</span> <span class="i">@CmpdLines</span><span class="cm">,</span> <span class="i">$Line</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
463 450 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
464 451
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
465 452 <span class="k">return</span> <span class="k">join</span> <span class="q">&quot;\n&quot;</span><span class="cm">,</span> <span class="i">@CmpdLines</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
466 453 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
467 454
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
468 455 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
469 456 <span class="c"># Using bond blocks, figure out the number of disconnected fragments and</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
470 457 <span class="c"># return their values along with the atom numbers in a string delimited by new</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
471 458 <span class="c"># line character.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
472 459 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
473 <a name="GetCmpdFragments-"></a> 460 <span class="k">sub </span><span class="m">GetCmpdFragments</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
474 461 <span class="k">my</span><span class="s">(</span><span class="i">$CmpdLines</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
475 462 <span class="k">my</span><span class="s">(</span><span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="cm">,</span> <span class="i">$FirstAtomNum</span><span class="cm">,</span> <span class="i">$SecondAtomNum</span><span class="cm">,</span> <span class="i">@AtomConnections</span><span class="cm">,</span> <span class="i">$BondType</span><span class="cm">,</span> <span class="i">$FragmentString</span><span class="cm">,</span> <span class="i">$FragmentCount</span><span class="cm">,</span> <span class="i">$LineIndex</span><span class="cm">,</span> <span class="i">$Index</span><span class="cm">,</span> <span class="i">$AtomNum</span><span class="cm">,</span> <span class="i">$NbrAtomNum</span><span class="cm">,</span> <span class="i">@ProcessedAtoms</span><span class="cm">,</span> <span class="i">$ProcessedAtomCount</span><span class="cm">,</span> <span class="i">$ProcessAtomNum</span><span class="cm">,</span> <span class="i">@ProcessingAtoms</span><span class="cm">,</span> <span class="i">@ConnectedAtoms</span><span class="cm">,</span> <span class="i">%Fragments</span><span class="cm">,</span> <span class="i">$FragmentNum</span><span class="cm">,</span> <span class="i">$AFragmentString</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
476 463
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
477 464 <span class="c"># Setup the connection table for each atom...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
478 465 <span class="i">@AtomConnections</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
479 466 <span class="s">(</span><span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="s">)</span> = <span class="i">ParseCmpdCountsLine</span><span class="s">(</span><span class="i">@$CmpdLines</span>[<span class="n">3</span>]<span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
480 467 <span class="k">for</span> <span class="i">$AtomNum</span> <span class="s">(</span><span class="n">1</span> .. <span class="i">$AtomCount</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
481 468 <span class="i">%</span>{<span class="i">$AtomConnections</span>[<span class="i">$AtomNum</span>]} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
482 469 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
483 470 <span class="k">for</span> <span class="s">(</span><span class="i">$LineIndex</span> = <span class="n">4</span> + <span class="i">$AtomCount</span><span class="sc">;</span> <span class="i">$LineIndex</span> &lt; <span class="s">(</span><span class="n">4</span> + <span class="i">$AtomCount</span> + <span class="i">$BondCount</span><span class="s">)</span><span class="sc">;</span> <span class="i">$LineIndex</span>++<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
484 471 <span class="s">(</span><span class="i">$FirstAtomNum</span><span class="cm">,</span> <span class="i">$SecondAtomNum</span><span class="cm">,</span> <span class="i">$BondType</span><span class="s">)</span> = <span class="i">ParseCmpdBondLine</span><span class="s">(</span><span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>]<span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
485 472 <span class="k">if</span> <span class="s">(</span>!<span class="i">$AtomConnections</span>[<span class="i">$FirstAtomNum</span>]{<span class="i">$SecondAtomNum</span>}<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
486 473 <span class="i">$AtomConnections</span>[<span class="i">$FirstAtomNum</span>]{<span class="i">$SecondAtomNum</span>} = <span class="i">$BondType</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
487 474 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
488 475 <span class="k">if</span> <span class="s">(</span>!<span class="i">$AtomConnections</span>[<span class="i">$SecondAtomNum</span>]{<span class="i">$FirstAtomNum</span>}<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
489 476 <span class="i">$AtomConnections</span>[<span class="i">$SecondAtomNum</span>]{<span class="i">$FirstAtomNum</span>} = <span class="i">$BondType</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
490 477 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
491 478 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
492 479
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
493 480 <span class="c">#Get set to count fragments...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
494 481 <span class="i">$ProcessedAtomCount</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
495 482 <span class="i">$FragmentNum</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
496 483 <span class="i">%Fragments</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
497 484 <span class="i">@ProcessedAtoms</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
498 485 <span class="k">for</span> <span class="i">$AtomNum</span> <span class="s">(</span><span class="n">1</span> .. <span class="i">$AtomCount</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
499 486 <span class="i">$ProcessedAtoms</span>[<span class="i">$AtomNum</span>] = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
500 487 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
501 488 <span class="k">while</span> <span class="s">(</span><span class="i">$ProcessedAtomCount</span> &lt; <span class="i">$AtomCount</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
502 489 <span class="i">@ProcessingAtoms</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
503 490 <span class="i">@ConnectedAtoms</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
504 491 <span class="j">ATOMNUM:</span> <span class="k">for</span> <span class="i">$AtomNum</span> <span class="s">(</span><span class="n">1</span> .. <span class="i">$AtomCount</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
505 492 <span class="k">if</span> <span class="s">(</span>!<span class="i">$ProcessedAtoms</span>[<span class="i">$AtomNum</span>]<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
506 493 <span class="i">$ProcessedAtomCount</span>++<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
507 494 <span class="i">$ProcessedAtoms</span>[<span class="i">$AtomNum</span>] = <span class="n">1</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
508 495 <span class="k">push</span> <span class="i">@ProcessingAtoms</span><span class="cm">,</span> <span class="i">$AtomNum</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
509 496 <span class="i">$FragmentNum</span>++<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
510 497 <span class="i">@</span>{<span class="i">$Fragments</span>{<span class="i">$FragmentNum</span>} } = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
511 498 <span class="k">push</span> <span class="i">@</span>{<span class="i">$Fragments</span>{<span class="i">$FragmentNum</span>} }<span class="cm">,</span> <span class="i">$AtomNum</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
512 499 <span class="k">last</span> <span class="j">ATOMNUM</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
513 500 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
514 501 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
515 502
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
516 503 <span class="c"># Go over the neighbors and follow the connection trail while collecting the</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
517 504 <span class="c"># atoms numbers present in the connected fragment...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
518 505 <span class="k">while</span> <span class="s">(</span><span class="i">@ProcessingAtoms</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
519 506 <span class="k">for</span> <span class="s">(</span><span class="i">$Index</span> = <span class="n">0</span><span class="sc">;</span> <span class="i">$Index</span> &lt; <span class="i">@ProcessingAtoms</span><span class="sc">;</span> <span class="i">$Index</span>++<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
520 507 <span class="i">$ProcessAtomNum</span> = <span class="i">$ProcessingAtoms</span>[<span class="i">$Index</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
521 508 <span class="k">for</span> <span class="i">$NbrAtomNum</span> <span class="s">(</span><span class="k">keys</span> <span class="i">%</span>{<span class="i">$AtomConnections</span>[<span class="i">$ProcessAtomNum</span>]}<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
522 509 <span class="k">if</span> <span class="s">(</span>!<span class="i">$ProcessedAtoms</span>[<span class="i">$NbrAtomNum</span>]<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
523 510 <span class="i">$ProcessedAtomCount</span>++<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
524 511 <span class="i">$ProcessedAtoms</span>[<span class="i">$NbrAtomNum</span>] = <span class="n">1</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
525 512 <span class="k">push</span> <span class="i">@ConnectedAtoms</span><span class="cm">,</span> <span class="i">$NbrAtomNum</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
526 513 <span class="k">push</span> <span class="i">@</span>{ <span class="i">$Fragments</span>{<span class="i">$FragmentNum</span>} }<span class="cm">,</span> <span class="i">$NbrAtomNum</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
527 514 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
528 515 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
529 516 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
530 517 <span class="i">@ProcessingAtoms</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
531 518 <span class="i">@ProcessingAtoms</span> = <span class="i">@ConnectedAtoms</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
532 519 <span class="i">@ConnectedAtoms</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
533 520 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
534 521 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
535 522 <span class="i">$FragmentCount</span> = <span class="i">$FragmentNum</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
536 523 <span class="i">$FragmentString</span> = <span class="q">&quot;&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
537 524
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
538 525 <span class="c"># Sort out the fragments by size...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
539 526 <span class="k">for</span> <span class="i">$FragmentNum</span> <span class="s">(</span><span class="k">sort</span> <span class="s">{</span> <span class="i">@</span>{<span class="i">$Fragments</span>{<span class="i">$b</span>}} &lt;=&gt; <span class="i">@</span>{<span class="i">$Fragments</span>{<span class="i">$a</span>}} <span class="s">}</span> <span class="k">keys</span> <span class="i">%Fragments</span> <span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
540 527 <span class="c"># Sort the atoms in a fragment by their numbers...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
541 528 <span class="i">$AFragmentString</span> = <span class="k">join</span> <span class="q">&quot; &quot;</span><span class="cm">,</span> <span class="k">sort</span> <span class="s">{</span> <span class="i">$a</span> &lt;=&gt; <span class="i">$b</span> <span class="s">}</span> <span class="i">@</span>{ <span class="i">$Fragments</span>{<span class="i">$FragmentNum</span>} }<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
542 529 <span class="k">if</span> <span class="s">(</span><span class="i">$FragmentString</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
543 530 <span class="i">$FragmentString</span> .= <span class="q">&quot;\n&quot;</span> . <span class="i">$AFragmentString</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
544 531 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
545 532 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
546 533 <span class="i">$FragmentString</span> = <span class="i">$AFragmentString</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
547 534 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
548 535 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
549 536 <span class="k">return</span> <span class="s">(</span><span class="i">$FragmentCount</span><span class="cm">,</span> <span class="i">$FragmentString</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
550 537 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
551 538
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
552 539 <span class="c"># Count number of lines present in between 4th and line containg &quot;M END&quot;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
553 <a name="GetCtabLinesCount-"></a> 540 <span class="k">sub </span><span class="m">GetCtabLinesCount</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
554 541 <span class="k">my</span><span class="s">(</span><span class="i">$CmpdLines</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
555 542 <span class="k">my</span><span class="s">(</span><span class="i">$LineIndex</span><span class="cm">,</span> <span class="i">$CtabLinesCount</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
556 543
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
557 544 <span class="i">$CtabLinesCount</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
558 545 <span class="j">LINE:</span> <span class="k">for</span> <span class="s">(</span><span class="i">$LineIndex</span> = <span class="n">4</span><span class="sc">;</span> <span class="i">$LineIndex</span> &lt; <span class="i">@$CmpdLines</span><span class="sc">;</span> <span class="i">$LineIndex</span>++<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
559 546 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
560 547 <span class="c"># Any line after atom and bond data starting with anything other than space or</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
561 548 <span class="c"># a digit indicates end of Ctab atom/bond data block...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
562 549 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
563 550 <span class="k">if</span> <span class="s">(</span><span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>] !~ <span class="q">/^[0-9 ]/</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
564 551 <span class="i">$CtabLinesCount</span> = <span class="i">$LineIndex</span> - <span class="n">4</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
565 552 <span class="k">last</span> <span class="j">LINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
566 553 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
567 554 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
568 555 <span class="k">return</span> <span class="i">$CtabLinesCount</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
569 556 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
570 557
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
571 558 <span class="c"># Using atom blocks, count the number of atoms which contain special element</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
572 559 <span class="c"># symbols not present in the periodic table.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
573 <a name="GetUnknownAtoms-"></a> 560 <span class="k">sub </span><span class="m">GetUnknownAtoms</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
574 561 <span class="k">my</span><span class="s">(</span><span class="i">$CmpdLines</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
575 562 <span class="k">my</span><span class="s">(</span><span class="i">$UnknownAtomCount</span><span class="cm">,</span> <span class="i">$UnknownAtoms</span><span class="cm">,</span> <span class="i">$UnknownAtomLines</span><span class="cm">,</span> <span class="i">$LineIndex</span><span class="cm">,</span> <span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$AtomSymbol</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
576 563
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
577 564 <span class="i">$UnknownAtomCount</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
578 565 <span class="i">$UnknownAtoms</span> = <span class="q">&quot;&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
579 566 <span class="i">$UnknownAtomLines</span> = <span class="q">&quot;&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
580 567 <span class="s">(</span><span class="i">$AtomCount</span><span class="s">)</span> = <span class="i">ParseCmpdCountsLine</span><span class="s">(</span><span class="i">@$CmpdLines</span>[<span class="n">3</span>]<span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
581 568 <span class="k">for</span> <span class="s">(</span><span class="i">$LineIndex</span> = <span class="n">4</span><span class="sc">;</span> <span class="i">$LineIndex</span> &lt; <span class="s">(</span><span class="n">4</span> + <span class="i">$AtomCount</span><span class="s">)</span><span class="sc">;</span> <span class="i">$LineIndex</span>++<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
582 569 <span class="s">(</span><span class="i">$AtomSymbol</span><span class="s">)</span> = <span class="i">ParseCmpdAtomLine</span><span class="s">(</span><span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>]<span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
583 570 <span class="k">if</span> <span class="s">(</span>!<span class="i">IsElement</span><span class="s">(</span><span class="i">$AtomSymbol</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
584 571 <span class="i">$UnknownAtomCount</span>++<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
585 572 <span class="i">$UnknownAtoms</span> .= <span class="q">&quot; $AtomSymbol&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
586 573 <span class="k">if</span> <span class="s">(</span><span class="i">$UnknownAtomLines</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
587 574 <span class="i">$UnknownAtomLines</span> .= <span class="q">&quot;\n&quot;</span> . <span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
588 575 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
589 576 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
590 577 <span class="i">$UnknownAtomLines</span> = <span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
591 578 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
592 579 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
593 580 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
594 581 <span class="k">return</span> <span class="s">(</span><span class="i">$UnknownAtomCount</span><span class="cm">,</span> <span class="i">$UnknownAtoms</span><span class="cm">,</span> <span class="i">$UnknownAtomLines</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
595 582 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
596 583
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
597 584 <span class="c"># Check z coordinates of all atoms to see whether any of them is non-zero</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
598 585 <span class="c"># which makes the compound geometry three dimensional...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
599 586 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
600 <a name="IsCmpd3D-"></a> 587 <span class="k">sub </span><span class="m">IsCmpd3D</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
601 588 <span class="k">my</span><span class="s">(</span><span class="i">$CmpdLines</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
602 589 <span class="k">my</span><span class="s">(</span><span class="i">$LineIndex</span><span class="cm">,</span> <span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$AtomSymbol</span><span class="cm">,</span> <span class="i">$AtomX</span><span class="cm">,</span> <span class="i">$AtomY</span><span class="cm">,</span> <span class="i">$AtomZ</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
603 590
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
604 591 <span class="s">(</span><span class="i">$AtomCount</span><span class="s">)</span> = <span class="i">ParseCmpdCountsLine</span><span class="s">(</span><span class="i">@$CmpdLines</span>[<span class="n">3</span>]<span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
605 592 <span class="k">for</span> <span class="s">(</span><span class="i">$LineIndex</span> = <span class="n">4</span><span class="sc">;</span> <span class="i">$LineIndex</span> &lt; <span class="s">(</span><span class="n">4</span> + <span class="i">$AtomCount</span><span class="s">)</span><span class="sc">;</span> <span class="i">$LineIndex</span>++<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
606 593 <span class="s">(</span><span class="i">$AtomSymbol</span><span class="cm">,</span> <span class="i">$AtomX</span><span class="cm">,</span> <span class="i">$AtomY</span><span class="cm">,</span> <span class="i">$AtomZ</span><span class="s">)</span> = <span class="i">ParseCmpdAtomLine</span><span class="s">(</span><span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>]<span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
607 594 <span class="k">if</span> <span class="s">(</span><span class="i">$AtomZ</span> != <span class="n">0</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
608 595 <span class="k">return</span> <span class="n">1</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
609 596 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
610 597 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
611 598 <span class="k">return</span> <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
612 599 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
613 600
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
614 601 <span class="c"># Check whether it&#39;s a 2D compound...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
615 602 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
616 <a name="IsCmpd2D-"></a> 603 <span class="k">sub </span><span class="m">IsCmpd2D</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
617 604 <span class="k">my</span><span class="s">(</span><span class="i">$CmpdLines</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
618 605
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
619 606 <span class="k">return</span> <span class="i">IsCmpd3D</span><span class="s">(</span><span class="i">$CmpdLines</span><span class="s">)</span> ? <span class="n">0</span> <span class="co">:</span> <span class="n">1</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
620 607 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
621 608
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
622 609 <span class="c"># Using bond blocks, count the number of bond lines which contain atom numbers</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
623 610 <span class="c"># greater than atom count specified in compound count line...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
624 611 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
625 <a name="GetInvalidAtomNumbers-"></a> 612 <span class="k">sub </span><span class="m">GetInvalidAtomNumbers</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
626 613 <span class="k">my</span><span class="s">(</span><span class="i">$CmpdLines</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
627 614 <span class="k">my</span><span class="s">(</span><span class="i">$LineIndex</span><span class="cm">,</span> <span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="cm">,</span> <span class="i">$FirstAtomNum</span><span class="cm">,</span> <span class="i">$SecondAtomNum</span><span class="cm">,</span> <span class="i">$InvalidAtomNumbersCount</span><span class="cm">,</span> <span class="i">$InvalidAtomNumbers</span><span class="cm">,</span> <span class="i">$InvalidAtomNumberLines</span><span class="cm">,</span> <span class="i">$Line</span><span class="cm">,</span> <span class="i">$InvalidAtomPropertyLine</span><span class="cm">,</span> <span class="i">$ValuePairIndex</span><span class="cm">,</span> <span class="i">$AtomNum</span><span class="cm">,</span> <span class="i">$Value</span><span class="cm">,</span> <span class="i">@ValuePairs</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
628 615
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
629 616 <span class="s">(</span><span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="s">)</span> = <span class="i">ParseCmpdCountsLine</span><span class="s">(</span><span class="i">@$CmpdLines</span>[<span class="n">3</span>]<span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
630 617
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
631 618 <span class="i">$InvalidAtomNumbersCount</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
632 619 <span class="i">$InvalidAtomNumbers</span> = <span class="q">&quot;&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
633 620 <span class="i">$InvalidAtomNumberLines</span> = <span class="q">&quot;&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
634 621
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
635 622 <span class="c"># Go over bond block lines...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
636 623 <span class="j">LINE:</span> <span class="k">for</span> <span class="s">(</span><span class="i">$LineIndex</span> = <span class="n">4</span> + <span class="i">$AtomCount</span><span class="sc">;</span> <span class="i">$LineIndex</span> &lt; <span class="s">(</span><span class="n">4</span> + <span class="i">$AtomCount</span> + <span class="i">$BondCount</span><span class="s">)</span><span class="sc">;</span> <span class="i">$LineIndex</span>++<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
637 624 <span class="s">(</span><span class="i">$FirstAtomNum</span><span class="cm">,</span> <span class="i">$SecondAtomNum</span><span class="s">)</span> = <span class="i">ParseCmpdBondLine</span><span class="s">(</span><span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>]<span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
638 625 <span class="k">if</span> <span class="s">(</span><span class="i">$FirstAtomNum</span> &lt;= <span class="i">$AtomCount</span> &amp;&amp; <span class="i">$SecondAtomNum</span> &lt;= <span class="i">$AtomCount</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
639 626 <span class="k">next</span> <span class="j">LINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
640 627 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
641 628 <span class="k">if</span> <span class="s">(</span><span class="i">$FirstAtomNum</span> &gt; <span class="i">$AtomCount</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
642 629 <span class="i">$InvalidAtomNumbersCount</span>++<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
643 630 <span class="i">$InvalidAtomNumbers</span> .= <span class="q">&quot; $FirstAtomNum&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
644 631 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
645 632 <span class="k">if</span> <span class="s">(</span><span class="i">$SecondAtomNum</span> &gt; <span class="i">$AtomCount</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
646 633 <span class="i">$InvalidAtomNumbersCount</span>++<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
647 634 <span class="i">$InvalidAtomNumbers</span> .= <span class="q">&quot; $SecondAtomNum&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
648 635 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
649 636 <span class="k">if</span> <span class="s">(</span><span class="i">$InvalidAtomNumberLines</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
650 637 <span class="i">$InvalidAtomNumberLines</span> .= <span class="q">&quot;\n&quot;</span> . <span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
651 638 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
652 639 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
653 640 <span class="i">$InvalidAtomNumberLines</span> = <span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
654 641 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
655 642 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
656 643 <span class="c"># Go over property lines before M END...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
657 644 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
658 645 <span class="j">LINE:</span> <span class="k">for</span> <span class="s">(</span><span class="i">$LineIndex</span> = <span class="s">(</span><span class="n">4</span> + <span class="i">$AtomCount</span> + <span class="i">$BondCount</span><span class="s">)</span><span class="sc">;</span> <span class="i">$LineIndex</span> &lt; <span class="i">@$CmpdLines</span><span class="sc">;</span> <span class="i">$LineIndex</span>++<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
659 646 <span class="i">$Line</span> = <span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
660 647 <span class="i">@ValuePairs</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
661 648 <span class="k">if</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^M END/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
662 649 <span class="k">last</span> <span class="j">LINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
663 650 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
664 651 <span class="i">@ValuePairs</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
665 652 <span class="k">if</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^M CHG/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
666 653 <span class="i">@ValuePairs</span> = <span class="i">ParseCmpdChargePropertyLine</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
667 654 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
668 655 <span class="k">elsif</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^M RAD/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
669 656 <span class="i">@ValuePairs</span> = <span class="i">ParseCmpdRadicalPropertyLine</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
670 657 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
671 658 <span class="k">elsif</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^M ISO/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
672 659 <span class="i">@ValuePairs</span> = <span class="i">ParseCmpdIsotopePropertyLine</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
673 660 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
674 661 <span class="k">elsif</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^A /i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
675 662 <span class="k">my</span><span class="s">(</span><span class="i">$NextLine</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
676 663 <span class="i">$LineIndex</span>++<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
677 664 <span class="i">$NextLine</span> = <span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
678 665 <span class="i">@ValuePairs</span> = <span class="i">ParseCmpdAtomAliasPropertyLine</span><span class="s">(</span><span class="i">$Line</span><span class="cm">,</span> <span class="i">$NextLine</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
679 666 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
680 667 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
681 668 <span class="k">next</span> <span class="j">LINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
682 669 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
683 670
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
684 671 <span class="i">$InvalidAtomPropertyLine</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
685 672 <span class="k">for</span> <span class="s">(</span><span class="i">$ValuePairIndex</span> = <span class="n">0</span><span class="sc">;</span> <span class="i">$ValuePairIndex</span> &lt; <span class="i">$#ValuePairs</span><span class="sc">;</span> <span class="i">$ValuePairIndex</span> += <span class="n">2</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
686 673 <span class="i">$AtomNum</span> = <span class="i">$ValuePairs</span>[<span class="i">$ValuePairIndex</span>]<span class="sc">;</span> <span class="i">$Value</span> = <span class="i">$ValuePairs</span>[<span class="i">$ValuePairIndex</span> + <span class="n">1</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
687 674 <span class="k">if</span> <span class="s">(</span><span class="i">$AtomNum</span> &gt; <span class="i">$AtomCount</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
688 675 <span class="i">$InvalidAtomPropertyLine</span> = <span class="n">1</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
689 676 <span class="i">$InvalidAtomNumbersCount</span>++<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
690 677 <span class="i">$InvalidAtomNumbers</span> .= <span class="q">&quot; $AtomNum&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
691 678 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
692 679 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
693 680 <span class="k">if</span> <span class="s">(</span><span class="i">$InvalidAtomPropertyLine</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
694 681 <span class="k">if</span> <span class="s">(</span><span class="i">$InvalidAtomNumberLines</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
695 682 <span class="i">$InvalidAtomNumberLines</span> .= <span class="q">&quot;\n&quot;</span> . <span class="i">$Line</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
696 683 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
697 684 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
698 685 <span class="i">$InvalidAtomNumberLines</span> = <span class="i">$Line</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
699 686 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
700 687 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
701 688 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
702 689
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
703 690 <span class="k">return</span> <span class="s">(</span><span class="i">$InvalidAtomNumbersCount</span><span class="cm">,</span> <span class="i">$InvalidAtomNumbers</span><span class="cm">,</span> <span class="i">$InvalidAtomNumberLines</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
704 691 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
705 692
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
706 693 <span class="c"># Ctab lines: Atom block</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
707 694 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
708 695 <span class="c"># Format: xxxxx.xxxxyyyyy.yyyyzzzzz.zzzz aaaddcccssshhhbbbvvvHHHrrriiimmmnnneee</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
709 696 <span class="c"># A10 A10 A10 xA3 A2A3 A3 A3 A3 A3 A3 A3 A3 A3 A3 A3</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
710 697 <span class="c"># x,y,z: Atom coordinates</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
711 698 <span class="c"># aaa: Atom symbol. Entry in periodic table or L for atom list, A, Q, * for unspecified</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
712 699 <span class="c"># atom, and LP for lone pair, or R# for Rgroup label</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
713 700 <span class="c"># dd: Mass difference. -3, -2, -1, 0, 1, 2, 3, 4 (0 for value beyond these limits)</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
714 701 <span class="c"># ccc: Charge. 0 = uncharged or value other than these, 1 = +3, 2 = +2, 3 = +1,</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
715 702 <span class="c"># 4 = doublet radical, 5 = -1, 6 = -2, 7 = -3</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
716 703 <span class="c"># sss: Atom stereo parity. 0 = not stereo, 1 = odd, 2 = even, 3 = either or unmarked stereo center</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
717 704 <span class="c"># hhh: Hydrogen count + 1. 1 = H0, 2 = H1, 3 = H2, 4 = H3, 5 = H4</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
718 705 <span class="c"># bbb: Stereo care box. 0 = ignore stereo configuration of this double bond atom, 1 = stereo</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
719 706 <span class="c"># configuration of double bond atom must match</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
720 707 <span class="c"># vvv: Valence. 0 = no marking (default)(1 to 14) = (1 to 14) 15 = zero valence</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
721 708 <span class="c"># HHH: H0 designator. 0 = not specified, 1 = no H atoms allowed (redundant due to hhh)</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
722 709 <span class="c"># rrr: Not used</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
723 710 <span class="c"># iii: Not used</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
724 711 <span class="c"># mmm: Atom-atom mapping number. 1 - number of atoms</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
725 712 <span class="c"># nnn: Inversion/retention flag. 0 = property not applied, 1 = configuration is inverted,</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
726 713 <span class="c"># 2 = configuration is retained.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
727 714 <span class="c"># eee: Exact change flag. 0 = property not applied, 1 = change on atom must be</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
728 715 <span class="c"># exactly as shown</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
729 716 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
730 717 <span class="c"># Notes:</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
731 718 <span class="c"># . StereoParity: 1 - ClockwiseStereo, 2 - AntiClockwiseStereo; 3 - Either; 0 - none. These</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
732 719 <span class="c"># values determine chirailty around the chiral center; a non zero value indicates atom</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
733 720 <span class="c"># has been marked as chiral center.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
734 721 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
735 <a name="ParseCmpdAtomLine-"></a> 722 <span class="k">sub </span><span class="m">ParseCmpdAtomLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
736 723 <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
737 724 <span class="k">my</span> <span class="s">(</span><span class="i">$LineIndex</span><span class="cm">,</span> <span class="i">$AtomX</span><span class="cm">,</span> <span class="i">$AtomY</span><span class="cm">,</span> <span class="i">$AtomZ</span><span class="cm">,</span> <span class="i">$AtomSymbol</span><span class="cm">,</span> <span class="i">$MassDifference</span><span class="cm">,</span> <span class="i">$Charge</span><span class="cm">,</span> <span class="i">$StereoParity</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
738 725
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
739 726 <span class="s">(</span><span class="i">$AtomX</span><span class="cm">,</span> <span class="i">$AtomY</span><span class="cm">,</span> <span class="i">$AtomZ</span><span class="cm">,</span> <span class="i">$AtomSymbol</span><span class="cm">,</span> <span class="i">$MassDifference</span><span class="cm">,</span> <span class="i">$Charge</span><span class="cm">,</span> <span class="i">$StereoParity</span><span class="s">)</span> = <span class="s">(</span><span class="q">&#39;&#39;</span><span class="s">)</span> x <span class="n">7</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
740 727 <span class="k">if</span> <span class="s">(</span><span class="k">length</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span> &gt; <span class="n">31</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
741 728 <span class="s">(</span><span class="i">$AtomX</span><span class="cm">,</span> <span class="i">$AtomY</span><span class="cm">,</span> <span class="i">$AtomZ</span><span class="cm">,</span> <span class="i">$AtomSymbol</span><span class="cm">,</span> <span class="i">$MassDifference</span><span class="cm">,</span> <span class="i">$Charge</span><span class="cm">,</span> <span class="i">$StereoParity</span><span class="s">)</span> = <span class="k">unpack</span><span class="s">(</span><span class="q">&quot;A10A10A10xA3A2A3A3&quot;</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
742 729 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
743 730 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
744 731 <span class="s">(</span><span class="i">$AtomX</span><span class="cm">,</span> <span class="i">$AtomY</span><span class="cm">,</span> <span class="i">$AtomZ</span><span class="cm">,</span> <span class="i">$AtomSymbol</span><span class="s">)</span> = <span class="k">unpack</span><span class="s">(</span><span class="q">&quot;A10A10A10&quot;</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
745 732 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
746 733 <span class="k">return</span> <span class="s">(</span><span class="i">$AtomSymbol</span><span class="cm">,</span> <span class="i">$AtomX</span><span class="cm">,</span> <span class="i">$AtomY</span><span class="cm">,</span> <span class="i">$AtomZ</span><span class="cm">,</span> <span class="i">$MassDifference</span><span class="cm">,</span> <span class="i">$Charge</span><span class="cm">,</span> <span class="i">$StereoParity</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
747 734 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
748 735
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
749 736 <span class="c"># Map MDL charge value used in SD and MOL files to internal charge used by MayaChemTools.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
750 737 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
751 <a name="MDLChargeToInternalCharge-"></a> 738 <span class="k">sub </span><span class="m">MDLChargeToInternalCharge</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
752 739 <span class="k">my</span><span class="s">(</span><span class="i">$MDLCharge</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
753 740 <span class="k">my</span><span class="s">(</span><span class="i">$InternalCharge</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
754 741
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
755 742 <span class="j">CHARGE:</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
756 743 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLCharge</span> == <span class="n">0</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalCharge</span> = <span class="n">0</span><span class="sc">;</span> <span class="k">last</span> <span class="j">CHARGE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
757 744 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLCharge</span> == <span class="n">1</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalCharge</span> = <span class="n">3</span><span class="sc">;</span> <span class="k">last</span> <span class="j">CHARGE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
758 745 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLCharge</span> == <span class="n">2</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalCharge</span> = <span class="n">2</span><span class="sc">;</span> <span class="k">last</span> <span class="j">CHARGE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
759 746 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLCharge</span> == <span class="n">3</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalCharge</span> = <span class="n">1</span><span class="sc">;</span> <span class="k">last</span> <span class="j">CHARGE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
760 747 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLCharge</span> == <span class="n">5</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalCharge</span> = <span class="n">-1</span><span class="sc">;</span> <span class="k">last</span> <span class="j">CHARGE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
761 748 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLCharge</span> == <span class="n">6</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalCharge</span> = <span class="n">-2</span><span class="sc">;</span> <span class="k">last</span> <span class="j">CHARGE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
762 749 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLCharge</span> == <span class="n">7</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalCharge</span> = <span class="n">-3</span><span class="sc">;</span> <span class="k">last</span> <span class="j">CHARGE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
763 750 <span class="c"># All other MDL charge values, including 4 corresponding to &quot;doublet radical&quot;,</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
764 751 <span class="c"># are assigned internal value of 0.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
765 752 <span class="i">$InternalCharge</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
766 753 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLCharge</span> != <span class="n">4</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
767 754 <span class="w">carp</span> <span class="q">&quot;Warning: MDLChargeToInternalCharge: MDL charge value, $MDLCharge, is not supported: An internal charge value, 0, has been assigned...&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
768 755 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
769 756 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
770 757 <span class="k">return</span> <span class="i">$InternalCharge</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
771 758 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
772 759
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
773 760 <span class="c"># Map internal charge used by MayaChemTools to MDL charge value used in SD and MOL files.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
774 761 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
775 <a name="InternalChargeToMDLCharge-"></a> 762 <span class="k">sub </span><span class="m">InternalChargeToMDLCharge</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
776 763 <span class="k">my</span><span class="s">(</span><span class="i">$InternalCharge</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
777 764 <span class="k">my</span><span class="s">(</span><span class="i">$MDLCharge</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
778 765
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
779 766 <span class="j">CHARGE:</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
780 767 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalCharge</span> == <span class="n">3</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLCharge</span> = <span class="n">1</span><span class="sc">;</span> <span class="k">last</span> <span class="j">CHARGE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
781 768 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalCharge</span> == <span class="n">2</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLCharge</span> = <span class="n">2</span><span class="sc">;</span> <span class="k">last</span> <span class="j">CHARGE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
782 769 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalCharge</span> == <span class="n">1</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLCharge</span> = <span class="n">3</span><span class="sc">;</span> <span class="k">last</span> <span class="j">CHARGE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
783 770 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalCharge</span> == <span class="n">-1</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLCharge</span> = <span class="n">5</span><span class="sc">;</span> <span class="k">last</span> <span class="j">CHARGE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
784 771 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalCharge</span> == <span class="n">-2</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLCharge</span> = <span class="n">6</span><span class="sc">;</span> <span class="k">last</span> <span class="j">CHARGE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
785 772 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalCharge</span> == <span class="n">-3</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLCharge</span> = <span class="n">7</span><span class="sc">;</span> <span class="k">last</span> <span class="j">CHARGE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
786 773 <span class="c"># All other MDL charge values, including 4 corresponding to &quot;doublet radical&quot;,</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
787 774 <span class="c"># are assigned internal value of 0.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
788 775 <span class="i">$MDLCharge</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
789 776 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
790 777 <span class="k">return</span> <span class="i">$MDLCharge</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
791 778 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
792 779
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
793 780 <span class="c"># Ctab lines: Bond block</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
794 781 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
795 782 <span class="c"># Format: 111222tttsssxxxrrrccc</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
796 783 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
797 784 <span class="c"># 111: First atom number.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
798 785 <span class="c"># 222: Second atom number.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
799 786 <span class="c"># ttt: Bond type. 1 = Single, 2 = Double, 3 = Triple, 4 = Aromatic, 5 = Single or Double,</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
800 787 <span class="c"># 6 = Single or Aromatic, 7 = Double or Aromatic, 8 = Any</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
801 788 <span class="c"># sss: Bond stereo. Single bonds: 0 = not stereo, 1 = Up, 4 = Either, 6 = Down,</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
802 789 <span class="c"># Double bonds: 0 = Use x-, y-, z-coords from atom block to determine cis or trans,</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
803 790 <span class="c"># 3 = Cis or trans (either) double bond</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
804 791 <span class="c"># xxx: Not used</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
805 792 <span class="c"># rrr: Bond topology. 0 = Either, 1 = Ring, 2 = Chain</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
806 793 <span class="c"># ccc: Reacting center status. 0 = unmarked, 1 = a center, -1 = not a center,</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
807 794 <span class="c"># Additional: 2 = no change,4 = bond made/broken, 8 = bond order changes 12 = 4+8</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
808 795 <span class="c"># (both made/broken and changes); 5 = (4 + 1), 9 = (8 + 1), and 13 = (12 + 1) are also possible</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
809 796 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
810 <a name="ParseCmpdBondLine-"></a> 797 <span class="k">sub </span><span class="m">ParseCmpdBondLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
811 798 <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
812 799 <span class="k">my</span><span class="s">(</span><span class="i">$FirstAtomNum</span><span class="cm">,</span> <span class="i">$SecondAtomNum</span><span class="cm">,</span> <span class="i">$BondType</span><span class="cm">,</span> <span class="i">$BondStereo</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
813 800
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
814 801 <span class="s">(</span><span class="i">$FirstAtomNum</span><span class="cm">,</span> <span class="i">$SecondAtomNum</span><span class="cm">,</span> <span class="i">$BondType</span><span class="cm">,</span> <span class="i">$BondStereo</span><span class="s">)</span> = <span class="k">map</span> <span class="s">{</span><span class="q">s/ //g</span><span class="sc">;</span> <span class="i">$_</span><span class="s">}</span> <span class="k">unpack</span><span class="s">(</span><span class="q">&quot;A3A3A3A3&quot;</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
815 802 <span class="k">return</span> <span class="s">(</span><span class="i">$FirstAtomNum</span><span class="cm">,</span> <span class="i">$SecondAtomNum</span><span class="cm">,</span> <span class="i">$BondType</span><span class="cm">,</span> <span class="i">$BondStereo</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
816 803 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
817 804
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
818 805 <span class="c"># Map MDL bond type value used in SD and MOL files to internal bond order and bond types</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
819 806 <span class="c"># values used by MayaChemTools...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
820 807 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
821 <a name="MDLBondTypeToInternalBondOrder-"></a> 808 <span class="k">sub </span><span class="m">MDLBondTypeToInternalBondOrder</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
822 809 <span class="k">my</span><span class="s">(</span><span class="i">$MDLBondType</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
823 810 <span class="k">my</span><span class="s">(</span><span class="i">$InternalBondOrder</span><span class="cm">,</span> <span class="i">$InternalBondType</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
824 811
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
825 812 <span class="i">$InternalBondType</span> = <span class="q">&#39;&#39;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
826 813
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
827 814 <span class="j">BONDTYPE:</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
828 815 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLBondType</span> == <span class="n">1</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalBondOrder</span> = <span class="n">1</span><span class="sc">;</span> <span class="i">$InternalBondType</span> = <span class="q">&#39;Single&#39;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDTYPE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
829 816 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLBondType</span> == <span class="n">2</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalBondOrder</span> = <span class="n">2</span><span class="sc">;</span> <span class="i">$InternalBondType</span> = <span class="q">&#39;Double&#39;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDTYPE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
830 817 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLBondType</span> == <span class="n">3</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalBondOrder</span> = <span class="n">3</span><span class="sc">;</span> <span class="i">$InternalBondType</span> = <span class="q">&#39;Triple&#39;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDTYPE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
831 818 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLBondType</span> == <span class="n">4</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalBondOrder</span> = <span class="n">1.5</span><span class="sc">;</span> <span class="i">$InternalBondType</span> = <span class="q">&#39;Aromatic&#39;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDTYPE</span><span class="sc">;</span><span class="s">}</span> <span class="c"># Aromatic</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
832 819 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLBondType</span> == <span class="n">5</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalBondOrder</span> = <span class="n">1</span><span class="sc">;</span> <span class="i">$InternalBondType</span> = <span class="q">&#39;SingleOrDouble&#39;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDTYPE</span><span class="sc">;</span><span class="s">}</span> <span class="c"># Aromatic</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
833 820 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLBondType</span> == <span class="n">6</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalBondOrder</span> = <span class="n">1</span><span class="sc">;</span> <span class="i">$InternalBondType</span> = <span class="q">&#39;SingleOrAromatic&#39;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDTYPE</span><span class="sc">;</span><span class="s">}</span> <span class="c"># Aromatic</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
834 821 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLBondType</span> == <span class="n">7</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalBondOrder</span> = <span class="n">2</span><span class="sc">;</span> <span class="i">$InternalBondType</span> = <span class="q">&#39;DoubleOrAromatic&#39;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDTYPE</span><span class="sc">;</span><span class="s">}</span> <span class="c"># Aromatic</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
835 822 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLBondType</span> == <span class="n">8</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalBondOrder</span> = <span class="n">1</span><span class="sc">;</span> <span class="i">$InternalBondType</span> = <span class="q">&#39;Any&#39;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDTYPE</span><span class="sc">;</span><span class="s">}</span> <span class="c"># Aromatic</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
836 823 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
837 824 <span class="c"># Although MDL aromatic bond values are used for query only and explicit Kekule bond order</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
838 825 <span class="c"># values must be assigned, internal value of 1.5 is allowed to indicate aromatic bond orders.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
839 826 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
840 827 <span class="c"># All other MDL bond type values - 5 = Single or Double, 6 = Single or Aromatic, 7 = Double or Aromatic,</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
841 828 <span class="c"># 8 = Any - are also assigned appropriate internal value of 1: These are meant to be used for</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
842 829 <span class="c"># structure queries by MDL products.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
843 830 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
844 831 <span class="i">$InternalBondOrder</span> = <span class="n">1</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
845 832 <span class="i">$InternalBondType</span> = <span class="q">&#39;Single&#39;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
846 833
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
847 834 <span class="w">carp</span> <span class="q">&quot;Warning: MDLBondTypeToInternalBondOrder: MDL bond type value, $MDLBondType, is not supported: An internal bond order value, 0, has been assigned...&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
848 835 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
849 836 <span class="k">return</span> <span class="s">(</span><span class="i">$InternalBondOrder</span><span class="cm">,</span> <span class="i">$InternalBondType</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
850 837 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
851 838
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
852 839 <span class="c"># Map internal bond order and bond type values used by MayaChemTools to MDL bond type value used</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
853 840 <span class="c"># in SD and MOL files...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
854 841 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
855 <a name="InternalBondOrderToMDLBondType-"></a> 842 <span class="k">sub </span><span class="m">InternalBondOrderToMDLBondType</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
856 843 <span class="k">my</span><span class="s">(</span><span class="i">$InternalBondOrder</span><span class="cm">,</span> <span class="i">$InternalBondType</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
857 844 <span class="k">my</span><span class="s">(</span><span class="i">$MDLBondType</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
858 845
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
859 846 <span class="j">BONDTYPE:</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
860 847 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalBondOrder</span> == <span class="n">1</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
861 848 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalBondType</span> =~ <span class="q">/^SingleOrDouble$/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
862 849 <span class="i">$MDLBondType</span> = <span class="n">5</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
863 850 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
864 851 <span class="k">elsif</span> <span class="s">(</span><span class="i">$InternalBondType</span> =~ <span class="q">/^SingleOrAromatic$/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
865 852 <span class="i">$MDLBondType</span> = <span class="n">6</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
866 853 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
867 854 <span class="k">elsif</span> <span class="s">(</span><span class="i">$InternalBondType</span> =~ <span class="q">/^Any$/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
868 855 <span class="i">$MDLBondType</span> = <span class="n">8</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
869 856 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
870 857 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
871 858 <span class="i">$MDLBondType</span> = <span class="n">1</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
872 859 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
873 860 <span class="i">$MDLBondType</span> = <span class="n">1</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
874 861 <span class="k">last</span> <span class="j">BONDTYPE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
875 862 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
876 863 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalBondOrder</span> == <span class="n">2</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
877 864 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalBondType</span> =~ <span class="q">/^DoubleOrAromatic$/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
878 865 <span class="i">$MDLBondType</span> = <span class="n">7</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
879 866 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
880 867 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
881 868 <span class="i">$MDLBondType</span> = <span class="n">2</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
882 869 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
883 870 <span class="k">last</span> <span class="j">BONDTYPE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
884 871 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
885 872 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalBondOrder</span> == <span class="n">3</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLBondType</span> = <span class="n">3</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDTYPE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
886 873 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalBondOrder</span> == <span class="n">1.5</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLBondType</span> = <span class="n">4</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDTYPE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
887 874 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalBondType</span> =~ <span class="q">/^Any$/i</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLBondType</span> = <span class="n">8</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDTYPE</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
888 875
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
889 876 <span class="i">$MDLBondType</span> = <span class="n">1</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
890 877
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
891 878 <span class="w">carp</span> <span class="q">&quot;Warning: InternalBondOrderToMDLBondType: Internal bond order and type values, $InternalBondOrder and $InternalBondType, don&#39;t match any valid MDL bond type: MDL bond type value, 1, has been assigned...&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
892 879 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
893 880 <span class="k">return</span> <span class="i">$MDLBondType</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
894 881 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
895 882
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
896 883 <span class="c"># Third line: Comments - A blank line is also allowed.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
897 <a name="ParseCmpdCommentsLine-"></a> 884 <span class="k">sub </span><span class="m">ParseCmpdCommentsLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
898 885 <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
899 886 <span class="k">my</span><span class="s">(</span><span class="i">$Comments</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
900 887
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
901 888 <span class="i">$Comments</span> = <span class="k">unpack</span><span class="s">(</span><span class="q">&quot;A80&quot;</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
902 889
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
903 890 <span class="k">return</span> <span class="s">(</span><span class="i">$Comments</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
904 891 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
905 892
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
906 893 <span class="c"># Map MDL bond stereo value used in SD and MOL files to internal bond stereochemistry values used by MayaChemTools...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
907 894 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
908 <a name="MDLBondStereoToInternalBondStereochemistry-"></a> 895 <span class="k">sub </span><span class="m">MDLBondStereoToInternalBondStereochemistry</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
909 896 <span class="k">my</span><span class="s">(</span><span class="i">$MDLBondStereo</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
910 897 <span class="k">my</span><span class="s">(</span><span class="i">$InternalBondStereo</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
911 898
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
912 899 <span class="i">$InternalBondStereo</span> = <span class="q">&#39;&#39;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
913 900
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
914 901 <span class="j">BONDSTEREO:</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
915 902 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLBondStereo</span> == <span class="n">1</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalBondStereo</span> = <span class="q">&#39;Up&#39;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDSTEREO</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
916 903 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLBondStereo</span> == <span class="n">4</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalBondStereo</span> = <span class="q">&#39;UpOrDown&#39;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDSTEREO</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
917 904 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLBondStereo</span> == <span class="n">6</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalBondStereo</span> = <span class="q">&#39;Down&#39;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDSTEREO</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
918 905 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLBondStereo</span> == <span class="n">3</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalBondStereo</span> = <span class="q">&#39;CisOrTrans&#39;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDSTEREO</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
919 906 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLBondStereo</span> == <span class="n">0</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalBondStereo</span> = <span class="q">&#39;None&#39;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDSTEREO</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
920 907
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
921 908 <span class="i">$InternalBondStereo</span> = <span class="q">&#39;&#39;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
922 909 <span class="w">carp</span> <span class="q">&quot;Warning: MDLBondStereoToInternalBondType: MDL bond stereo value, $MDLBondStereo, is not supported: It has been ignored and bond order would be used to determine bond type...&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
923 910 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
924 911 <span class="k">return</span> <span class="i">$InternalBondStereo</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
925 912 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
926 913
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
927 914 <span class="c"># Map internal bond stereochemistry values used by MayaChemTools to MDL bond stereo value used in SD and MOL files...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
928 915 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
929 <a name="InternalBondStereochemistryToMDLBondStereo-"></a> 916 <span class="k">sub </span><span class="m">InternalBondStereochemistryToMDLBondStereo</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
930 917 <span class="k">my</span><span class="s">(</span><span class="i">$InternalBondStereo</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
931 918 <span class="k">my</span><span class="s">(</span><span class="i">$MDLBondStereo</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
932 919
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
933 920 <span class="i">$MDLBondStereo</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
934 921
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
935 922 <span class="j">BONDSTEREO:</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
936 923 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalBondStereo</span> =~ <span class="q">/^Up$/i</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLBondStereo</span> = <span class="n">1</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDSTEREO</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
937 924 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalBondStereo</span> =~ <span class="q">/^UpOrDown$/i</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLBondStereo</span> = <span class="n">4</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDSTEREO</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
938 925 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalBondStereo</span> =~ <span class="q">/^Down$/</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLBondStereo</span> = <span class="n">6</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDSTEREO</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
939 926 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalBondStereo</span> =~ <span class="q">/^CisOrTrans$/</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLBondStereo</span> = <span class="n">3</span><span class="sc">;</span> <span class="k">last</span> <span class="j">BONDSTEREO</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
940 927
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
941 928 <span class="i">$MDLBondStereo</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
942 929 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
943 930 <span class="k">return</span> <span class="i">$MDLBondStereo</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
944 931 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
945 932
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
946 933 <span class="c"># Fourth line: Counts</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
947 934 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
948 935 <span class="c"># Format: aaabbblllfffcccsssxxxrrrpppiiimmmvvvvvv</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
949 936 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
950 937 <span class="c"># aaa: number of atoms; bbb: number of bonds; lll: number of atom lists; fff: (obsolete)</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
951 938 <span class="c"># ccc: chiral flag: 0=not chiral, 1=chiral; sss: number of stext entries; xxx,rrr,ppp,iii:</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
952 939 <span class="c"># (obsolete); mmm: number of lines of additional properties, including the M END line, No</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
953 940 <span class="c"># longer supported, default is set to 999; vvvvvv: version</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
954 941
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
955 <a name="ParseCmpdCountsLine-"></a> 942 <span class="k">sub </span><span class="m">ParseCmpdCountsLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
956 943 <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
957 944 <span class="k">my</span><span class="s">(</span><span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="cm">,</span> <span class="i">$ChiralFlag</span><span class="cm">,</span> <span class="i">$PropertyCount</span><span class="cm">,</span> <span class="i">$Version</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
958 945
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
959 946 <span class="k">if</span> <span class="s">(</span><span class="k">length</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span> &gt;= <span class="n">39</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
960 947 <span class="s">(</span><span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="cm">,</span> <span class="i">$ChiralFlag</span><span class="cm">,</span> <span class="i">$PropertyCount</span><span class="cm">,</span> <span class="i">$Version</span><span class="s">)</span> = <span class="k">unpack</span><span class="s">(</span><span class="q">&quot;A3A3x3x3A3x3x3x3x3x3A3A6&quot;</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
961 948 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
962 949 <span class="k">elsif</span> <span class="s">(</span><span class="k">length</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span> &gt;= <span class="n">15</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
963 950 <span class="s">(</span><span class="i">$PropertyCount</span><span class="cm">,</span> <span class="i">$Version</span><span class="s">)</span> = <span class="s">(</span><span class="q">&quot;999&quot;</span><span class="cm">,</span> <span class="q">&quot;v2000&quot;</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
964 951 <span class="s">(</span><span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="cm">,</span> <span class="i">$ChiralFlag</span><span class="s">)</span> = <span class="k">unpack</span><span class="s">(</span><span class="q">&quot;A3A3x3x3A3&quot;</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
965 952 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
966 953 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
967 954 <span class="s">(</span><span class="i">$ChiralFlag</span><span class="cm">,</span> <span class="i">$PropertyCount</span><span class="cm">,</span> <span class="i">$Version</span><span class="s">)</span> = <span class="s">(</span><span class="q">&quot;0&quot;</span><span class="cm">,</span> <span class="q">&quot;999&quot;</span><span class="cm">,</span> <span class="q">&quot;v2000&quot;</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
968 955 <span class="s">(</span><span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="s">)</span> = <span class="k">unpack</span><span class="s">(</span><span class="q">&quot;A3A3&quot;</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
969 956 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
970 957
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
971 958 <span class="k">if</span> <span class="s">(</span><span class="i">$Version</span> =~ <span class="q">/V3000/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
972 959 <span class="c"># Current version of MayaChemTools modules and classes for processing MDL MOL and SD don&#39;t support</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
973 960 <span class="c"># V3000. So instead of relying on callers, just exit with an error to disable any processing of V3000</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
974 961 <span class="c"># format.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
975 962 <span class="w">croak</span> <span class="q">&quot;Error: SDFileUtil::ParseCmpdCountsLine: The Extended Connection Table (V3000) format in MDL MOL and SD files is not supported by the current release of MayaChemTools...&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
976 963 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
977 964
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
978 965 <span class="k">return</span> <span class="s">(</span><span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="cm">,</span> <span class="i">$ChiralFlag</span><span class="cm">,</span> <span class="i">$PropertyCount</span><span class="cm">,</span> <span class="i">$Version</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
979 966 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
980 967
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
981 968 <span class="c"># Second line: Misc info</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
982 969 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
983 970 <span class="c"># Format: IIPPPPPPPPMMDDYYHHmmddSSssssssssssEEEEEEEEEEEERRRRRR</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
984 971 <span class="c"># A2A8 A10 A2I2A10 A12 A6</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
985 972 <span class="c"># User&#39;s first and last initials (I), program name (P), date/time (M/D/Y,H:m),</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
986 973 <span class="c"># dimensional codes - 2D or 3D (d),scaling factors (S, s), energy (E) if modeling program input,</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
987 974 <span class="c"># internal registry number (R) if input through MDL form. A blank line is also allowed.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
988 <a name="ParseCmpdMiscInfoLine-"></a> 975 <span class="k">sub </span><span class="m">ParseCmpdMiscInfoLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
989 976 <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
990 977 <span class="k">my</span><span class="s">(</span><span class="i">$UserInitial</span><span class="cm">,</span> <span class="i">$ProgramName</span><span class="cm">,</span> <span class="i">$Date</span><span class="cm">,</span> <span class="i">$Code</span><span class="cm">,</span> <span class="i">$ScalingFactor1</span><span class="cm">,</span> <span class="i">$ScalingFactor2</span><span class="cm">,</span> <span class="i">$Energy</span><span class="cm">,</span> <span class="i">$RegistryNum</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
991 978
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
992 979 <span class="s">(</span><span class="i">$UserInitial</span><span class="cm">,</span> <span class="i">$ProgramName</span><span class="cm">,</span> <span class="i">$Date</span><span class="cm">,</span> <span class="i">$Code</span><span class="cm">,</span> <span class="i">$ScalingFactor1</span><span class="cm">,</span> <span class="i">$ScalingFactor2</span><span class="cm">,</span> <span class="i">$Energy</span><span class="cm">,</span> <span class="i">$RegistryNum</span><span class="s">)</span> = <span class="k">unpack</span><span class="s">(</span><span class="q">&quot;A2A8A10A2A2A10A12A6&quot;</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
993 980 <span class="k">return</span> <span class="s">(</span><span class="i">$UserInitial</span><span class="cm">,</span> <span class="i">$ProgramName</span><span class="cm">,</span> <span class="i">$Date</span><span class="cm">,</span> <span class="i">$Code</span><span class="cm">,</span> <span class="i">$ScalingFactor1</span><span class="cm">,</span> <span class="i">$ScalingFactor2</span><span class="cm">,</span> <span class="i">$Energy</span><span class="cm">,</span> <span class="i">$RegistryNum</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
994 981 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
995 982
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
996 983 <span class="c"># First line: Molecule name. This line is unformatted, but like all other lines in a</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
997 984 <span class="c"># molfile may not extend beyond column 80. A blank line is also allowed.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
998 <a name="ParseCmpdMolNameLine-"></a> 985 <span class="k">sub </span><span class="m">ParseCmpdMolNameLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
999 986 <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1000 987 <span class="k">my</span><span class="s">(</span><span class="i">$MolName</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1001 988
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1002 989 <span class="i">$MolName</span> = <span class="k">unpack</span><span class="s">(</span><span class="q">&quot;A80&quot;</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1003 990
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1004 991 <span class="k">return</span> <span class="s">(</span><span class="i">$MolName</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1005 992 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1006 993
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1007 994 <span class="c"># Parse atom alias property line in CTAB generic properties block.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1008 995 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1009 996 <span class="c"># Atom alias property line format:</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1010 997 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1011 998 <span class="c"># A aaa</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1012 999 <span class="c"># x...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1013 1000 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1014 1001 <span class="c"># aaa: Atom number</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1015 1002 <span class="c"># x: Atom alias in next line</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1016 1003 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1017 <a name="ParseCmpdAtomAliasPropertyLine-"></a>1004 <span class="k">sub </span><span class="m">ParseCmpdAtomAliasPropertyLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1018 1005 <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="cm">,</span> <span class="i">$NextLine</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1019 1006 <span class="k">my</span><span class="s">(</span><span class="i">$Label</span><span class="cm">,</span> <span class="i">$AtomNumber</span><span class="cm">,</span> <span class="i">$AtomAlias</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1020 1007
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1021 1008 <span class="s">(</span><span class="i">$Label</span><span class="cm">,</span> <span class="i">$AtomNumber</span><span class="s">)</span> = <span class="k">split</span><span class="s">(</span><span class="q">&#39; &#39;</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1022 1009 <span class="i">$AtomAlias</span> = <span class="i">$NextLine</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1023 1010
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1024 1011 <span class="k">if</span> <span class="s">(</span>!<span class="i">$AtomAlias</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1025 1012 <span class="w">carp</span> <span class="q">&quot;Warning: _ParseCmpdAtomAliasPropertyLine: No atom alias value specified on the line following atom alias property line...&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1026 1013 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1027 1014
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1028 1015 <span class="k">return</span> <span class="s">(</span><span class="i">$AtomNumber</span><span class="cm">,</span> <span class="i">$AtomAlias</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1029 1016 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1030 1017
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1031 1018 <span class="c"># Parse charge property line in CTAB generic properties block.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1032 1019 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1033 1020 <span class="c"># Charge property line format:</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1034 1021 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1035 1022 <span class="c"># M CHGnn8 aaa vvv ...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1036 1023 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1037 1024 <span class="c"># nn8: Number of value pairs. Maximum of 8 pairs allowed.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1038 1025 <span class="c"># aaa: Atom number</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1039 1026 <span class="c"># vvv: -15 to +15. Default of 0 = uncharged atom. When present, this property supersedes</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1040 1027 <span class="c"># all charge and radical values in the atom block, forcing a 0 charge on all atoms not</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1041 1028 <span class="c"># listed in an M CHG or M RAD line.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1042 1029 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1043 <a name="ParseCmpdChargePropertyLine-"></a>1030 <span class="k">sub </span><span class="m">ParseCmpdChargePropertyLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1044 1031 <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1045 1032
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1046 1033 <span class="k">return</span> <span class="i">_ParseCmpdGenericPropertyLine</span><span class="s">(</span><span class="q">&#39;Charge&#39;</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1047 1034 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1048 1035
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1049 1036
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1050 1037 <span class="c"># Parse isotope property line in CTAB generic properties block.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1051 1038 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1052 1039 <span class="c"># Isoptope property line format:</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1053 1040 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1054 1041 <span class="c"># M ISOnn8 aaa vvv ...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1055 1042 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1056 1043 <span class="c"># nn8: Number of value paris. Maximum of 8 pairs allowed.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1057 1044 <span class="c"># aaa: Atom number</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1058 1045 <span class="c"># vvv: Absolute mass of the atom isotope as a positive integer. When present, this property</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1059 1046 <span class="c"># supersedes all isotope values in the atom block. Default (no entry) means natural</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1060 1047 <span class="c"># abundance. The difference between this absolute mass value and the natural</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1061 1048 <span class="c"># abundance value specified in the PTABLE.DAT file must be within the range of -18</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1062 1049 <span class="c"># to +12</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1063 1050 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1064 1051 <span class="c"># Notes:</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1065 1052 <span class="c"># . Values correspond to mass numbers...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1066 1053 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1067 <a name="ParseCmpdIsotopePropertyLine-"></a>1054 <span class="k">sub </span><span class="m">ParseCmpdIsotopePropertyLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1068 1055 <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1069 1056
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1070 1057 <span class="k">return</span> <span class="i">_ParseCmpdGenericPropertyLine</span><span class="s">(</span><span class="q">&#39;Isotope&#39;</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1071 1058 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1072 1059
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1073 1060 <span class="c"># Parse radical property line in CTAB generic properties block.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1074 1061 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1075 1062 <span class="c"># Radical property line format:</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1076 1063 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1077 1064 <span class="c"># M RADnn8 aaa vvv ...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1078 1065 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1079 1066 <span class="c"># nn8: Number of value paris. Maximum of 8 pairs allowed.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1080 1067 <span class="c"># aaa: Atom number</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1081 1068 <span class="c"># vvv: Default of 0 = no radical, 1 = singlet, 2 = doublet, 3 = triplet . When</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1082 1069 <span class="c"># present, this property supersedes all charge and radical values in the atom block,</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1083 1070 <span class="c"># forcing a 0 (zero) charge and radical on all atoms not listed in an M CHG or</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1084 1071 <span class="c"># M RAD line.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1085 1072 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1086 <a name="ParseCmpdRadicalPropertyLine-"></a>1073 <span class="k">sub </span><span class="m">ParseCmpdRadicalPropertyLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1087 1074 <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1088 1075
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1089 1076 <span class="k">return</span> <span class="i">_ParseCmpdGenericPropertyLine</span><span class="s">(</span><span class="q">&#39;Radical&#39;</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1090 1077 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1091 1078
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1092 1079 <span class="c"># Map MDL radical stereo value used in SD and MOL files to internal spin multiplicity values used by MayaChemTools...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1093 1080 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1094 <a name="MDLRadicalToInternalSpinMultiplicity-"></a>1081 <span class="k">sub </span><span class="m">MDLRadicalToInternalSpinMultiplicity</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1095 1082 <span class="k">my</span><span class="s">(</span><span class="i">$MDLRadical</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1096 1083 <span class="k">my</span><span class="s">(</span><span class="i">$InternalSpinMultiplicity</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1097 1084
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1098 1085 <span class="i">$InternalSpinMultiplicity</span> = <span class="q">&#39;&#39;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1099 1086
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1100 1087 <span class="j">SPINMULTIPLICITY:</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1101 1088 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLRadical</span> == <span class="n">0</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalSpinMultiplicity</span> = <span class="n">0</span><span class="sc">;</span> <span class="k">last</span> <span class="j">SPINMULTIPLICITY</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1102 1089 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLRadical</span> == <span class="n">1</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalSpinMultiplicity</span> = <span class="n">1</span><span class="sc">;</span> <span class="k">last</span> <span class="j">SPINMULTIPLICITY</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1103 1090 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLRadical</span> == <span class="n">2</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalSpinMultiplicity</span> = <span class="n">2</span><span class="sc">;</span> <span class="k">last</span> <span class="j">SPINMULTIPLICITY</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1104 1091 <span class="k">if</span> <span class="s">(</span><span class="i">$MDLRadical</span> == <span class="n">3</span><span class="s">)</span> <span class="s">{</span> <span class="i">$InternalSpinMultiplicity</span> = <span class="n">3</span><span class="sc">;</span> <span class="k">last</span> <span class="j">SPINMULTIPLICITY</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1105 1092 <span class="i">$InternalSpinMultiplicity</span> = <span class="q">&#39;&#39;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1106 1093 <span class="w">carp</span> <span class="q">&quot;Warning: MDLRadicalToInternalSpinMultiplicity: MDL radical value, $MDLRadical, specifed on line M RAD is not supported...&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1107 1094 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1108 1095 <span class="k">return</span> <span class="i">$InternalSpinMultiplicity</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1109 1096 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1110 1097
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1111 1098 <span class="c"># Map internal spin multiplicity values used by MayaChemTools to MDL radical stereo value used in SD and MOL files...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1112 1099 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1113 <a name="InternalSpinMultiplicityToMDLRadical-"></a>1100 <span class="k">sub </span><span class="m">InternalSpinMultiplicityToMDLRadical</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1114 1101 <span class="k">my</span><span class="s">(</span><span class="i">$InternalSpinMultiplicity</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1115 1102 <span class="k">my</span><span class="s">(</span><span class="i">$MDLRadical</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1116 1103
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1117 1104 <span class="i">$MDLRadical</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1118 1105
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1119 1106 <span class="j">SPINMULTIPLICITY:</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1120 1107 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalSpinMultiplicity</span> == <span class="n">1</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLRadical</span> = <span class="n">1</span><span class="sc">;</span> <span class="k">last</span> <span class="j">SPINMULTIPLICITY</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1121 1108 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalSpinMultiplicity</span> == <span class="n">2</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLRadical</span> = <span class="n">2</span><span class="sc">;</span> <span class="k">last</span> <span class="j">SPINMULTIPLICITY</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1122 1109 <span class="k">if</span> <span class="s">(</span><span class="i">$InternalSpinMultiplicity</span> == <span class="n">3</span><span class="s">)</span> <span class="s">{</span> <span class="i">$MDLRadical</span> = <span class="n">3</span><span class="sc">;</span> <span class="k">last</span> <span class="j">SPINMULTIPLICITY</span><span class="sc">;</span><span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1123 1110 <span class="i">$MDLRadical</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1124 1111 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1125 1112 <span class="k">return</span> <span class="i">$MDLRadical</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1126 1113 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1127 1114
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1128 1115 <span class="c"># Process generic CTAB property line...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1129 <a name="_ParseCmpdGenericPropertyLine-"></a>1116 <span class="k">sub </span><span class="m">_ParseCmpdGenericPropertyLine</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1130 1117 <span class="k">my</span><span class="s">(</span><span class="i">$PropertyName</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1131 1118
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1132 1119 <span class="k">my</span><span class="s">(</span><span class="i">$Label</span><span class="cm">,</span> <span class="i">$PropertyLabel</span><span class="cm">,</span> <span class="i">$ValuesCount</span><span class="cm">,</span> <span class="i">$ValuePairsCount</span><span class="cm">,</span> <span class="i">@ValuePairs</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1133 1120
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1134 1121 <span class="i">@ValuePairs</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1135 1122 <span class="s">(</span><span class="i">$Label</span><span class="cm">,</span> <span class="i">$PropertyLabel</span><span class="cm">,</span> <span class="i">$ValuesCount</span><span class="cm">,</span> <span class="i">@ValuePairs</span><span class="s">)</span> = <span class="k">split</span><span class="s">(</span><span class="q">&#39; &#39;</span><span class="cm">,</span> <span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1136 1123 <span class="i">$ValuePairsCount</span> = <span class="s">(</span><span class="k">scalar</span> <span class="i">@ValuePairs</span><span class="s">)</span>/<span class="n">2</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1137 1124 <span class="k">if</span> <span class="s">(</span><span class="i">$ValuesCount</span> != <span class="i">$ValuePairsCount</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1138 1125 <span class="w">carp</span> <span class="q">&quot;Warning: _ParseCmpdGenericPropertyLine: Number of atom number and $PropertyName value paris specified on $Label $PropertyLabel property line, $ValuePairsCount, does not match expected value of $ValuesCount...&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1139 1126 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1140 1127
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1141 1128 <span class="k">return</span> <span class="s">(</span><span class="i">@ValuePairs</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1142 1129 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1143 1130
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1144 1131 <span class="c"># Generic CTAB property lines for charge, istope and radical properties...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1145 1132 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1146 <a name="_GenerateCmpdGenericPropertyLines-"></a>1133 <span class="k">sub </span><span class="m">_GenerateCmpdGenericPropertyLines</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1147 1134 <span class="k">my</span><span class="s">(</span><span class="i">$PropertyName</span><span class="cm">,</span> <span class="i">$PropertyValuePairsRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1148 1135 <span class="k">my</span><span class="s">(</span><span class="i">$Index</span><span class="cm">,</span> <span class="i">$PropertyLabel</span><span class="cm">,</span> <span class="i">$Line</span><span class="cm">,</span> <span class="i">$PropertyCount</span><span class="cm">,</span> <span class="i">$AtomNum</span><span class="cm">,</span> <span class="i">$PropertyValue</span><span class="cm">,</span> <span class="i">@PropertyLines</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1149 1136
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1150 1137 <span class="i">@PropertyLines</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1151 1138 <span class="j">NAME:</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1152 1139 <span class="k">if</span> <span class="s">(</span><span class="i">$PropertyName</span> =~ <span class="q">/^Charge$/i</span><span class="s">)</span> <span class="s">{</span> <span class="i">$PropertyLabel</span> = <span class="q">&quot;M CHG&quot;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">NAME</span><span class="sc">;</span> <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1153 1140 <span class="k">if</span> <span class="s">(</span><span class="i">$PropertyName</span> =~ <span class="q">/^Isotope$/i</span><span class="s">)</span> <span class="s">{</span> <span class="i">$PropertyLabel</span> = <span class="q">&quot;M ISO&quot;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">NAME</span><span class="sc">;</span> <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1154 1141 <span class="k">if</span> <span class="s">(</span><span class="i">$PropertyName</span> =~ <span class="q">/^Radical$/i</span><span class="s">)</span> <span class="s">{</span> <span class="i">$PropertyLabel</span> = <span class="q">&quot;M RAD&quot;</span><span class="sc">;</span> <span class="k">last</span> <span class="j">NAME</span><span class="sc">;</span> <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1155 1142 <span class="w">carp</span> <span class="q">&quot;Warning: _GenerateCmpdGenericPropertyLines: Unknown property name, $PropertyName, specified...&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1156 1143 <span class="k">return</span> <span class="i">@PropertyLines</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1157 1144 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1158 1145
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1159 1146 <span class="c"># A maximum of 8 property pair values allowed per line...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1160 1147 <span class="i">$PropertyCount</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1161 1148 <span class="i">$Line</span> = <span class="q">&#39;&#39;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1162 1149 <span class="k">for</span> <span class="s">(</span><span class="i">$Index</span> = <span class="n">0</span><span class="sc">;</span> <span class="i">$Index</span> &lt; <span class="i">$#</span>{<span class="i">$PropertyValuePairsRef</span>}<span class="sc">;</span> <span class="i">$Index</span> += <span class="n">2</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1163 1150 <span class="k">if</span> <span class="s">(</span><span class="i">$PropertyCount</span> &gt; <span class="n">8</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1164 1151 <span class="c"># Setup property line...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1165 1152 <span class="i">$Line</span> = <span class="q">&quot;${PropertyLabel} 8${Line}&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1166 1153 <span class="k">push</span> <span class="i">@PropertyLines</span><span class="cm">,</span> <span class="i">$Line</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1167 1154
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1168 1155 <span class="i">$PropertyCount</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1169 1156 <span class="i">$Line</span> = <span class="q">&#39;&#39;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1170 1157 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1171 1158 <span class="i">$PropertyCount</span>++<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1172 1159 <span class="i">$AtomNum</span> = <span class="i">$PropertyValuePairsRef</span>-&gt;[<span class="i">$Index</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1173 1160 <span class="i">$PropertyValue</span> = <span class="i">$PropertyValuePairsRef</span>-&gt;[<span class="i">$Index</span> + <span class="n">1</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1174 1161 <span class="i">$Line</span> .= <span class="k">sprintf</span> <span class="q">&quot; %3i %3i&quot;</span><span class="cm">,</span> <span class="i">$AtomNum</span><span class="cm">,</span> <span class="i">$PropertyValue</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1175 1162 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1176 1163 <span class="k">if</span> <span class="s">(</span><span class="i">$Line</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1177 1164 <span class="i">$Line</span> = <span class="q">&quot;${PropertyLabel} ${PropertyCount}${Line}&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1178 1165 <span class="k">push</span> <span class="i">@PropertyLines</span><span class="cm">,</span> <span class="i">$Line</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1179 1166 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1180 1167 <span class="k">return</span> <span class="i">@PropertyLines</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1181 1168 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1182 1169
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1183 1170 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1184 1171 <span class="c"># Read compound data into a string and return its value</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1185 <a name="ReadCmpdString-"></a>1172 <span class="k">sub </span><span class="m">ReadCmpdString</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1186 1173 <span class="k">my</span><span class="s">(</span><span class="i">$SDFileRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1187 1174 <span class="k">my</span><span class="s">(</span><span class="i">$CmpdString</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1188 1175
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1189 1176 <span class="i">$CmpdString</span> = <span class="q">&quot;&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1190 1177 <span class="j">LINE:</span> <span class="k">while</span> <span class="s">(</span><span class="k">defined</span><span class="s">(</span><span class="i">$_</span> = <span class="q">&lt;$SDFileRef&gt;</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1191 1178 <span class="c"># Change Windows and Mac new line char to UNIX...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1192 1179 <span class="q">s/(\r\n)|(\r)/\n/g</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1193 1180
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1194 1181 <span class="k">if</span> <span class="s">(</span><span class="q">/^\$\$\$\$/</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1195 1182 <span class="c"># Take out any new line char at the end by explicitly removing it instead of using</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1196 1183 <span class="c"># chomp, which might not always work correctly on files generated on a system</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1197 1184 <span class="c"># with a value of input line separator different from the current system...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1198 1185 <span class="q">s/\n$//g</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1199 1186
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1200 1187 <span class="c"># Doesn&#39;t hurt to chomp...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1201 1188 <span class="k">chomp</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1202 1189
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1203 1190 <span class="i">$CmpdString</span> .= <span class="i">$_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1204 1191 <span class="k">last</span> <span class="j">LINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1205 1192 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1206 1193 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1207 1194 <span class="i">$CmpdString</span> .= <span class="i">$_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1208 1195 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1209 1196 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1210 1197 <span class="k">return</span> <span class="i">$CmpdString</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1211 1198 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1212 1199
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1213 1200 <span class="c"># Find out the number of fragements in the compounds. And for the compound with</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1214 1201 <span class="c"># more than one fragment, remove all the others besides the largest one.</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1215 <a name="WashCmpd-"></a>1202 <span class="k">sub </span><span class="m">WashCmpd</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1216 1203 <span class="k">my</span><span class="s">(</span><span class="i">$CmpdLines</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1217 1204 <span class="k">my</span><span class="s">(</span><span class="i">$WashedCmpdString</span><span class="cm">,</span> <span class="i">$FragmentCount</span><span class="cm">,</span> <span class="i">$Fragments</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1218 1205
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1219 1206 <span class="i">$WashedCmpdString</span> = <span class="q">&quot;&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1220 1207 <span class="s">(</span><span class="i">$FragmentCount</span><span class="cm">,</span> <span class="i">$Fragments</span><span class="s">)</span> = <span class="i">GetCmpdFragments</span><span class="s">(</span><span class="i">$CmpdLines</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1221 1208 <span class="k">if</span> <span class="s">(</span><span class="i">$FragmentCount</span> &gt; <span class="n">1</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1222 1209 <span class="c"># Go over the compound data for the largest fragment including property</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1223 1210 <span class="c"># data...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1224 1211 <span class="k">my</span> <span class="s">(</span><span class="i">@AllFragments</span><span class="cm">,</span> <span class="i">@LargestFragment</span><span class="cm">,</span> <span class="i">%LargestFragmentAtoms</span><span class="cm">,</span> <span class="i">@WashedCmpdLines</span><span class="cm">,</span> <span class="i">$Index</span><span class="cm">,</span> <span class="i">$LineIndex</span><span class="cm">,</span> <span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="cm">,</span> <span class="i">$NewAtomCount</span><span class="cm">,</span> <span class="i">$NewBondCount</span><span class="cm">,</span> <span class="i">$FirstAtomNum</span><span class="cm">,</span> <span class="i">$SecondAtomNum</span><span class="cm">,</span> <span class="i">$BondType</span><span class="cm">,</span> <span class="i">$BondStereo</span><span class="cm">,</span> <span class="i">$FirstNewAtomNum</span><span class="cm">,</span> <span class="i">$SecondNewAtomNum</span><span class="cm">,</span> <span class="i">$AtomNum</span><span class="cm">,</span> <span class="i">$ChiralFlag</span><span class="cm">,</span> <span class="i">$BondLine</span><span class="cm">,</span> <span class="i">$MENDLineIndex</span><span class="cm">,</span> <span class="i">$Line</span><span class="cm">,</span> <span class="i">$Value</span><span class="cm">,</span> <span class="i">@ValuePairs</span><span class="cm">,</span> <span class="i">@NewValuePairs</span><span class="cm">,</span> <span class="i">$ValuePairIndex</span><span class="cm">,</span> <span class="i">$NewAtomNum</span><span class="cm">,</span> <span class="i">@NewPropertyLines</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1225 1212
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1226 1213 <span class="i">@AllFragments</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> <span class="i">@LargestFragment</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1227 1214 <span class="i">%LargestFragmentAtoms</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1228 1215 <span class="i">@AllFragments</span> = <span class="k">split</span> <span class="q">&quot;\n&quot;</span><span class="cm">,</span> <span class="i">$Fragments</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1229 1216 <span class="i">@LargestFragment</span> = <span class="k">split</span> <span class="q">&quot; &quot;</span><span class="cm">,</span> <span class="i">$AllFragments</span>[<span class="n">0</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1230 1217 <span class="k">for</span> <span class="i">$Index</span> <span class="s">(</span><span class="n">0</span> .. <span class="i">$#LargestFragment</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1231 1218 <span class="c"># Map old atom numbers to new atom numbers as the fragment atom numbers are sorted</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1232 1219 <span class="c"># from lowest to highest old atom numbers...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1233 1220 <span class="i">$LargestFragmentAtoms</span>{<span class="i">$LargestFragment</span>[<span class="i">$Index</span>]} = <span class="i">$Index</span> + <span class="n">1</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1234 1221 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1235 1222 <span class="i">@WashedCmpdLines</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1236 1223 <span class="k">push</span> <span class="i">@WashedCmpdLines</span><span class="cm">,</span> <span class="i">@$CmpdLines</span>[<span class="n">0</span>]<span class="cm">,</span> <span class="i">@$CmpdLines</span>[<span class="n">1</span>]<span class="cm">,</span> <span class="i">@$CmpdLines</span>[<span class="n">2</span>]<span class="cm">,</span> <span class="i">@$CmpdLines</span>[<span class="n">3</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1237 1224 <span class="s">(</span><span class="i">$AtomCount</span><span class="cm">,</span> <span class="i">$BondCount</span><span class="cm">,</span> <span class="i">$ChiralFlag</span><span class="s">)</span> = <span class="i">ParseCmpdCountsLine</span><span class="s">(</span><span class="i">@$CmpdLines</span>[<span class="n">3</span>]<span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1238 1225 <span class="i">$NewAtomCount</span> = <span class="i">@LargestFragment</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1239 1226 <span class="i">$NewBondCount</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1240 1227 <span class="i">$AtomNum</span> = <span class="n">0</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1241 1228 <span class="c"># Retrieve the largest fragment atom lines...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1242 1229 <span class="k">for</span> <span class="s">(</span><span class="i">$LineIndex</span> = <span class="n">4</span><span class="sc">;</span> <span class="i">$LineIndex</span> &lt; <span class="s">(</span><span class="n">4</span> + <span class="i">$AtomCount</span><span class="s">)</span><span class="sc">;</span> <span class="i">$LineIndex</span>++<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1243 1230 <span class="i">$AtomNum</span>++<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1244 1231 <span class="k">if</span> <span class="s">(</span><span class="i">$LargestFragmentAtoms</span>{<span class="i">$AtomNum</span>}<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1245 1232 <span class="k">push</span> <span class="i">@WashedCmpdLines</span><span class="cm">,</span> <span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1246 1233 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1247 1234 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1248 1235 <span class="c"># Retrieve the largest fragment bond lines...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1249 1236 <span class="k">for</span> <span class="s">(</span><span class="i">$LineIndex</span> = <span class="n">4</span> + <span class="i">$AtomCount</span><span class="sc">;</span> <span class="i">$LineIndex</span> &lt; <span class="s">(</span><span class="n">4</span> + <span class="i">$AtomCount</span> + <span class="i">$BondCount</span><span class="s">)</span><span class="sc">;</span> <span class="i">$LineIndex</span>++<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1250 1237 <span class="s">(</span><span class="i">$FirstAtomNum</span><span class="cm">,</span> <span class="i">$SecondAtomNum</span><span class="cm">,</span> <span class="i">$BondType</span><span class="cm">,</span> <span class="i">$BondStereo</span><span class="s">)</span> = <span class="i">ParseCmpdBondLine</span><span class="s">(</span><span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>]<span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1251 1238 <span class="k">if</span> <span class="s">(</span><span class="i">$LargestFragmentAtoms</span>{<span class="i">$FirstAtomNum</span>} &amp;&amp; <span class="i">$LargestFragmentAtoms</span>{<span class="i">$SecondAtomNum</span>}<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1252 1239 <span class="i">$NewBondCount</span>++<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1253 1240 <span class="c"># Set up bond line with new atom number mapping...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1254 1241 <span class="i">$FirstNewAtomNum</span> = <span class="i">$LargestFragmentAtoms</span>{<span class="i">$FirstAtomNum</span>}<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1255 1242 <span class="i">$SecondNewAtomNum</span> = <span class="i">$LargestFragmentAtoms</span>{<span class="i">$SecondAtomNum</span>}<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1256 1243 <span class="i">$BondLine</span> = <span class="i">GenerateCmpdBondLine</span><span class="s">(</span><span class="i">$FirstNewAtomNum</span><span class="cm">,</span> <span class="i">$SecondNewAtomNum</span><span class="cm">,</span> <span class="i">$BondType</span><span class="cm">,</span> <span class="i">$BondStereo</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1257 1244 <span class="k">push</span> <span class="i">@WashedCmpdLines</span><span class="cm">,</span> <span class="i">$BondLine</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1258 1245 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1259 1246 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1260 1247 <span class="c"># Get property lines for CHG, ISO and RAD label and map the old atom numbers to new</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1261 1248 <span class="c"># atom numners; Others, property lines before M END line are skipped as atom numbers for</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1262 1249 <span class="c"># other properties might not valid anymore...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1263 1250 <span class="c">#</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1264 1251 <span class="i">$MENDLineIndex</span> = <span class="i">$LineIndex</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1265 1252 <span class="j">LINE:</span> <span class="k">for</span> <span class="s">(</span><span class="i">$LineIndex</span> = <span class="s">(</span><span class="n">4</span> + <span class="i">$AtomCount</span> + <span class="i">$BondCount</span><span class="s">)</span><span class="sc">;</span> <span class="i">$LineIndex</span> &lt; <span class="i">@$CmpdLines</span><span class="sc">;</span> <span class="i">$LineIndex</span>++<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1266 1253 <span class="i">$Line</span> = <span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1267 1254 <span class="k">if</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^M END/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1268 1255 <span class="k">push</span> <span class="i">@WashedCmpdLines</span><span class="cm">,</span> <span class="q">&quot;M END&quot;</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1269 1256 <span class="i">$MENDLineIndex</span> = <span class="i">$LineIndex</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1270 1257 <span class="k">last</span> <span class="j">LINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1271 1258 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1272 1259
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1273 1260 <span class="i">@ValuePairs</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1274 1261 <span class="k">if</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^M CHG/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1275 1262 <span class="i">@ValuePairs</span> = <span class="i">ParseCmpdChargePropertyLine</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1276 1263 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1277 1264 <span class="k">elsif</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^M RAD/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1278 1265 <span class="i">@ValuePairs</span> = <span class="i">ParseCmpdRadicalPropertyLine</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1279 1266 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1280 1267 <span class="k">elsif</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^M ISO/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1281 1268 <span class="i">@ValuePairs</span> = <span class="i">ParseCmpdIsotopePropertyLine</span><span class="s">(</span><span class="i">$Line</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1282 1269 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1283 1270 <span class="k">elsif</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^A /i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1284 1271 <span class="k">my</span><span class="s">(</span><span class="i">$NextLine</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1285 1272 <span class="i">$LineIndex</span>++<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1286 1273 <span class="i">$NextLine</span> = <span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1287 1274 <span class="i">@ValuePairs</span> = <span class="i">ParseCmpdAtomAliasPropertyLine</span><span class="s">(</span><span class="i">$Line</span><span class="cm">,</span> <span class="i">$NextLine</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1288 1275 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1289 1276 <span class="k">else</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1290 1277 <span class="k">next</span> <span class="j">LINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1291 1278 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1292 1279
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1293 1280 <span class="k">if</span> <span class="s">(</span>!<span class="i">@ValuePairs</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1294 1281 <span class="k">next</span> <span class="j">LINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1295 1282 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1296 1283
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1297 1284 <span class="c"># Collect values for valid atom numbers with mapping to new atom numbers...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1298 1285 <span class="i">@NewValuePairs</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1299 1286 <span class="j">VALUEINDEX:</span> <span class="k">for</span> <span class="s">(</span><span class="i">$ValuePairIndex</span> = <span class="n">0</span><span class="sc">;</span> <span class="i">$ValuePairIndex</span> &lt; <span class="i">$#ValuePairs</span><span class="sc">;</span> <span class="i">$ValuePairIndex</span> += <span class="n">2</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1300 1287 <span class="i">$AtomNum</span> = <span class="i">$ValuePairs</span>[<span class="i">$ValuePairIndex</span>]<span class="sc">;</span> <span class="i">$Value</span> = <span class="i">$ValuePairs</span>[<span class="i">$ValuePairIndex</span> + <span class="n">1</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1301 1288 <span class="k">if</span> <span class="s">(</span>!<span class="k">exists</span> <span class="i">$LargestFragmentAtoms</span>{<span class="i">$AtomNum</span>}<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1302 1289 <span class="k">next</span> <span class="j">VALUEINDEX</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1303 1290 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1304 1291 <span class="i">$NewAtomNum</span> = <span class="i">$LargestFragmentAtoms</span>{<span class="i">$AtomNum</span>}<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1305 1292 <span class="k">push</span> <span class="i">@NewValuePairs</span><span class="cm">,</span> <span class="s">(</span><span class="i">$NewAtomNum</span><span class="cm">,</span> <span class="i">$Value</span><span class="s">)</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1306 1293 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1307 1294 <span class="k">if</span> <span class="s">(</span>!<span class="i">@NewValuePairs</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1308 1295 <span class="k">next</span> <span class="j">LINE</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1309 1296 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1310 1297 <span class="i">@NewPropertyLines</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1311 1298 <span class="k">if</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^M CHG/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1312 1299 <span class="i">@NewPropertyLines</span> = <span class="i">GenerateCmpdChargePropertyLines</span><span class="s">(</span>\<span class="i">@NewValuePairs</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1313 1300 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1314 1301 <span class="k">elsif</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^M RAD/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1315 1302 <span class="i">@NewPropertyLines</span> = <span class="i">GenerateCmpdRadicalPropertyLines</span><span class="s">(</span>\<span class="i">@NewValuePairs</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1316 1303 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1317 1304 <span class="k">elsif</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^M ISO/i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1318 1305 <span class="i">@NewPropertyLines</span> = <span class="i">GenerateCmpdIsotopePropertyLines</span><span class="s">(</span>\<span class="i">@NewValuePairs</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1319 1306 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1320 1307 <span class="k">elsif</span> <span class="s">(</span><span class="i">$Line</span> =~ <span class="q">/^A /i</span><span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1321 1308 <span class="i">@NewPropertyLines</span> = <span class="i">GenerateCmpdAtomAliasPropertyLines</span><span class="s">(</span>\<span class="i">@NewValuePairs</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1322 1309 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1323 1310 <span class="k">push</span> <span class="i">@WashedCmpdLines</span><span class="cm">,</span> <span class="i">@NewPropertyLines</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1324 1311 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1325 1312
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1326 1313 <span class="c"># Retrieve rest of the data label and value property data...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1327 1314 <span class="k">for</span> <span class="s">(</span><span class="i">$LineIndex</span> = <span class="s">(</span><span class="n">1</span> + <span class="i">$MENDLineIndex</span><span class="s">)</span><span class="sc">;</span> <span class="i">$LineIndex</span> &lt; <span class="i">@$CmpdLines</span><span class="sc">;</span> <span class="i">$LineIndex</span>++<span class="s">)</span> <span class="s">{</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1328 1315 <span class="k">push</span> <span class="i">@WashedCmpdLines</span><span class="cm">,</span> <span class="i">@$CmpdLines</span>[<span class="i">$LineIndex</span>]<span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1329 1316 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1330 1317 <span class="c"># Update atom and bond count line...</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1331 1318 <span class="i">$WashedCmpdLines</span>[<span class="n">3</span>] = <span class="i">GenerateCmpdCountsLine</span><span class="s">(</span><span class="i">$NewAtomCount</span><span class="cm">,</span> <span class="i">$NewBondCount</span><span class="cm">,</span> <span class="i">$ChiralFlag</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1332 1319
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1333 1320 <span class="i">$WashedCmpdString</span> = <span class="k">join</span> <span class="q">&quot;\n&quot;</span><span class="cm">,</span> <span class="i">@WashedCmpdLines</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1334 1321 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1335 1322 <span class="k">return</span> <span class="s">(</span><span class="i">$FragmentCount</span><span class="cm">,</span> <span class="i">$Fragments</span><span class="cm">,</span> <span class="i">$WashedCmpdString</span><span class="s">)</span><span class="sc">;</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1336 1323 <span class="s">}</span>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1337 1324
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1338 <a name="EOF-"></a></pre>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1339 <p>&nbsp;</p>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1340 <br />
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1341 <center>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1342 <img src="../../../images/h2o2.png">
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1343 </center>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1344 </body>
73ae111cf86f Uploaded
deepakjadmin
parents:
diff changeset
1345 </html>