<!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/atmdat_dielrec_fe.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">atmdat_dielrec_fe.cpp</div> </div> </div><!--header--> <div class="contents"> <a href="atmdat__dielrec__fe_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">/*atmdat_dielrec_fe Dielectronic recombination rates for Fe from Arnaud & Raymond 1992 */</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="atmdat_8h.html">atmdat.h</a>"</span></div> <div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="preprocessor">#include "<a class="code" href="physconst_8h.html">physconst.h</a>"</span></div> <div class="line"><a name="l00007"></a><span class="lineno"> 7</span> </div> <div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="comment">/*atmdat_dielrec_fe Dielectronic recombination rates for Fe from Arnaud & Raymond 1992 */</span></div> <div class="line"><a name="l00009"></a><span class="lineno"><a class="line" href="atmdat__dielrec__fe_8cpp.html#a3283acfcfc73aca51218308909b4091c"> 9</a></span> <span class="keywordtype">double</span> <a class="code" href="atmdat_8h.html#a3283acfcfc73aca51218308909b4091c">atmdat_dielrec_fe</a>(</div> <div class="line"><a name="l00010"></a><span class="lineno"> 10</span>  <span class="comment">// ion on physics scale</span></div> <div class="line"><a name="l00011"></a><span class="lineno"> 11</span>  <span class="keywordtype">long</span> <span class="keywordtype">int</span> ion,</div> <div class="line"><a name="l00012"></a><span class="lineno"> 12</span>  <span class="comment">// electron temperture</span></div> <div class="line"><a name="l00013"></a><span class="lineno"> 13</span>  <span class="keywordtype">double</span> t)</div> <div class="line"><a name="l00014"></a><span class="lineno"> 14</span> {</div> <div class="line"><a name="l00015"></a><span class="lineno"> 15</span>  <span class="keyword">static</span> <span class="keyword">const</span> <span class="keywordtype">double</span> dfe[25][8] = {</div> <div class="line"><a name="l00016"></a><span class="lineno"> 16</span>  {5.120e00,1.29e01,0.00e00,0.00e00,2.20e-04,1.00e-04,0.00e00,0.00e00},</div> <div class="line"><a name="l00017"></a><span class="lineno"> 17</span>  {1.670e01,3.14e01,0.00e00,0.00e00,2.30e-03,2.70e-03,0.00e00,0.00e00},</div> <div class="line"><a name="l00018"></a><span class="lineno"> 18</span>  {2.860e01,5.21e01,0.00e00,0.00e00,1.50e-02,4.70e-03,0.00e00,0.00e00},</div> <div class="line"><a name="l00019"></a><span class="lineno"> 19</span>  {3.730e01,6.74e01,0.00e00,0.00e00,3.80e-02,1.60e-02,0.00e00,0.00e00},</div> <div class="line"><a name="l00020"></a><span class="lineno"> 20</span>  {5.420e01,1.00e02,0.00e00,0.00e00,8.00e-02,2.40e-02,0.00e00,0.00e00},</div> <div class="line"><a name="l00021"></a><span class="lineno"> 21</span>  {4.550e01,3.60e02,0.00e00,0.00e00,9.20e-02,4.10e-02,0.00e00,0.00e00},</div> <div class="line"><a name="l00022"></a><span class="lineno"> 22</span>  {6.670e01,1.23e02,0.00e00,0.00e00,1.60e-01,3.60e-02,0.00e00,0.00e00},</div> <div class="line"><a name="l00023"></a><span class="lineno"> 23</span>  {6.610e01,1.29e02,0.00e00,0.00e00,1.80e-01,7.00e-02,0.00e00,0.00e00},</div> <div class="line"><a name="l00024"></a><span class="lineno"> 24</span>  {2.160e01,1.36e02,0.00e00,0.00e00,1.40e-01,2.60e-01,0.00e00,0.00e00},</div> <div class="line"><a name="l00025"></a><span class="lineno"> 25</span>  {2.220e01,1.44e02,0.00e00,0.00e00,1.00e-01,2.80e-01,0.00e00,0.00e00},</div> <div class="line"><a name="l00026"></a><span class="lineno"> 26</span>  {5.960e01,3.62e02,0.00e00,0.00e00,2.25e-01,2.31e-01,0.00e00,0.00e00},</div> <div class="line"><a name="l00027"></a><span class="lineno"> 27</span>  {7.500e01,2.05e02,0.00e00,0.00e00,2.40e-01,1.70e-01,0.00e00,0.00e00},</div> <div class="line"><a name="l00028"></a><span class="lineno"> 28</span>  {3.630e01,1.93e02,0.00e00,0.00e00,2.60e-01,1.60e-01,0.00e00,0.00e00},</div> <div class="line"><a name="l00029"></a><span class="lineno"> 29</span>  {3.940e01,1.98e02,0.00e00,0.00e00,1.90e-01,9.00e-02,0.00e00,0.00e00},</div> <div class="line"><a name="l00030"></a><span class="lineno"> 30</span>  {2.460e01,2.48e02,5.60e02,0.00e00,1.20e-01,1.20e-01,6.00e-01,0.00e00},</div> <div class="line"><a name="l00031"></a><span class="lineno"> 31</span>  {5.600e02,0.00e00,0.00e00,0.00e00,1.23e00,0.00e00,0.00e00,0.00e00},</div> <div class="line"><a name="l00032"></a><span class="lineno"> 32</span>  {2.250e01,1.17e02,3.41e02,6.83e02,2.53e-03,3.36e-02,1.81e-01,1.92e00},</div> <div class="line"><a name="l00033"></a><span class="lineno"> 33</span>  {1.620e01,9.60e01,3.30e02,7.29e02,5.67e-03,7.82e-02,3.18e-02,1.26e00},</div> <div class="line"><a name="l00034"></a><span class="lineno"> 34</span>  {2.370e01,8.51e01,3.29e02,7.87e02,1.60e-02,7.17e-02,9.06e-02,7.39e-01},</div> <div class="line"><a name="l00035"></a><span class="lineno"> 35</span>  {1.320e01,6.66e01,2.97e02,7.14e02,1.85e-02,9.53e-02,7.90e-02,1.23e00},</div> <div class="line"><a name="l00036"></a><span class="lineno"> 36</span>  {3.910e01,8.03e01,3.92e02,9.19e02,9.20e-04,1.29e-01,1.92e-01,9.12e-01},</div> <div class="line"><a name="l00037"></a><span class="lineno"> 37</span>  {7.320e01,3.16e02,8.77e02,0.00e00,1.31e-01,8.49e-02,6.13e-01,0.00e00},</div> <div class="line"><a name="l00038"></a><span class="lineno"> 38</span>  {1.000e-01,3.62e01,3.06e02,9.28e02,1.10e-02,4.88e-02,8.01e-02,5.29e-01},</div> <div class="line"><a name="l00039"></a><span class="lineno"> 39</span>  {4.625e03,6.00e03,0.00e00,0.00e00,2.56e-01,4.52e-01,0.00e00,0.00e00},</div> <div class="line"><a name="l00040"></a><span class="lineno"> 40</span>  {5.300e03,0.00e00,0.00e00,0.00e00,4.30e-01,0.00e00,0.00e00,0.00e00}</div> <div class="line"><a name="l00041"></a><span class="lineno"> 41</span>  };</div> <div class="line"><a name="l00042"></a><span class="lineno"> 42</span> </div> <div class="line"><a name="l00043"></a><span class="lineno"> 43</span>  <span class="keywordtype">double</span> rate, te;</div> <div class="line"><a name="l00044"></a><span class="lineno"> 44</span> </div> <div class="line"><a name="l00045"></a><span class="lineno"> 45</span>  <a class="code" href="cddefines_8h.html#ada0cfa6927ba3683adc971d17643f2e5">DEBUG_ENTRY</a>( <span class="stringliteral">"atmdat_dielrec_fe()"</span> );</div> <div class="line"><a name="l00046"></a><span class="lineno"> 46</span>  <span class="comment">/*Dielectronic recombination rates for Fe from </span></div> <div class="line"><a name="l00047"></a><span class="lineno"> 47</span> <span class="comment"> * >>refer Fe rec Arnaud, M. & Raymond, J. 1992, ApJ, 398, 394 */</span></div> <div class="line"><a name="l00048"></a><span class="lineno"> 48</span> </div> <div class="line"><a name="l00049"></a><span class="lineno"> 49</span>  <span class="comment">/* ion - spectroscopic symbol of final ion</span></div> <div class="line"><a name="l00050"></a><span class="lineno"> 50</span> <span class="comment"> * t - temperature, K</span></div> <div class="line"><a name="l00051"></a><span class="lineno"> 51</span> <span class="comment"> * d - rate coefficient, cm^3 s^-1 */</span></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>  <span class="keywordflow">if</span>( ion > 26 )</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>  fprintf( <a class="code" href="cddefines_8cpp.html#a94bc49f6569c17375ac7fbf6a66f6dd0">ioQQQ</a>, <span class="stringliteral">" atmdat_dielrec_fe invalid ion%10ld\n"</span>, ion );</div> <div class="line"><a name="l00056"></a><span class="lineno"> 56</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="l00057"></a><span class="lineno"> 57</span>  }</div> <div class="line"><a name="l00058"></a><span class="lineno"> 58</span> </div> <div class="line"><a name="l00059"></a><span class="lineno"> 59</span>  <span class="keywordflow">else</span> <span class="keywordflow">if</span>( ion == 26 )</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>  <span class="comment">/* d is the rate */</span></div> <div class="line"><a name="l00062"></a><span class="lineno"> 62</span>  rate = 0.0;</div> <div class="line"><a name="l00063"></a><span class="lineno"> 63</span>  }</div> <div class="line"><a name="l00064"></a><span class="lineno"> 64</span> </div> <div class="line"><a name="l00065"></a><span class="lineno"> 65</span>  <span class="keywordflow">else</span></div> <div class="line"><a name="l00066"></a><span class="lineno"> 66</span>  {</div> <div class="line"><a name="l00067"></a><span class="lineno"> 67</span>  te = t*<a class="code" href="physconst_8h.html#ab1891e433b4c259d18bc1c800ea15316">EVRYD</a>/<a class="code" href="physconst_8h.html#a0db3a0b9c80b92a0efc50282f0a813a0">TE1RYD</a>;</div> <div class="line"><a name="l00068"></a><span class="lineno"> 68</span>  rate = 0.0;</div> <div class="line"><a name="l00069"></a><span class="lineno"> 69</span>  <span class="keywordflow">for</span>( <span class="keywordtype">int</span> j=0; j < 4; j++ )</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="keywordtype">int</span> k = j + 4;</div> <div class="line"><a name="l00072"></a><span class="lineno"> 72</span>  rate += (dfe[ion-1][k]*<a class="code" href="cddefines_8h.html#a0949cd6050ce9bd2506528396ef7424f">sexp</a>(dfe[ion-1][j]/te));</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>  rate = rate/pow(t,1.5);</div> <div class="line"><a name="l00075"></a><span class="lineno"> 75</span>  }</div> <div class="line"><a name="l00076"></a><span class="lineno"> 76</span>  <span class="keywordflow">return</span> rate;</div> <div class="line"><a name="l00077"></a><span class="lineno"> 77</span> }</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="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="physconst_8h_html_ab1891e433b4c259d18bc1c800ea15316"><div class="ttname"><a href="physconst_8h.html#ab1891e433b4c259d18bc1c800ea15316">EVRYD</a></div><div class="ttdeci">const double EVRYD</div><div class="ttdef"><b>Definition:</b> <a href="physconst_8h_source.html#l00189">physconst.h:189</a></div></div> <div class="ttc" id="atmdat_8h_html"><div class="ttname"><a href="atmdat_8h.html">atmdat.h</a></div></div> <div class="ttc" id="atmdat_8h_html_a3283acfcfc73aca51218308909b4091c"><div class="ttname"><a href="atmdat_8h.html#a3283acfcfc73aca51218308909b4091c">atmdat_dielrec_fe</a></div><div class="ttdeci">double atmdat_dielrec_fe(long int ion, double t)</div><div class="ttdef"><b>Definition:</b> <a href="atmdat__dielrec__fe_8cpp_source.html#l00009">atmdat_dielrec_fe.cpp:9</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="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="physconst_8h_html"><div class="ttname"><a href="physconst_8h.html">physconst.h</a></div></div> <div class="ttc" id="physconst_8h_html_a0db3a0b9c80b92a0efc50282f0a813a0"><div class="ttname"><a href="physconst_8h.html#a0db3a0b9c80b92a0efc50282f0a813a0">TE1RYD</a></div><div class="ttdeci">const double TE1RYD</div><div class="ttdef"><b>Definition:</b> <a href="physconst_8h_source.html#l00183">physconst.h:183</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>