<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<meta name="generator" content="Doxygen 1.8.5"/>
<title>cloudy: /home/iurt/rpmbuild/BUILD/c13.01/source/atom_pop5.cpp Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<link href="search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="search/search.js"></script>
<script type="text/javascript">
$(document).ready(function() { searchBox.OnSelectItem(0); });
</script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
<tbody>
<tr style="height: 56px;">
<td style="padding-left: 0.5em;">
<div id="projectname">cloudy
 <span id="projectnumber">trunk</span>
</div>
</td>
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.5 -->
<script type="text/javascript">
var searchBox = new SearchBox("searchBox", "search",false,'Search');
</script>
<div id="navrow1" class="tabs">
<ul class="tablist">
<li><a href="index.html"><span>Main Page</span></a></li>
<li><a href="pages.html"><span>Related Pages</span></a></li>
<li><a href="namespaces.html"><span>Namespaces</span></a></li>
<li><a href="annotated.html"><span>Data Structures</span></a></li>
<li class="current"><a href="files.html"><span>Files</span></a></li>
<li>
<div id="MSearchBox" class="MSearchBoxInactive">
<span class="left">
<img id="MSearchSelect" src="search/mag_sel.png"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
alt=""/>
<input type="text" id="MSearchField" value="Search" accesskey="S"
onfocus="searchBox.OnSearchFieldFocus(true)"
onblur="searchBox.OnSearchFieldFocus(false)"
onkeyup="searchBox.OnSearchFieldChange(event)"/>
</span><span class="right">
<a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="search/close.png" alt=""/></a>
</span>
</div>
</li>
</ul>
</div>
<div id="navrow2" class="tabs2">
<ul class="tablist">
<li><a href="files.html"><span>File List</span></a></li>
<li><a href="globals.html"><span>Globals</span></a></li>
</ul>
</div>
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
onkeydown="return searchBox.OnSearchSelectKey(event)">
<a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(0)"><span class="SelectionMark"> </span>All</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(1)"><span class="SelectionMark"> </span>Data Structures</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(2)"><span class="SelectionMark"> </span>Namespaces</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(3)"><span class="SelectionMark"> </span>Files</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(4)"><span class="SelectionMark"> </span>Functions</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(5)"><span class="SelectionMark"> </span>Variables</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(6)"><span class="SelectionMark"> </span>Typedefs</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(7)"><span class="SelectionMark"> </span>Enumerations</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(8)"><span class="SelectionMark"> </span>Enumerator</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(9)"><span class="SelectionMark"> </span>Friends</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(10)"><span class="SelectionMark"> </span>Macros</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(11)"><span class="SelectionMark"> </span>Pages</a></div>
<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="javascript:void(0)" frameborder="0"
name="MSearchResults" id="MSearchResults">
</iframe>
</div>
<div id="nav-path" class="navpath">
<ul>
<li class="navelem"><a class="el" href="dir_b2f33c71d4aa5e7af42a1ca61ff5af1b.html">source</a></li> </ul>
</div>
</div><!-- top -->
<div class="header">
<div class="headertitle">
<div class="title">atom_pop5.cpp</div> </div>
</div><!--header-->
<div class="contents">
<a href="atom__pop5_8cpp.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span> <span class="comment">/* This file is part of Cloudy and is copyright (C)1978-2013 by Gary J. Ferland and</span></div>
<div class="line"><a name="l00002"></a><span class="lineno"> 2</span> <span class="comment"> * others. For conditions of distribution and use see copyright notice in license.txt */</span></div>
<div class="line"><a name="l00003"></a><span class="lineno"> 3</span> <span class="comment">/*atom_pop5 do five level atom population and cooling */</span></div>
<div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="preprocessor">#include "<a class="code" href="cddefines_8h.html">cddefines.h</a>"</span></div>
<div class="line"><a name="l00005"></a><span class="lineno"> 5</span> <span class="preprocessor">#include "<a class="code" href="physconst_8h.html">physconst.h</a>"</span></div>
<div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="preprocessor">#include "<a class="code" href="phycon_8h.html">phycon.h</a>"</span></div>
<div class="line"><a name="l00007"></a><span class="lineno"> 7</span> <span class="preprocessor">#include "<a class="code" href="thermal_8h.html">thermal.h</a>"</span></div>
<div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="preprocessor">#include "<a class="code" href="dense_8h.html">dense.h</a>"</span></div>
<div class="line"><a name="l00009"></a><span class="lineno"> 9</span> <span class="preprocessor">#include "<a class="code" href="thirdparty_8h.html">thirdparty.h</a>"</span></div>
<div class="line"><a name="l00010"></a><span class="lineno"> 10</span> <span class="preprocessor">#include "<a class="code" href="atoms_8h.html">atoms.h</a>"</span></div>
<div class="line"><a name="l00011"></a><span class="lineno"> 11</span> </div>
<div class="line"><a name="l00012"></a><span class="lineno"> 12</span> <span class="comment">/*atom_pop5 do five level atom population and cooling */</span></div>
<div class="line"><a name="l00013"></a><span class="lineno"><a class="line" href="atoms_8h.html#a0a09c35ed23f91dcf17edac9bbfb0457"> 13</a></span> <span class="keywordtype">void</span> <a class="code" href="atom__pop5_8cpp.html#ae6b795511d72534d617cf74989eb6a7d">atom_pop5</a>(</div>
<div class="line"><a name="l00014"></a><span class="lineno"> 14</span>  <span class="comment">/* vector giving statistical weights on the five levels */</span></div>
<div class="line"><a name="l00015"></a><span class="lineno"> 15</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> <a class="code" href="species2_8cpp.html#a6a53680399f0706f259b732761169657">g</a>[], </div>
<div class="line"><a name="l00016"></a><span class="lineno"> 16</span>  <span class="comment">/* vector giving the excitation energy differences of the 5 levels. The energies</span></div>
<div class="line"><a name="l00017"></a><span class="lineno"> 17</span> <span class="comment"> * are the energy in wavenumbers between adjacent levels. So EnerWN[0] is the energy</span></div>
<div class="line"><a name="l00018"></a><span class="lineno"> 18</span> <span class="comment"> * 1-2, EnerWN[1] is the energy between 2-3, etc */</span></div>
<div class="line"><a name="l00019"></a><span class="lineno"> 19</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> EnerWN[], </div>
<div class="line"><a name="l00020"></a><span class="lineno"> 20</span>  <span class="comment">/* the collision strengths for the levels */</span></div>
<div class="line"><a name="l00021"></a><span class="lineno"> 21</span>  <span class="keywordtype">double</span> cs12, </div>
<div class="line"><a name="l00022"></a><span class="lineno"> 22</span>  <span class="keywordtype">double</span> cs13, </div>
<div class="line"><a name="l00023"></a><span class="lineno"> 23</span>  <span class="keywordtype">double</span> cs14, </div>
<div class="line"><a name="l00024"></a><span class="lineno"> 24</span>  <span class="keywordtype">double</span> cs15, </div>
<div class="line"><a name="l00025"></a><span class="lineno"> 25</span>  <span class="keywordtype">double</span> cs23, </div>
<div class="line"><a name="l00026"></a><span class="lineno"> 26</span>  <span class="keywordtype">double</span> cs24, </div>
<div class="line"><a name="l00027"></a><span class="lineno"> 27</span>  <span class="keywordtype">double</span> cs25, </div>
<div class="line"><a name="l00028"></a><span class="lineno"> 28</span>  <span class="keywordtype">double</span> cs34, </div>
<div class="line"><a name="l00029"></a><span class="lineno"> 29</span>  <span class="keywordtype">double</span> cs35, </div>
<div class="line"><a name="l00030"></a><span class="lineno"> 30</span>  <span class="keywordtype">double</span> cs45, </div>
<div class="line"><a name="l00031"></a><span class="lineno"> 31</span>  <span class="comment">/* the transition probabilities between the various levels */</span></div>
<div class="line"><a name="l00032"></a><span class="lineno"> 32</span>  <span class="keywordtype">double</span> a21, </div>
<div class="line"><a name="l00033"></a><span class="lineno"> 33</span>  <span class="keywordtype">double</span> a31, </div>
<div class="line"><a name="l00034"></a><span class="lineno"> 34</span>  <span class="keywordtype">double</span> a41, </div>
<div class="line"><a name="l00035"></a><span class="lineno"> 35</span>  <span class="keywordtype">double</span> a51, </div>
<div class="line"><a name="l00036"></a><span class="lineno"> 36</span>  <span class="keywordtype">double</span> a32, </div>
<div class="line"><a name="l00037"></a><span class="lineno"> 37</span>  <span class="keywordtype">double</span> a42, </div>
<div class="line"><a name="l00038"></a><span class="lineno"> 38</span>  <span class="keywordtype">double</span> a52, </div>
<div class="line"><a name="l00039"></a><span class="lineno"> 39</span>  <span class="keywordtype">double</span> a43, </div>
<div class="line"><a name="l00040"></a><span class="lineno"> 40</span>  <span class="keywordtype">double</span> a53, </div>
<div class="line"><a name="l00041"></a><span class="lineno"> 41</span>  <span class="keywordtype">double</span> a54, </div>
<div class="line"><a name="l00042"></a><span class="lineno"> 42</span>  <span class="comment">/* the destroyed level populations (units cm-3) for the five levels */</span></div>
<div class="line"><a name="l00043"></a><span class="lineno"> 43</span>  <span class="keywordtype">double</span> p[], </div>
<div class="line"><a name="l00044"></a><span class="lineno"> 44</span>  <span class="comment">/* the total density of this ion */</span></div>
<div class="line"><a name="l00045"></a><span class="lineno"> 45</span>  <a class="code" href="cddefines_8h.html#a9aeda5ccf905f7d5241c0fe35434c820">realnum</a> <a class="code" href="abund_8cpp.html#ad2cccc5ee3b32558e3e5c1122823354f">abund</a>,</div>
<div class="line"><a name="l00046"></a><span class="lineno"> 46</span>  <span class="keywordtype">double</span> *Cooling, </div>
<div class="line"><a name="l00047"></a><span class="lineno"> 47</span>  <span class="keywordtype">double</span> *CoolingDeriv,</div>
<div class="line"><a name="l00048"></a><span class="lineno"> 48</span>  <span class="keywordtype">double</span> pump01,</div>
<div class="line"><a name="l00049"></a><span class="lineno"> 49</span>  <span class="keywordtype">double</span> pump02,</div>
<div class="line"><a name="l00050"></a><span class="lineno"> 50</span>  <span class="keywordtype">double</span> pump03,</div>
<div class="line"><a name="l00051"></a><span class="lineno"> 51</span>  <span class="keywordtype">double</span> pump04</div>
<div class="line"><a name="l00052"></a><span class="lineno"> 52</span>  )</div>
<div class="line"><a name="l00053"></a><span class="lineno"> 53</span> {</div>
<div class="line"><a name="l00054"></a><span class="lineno"> 54</span> </div>
<div class="line"><a name="l00055"></a><span class="lineno"> 55</span>  <a class="code" href="cddefines_8h.html#ada0cfa6927ba3683adc971d17643f2e5">DEBUG_ENTRY</a>( <span class="stringliteral">"atom_pop5()"</span> );</div>
<div class="line"><a name="l00056"></a><span class="lineno"> 56</span> </div>
<div class="line"><a name="l00057"></a><span class="lineno"> 57</span>  <span class="comment">/* quit if no species present */</span></div>
<div class="line"><a name="l00058"></a><span class="lineno"> 58</span>  <a class="code" href="cddefines_8h.html#a092bd7bb2cb7fd9f483b6995cee61bc0">ASSERT</a>( abund>=0. );</div>
<div class="line"><a name="l00059"></a><span class="lineno"> 59</span>  <span class="keywordflow">if</span>( abund == 0. )</div>
<div class="line"><a name="l00060"></a><span class="lineno"> 60</span>  {</div>
<div class="line"><a name="l00061"></a><span class="lineno"> 61</span>  p[0] = 0.;</div>
<div class="line"><a name="l00062"></a><span class="lineno"> 62</span>  p[1] = 0.;</div>
<div class="line"><a name="l00063"></a><span class="lineno"> 63</span>  p[2] = 0.;</div>
<div class="line"><a name="l00064"></a><span class="lineno"> 64</span>  p[3] = 0.;</div>
<div class="line"><a name="l00065"></a><span class="lineno"> 65</span>  p[4] = 0.;</div>
<div class="line"><a name="l00066"></a><span class="lineno"> 66</span>  *Cooling = 0.;</div>
<div class="line"><a name="l00067"></a><span class="lineno"> 67</span>  *CoolingDeriv = 0.;</div>
<div class="line"><a name="l00068"></a><span class="lineno"> 68</span>  <span class="keywordflow">return</span>;</div>
<div class="line"><a name="l00069"></a><span class="lineno"> 69</span>  }</div>
<div class="line"><a name="l00070"></a><span class="lineno"> 70</span> </div>
<div class="line"><a name="l00071"></a><span class="lineno"> 71</span>  <span class="comment">// tf = 1.438 / te, converts energy in wavenumbers into Boltzmann factor</span></div>
<div class="line"><a name="l00072"></a><span class="lineno"> 72</span>  <span class="keywordtype">double</span> tf = <a class="code" href="physconst_8h.html#a93f822e34d056465aab7d7a6d115709c">T1CM</a>/<a class="code" href="phycon_8cpp.html#a58e7356b28193ac2b0ee8b08f06939ba">phycon</a>.<a class="code" href="structt__phycon.html#a38484680d6075e5c0ba33fab1e514fee">te</a>;</div>
<div class="line"><a name="l00073"></a><span class="lineno"> 73</span> </div>
<div class="line"><a name="l00074"></a><span class="lineno"> 74</span>  <span class="comment">// define Boltzmann factors</span></div>
<div class="line"><a name="l00075"></a><span class="lineno"> 75</span>  <span class="keywordtype">double</span> BoltzFac[5][5];</div>
<div class="line"><a name="l00076"></a><span class="lineno"> 76</span> </div>
<div class="line"><a name="l00077"></a><span class="lineno"> 77</span>  BoltzFac[0][1] = <a class="code" href="cddefines_8h.html#a0949cd6050ce9bd2506528396ef7424f">sexp</a>(EnerWN[0]*tf);</div>
<div class="line"><a name="l00078"></a><span class="lineno"> 78</span>  BoltzFac[1][2] = <a class="code" href="cddefines_8h.html#a0949cd6050ce9bd2506528396ef7424f">sexp</a>(EnerWN[1]*tf);</div>
<div class="line"><a name="l00079"></a><span class="lineno"> 79</span>  BoltzFac[2][3] = <a class="code" href="cddefines_8h.html#a0949cd6050ce9bd2506528396ef7424f">sexp</a>(EnerWN[2]*tf);</div>
<div class="line"><a name="l00080"></a><span class="lineno"> 80</span>  BoltzFac[3][4] = <a class="code" href="cddefines_8h.html#a0949cd6050ce9bd2506528396ef7424f">sexp</a>(EnerWN[3]*tf);</div>
<div class="line"><a name="l00081"></a><span class="lineno"> 81</span>  BoltzFac[0][2] = BoltzFac[0][1]*BoltzFac[1][2];</div>
<div class="line"><a name="l00082"></a><span class="lineno"> 82</span>  BoltzFac[0][3] = BoltzFac[0][2]*BoltzFac[2][3];</div>
<div class="line"><a name="l00083"></a><span class="lineno"> 83</span>  BoltzFac[0][4] = BoltzFac[0][3]*BoltzFac[3][4];</div>
<div class="line"><a name="l00084"></a><span class="lineno"> 84</span>  BoltzFac[1][3] = BoltzFac[1][2]*BoltzFac[2][3];</div>
<div class="line"><a name="l00085"></a><span class="lineno"> 85</span>  BoltzFac[1][4] = BoltzFac[1][3]*BoltzFac[3][4];</div>
<div class="line"><a name="l00086"></a><span class="lineno"> 86</span>  BoltzFac[2][4] = BoltzFac[2][3]*BoltzFac[3][4];</div>
<div class="line"><a name="l00087"></a><span class="lineno"> 87</span> </div>
<div class="line"><a name="l00088"></a><span class="lineno"> 88</span>  <span class="comment">/* quit it highest level Boltzmann factor too large */</span></div>
<div class="line"><a name="l00089"></a><span class="lineno"> 89</span>  <span class="keywordflow">if</span>( (BoltzFac[0][4]+pump04) == 0. )</div>
<div class="line"><a name="l00090"></a><span class="lineno"> 90</span>  {</div>
<div class="line"><a name="l00091"></a><span class="lineno"> 91</span>  p[0] = 0.;</div>
<div class="line"><a name="l00092"></a><span class="lineno"> 92</span>  p[1] = 0.;</div>
<div class="line"><a name="l00093"></a><span class="lineno"> 93</span>  p[2] = 0.;</div>
<div class="line"><a name="l00094"></a><span class="lineno"> 94</span>  p[3] = 0.;</div>
<div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  p[4] = 0.;</div>
<div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  *Cooling = 0.;</div>
<div class="line"><a name="l00097"></a><span class="lineno"> 97</span>  *CoolingDeriv = 0.;</div>
<div class="line"><a name="l00098"></a><span class="lineno"> 98</span>  <span class="keywordflow">return</span>;</div>
<div class="line"><a name="l00099"></a><span class="lineno"> 99</span>  }</div>
<div class="line"><a name="l00100"></a><span class="lineno"> 100</span> </div>
<div class="line"><a name="l00101"></a><span class="lineno"> 101</span>  <span class="comment">// get collision rates, dense.cdsqte is 8.629e-6 / sqrte * eden */</span></div>
<div class="line"><a name="l00102"></a><span class="lineno"> 102</span>  <span class="comment">// rates units are s^-1</span></div>
<div class="line"><a name="l00103"></a><span class="lineno"> 103</span>  <span class="keywordtype">double</span> col[5][5];</div>
<div class="line"><a name="l00104"></a><span class="lineno"> 104</span>  col[1][0] = <a class="code" href="dense_8cpp.html#a61a760c32a26242e4eb24e230b3d5c57">dense</a>.<a class="code" href="structt__dense.html#a53085031083ac598a3b49661acd03792">cdsqte</a>*cs12/g[1];</div>
<div class="line"><a name="l00105"></a><span class="lineno"> 105</span>  col[0][1] = col[1][0]*g[1]/g[0]*BoltzFac[0][1];</div>
<div class="line"><a name="l00106"></a><span class="lineno"> 106</span> </div>
<div class="line"><a name="l00107"></a><span class="lineno"> 107</span>  col[2][0] = <a class="code" href="dense_8cpp.html#a61a760c32a26242e4eb24e230b3d5c57">dense</a>.<a class="code" href="structt__dense.html#a53085031083ac598a3b49661acd03792">cdsqte</a>*cs13/g[2];</div>
<div class="line"><a name="l00108"></a><span class="lineno"> 108</span>  col[0][2] = col[2][0]*g[2]/g[0]*BoltzFac[0][2];</div>
<div class="line"><a name="l00109"></a><span class="lineno"> 109</span> </div>
<div class="line"><a name="l00110"></a><span class="lineno"> 110</span>  col[3][0] = <a class="code" href="dense_8cpp.html#a61a760c32a26242e4eb24e230b3d5c57">dense</a>.<a class="code" href="structt__dense.html#a53085031083ac598a3b49661acd03792">cdsqte</a>*cs14/g[3];</div>
<div class="line"><a name="l00111"></a><span class="lineno"> 111</span>  col[0][3] = col[3][0]*g[3]/g[0]*BoltzFac[0][3];</div>
<div class="line"><a name="l00112"></a><span class="lineno"> 112</span> </div>
<div class="line"><a name="l00113"></a><span class="lineno"> 113</span>  col[4][0] = <a class="code" href="dense_8cpp.html#a61a760c32a26242e4eb24e230b3d5c57">dense</a>.<a class="code" href="structt__dense.html#a53085031083ac598a3b49661acd03792">cdsqte</a>*cs15/g[4];</div>
<div class="line"><a name="l00114"></a><span class="lineno"> 114</span>  col[0][4] = col[4][0]*g[4]/g[0]*BoltzFac[0][4];</div>
<div class="line"><a name="l00115"></a><span class="lineno"> 115</span> </div>
<div class="line"><a name="l00116"></a><span class="lineno"> 116</span>  col[2][1] = <a class="code" href="dense_8cpp.html#a61a760c32a26242e4eb24e230b3d5c57">dense</a>.<a class="code" href="structt__dense.html#a53085031083ac598a3b49661acd03792">cdsqte</a>*cs23/g[2];</div>
<div class="line"><a name="l00117"></a><span class="lineno"> 117</span>  col[1][2] = col[2][1]*g[2]/g[1]*BoltzFac[1][2];</div>
<div class="line"><a name="l00118"></a><span class="lineno"> 118</span> </div>
<div class="line"><a name="l00119"></a><span class="lineno"> 119</span>  col[3][1] = <a class="code" href="dense_8cpp.html#a61a760c32a26242e4eb24e230b3d5c57">dense</a>.<a class="code" href="structt__dense.html#a53085031083ac598a3b49661acd03792">cdsqte</a>*cs24/g[3];</div>
<div class="line"><a name="l00120"></a><span class="lineno"> 120</span>  col[1][3] = col[3][1]*g[3]/g[1]*BoltzFac[1][3];</div>
<div class="line"><a name="l00121"></a><span class="lineno"> 121</span> </div>
<div class="line"><a name="l00122"></a><span class="lineno"> 122</span>  col[4][1] = <a class="code" href="dense_8cpp.html#a61a760c32a26242e4eb24e230b3d5c57">dense</a>.<a class="code" href="structt__dense.html#a53085031083ac598a3b49661acd03792">cdsqte</a>*cs25/g[4];</div>
<div class="line"><a name="l00123"></a><span class="lineno"> 123</span>  col[1][4] = col[4][1]*g[4]/g[1]*BoltzFac[1][4];</div>
<div class="line"><a name="l00124"></a><span class="lineno"> 124</span> </div>
<div class="line"><a name="l00125"></a><span class="lineno"> 125</span>  col[3][2] = <a class="code" href="dense_8cpp.html#a61a760c32a26242e4eb24e230b3d5c57">dense</a>.<a class="code" href="structt__dense.html#a53085031083ac598a3b49661acd03792">cdsqte</a>*cs34/g[3];</div>
<div class="line"><a name="l00126"></a><span class="lineno"> 126</span>  col[2][3] = col[3][2]*g[3]/g[2]*BoltzFac[2][3];</div>
<div class="line"><a name="l00127"></a><span class="lineno"> 127</span> </div>
<div class="line"><a name="l00128"></a><span class="lineno"> 128</span>  col[4][2] = <a class="code" href="dense_8cpp.html#a61a760c32a26242e4eb24e230b3d5c57">dense</a>.<a class="code" href="structt__dense.html#a53085031083ac598a3b49661acd03792">cdsqte</a>*cs35/g[4];</div>
<div class="line"><a name="l00129"></a><span class="lineno"> 129</span>  col[2][4] = col[4][2]*g[4]/g[2]*BoltzFac[2][4];</div>
<div class="line"><a name="l00130"></a><span class="lineno"> 130</span> </div>
<div class="line"><a name="l00131"></a><span class="lineno"> 131</span>  col[4][3] = <a class="code" href="dense_8cpp.html#a61a760c32a26242e4eb24e230b3d5c57">dense</a>.<a class="code" href="structt__dense.html#a53085031083ac598a3b49661acd03792">cdsqte</a>*cs45/g[4];</div>
<div class="line"><a name="l00132"></a><span class="lineno"> 132</span>  col[3][4] = col[4][3]*g[4]/g[3]*BoltzFac[3][4];</div>
<div class="line"><a name="l00133"></a><span class="lineno"> 133</span> </div>
<div class="line"><a name="l00134"></a><span class="lineno"> 134</span>  <span class="keywordtype">double</span> <a class="code" href="atom__feii_8cpp.html#a48d7a59fc34dbf78eb423637699dd593">amat</a>[5][5], bvec[5];</div>
<div class="line"><a name="l00135"></a><span class="lineno"> 135</span>  <span class="comment">// homogeneous matrix - no source or sink - use conservation at 5th equation</span></div>
<div class="line"><a name="l00136"></a><span class="lineno"> 136</span>  <span class="keywordflow">for</span>( <span class="keywordtype">long</span> i=0; i<5; ++i )</div>
<div class="line"><a name="l00137"></a><span class="lineno"> 137</span>  {</div>
<div class="line"><a name="l00138"></a><span class="lineno"> 138</span>  amat[i][4] = 1.;</div>
<div class="line"><a name="l00139"></a><span class="lineno"> 139</span>  bvec[i] = 0.;</div>
<div class="line"><a name="l00140"></a><span class="lineno"> 140</span>  }</div>
<div class="line"><a name="l00141"></a><span class="lineno"> 141</span>  bvec[4] = <a class="code" href="abund_8cpp.html#ad2cccc5ee3b32558e3e5c1122823354f">abund</a>;</div>
<div class="line"><a name="l00142"></a><span class="lineno"> 142</span> </div>
<div class="line"><a name="l00143"></a><span class="lineno"> 143</span>  <span class="comment">/* level one balance equation */</span></div>
<div class="line"><a name="l00144"></a><span class="lineno"> 144</span>  amat[0][0] = col[0][1] + col[0][2] + col[0][3] + col[0][4] +pump01+pump02+pump03+pump04;</div>
<div class="line"><a name="l00145"></a><span class="lineno"> 145</span>  amat[1][0] = -a21 - col[1][0];</div>
<div class="line"><a name="l00146"></a><span class="lineno"> 146</span>  amat[2][0] = -a31 - col[2][0];</div>
<div class="line"><a name="l00147"></a><span class="lineno"> 147</span>  amat[3][0] = -a41 - col[3][0];</div>
<div class="line"><a name="l00148"></a><span class="lineno"> 148</span>  amat[4][0] = -a51 - col[4][0];</div>
<div class="line"><a name="l00149"></a><span class="lineno"> 149</span> </div>
<div class="line"><a name="l00150"></a><span class="lineno"> 150</span>  <span class="comment">/* level two balance equation */</span></div>
<div class="line"><a name="l00151"></a><span class="lineno"> 151</span>  amat[0][1] = -col[0][1] - pump01;</div>
<div class="line"><a name="l00152"></a><span class="lineno"> 152</span>  amat[1][1] = col[1][0] + a21 + col[1][2] + col[1][3] + col[1][4];</div>
<div class="line"><a name="l00153"></a><span class="lineno"> 153</span>  amat[2][1] = -col[2][1] - a32;</div>
<div class="line"><a name="l00154"></a><span class="lineno"> 154</span>  amat[3][1] = -col[3][1] - a42;</div>
<div class="line"><a name="l00155"></a><span class="lineno"> 155</span>  amat[4][1] = -col[4][1] - a52;</div>
<div class="line"><a name="l00156"></a><span class="lineno"> 156</span> </div>
<div class="line"><a name="l00157"></a><span class="lineno"> 157</span>  <span class="comment">/* level three balance equation */</span></div>
<div class="line"><a name="l00158"></a><span class="lineno"> 158</span>  amat[0][2] = -col[0][2] - pump02;</div>
<div class="line"><a name="l00159"></a><span class="lineno"> 159</span>  amat[1][2] = -col[1][2];</div>
<div class="line"><a name="l00160"></a><span class="lineno"> 160</span>  amat[2][2] = a31 + a32 + col[2][0] + col[2][1] + col[2][3] + col[2][4];</div>
<div class="line"><a name="l00161"></a><span class="lineno"> 161</span>  amat[3][2] = -col[3][2] - a43;</div>
<div class="line"><a name="l00162"></a><span class="lineno"> 162</span>  amat[4][2] = -col[4][2] - a53;</div>
<div class="line"><a name="l00163"></a><span class="lineno"> 163</span> </div>
<div class="line"><a name="l00164"></a><span class="lineno"> 164</span>  <span class="comment">/* level four balance equation */</span></div>
<div class="line"><a name="l00165"></a><span class="lineno"> 165</span>  amat[0][3] = -col[0][3] - pump03;</div>
<div class="line"><a name="l00166"></a><span class="lineno"> 166</span>  amat[1][3] = -col[1][3];</div>
<div class="line"><a name="l00167"></a><span class="lineno"> 167</span>  amat[2][3] = -col[2][3];</div>
<div class="line"><a name="l00168"></a><span class="lineno"> 168</span>  amat[3][3] = a41 + col[3][0] + a42 + col[3][1] + a43 + col[3][2] + col[3][4];</div>
<div class="line"><a name="l00169"></a><span class="lineno"> 169</span>  amat[4][3] = -col[4][3] - a54;</div>
<div class="line"><a name="l00170"></a><span class="lineno"> 170</span> </div>
<div class="line"><a name="l00171"></a><span class="lineno"> 171</span> <span class="preprocessor"># if 0</span></div>
<div class="line"><a name="l00172"></a><span class="lineno"> 172</span> <span class="preprocessor"></span> <span class="comment">// is it necessary to precondition the vars?</span></div>
<div class="line"><a name="l00173"></a><span class="lineno"> 173</span>  <span class="comment">/* divide both sides of equation by largest number to stop overflow */</span></div>
<div class="line"><a name="l00174"></a><span class="lineno"> 174</span>  <span class="keywordtype">double</span> dmax = -1e0;</div>
<div class="line"><a name="l00175"></a><span class="lineno"> 175</span>  <span class="keywordflow">for</span>( i=0; i < 6; i++ )</div>
<div class="line"><a name="l00176"></a><span class="lineno"> 176</span>  <span class="keywordflow">for</span>( j=0; j < 5; j++ )</div>
<div class="line"><a name="l00177"></a><span class="lineno"> 177</span>  dmax = <a class="code" href="cddefines_8h.html#aa0c19e95ba3c1d759e15a79a136ce3c6">MAX2</a>(zz[i][j],dmax);</div>
<div class="line"><a name="l00178"></a><span class="lineno"> 178</span> </div>
<div class="line"><a name="l00179"></a><span class="lineno"> 179</span>  <span class="keywordflow">for</span>( i=0; i < 6; i++ )</div>
<div class="line"><a name="l00180"></a><span class="lineno"> 180</span>  <span class="keywordflow">for</span>( j=0; j < 5; j++ )</div>
<div class="line"><a name="l00181"></a><span class="lineno"> 181</span>  zz[i][j] /= dmax;</div>
<div class="line"><a name="l00182"></a><span class="lineno"> 182</span> <span class="preprocessor"># endif</span></div>
<div class="line"><a name="l00183"></a><span class="lineno"> 183</span> <span class="preprocessor"></span></div>
<div class="line"><a name="l00184"></a><span class="lineno"> 184</span>  int32 ipiv[5], ner=0;</div>
<div class="line"><a name="l00185"></a><span class="lineno"> 185</span>  <span class="comment">/* solve matrix */</span></div>
<div class="line"><a name="l00186"></a><span class="lineno"> 186</span>  <a class="code" href="thirdparty_8h.html#adc100df3d6f1a56dade362953e84e726">getrf_wrapper</a>(5,5,(<span class="keywordtype">double</span>*)amat,5,ipiv,&ner);</div>
<div class="line"><a name="l00187"></a><span class="lineno"> 187</span>  <a class="code" href="thirdparty_8h.html#a539cec613446b3d750ce185b13b5738a">getrs_wrapper</a>(<span class="charliteral">'N'</span>,5,1,(<span class="keywordtype">double</span>*)amat,5,ipiv,bvec,5,&ner);</div>
<div class="line"><a name="l00188"></a><span class="lineno"> 188</span> </div>
<div class="line"><a name="l00189"></a><span class="lineno"> 189</span>  <span class="keywordflow">if</span>( ner != 0 )</div>
<div class="line"><a name="l00190"></a><span class="lineno"> 190</span>  {</div>
<div class="line"><a name="l00191"></a><span class="lineno"> 191</span>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">"DISASTER PROBLEM atom_pop5: dgetrs finds singular or ill-conditioned matrix\n"</span> );</div>
<div class="line"><a name="l00192"></a><span class="lineno"> 192</span>  <a class="code" href="cddefines_8h.html#a937636e18a3a517d1aeb2c111adb5160">cdEXIT</a>(<a class="code" href="cddefines_8h.html#a73efe787c131b385070f25d18b7c9aa4">EXIT_FAILURE</a>);</div>
<div class="line"><a name="l00193"></a><span class="lineno"> 193</span>  }</div>
<div class="line"><a name="l00194"></a><span class="lineno"> 194</span> </div>
<div class="line"><a name="l00195"></a><span class="lineno"> 195</span>  <span class="comment">/* p(5) was very slightly negative (1e-40) for SII in dqher.in - highest</span></div>
<div class="line"><a name="l00196"></a><span class="lineno"> 196</span> <span class="comment"> * level has smallest excitation rates may be closest to ill conditioned*/</span></div>
<div class="line"><a name="l00197"></a><span class="lineno"> 197</span>  p[1] = <a class="code" href="cddefines_8h.html#aa0c19e95ba3c1d759e15a79a136ce3c6">MAX2</a>(0.e0,bvec[1]);</div>
<div class="line"><a name="l00198"></a><span class="lineno"> 198</span>  p[2] = <a class="code" href="cddefines_8h.html#aa0c19e95ba3c1d759e15a79a136ce3c6">MAX2</a>(0.e0,bvec[2]);</div>
<div class="line"><a name="l00199"></a><span class="lineno"> 199</span>  p[3] = <a class="code" href="cddefines_8h.html#aa0c19e95ba3c1d759e15a79a136ce3c6">MAX2</a>(0.e0,bvec[3]);</div>
<div class="line"><a name="l00200"></a><span class="lineno"> 200</span>  p[4] = <a class="code" href="cddefines_8h.html#aa0c19e95ba3c1d759e15a79a136ce3c6">MAX2</a>(0.e0,bvec[4]);</div>
<div class="line"><a name="l00201"></a><span class="lineno"> 201</span>  <span class="comment">// this should be majority and so numerically more</span></div>
<div class="line"><a name="l00202"></a><span class="lineno"> 202</span>  p[0] = abund - p[1] - p[2] - p[3] - p[4];</div>
<div class="line"><a name="l00203"></a><span class="lineno"> 203</span> </div>
<div class="line"><a name="l00204"></a><span class="lineno"> 204</span>  <span class="comment">// level energies in ergs, needed for energy exchange</span></div>
<div class="line"><a name="l00205"></a><span class="lineno"> 205</span>  <span class="keywordtype">double</span> Erg[5] , EnergyKelvin[5];</div>
<div class="line"><a name="l00206"></a><span class="lineno"> 206</span>  Erg[0] = 0.;</div>
<div class="line"><a name="l00207"></a><span class="lineno"> 207</span>  EnergyKelvin[0] = 0.;</div>
<div class="line"><a name="l00208"></a><span class="lineno"> 208</span>  <span class="keywordflow">for</span>( <span class="keywordtype">long</span> i=1; i<5; ++i )</div>
<div class="line"><a name="l00209"></a><span class="lineno"> 209</span>  {</div>
<div class="line"><a name="l00210"></a><span class="lineno"> 210</span>  Erg[i] = Erg[i-1] + EnerWN[i-1]*<a class="code" href="physconst_8h.html#a4afe06a73410b1bc5c10c8431e8df368">ERG1CM</a>;</div>
<div class="line"><a name="l00211"></a><span class="lineno"> 211</span>  EnergyKelvin[i] = EnergyKelvin[i-1] + EnerWN[i-1]*<a class="code" href="physconst_8h.html#a93f822e34d056465aab7d7a6d115709c">T1CM</a>;</div>
<div class="line"><a name="l00212"></a><span class="lineno"> 212</span>  }</div>
<div class="line"><a name="l00213"></a><span class="lineno"> 213</span> </div>
<div class="line"><a name="l00214"></a><span class="lineno"> 214</span>  *Cooling = 0.;</div>
<div class="line"><a name="l00215"></a><span class="lineno"> 215</span>  *CoolingDeriv = 0.;</div>
<div class="line"><a name="l00216"></a><span class="lineno"> 216</span>  <span class="keywordflow">for</span>( <span class="keywordtype">long</span> ihi=1; ihi<5; ++ihi )</div>
<div class="line"><a name="l00217"></a><span class="lineno"> 217</span>  {</div>
<div class="line"><a name="l00218"></a><span class="lineno"> 218</span>  <span class="keywordflow">for</span>( <span class="keywordtype">long</span> ilo=0; ilo<ihi; ++ilo )</div>
<div class="line"><a name="l00219"></a><span class="lineno"> 219</span>  {</div>
<div class="line"><a name="l00220"></a><span class="lineno"> 220</span>  <span class="keywordtype">double</span> CoolOne = (p[ilo]*col[ilo][ihi] - p[ihi]*col[ihi][ilo])* </div>
<div class="line"><a name="l00221"></a><span class="lineno"> 221</span>  (Erg[ihi]-Erg[ilo]);</div>
<div class="line"><a name="l00222"></a><span class="lineno"> 222</span> </div>
<div class="line"><a name="l00223"></a><span class="lineno"> 223</span>  *Cooling += CoolOne;</div>
<div class="line"><a name="l00224"></a><span class="lineno"> 224</span>  *CoolingDeriv += CoolOne*( EnergyKelvin[ihi]*<a class="code" href="thermal_8cpp.html#aa9fb6393d963cd8537ddeb8735d4eb6a">thermal</a>.<a class="code" href="structt__thermal.html#a49b1a07a4f89e537b0b67aa210bd742a">tsq1</a> - <a class="code" href="thermal_8cpp.html#aa9fb6393d963cd8537ddeb8735d4eb6a">thermal</a>.<a class="code" href="structt__thermal.html#ae2dbe2550225f3c34a52c8d4b00f1a02">halfte</a> );</div>
<div class="line"><a name="l00225"></a><span class="lineno"> 225</span>  }</div>
<div class="line"><a name="l00226"></a><span class="lineno"> 226</span>  }</div>
<div class="line"><a name="l00227"></a><span class="lineno"> 227</span>  </div>
<div class="line"><a name="l00228"></a><span class="lineno"> 228</span>  <span class="keywordflow">return</span>;</div>
<div class="line"><a name="l00229"></a><span class="lineno"> 229</span> }</div>
<div class="ttc" id="structt__dense_html_a53085031083ac598a3b49661acd03792"><div class="ttname"><a href="structt__dense.html#a53085031083ac598a3b49661acd03792">t_dense::cdsqte</a></div><div class="ttdeci">double cdsqte</div><div class="ttdef"><b>Definition:</b> <a href="dense_8h_source.html#l00235">dense.h:235</a></div></div>
<div class="ttc" id="atoms_8h_html"><div class="ttname"><a href="atoms_8h.html">atoms.h</a></div></div>
<div class="ttc" id="thermal_8cpp_html_aa9fb6393d963cd8537ddeb8735d4eb6a"><div class="ttname"><a href="thermal_8cpp.html#aa9fb6393d963cd8537ddeb8735d4eb6a">thermal</a></div><div class="ttdeci">t_thermal thermal</div><div class="ttdef"><b>Definition:</b> <a href="thermal_8cpp_source.html#l00005">thermal.cpp:5</a></div></div>
<div class="ttc" id="thirdparty_8h_html_a539cec613446b3d750ce185b13b5738a"><div class="ttname"><a href="thirdparty_8h.html#a539cec613446b3d750ce185b13b5738a">getrs_wrapper</a></div><div class="ttdeci">void getrs_wrapper(char trans, long N, long nrhs, double *A, long lda, int32 *ipiv, double *B, long ldb, int32 *info)</div><div class="ttdef"><b>Definition:</b> <a href="thirdparty__lapack_8cpp_source.html#l00069">thirdparty_lapack.cpp:69</a></div></div>
<div class="ttc" id="cddefines_8h_html_aa0c19e95ba3c1d759e15a79a136ce3c6"><div class="ttname"><a href="cddefines_8h.html#aa0c19e95ba3c1d759e15a79a136ce3c6">MAX2</a></div><div class="ttdeci">#define MAX2</div><div class="ttdef"><b>Definition:</b> <a href="cddefines_8h_source.html#l00786">cddefines.h:786</a></div></div>
<div class="ttc" id="phycon_8cpp_html_a58e7356b28193ac2b0ee8b08f06939ba"><div class="ttname"><a href="phycon_8cpp.html#a58e7356b28193ac2b0ee8b08f06939ba">phycon</a></div><div class="ttdeci">t_phycon phycon</div><div class="ttdef"><b>Definition:</b> <a href="phycon_8cpp_source.html#l00006">phycon.cpp:6</a></div></div>
<div class="ttc" id="dense_8cpp_html_a61a760c32a26242e4eb24e230b3d5c57"><div class="ttname"><a href="dense_8cpp.html#a61a760c32a26242e4eb24e230b3d5c57">dense</a></div><div class="ttdeci">t_dense dense</div><div class="ttdef"><b>Definition:</b> <a href="dense_8cpp_source.html#l00024">dense.cpp:24</a></div></div>
<div class="ttc" id="phycon_8h_html"><div class="ttname"><a href="phycon_8h.html">phycon.h</a></div></div>
<div class="ttc" id="thermal_8h_html"><div class="ttname"><a href="thermal_8h.html">thermal.h</a></div></div>
<div class="ttc" id="cddefines_8h_html_a0949cd6050ce9bd2506528396ef7424f"><div class="ttname"><a href="cddefines_8h.html#a0949cd6050ce9bd2506528396ef7424f">sexp</a></div><div class="ttdeci">sys_float sexp(sys_float x)</div><div class="ttdef"><b>Definition:</b> <a href="service_8cpp_source.html#l00914">service.cpp:914</a></div></div>
<div class="ttc" id="species2_8cpp_html_a6a53680399f0706f259b732761169657"><div class="ttname"><a href="species2_8cpp.html#a6a53680399f0706f259b732761169657">g</a></div><div class="ttdeci">static double * g</div><div class="ttdef"><b>Definition:</b> <a href="species2_8cpp_source.html#l00028">species2.cpp:28</a></div></div>
<div class="ttc" id="cddefines_8cpp_html_a94bc49f6569c17375ac7fbf6a66f6dd0"><div class="ttname"><a href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a></div><div class="ttdeci">FILE * ioQQQ</div><div class="ttdef"><b>Definition:</b> <a href="cddefines_8cpp_source.html#l00007">cddefines.cpp:7</a></div></div>
<div class="ttc" id="dense_8h_html"><div class="ttname"><a href="dense_8h.html">dense.h</a></div></div>
<div class="ttc" id="abund_8cpp_html_ad2cccc5ee3b32558e3e5c1122823354f"><div class="ttname"><a href="abund_8cpp.html#ad2cccc5ee3b32558e3e5c1122823354f">abund</a></div><div class="ttdeci">t_abund abund</div><div class="ttdef"><b>Definition:</b> <a href="abund_8cpp_source.html#l00005">abund.cpp:5</a></div></div>
<div class="ttc" id="thirdparty_8h_html_adc100df3d6f1a56dade362953e84e726"><div class="ttname"><a href="thirdparty_8h.html#adc100df3d6f1a56dade362953e84e726">getrf_wrapper</a></div><div class="ttdeci">void getrf_wrapper(long M, long N, double *A, long lda, int32 *ipiv, int32 *info)</div><div class="ttdef"><b>Definition:</b> <a href="thirdparty__lapack_8cpp_source.html#l00047">thirdparty_lapack.cpp:47</a></div></div>
<div class="ttc" id="physconst_8h_html_a4afe06a73410b1bc5c10c8431e8df368"><div class="ttname"><a href="physconst_8h.html#a4afe06a73410b1bc5c10c8431e8df368">ERG1CM</a></div><div class="ttdeci">const double ERG1CM</div><div class="ttdef"><b>Definition:</b> <a href="physconst_8h_source.html#l00164">physconst.h:164</a></div></div>
<div class="ttc" id="atom__pop5_8cpp_html_ae6b795511d72534d617cf74989eb6a7d"><div class="ttname"><a href="atom__pop5_8cpp.html#ae6b795511d72534d617cf74989eb6a7d">atom_pop5</a></div><div class="ttdeci">void atom_pop5(const double g[], const double EnerWN[], double cs12, double cs13, double cs14, double cs15, double cs23, double cs24, double cs25, double cs34, double cs35, double cs45, double a21, double a31, double a41, double a51, double a32, double a42, double a52, double a43, double a53, double a54, double p[], realnum abund, double *Cooling, double *CoolingDeriv, double pump01, double pump02, double pump03, double pump04)</div><div class="ttdef"><b>Definition:</b> <a href="atom__pop5_8cpp_source.html#l00013">atom_pop5.cpp:13</a></div></div>
<div class="ttc" id="physconst_8h_html_a93f822e34d056465aab7d7a6d115709c"><div class="ttname"><a href="physconst_8h.html#a93f822e34d056465aab7d7a6d115709c">T1CM</a></div><div class="ttdeci">const double T1CM</div><div class="ttdef"><b>Definition:</b> <a href="physconst_8h_source.html#l00167">physconst.h:167</a></div></div>
<div class="ttc" id="cddefines_8h_html_a9aeda5ccf905f7d5241c0fe35434c820"><div class="ttname"><a href="cddefines_8h.html#a9aeda5ccf905f7d5241c0fe35434c820">realnum</a></div><div class="ttdeci">float realnum</div><div class="ttdef"><b>Definition:</b> <a href="cddefines_8h_source.html#l00107">cddefines.h:107</a></div></div>
<div class="ttc" id="cddefines_8h_html_a73efe787c131b385070f25d18b7c9aa4"><div class="ttname"><a href="cddefines_8h.html#a73efe787c131b385070f25d18b7c9aa4">EXIT_FAILURE</a></div><div class="ttdeci">#define EXIT_FAILURE</div><div class="ttdef"><b>Definition:</b> <a href="cddefines_8h_source.html#l00144">cddefines.h:144</a></div></div>
<div class="ttc" id="cddefines_8h_html_a937636e18a3a517d1aeb2c111adb5160"><div class="ttname"><a href="cddefines_8h.html#a937636e18a3a517d1aeb2c111adb5160">cdEXIT</a></div><div class="ttdeci">#define cdEXIT(FAIL)</div><div class="ttdef"><b>Definition:</b> <a href="cddefines_8h_source.html#l00438">cddefines.h:438</a></div></div>
<div class="ttc" id="atom__feii_8cpp_html_a48d7a59fc34dbf78eb423637699dd593"><div class="ttname"><a href="atom__feii_8cpp.html#a48d7a59fc34dbf78eb423637699dd593">amat</a></div><div class="ttdeci">static double * amat</div><div class="ttdef"><b>Definition:</b> <a href="atom__feii_8cpp_source.html#l00167">atom_feii.cpp:167</a></div></div>
<div class="ttc" id="thirdparty_8h_html"><div class="ttname"><a href="thirdparty_8h.html">thirdparty.h</a></div></div>
<div class="ttc" id="cddefines_8h_html_a092bd7bb2cb7fd9f483b6995cee61bc0"><div class="ttname"><a href="cddefines_8h.html#a092bd7bb2cb7fd9f483b6995cee61bc0">ASSERT</a></div><div class="ttdeci">#define ASSERT(exp)</div><div class="ttdef"><b>Definition:</b> <a href="cddefines_8h_source.html#l00582">cddefines.h:582</a></div></div>
<div class="ttc" id="cddefines_8h_html"><div class="ttname"><a href="cddefines_8h.html">cddefines.h</a></div></div>
<div class="ttc" id="cddefines_8h_html_ada0cfa6927ba3683adc971d17643f2e5"><div class="ttname"><a href="cddefines_8h.html#ada0cfa6927ba3683adc971d17643f2e5">DEBUG_ENTRY</a></div><div class="ttdeci">#define DEBUG_ENTRY(funcname)</div><div class="ttdef"><b>Definition:</b> <a href="cddefines_8h_source.html#l00688">cddefines.h:688</a></div></div>
<div class="ttc" id="structt__thermal_html_ae2dbe2550225f3c34a52c8d4b00f1a02"><div class="ttname"><a href="structt__thermal.html#ae2dbe2550225f3c34a52c8d4b00f1a02">t_thermal::halfte</a></div><div class="ttdeci">double halfte</div><div class="ttdef"><b>Definition:</b> <a href="thermal_8h_source.html#l00122">thermal.h:122</a></div></div>
<div class="ttc" id="physconst_8h_html"><div class="ttname"><a href="physconst_8h.html">physconst.h</a></div></div>
<div class="ttc" id="structt__phycon_html_a38484680d6075e5c0ba33fab1e514fee"><div class="ttname"><a href="structt__phycon.html#a38484680d6075e5c0ba33fab1e514fee">t_phycon::te</a></div><div class="ttdeci">double te</div><div class="ttdef"><b>Definition:</b> <a href="phycon_8h_source.html#l00011">phycon.h:11</a></div></div>
<div class="ttc" id="structt__thermal_html_a49b1a07a4f89e537b0b67aa210bd742a"><div class="ttname"><a href="structt__thermal.html#a49b1a07a4f89e537b0b67aa210bd742a">t_thermal::tsq1</a></div><div class="ttdeci">double tsq1</div><div class="ttdef"><b>Definition:</b> <a href="thermal_8h_source.html#l00122">thermal.h:122</a></div></div>
</div><!-- fragment --></div><!-- contents -->
<!-- start footer part -->
<hr class="footer"/><address class="footer"><small>
Generated on Sat Oct 19 2013 10:43:44 for cloudy by  <a href="http://www.doxygen.org/index.html">
<img class="footer" src="doxygen.png" alt="doxygen"/>
</a> 1.8.5
</small></address>
</body>
</html>
