<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "https://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.15"/>
<meta name="viewport" content="width=device-width, initial-scale=1"/>
<title>dfnTrans: Docs_PT/ParticleTracking_source/InitialPartPositions.c 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="navtree.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="resize.js"></script>
<script type="text/javascript" src="navtreedata.js"></script>
<script type="text/javascript" src="navtree.js"></script>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&amp;dn=gpl-2.0.txt GPL-v2 */
  $(document).ready(initResizable);
/* @license-end */</script>
<link href="search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="search/search.js"></script>
<link rel="search" href="search_opensearch.php?v=opensearch.xml" type="application/opensearchdescription+xml" title="dfnTrans"/>
<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 id="projectalign" style="padding-left: 0.5em;">
   <div id="projectname">dfnTrans
   </div>
   <div id="projectbrief">Code for Particle Tracking simulations in 3D DFN</div>
  </td>
 </tr>
 </tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.15 -->
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&amp;dn=gpl-2.0.txt GPL-v2 */
var searchBox = new SearchBox("searchBox", "search",false,'Search');
/* @license-end */
</script>
<script type="text/javascript" src="menudata.js"></script>
<script type="text/javascript" src="menu.js"></script>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&amp;dn=gpl-2.0.txt GPL-v2 */
$(function() {
  initMenu('',true,true,'search.php','Search');
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&amp;dn=gpl-2.0.txt GPL-v2 */
  $(document).ready(function() {
    if ($('.searchresults').length > 0) { searchBox.DOMSearchField().focus(); }
  });
});
/* @license-end */</script>
<div id="main-nav"></div>
</div><!-- top -->
<div id="side-nav" class="ui-resizable side-nav-resizable">
  <div id="nav-tree">
    <div id="nav-tree-contents">
      <div id="nav-sync" class="sync"></div>
    </div>
  </div>
  <div id="splitbar" style="-moz-user-select:none;" 
       class="ui-resizable-handle">
  </div>
</div>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&amp;dn=gpl-2.0.txt GPL-v2 */
$(document).ready(function(){initNavTree('_initial_part_positions_8c_source.html','');});
/* @license-end */
</script>
<div id="doc-content">
<div class="header">
  <div class="headertitle">
<div class="title">InitialPartPositions.c</div>  </div>
</div><!--header-->
<div class="contents">
<a href="_initial_part_positions_8c.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno">    1</span>&#160;<span class="preprocessor">#include &lt;stdio.h&gt;</span></div><div class="line"><a name="l00002"></a><span class="lineno">    2</span>&#160;<span class="preprocessor">#include &lt;search.h&gt;</span></div><div class="line"><a name="l00003"></a><span class="lineno">    3</span>&#160;<span class="preprocessor">#include &lt;stdlib.h&gt;</span></div><div class="line"><a name="l00004"></a><span class="lineno">    4</span>&#160;<span class="preprocessor">#include &lt;math.h&gt;</span></div><div class="line"><a name="l00005"></a><span class="lineno">    5</span>&#160;<span class="preprocessor">#include &lt;string.h&gt;</span></div><div class="line"><a name="l00006"></a><span class="lineno">    6</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="_func_def_8h.html">FuncDef.h</a>&quot;</span></div><div class="line"><a name="l00007"></a><span class="lineno">    7</span>&#160;<span class="preprocessor">#include &lt;unistd.h&gt;</span></div><div class="line"><a name="l00008"></a><span class="lineno">    8</span>&#160;<span class="preprocessor">#include &lt;time.h&gt;</span></div><div class="line"><a name="l00009"></a><span class="lineno">    9</span>&#160;</div><div class="line"><a name="l00010"></a><span class="lineno"><a class="line" href="structinpfile.html">   10</a></span>&#160;<span class="keyword">struct </span><a class="code" href="structinpfile.html">inpfile</a> {</div><div class="line"><a name="l00011"></a><span class="lineno"><a class="line" href="structinpfile.html#a66e8759aa38ee20350f04cdc58dee702">   11</a></span>&#160;    <span class="keywordtype">char</span> <a class="code" href="structinpfile.html#a66e8759aa38ee20350f04cdc58dee702">filename</a>[120];</div><div class="line"><a name="l00012"></a><span class="lineno"><a class="line" href="structinpfile.html#a23ae60bf94767194de9a09a4eee8cc1d">   12</a></span>&#160;    <span class="keywordtype">long</span> <span class="keywordtype">int</span> <a class="code" href="structinpfile.html#a23ae60bf94767194de9a09a4eee8cc1d">flag</a>;</div><div class="line"><a name="l00013"></a><span class="lineno"><a class="line" href="structinpfile.html#a57ff2cce78c35a8a128927da9eeccfbc">   13</a></span>&#160;    <span class="keywordtype">double</span> <a class="code" href="structinpfile.html#a57ff2cce78c35a8a128927da9eeccfbc">param</a>;</div><div class="line"><a name="l00014"></a><span class="lineno">   14</span>&#160;};</div><div class="line"><a name="l00015"></a><span class="lineno">   15</span>&#160;</div><div class="line"><a name="l00016"></a><span class="lineno"><a class="line" href="structposit3d.html">   16</a></span>&#160;<span class="keyword">struct </span><a class="code" href="structposit3d.html">posit3d</a> {</div><div class="line"><a name="l00017"></a><span class="lineno"><a class="line" href="structposit3d.html#aefc0a716be214a76cae4e2176134c28e">   17</a></span>&#160;    <span class="keywordtype">double</span> <a class="code" href="structposit3d.html#aefc0a716be214a76cae4e2176134c28e">cord3</a>[3];</div><div class="line"><a name="l00018"></a><span class="lineno">   18</span>&#160;};</div><div class="line"><a name="l00019"></a><span class="lineno">   19</span>&#160;</div><div class="line"><a name="l00020"></a><span class="lineno">   20</span>&#160;</div><div class="line"><a name="l00021"></a><span class="lineno"><a class="line" href="_initial_part_positions_8c.html#ae2df2bd90c8d4aa96e508002a36327c9">   21</a></span>&#160;<span class="keywordtype">int</span> <a class="code" href="_initial_part_positions_8c.html#ae2df2bd90c8d4aa96e508002a36327c9">InitPos</a>()</div><div class="line"><a name="l00023"></a><span class="lineno">   23</span>&#160;{</div><div class="line"><a name="l00024"></a><span class="lineno">   24</span>&#160;    </div><div class="line"><a name="l00025"></a><span class="lineno">   25</span>&#160;    <span class="keywordtype">int</span> i, k_current, k_new=0, numbf=1,  frc,firstn,lastn, flag_in=0,parts_fracture;</div><div class="line"><a name="l00026"></a><span class="lineno">   26</span>&#160;    <span class="keyword">struct </span><a class="code" href="structinpfile.html">inpfile</a> initfile;</div><div class="line"><a name="l00027"></a><span class="lineno">   27</span>&#160;    <span class="keywordtype">double</span> parts_dist=0;</div><div class="line"><a name="l00028"></a><span class="lineno">   28</span>&#160;    <span class="keywordtype">int</span>  zonenumb_in=0, first_ind=0, last_ind=0;</div><div class="line"><a name="l00029"></a><span class="lineno">   29</span>&#160;    <span class="keywordtype">double</span> ixmin=0, ixmax=0, iymin=0, iymax=0, izmin=0, izmax=0;</div><div class="line"><a name="l00030"></a><span class="lineno">   30</span>&#160;    <span class="keywordtype">double</span> px[4]={0.0, 0.0, 0.0, 0.0}, py[4]={0.0, 0.0, 0.0, 0.0};</div><div class="line"><a name="l00031"></a><span class="lineno">   31</span>&#160;    <span class="keywordtype">double</span> weight_p=0.0;</div><div class="line"><a name="l00032"></a><span class="lineno">   32</span>&#160;    </div><div class="line"><a name="l00033"></a><span class="lineno">   33</span>&#160;    <span class="comment">/* calculating number of fractures in in-flow boundary face of the domain ****/</span></div><div class="line"><a name="l00034"></a><span class="lineno">   34</span>&#160;    </div><div class="line"><a name="l00035"></a><span class="lineno">   35</span>&#160;    <span class="keywordflow">for</span> (i=0; i&lt;<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>; i++)</div><div class="line"><a name="l00036"></a><span class="lineno">   36</span>&#160;    {</div><div class="line"><a name="l00037"></a><span class="lineno">   37</span>&#160;        </div><div class="line"><a name="l00038"></a><span class="lineno">   38</span>&#160;        <span class="keywordflow">if</span> ((i&gt;0)&amp;&amp;(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[0]!=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i-1]-1].<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[0]))</div><div class="line"><a name="l00039"></a><span class="lineno">   39</span>&#160;        {</div><div class="line"><a name="l00040"></a><span class="lineno">   40</span>&#160;            <span class="keywordflow">if</span>  (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[1]!=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i-1]-1].<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[0])</div><div class="line"><a name="l00041"></a><span class="lineno">   41</span>&#160;                </div><div class="line"><a name="l00042"></a><span class="lineno">   42</span>&#160;                <span class="comment">//      printf(&quot;fract %d &quot;, node[nodezonein[i]-1].fracture[0]);</span></div><div class="line"><a name="l00043"></a><span class="lineno">   43</span>&#160;                numbf++;</div><div class="line"><a name="l00044"></a><span class="lineno">   44</span>&#160;            </div><div class="line"><a name="l00045"></a><span class="lineno">   45</span>&#160;            </div><div class="line"><a name="l00046"></a><span class="lineno">   46</span>&#160;        }</div><div class="line"><a name="l00047"></a><span class="lineno">   47</span>&#160;    }</div><div class="line"><a name="l00048"></a><span class="lineno">   48</span>&#160;    </div><div class="line"><a name="l00049"></a><span class="lineno">   49</span>&#160;    <span class="keywordflow">if</span> ((numbf==0)&amp;&amp;(<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>!=0))</div><div class="line"><a name="l00050"></a><span class="lineno">   50</span>&#160;        numbf=1;</div><div class="line"><a name="l00051"></a><span class="lineno">   51</span>&#160;    printf(<span class="stringliteral">&quot; %d fracture(s) in in-flow boundary zone \n&quot;</span>,numbf);</div><div class="line"><a name="l00052"></a><span class="lineno">   52</span>&#160;    </div><div class="line"><a name="l00053"></a><span class="lineno">   53</span>&#160;    <span class="comment">// define first and last node and their ID for each boundary fracture edge</span></div><div class="line"><a name="l00054"></a><span class="lineno">   54</span>&#160;    <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> firstind[numbf], lastind[numbf], firstnode[numbf], lastnode[numbf], numberf=1, curfr=0;</div><div class="line"><a name="l00055"></a><span class="lineno">   55</span>&#160;    firstnode[0]=<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[0];</div><div class="line"><a name="l00056"></a><span class="lineno">   56</span>&#160;    lastnode[0]=<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[0];</div><div class="line"><a name="l00057"></a><span class="lineno">   57</span>&#160;    firstind[0]=0;</div><div class="line"><a name="l00058"></a><span class="lineno">   58</span>&#160;    lastind[0]=0;</div><div class="line"><a name="l00059"></a><span class="lineno">   59</span>&#160;    curfr=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[0]-1].<a class="code" href="structvertex.html#aa2cd06573422bc415cb8f38a67232de7">fracture</a>[0];</div><div class="line"><a name="l00060"></a><span class="lineno">   60</span>&#160;    </div><div class="line"><a name="l00061"></a><span class="lineno">   61</span>&#160;    <span class="keywordflow">for</span> (i=1; i&lt;<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>; i++)</div><div class="line"><a name="l00062"></a><span class="lineno">   62</span>&#160;    {</div><div class="line"><a name="l00063"></a><span class="lineno">   63</span>&#160;        </div><div class="line"><a name="l00064"></a><span class="lineno">   64</span>&#160;        <span class="keywordflow">if</span> ((<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[0]!=curfr) &amp;&amp; (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[1]!=curfr))</div><div class="line"><a name="l00065"></a><span class="lineno">   65</span>&#160;        {</div><div class="line"><a name="l00066"></a><span class="lineno">   66</span>&#160;            lastnode[numberf-1]=<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i-1];</div><div class="line"><a name="l00067"></a><span class="lineno">   67</span>&#160;            lastind[numberf-1]=i-1;</div><div class="line"><a name="l00068"></a><span class="lineno">   68</span>&#160;            <span class="keywordflow">if</span> (i&lt;<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>-1)</div><div class="line"><a name="l00069"></a><span class="lineno">   69</span>&#160;            {</div><div class="line"><a name="l00070"></a><span class="lineno">   70</span>&#160;                numberf++;</div><div class="line"><a name="l00071"></a><span class="lineno">   71</span>&#160;                firstnode[numberf-1]=<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i];</div><div class="line"><a name="l00072"></a><span class="lineno">   72</span>&#160;                firstind[numberf-1]=i;</div><div class="line"><a name="l00073"></a><span class="lineno">   73</span>&#160;                curfr=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].<a class="code" href="structvertex.html#aa2cd06573422bc415cb8f38a67232de7">fracture</a>[0];</div><div class="line"><a name="l00074"></a><span class="lineno">   74</span>&#160;            }</div><div class="line"><a name="l00075"></a><span class="lineno">   75</span>&#160;            </div><div class="line"><a name="l00076"></a><span class="lineno">   76</span>&#160;        }</div><div class="line"><a name="l00077"></a><span class="lineno">   77</span>&#160;    }</div><div class="line"><a name="l00078"></a><span class="lineno">   78</span>&#160;    </div><div class="line"><a name="l00079"></a><span class="lineno">   79</span>&#160;    lastnode[numberf-1]=<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>-1];</div><div class="line"><a name="l00080"></a><span class="lineno">   80</span>&#160;    lastind[numberf-1]=<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>-1;</div><div class="line"><a name="l00081"></a><span class="lineno">   81</span>&#160;    </div><div class="line"><a name="l00082"></a><span class="lineno">   82</span>&#160;    <span class="keywordflow">if</span> (numberf&gt;numbf)</div><div class="line"><a name="l00083"></a><span class="lineno">   83</span>&#160;        printf(<span class="stringliteral">&quot;check the boundary zone file\n&quot;</span>);</div><div class="line"><a name="l00084"></a><span class="lineno">   84</span>&#160;    </div><div class="line"><a name="l00085"></a><span class="lineno">   85</span>&#160;    </div><div class="line"><a name="l00086"></a><span class="lineno">   86</span>&#160;    <span class="comment">//  check the input options of particles positions</span></div><div class="line"><a name="l00087"></a><span class="lineno">   87</span>&#160;    <span class="keywordtype">double</span> inter_p[numbf][4];</div><div class="line"><a name="l00088"></a><span class="lineno">   88</span>&#160;    <span class="keywordtype">int</span> inter_fr[numbf];</div><div class="line"><a name="l00089"></a><span class="lineno">   89</span>&#160;    <span class="keywordtype">int</span> res;</div><div class="line"><a name="l00090"></a><span class="lineno">   90</span>&#160;    initfile=<a class="code" href="_func_def_8h.html#a2afbea1a5d0b14d8a098f551fe055420">Control_File</a>(<span class="stringliteral">&quot;init_nf:&quot;</span>,8);</div><div class="line"><a name="l00091"></a><span class="lineno">   91</span>&#160;    res=strncmp(initfile.<a class="code" href="structinpfile.html#a66e8759aa38ee20350f04cdc58dee702">filename</a>,<span class="stringliteral">&quot;yes&quot;</span>,3);</div><div class="line"><a name="l00093"></a><span class="lineno">   93</span>&#160;    <span class="keywordflow">if</span> (res==0)</div><div class="line"><a name="l00094"></a><span class="lineno">   94</span>&#160;    {</div><div class="line"><a name="l00095"></a><span class="lineno">   95</span>&#160;        </div><div class="line"><a name="l00096"></a><span class="lineno">   96</span>&#160;        flag_in=1;</div><div class="line"><a name="l00097"></a><span class="lineno">   97</span>&#160;        <a class="code" href="_func_def_8h.html#adef0bd5e77533cdd4cdd0bcf0d5dc96b">flag_w</a>=1;</div><div class="line"><a name="l00098"></a><span class="lineno">   98</span>&#160;        initfile = <a class="code" href="_func_def_8h.html#a9bb4185d1b46d02dc42d5889e4d70ebf">Control_Data</a>(<span class="stringliteral">&quot;init_partn:&quot;</span>,11 );</div><div class="line"><a name="l00099"></a><span class="lineno">   99</span>&#160;        parts_fracture=initfile.<a class="code" href="structinpfile.html#a23ae60bf94767194de9a09a4eee8cc1d">flag</a>;</div><div class="line"><a name="l00100"></a><span class="lineno">  100</span>&#160;        printf(<span class="stringliteral">&quot;\n  %d  particles per boundary fracture edge \n&quot;</span>, parts_fracture);</div><div class="line"><a name="l00101"></a><span class="lineno">  101</span>&#160;        <a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>=(numbf+1)*parts_fracture;</div><div class="line"><a name="l00102"></a><span class="lineno">  102</span>&#160;        </div><div class="line"><a name="l00103"></a><span class="lineno">  103</span>&#160;        </div><div class="line"><a name="l00104"></a><span class="lineno">  104</span>&#160;        <span class="comment">/*  memory allocatation for particle  */</span></div><div class="line"><a name="l00105"></a><span class="lineno">  105</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>=(<span class="keyword">struct </span><a class="code" href="structcontam.html">contam</a>*) malloc (<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>*<span class="keyword">sizeof</span>(<span class="keyword">struct</span> <a class="code" href="structcontam.html">contam</a>));</div><div class="line"><a name="l00106"></a><span class="lineno">  106</span>&#160;    }</div><div class="line"><a name="l00107"></a><span class="lineno">  107</span>&#160;    <span class="keywordflow">else</span></div><div class="line"><a name="l00108"></a><span class="lineno">  108</span>&#160;    {</div><div class="line"><a name="l00109"></a><span class="lineno">  109</span>&#160;        initfile=<a class="code" href="_func_def_8h.html#a2afbea1a5d0b14d8a098f551fe055420">Control_File</a>(<span class="stringliteral">&quot;init_eqd:&quot;</span>,9);</div><div class="line"><a name="l00110"></a><span class="lineno">  110</span>&#160;        res=strncmp(initfile.<a class="code" href="structinpfile.html#a66e8759aa38ee20350f04cdc58dee702">filename</a>,<span class="stringliteral">&quot;yes&quot;</span>,3);</div><div class="line"><a name="l00111"></a><span class="lineno">  111</span>&#160;        <span class="keywordflow">if</span> (res==0)</div><div class="line"><a name="l00112"></a><span class="lineno">  112</span>&#160;        {</div><div class="line"><a name="l00113"></a><span class="lineno">  113</span>&#160;            flag_in=2;</div><div class="line"><a name="l00114"></a><span class="lineno">  114</span>&#160;            <a class="code" href="_func_def_8h.html#adef0bd5e77533cdd4cdd0bcf0d5dc96b">flag_w</a>=1;</div><div class="line"><a name="l00119"></a><span class="lineno">  119</span>&#160;            initfile = <a class="code" href="_func_def_8h.html#a9bb4185d1b46d02dc42d5889e4d70ebf">Control_Data</a>(<span class="stringliteral">&quot;init_npart:&quot;</span>,11 );</div><div class="line"><a name="l00120"></a><span class="lineno">  120</span>&#160;            parts_fracture=initfile.<a class="code" href="structinpfile.html#a23ae60bf94767194de9a09a4eee8cc1d">flag</a>;</div><div class="line"><a name="l00121"></a><span class="lineno">  121</span>&#160;            <a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>=(numbf)*parts_fracture*2;</div><div class="line"><a name="l00122"></a><span class="lineno">  122</span>&#160;            <span class="comment">/*  memory allocatation for particle structures */</span></div><div class="line"><a name="l00123"></a><span class="lineno">  123</span>&#160;            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>=(<span class="keyword">struct </span><a class="code" href="structcontam.html">contam</a>*) malloc ((<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>)*<span class="keyword">sizeof</span>(<span class="keyword">struct </span><a class="code" href="structcontam.html">contam</a>));</div><div class="line"><a name="l00124"></a><span class="lineno">  124</span>&#160;            <span class="comment">// define a distance between particles</span></div><div class="line"><a name="l00125"></a><span class="lineno">  125</span>&#160;            </div><div class="line"><a name="l00126"></a><span class="lineno">  126</span>&#160;            <span class="keywordtype">double</span> length2=0.0, t_length=0.0;</div><div class="line"><a name="l00127"></a><span class="lineno">  127</span>&#160;            <span class="keywordflow">for</span> (i=0; i&lt;numbf; i++)</div><div class="line"><a name="l00128"></a><span class="lineno">  128</span>&#160;            {</div><div class="line"><a name="l00129"></a><span class="lineno">  129</span>&#160;                </div><div class="line"><a name="l00130"></a><span class="lineno">  130</span>&#160;                length2=pow((<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstnode[i]-1].coord[0]-<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastnode[i]-1].coord[0]),2)+pow((<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstnode[i]-1].coord[1]-<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastnode[i]-1].coord[1]),2)+pow((<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstnode[i]-1].coord[2]-<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastnode[i]-1].coord[2]),2);</div><div class="line"><a name="l00131"></a><span class="lineno">  131</span>&#160;                t_length=t_length+sqrt(length2);</div><div class="line"><a name="l00132"></a><span class="lineno">  132</span>&#160;                </div><div class="line"><a name="l00133"></a><span class="lineno">  133</span>&#160;            }</div><div class="line"><a name="l00134"></a><span class="lineno">  134</span>&#160;            </div><div class="line"><a name="l00135"></a><span class="lineno">  135</span>&#160;            </div><div class="line"><a name="l00136"></a><span class="lineno">  136</span>&#160;            parts_dist=t_length/(parts_fracture*numbf);</div><div class="line"><a name="l00137"></a><span class="lineno">  137</span>&#160;            printf(<span class="stringliteral">&quot;\n  Particles placed on %f [m]  from each other. Total length %lf  \n&quot;</span>, parts_dist, t_length);</div><div class="line"><a name="l00138"></a><span class="lineno">  138</span>&#160;        }</div><div class="line"><a name="l00139"></a><span class="lineno">  139</span>&#160;        <span class="keywordflow">else</span></div><div class="line"><a name="l00140"></a><span class="lineno">  140</span>&#160;        {</div><div class="line"><a name="l00141"></a><span class="lineno">  141</span>&#160;            initfile=<a class="code" href="_func_def_8h.html#a2afbea1a5d0b14d8a098f551fe055420">Control_File</a>(<span class="stringliteral">&quot;init_oneregion:&quot;</span>,15);</div><div class="line"><a name="l00142"></a><span class="lineno">  142</span>&#160;            res=strncmp(initfile.<a class="code" href="structinpfile.html#a66e8759aa38ee20350f04cdc58dee702">filename</a>,<span class="stringliteral">&quot;yes&quot;</span>,3);</div><div class="line"><a name="l00143"></a><span class="lineno">  143</span>&#160;            <span class="keywordflow">if</span> (res==0)</div><div class="line"><a name="l00144"></a><span class="lineno">  144</span>&#160;            {</div><div class="line"><a name="l00145"></a><span class="lineno">  145</span>&#160;                flag_in=3;</div><div class="line"><a name="l00146"></a><span class="lineno">  146</span>&#160;                <a class="code" href="_func_def_8h.html#adef0bd5e77533cdd4cdd0bcf0d5dc96b">flag_w</a>=1;</div><div class="line"><a name="l00150"></a><span class="lineno">  150</span>&#160;                initfile = <a class="code" href="_func_def_8h.html#a9bb4185d1b46d02dc42d5889e4d70ebf">Control_Data</a>(<span class="stringliteral">&quot;in_partn:&quot;</span>,9 );</div><div class="line"><a name="l00151"></a><span class="lineno">  151</span>&#160;                <a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>=initfile.<a class="code" href="structinpfile.html#a23ae60bf94767194de9a09a4eee8cc1d">flag</a>;</div><div class="line"><a name="l00152"></a><span class="lineno">  152</span>&#160;                <span class="comment">/*  memory allocatation for particle structures */</span></div><div class="line"><a name="l00153"></a><span class="lineno">  153</span>&#160;                <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>=(<span class="keyword">struct </span><a class="code" href="structcontam.html">contam</a>*) malloc ((<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>*2)*<span class="keyword">sizeof</span>(<span class="keyword">struct </span><a class="code" href="structcontam.html">contam</a>));</div><div class="line"><a name="l00154"></a><span class="lineno">  154</span>&#160;                printf(<span class="stringliteral">&quot;\n Initially particles have the same starting region \n&quot;</span>);</div><div class="line"><a name="l00155"></a><span class="lineno">  155</span>&#160;                </div><div class="line"><a name="l00156"></a><span class="lineno">  156</span>&#160;                initfile = <a class="code" href="_func_def_8h.html#a922a07b425353141b29d8d57e830a7c1">Control_Param</a>(<span class="stringliteral">&quot;in_xmin:&quot;</span>,8 );</div><div class="line"><a name="l00157"></a><span class="lineno">  157</span>&#160;                ixmin=initfile.<a class="code" href="structinpfile.html#a57ff2cce78c35a8a128927da9eeccfbc">param</a>;</div><div class="line"><a name="l00158"></a><span class="lineno">  158</span>&#160;                initfile = <a class="code" href="_func_def_8h.html#a922a07b425353141b29d8d57e830a7c1">Control_Param</a>(<span class="stringliteral">&quot;in_xmax:&quot;</span>,8 );</div><div class="line"><a name="l00159"></a><span class="lineno">  159</span>&#160;                ixmax=initfile.<a class="code" href="structinpfile.html#a57ff2cce78c35a8a128927da9eeccfbc">param</a>;</div><div class="line"><a name="l00160"></a><span class="lineno">  160</span>&#160;                initfile = <a class="code" href="_func_def_8h.html#a922a07b425353141b29d8d57e830a7c1">Control_Param</a>(<span class="stringliteral">&quot;in_ymin:&quot;</span>,8 );</div><div class="line"><a name="l00161"></a><span class="lineno">  161</span>&#160;                iymin=initfile.<a class="code" href="structinpfile.html#a57ff2cce78c35a8a128927da9eeccfbc">param</a>;</div><div class="line"><a name="l00162"></a><span class="lineno">  162</span>&#160;                initfile = <a class="code" href="_func_def_8h.html#a922a07b425353141b29d8d57e830a7c1">Control_Param</a>(<span class="stringliteral">&quot;in_ymax:&quot;</span>,8 );</div><div class="line"><a name="l00163"></a><span class="lineno">  163</span>&#160;                iymax=initfile.<a class="code" href="structinpfile.html#a57ff2cce78c35a8a128927da9eeccfbc">param</a>;</div><div class="line"><a name="l00164"></a><span class="lineno">  164</span>&#160;                initfile = <a class="code" href="_func_def_8h.html#a922a07b425353141b29d8d57e830a7c1">Control_Param</a>(<span class="stringliteral">&quot;in_zmin:&quot;</span>,8 );</div><div class="line"><a name="l00165"></a><span class="lineno">  165</span>&#160;                izmin=initfile.<a class="code" href="structinpfile.html#a57ff2cce78c35a8a128927da9eeccfbc">param</a>;</div><div class="line"><a name="l00166"></a><span class="lineno">  166</span>&#160;                initfile = <a class="code" href="_func_def_8h.html#a922a07b425353141b29d8d57e830a7c1">Control_Param</a>(<span class="stringliteral">&quot;in_zmax:&quot;</span>,8 );</div><div class="line"><a name="l00167"></a><span class="lineno">  167</span>&#160;                izmax=initfile.<a class="code" href="structinpfile.html#a57ff2cce78c35a8a128927da9eeccfbc">param</a>;</div><div class="line"><a name="l00168"></a><span class="lineno">  168</span>&#160;                </div><div class="line"><a name="l00169"></a><span class="lineno">  169</span>&#160;                initfile = <a class="code" href="_func_def_8h.html#a9bb4185d1b46d02dc42d5889e4d70ebf">Control_Data</a>(<span class="stringliteral">&quot;in-flow-boundary:&quot;</span>,17 );</div><div class="line"><a name="l00170"></a><span class="lineno">  170</span>&#160;                zonenumb_in=initfile.<a class="code" href="structinpfile.html#a23ae60bf94767194de9a09a4eee8cc1d">flag</a>;</div><div class="line"><a name="l00171"></a><span class="lineno">  171</span>&#160;                <span class="comment">/* define coordinations of region according to in-flow zone */</span></div><div class="line"><a name="l00172"></a><span class="lineno">  172</span>&#160;                <span class="keywordflow">if</span> ((zonenumb_in==1)||(zonenumb_in==2))</div><div class="line"><a name="l00173"></a><span class="lineno">  173</span>&#160;                {</div><div class="line"><a name="l00174"></a><span class="lineno">  174</span>&#160;                    px[0]=ixmin;</div><div class="line"><a name="l00175"></a><span class="lineno">  175</span>&#160;                    py[0]=iymin;</div><div class="line"><a name="l00176"></a><span class="lineno">  176</span>&#160;                    px[1]=ixmin;</div><div class="line"><a name="l00177"></a><span class="lineno">  177</span>&#160;                    py[1]=iymax;</div><div class="line"><a name="l00178"></a><span class="lineno">  178</span>&#160;                    px[2]=ixmax;</div><div class="line"><a name="l00179"></a><span class="lineno">  179</span>&#160;                    py[2]=iymax;</div><div class="line"><a name="l00180"></a><span class="lineno">  180</span>&#160;                    px[3]=ixmax;</div><div class="line"><a name="l00181"></a><span class="lineno">  181</span>&#160;                    py[3]=iymin;</div><div class="line"><a name="l00182"></a><span class="lineno">  182</span>&#160;                }</div><div class="line"><a name="l00183"></a><span class="lineno">  183</span>&#160;                </div><div class="line"><a name="l00184"></a><span class="lineno">  184</span>&#160;                <span class="keywordflow">if</span> ((zonenumb_in==3)||(zonenumb_in==5))</div><div class="line"><a name="l00185"></a><span class="lineno">  185</span>&#160;                {</div><div class="line"><a name="l00186"></a><span class="lineno">  186</span>&#160;                    px[0]=izmin;</div><div class="line"><a name="l00187"></a><span class="lineno">  187</span>&#160;                    py[0]=iymin;</div><div class="line"><a name="l00188"></a><span class="lineno">  188</span>&#160;                    px[1]=izmin;</div><div class="line"><a name="l00189"></a><span class="lineno">  189</span>&#160;                    py[1]=iymax;</div><div class="line"><a name="l00190"></a><span class="lineno">  190</span>&#160;                    px[2]=izmax;</div><div class="line"><a name="l00191"></a><span class="lineno">  191</span>&#160;                    py[2]=iymax;</div><div class="line"><a name="l00192"></a><span class="lineno">  192</span>&#160;                    px[3]=izmax;</div><div class="line"><a name="l00193"></a><span class="lineno">  193</span>&#160;                    py[3]=iymin;</div><div class="line"><a name="l00194"></a><span class="lineno">  194</span>&#160;                }</div><div class="line"><a name="l00195"></a><span class="lineno">  195</span>&#160;                </div><div class="line"><a name="l00196"></a><span class="lineno">  196</span>&#160;                <span class="keywordflow">if</span> ((zonenumb_in==4)||(zonenumb_in==6))</div><div class="line"><a name="l00197"></a><span class="lineno">  197</span>&#160;                {</div><div class="line"><a name="l00198"></a><span class="lineno">  198</span>&#160;                    px[0]=ixmin;</div><div class="line"><a name="l00199"></a><span class="lineno">  199</span>&#160;                    py[0]=izmin;</div><div class="line"><a name="l00200"></a><span class="lineno">  200</span>&#160;                    px[1]=ixmin;</div><div class="line"><a name="l00201"></a><span class="lineno">  201</span>&#160;                    py[1]=izmax;</div><div class="line"><a name="l00202"></a><span class="lineno">  202</span>&#160;                    px[2]=ixmax;</div><div class="line"><a name="l00203"></a><span class="lineno">  203</span>&#160;                    py[2]=izmax;</div><div class="line"><a name="l00204"></a><span class="lineno">  204</span>&#160;                    px[3]=ixmax;</div><div class="line"><a name="l00205"></a><span class="lineno">  205</span>&#160;                    py[3]=izmin;</div><div class="line"><a name="l00206"></a><span class="lineno">  206</span>&#160;                }</div><div class="line"><a name="l00207"></a><span class="lineno">  207</span>&#160;            }</div><div class="line"><a name="l00208"></a><span class="lineno">  208</span>&#160;            <span class="keywordflow">else</span></div><div class="line"><a name="l00209"></a><span class="lineno">  209</span>&#160;            {</div><div class="line"><a name="l00211"></a><span class="lineno">  211</span>&#160;                initfile=<a class="code" href="_func_def_8h.html#a2afbea1a5d0b14d8a098f551fe055420">Control_File</a>(<span class="stringliteral">&quot;init_random:&quot;</span>,12);</div><div class="line"><a name="l00212"></a><span class="lineno">  212</span>&#160;                res=strncmp(initfile.<a class="code" href="structinpfile.html#a66e8759aa38ee20350f04cdc58dee702">filename</a>,<span class="stringliteral">&quot;yes&quot;</span>,3);</div><div class="line"><a name="l00213"></a><span class="lineno">  213</span>&#160;                <span class="keywordflow">if</span> (res==0)</div><div class="line"><a name="l00214"></a><span class="lineno">  214</span>&#160;                {</div><div class="line"><a name="l00215"></a><span class="lineno">  215</span>&#160;                    flag_in=4;</div><div class="line"><a name="l00216"></a><span class="lineno">  216</span>&#160;                    initfile = <a class="code" href="_func_def_8h.html#a9bb4185d1b46d02dc42d5889e4d70ebf">Control_Data</a>(<span class="stringliteral">&quot;in_randpart:&quot;</span>,12 );</div><div class="line"><a name="l00217"></a><span class="lineno">  217</span>&#160;                    <a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>=initfile.<a class="code" href="structinpfile.html#a23ae60bf94767194de9a09a4eee8cc1d">flag</a>;</div><div class="line"><a name="l00218"></a><span class="lineno">  218</span>&#160;                    <span class="comment">/*  memory allocatation for particle structures */</span></div><div class="line"><a name="l00219"></a><span class="lineno">  219</span>&#160;                    <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>=(<span class="keyword">struct </span><a class="code" href="structcontam.html">contam</a>*) malloc ((<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>+1)*<span class="keyword">sizeof</span>(<span class="keyword">struct </span><a class="code" href="structcontam.html">contam</a>));</div><div class="line"><a name="l00220"></a><span class="lineno">  220</span>&#160;                    printf(<span class="stringliteral">&quot;\n Initially particles will be distributed randomly over all fracture surfaces \n&quot;</span>);</div><div class="line"><a name="l00221"></a><span class="lineno">  221</span>&#160;                    </div><div class="line"><a name="l00222"></a><span class="lineno">  222</span>&#160;                    <span class="keywordtype">double</span> random_number=0, sum_aperture=0.0;</div><div class="line"><a name="l00223"></a><span class="lineno">  223</span>&#160;                    </div><div class="line"><a name="l00224"></a><span class="lineno">  224</span>&#160;                    <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> currentcell, k_curr=0;</div><div class="line"><a name="l00225"></a><span class="lineno">  225</span>&#160;                    </div><div class="line"><a name="l00226"></a><span class="lineno">  226</span>&#160;                    </div><div class="line"><a name="l00227"></a><span class="lineno">  227</span>&#160;                    <span class="keywordflow">do</span></div><div class="line"><a name="l00228"></a><span class="lineno">  228</span>&#160;                    {</div><div class="line"><a name="l00229"></a><span class="lineno">  229</span>&#160;                        </div><div class="line"><a name="l00230"></a><span class="lineno">  230</span>&#160;                        random_number=drand48();</div><div class="line"><a name="l00231"></a><span class="lineno">  231</span>&#160;                        </div><div class="line"><a name="l00232"></a><span class="lineno">  232</span>&#160;                        currentcell=random_number*<a class="code" href="_func_def_8h.html#a5225ceefbfdde0f9cb32cb84174a4863">ncells</a>;</div><div class="line"><a name="l00233"></a><span class="lineno">  233</span>&#160;                        </div><div class="line"><a name="l00234"></a><span class="lineno">  234</span>&#160;                        </div><div class="line"><a name="l00235"></a><span class="lineno">  235</span>&#160;                        <span class="keywordflow">if</span> ((currentcell!=0) &amp;&amp; (((<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#a26a18460ba7bc9d4e9ace81fda1819ee">cell</a>[currentcell-1].node_ind[0]-1].typeN&lt;200)||(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#a26a18460ba7bc9d4e9ace81fda1819ee">cell</a>[currentcell-1].node_ind[0]-1].typeN&gt;250))&amp;&amp; ((<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#a26a18460ba7bc9d4e9ace81fda1819ee">cell</a>[currentcell-1].node_ind[1]-1].typeN&lt;200)||(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#a26a18460ba7bc9d4e9ace81fda1819ee">cell</a>[currentcell-1].node_ind[1]-1].typeN&gt;250)) &amp;&amp; ((<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#a26a18460ba7bc9d4e9ace81fda1819ee">cell</a>[currentcell-1].node_ind[2]-1].typeN&lt;200)||(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#a26a18460ba7bc9d4e9ace81fda1819ee">cell</a>[currentcell-1].node_ind[2]-1].typeN&gt;250))))</div><div class="line"><a name="l00236"></a><span class="lineno">  236</span>&#160;                        {</div><div class="line"><a name="l00237"></a><span class="lineno">  237</span>&#160;                            </div><div class="line"><a name="l00238"></a><span class="lineno">  238</span>&#160;                            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_curr].<a class="code" href="structcontam.html#a70e44829116675e21bda7422f0c2b777">velocity</a>[0]=0.;</div><div class="line"><a name="l00239"></a><span class="lineno">  239</span>&#160;                            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_curr].<a class="code" href="structcontam.html#a70e44829116675e21bda7422f0c2b777">velocity</a>[1]=0.;</div><div class="line"><a name="l00240"></a><span class="lineno">  240</span>&#160;                            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_curr].<a class="code" href="structcontam.html#a0b9253ce268d781958145f9d110df2de">fracture</a>=<a class="code" href="_func_def_8h.html#a26a18460ba7bc9d4e9ace81fda1819ee">cell</a>[currentcell-1].<a class="code" href="structelement.html#af51b08146a1c495aadf45e7f085e1267">fracture</a>;</div><div class="line"><a name="l00241"></a><span class="lineno">  241</span>&#160;                            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_curr].<a class="code" href="structcontam.html#a858c05f3221a70f6e212c1e7e172252c">cell</a>=currentcell;</div><div class="line"><a name="l00242"></a><span class="lineno">  242</span>&#160;                            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_curr].<a class="code" href="structcontam.html#a11c5e87e6e62d258648d5a02472251d5">time</a>=0.0;</div><div class="line"><a name="l00243"></a><span class="lineno">  243</span>&#160;                            <a class="code" href="_func_def_8h.html#a69f9dfd9f5eba61f40439fc5dc92f3b6">Moving2Center</a> (k_curr, currentcell);</div><div class="line"><a name="l00244"></a><span class="lineno">  244</span>&#160;                            <span class="keywordtype">int</span> insc;</div><div class="line"><a name="l00245"></a><span class="lineno">  245</span>&#160;                            insc=<a class="code" href="_func_def_8h.html#a5c8fe8978ce1689f5dcf12bc7b56793d">InsideCell</a> (currentcell);</div><div class="line"><a name="l00246"></a><span class="lineno">  246</span>&#160;                            sum_aperture=sum_aperture+<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#a26a18460ba7bc9d4e9ace81fda1819ee">cell</a>[currentcell-1].<a class="code" href="structelement.html#aadc95dc4817009223151c40387f40e06">node_ind</a>[0]-1].<a class="code" href="structvertex.html#a7588fad77cf759d8db0957c9e957e23a">aperture</a>;</div><div class="line"><a name="l00247"></a><span class="lineno">  247</span>&#160;                            k_curr++;</div><div class="line"><a name="l00248"></a><span class="lineno">  248</span>&#160;                        }</div><div class="line"><a name="l00249"></a><span class="lineno">  249</span>&#160;                        </div><div class="line"><a name="l00250"></a><span class="lineno">  250</span>&#160;                    }</div><div class="line"><a name="l00251"></a><span class="lineno">  251</span>&#160;                    <span class="keywordflow">while</span>(k_curr!=<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>);</div><div class="line"><a name="l00252"></a><span class="lineno">  252</span>&#160;                    </div><div class="line"><a name="l00253"></a><span class="lineno">  253</span>&#160;                    </div><div class="line"><a name="l00254"></a><span class="lineno">  254</span>&#160;                    </div><div class="line"><a name="l00255"></a><span class="lineno">  255</span>&#160;                    k_new=k_curr;</div><div class="line"><a name="l00256"></a><span class="lineno">  256</span>&#160;                    </div><div class="line"><a name="l00257"></a><span class="lineno">  257</span>&#160;                    <span class="keywordflow">for</span> (i=0; i&lt;<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>; i++)</div><div class="line"><a name="l00258"></a><span class="lineno">  258</span>&#160;                    {</div><div class="line"><a name="l00259"></a><span class="lineno">  259</span>&#160;                        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[i].<a class="code" href="structcontam.html#ad50f20b01e37af1fedfd7bbe0cd1d7e3">fl_weight</a>=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#a26a18460ba7bc9d4e9ace81fda1819ee">cell</a>[<a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[i].<a class="code" href="structcontam.html#a858c05f3221a70f6e212c1e7e172252c">cell</a>-1].<a class="code" href="structelement.html#aadc95dc4817009223151c40387f40e06">node_ind</a>[0]-1].<a class="code" href="structvertex.html#a7588fad77cf759d8db0957c9e957e23a">aperture</a>/sum_aperture;</div><div class="line"><a name="l00260"></a><span class="lineno">  260</span>&#160;                    }</div><div class="line"><a name="l00261"></a><span class="lineno">  261</span>&#160;                } <span class="comment">//end if flag_in=4</span></div><div class="line"><a name="l00262"></a><span class="lineno">  262</span>&#160;                <span class="keywordflow">else</span></div><div class="line"><a name="l00263"></a><span class="lineno">  263</span>&#160;                {</div><div class="line"><a name="l00265"></a><span class="lineno">  265</span>&#160;                    initfile=<a class="code" href="_func_def_8h.html#a484100bd5902d4efcdd15c21d1cb9048">Control_File_Optional</a>(<span class="stringliteral">&quot;init_matrix:&quot;</span>,12);</div><div class="line"><a name="l00266"></a><span class="lineno">  266</span>&#160;                    <span class="keywordflow">if</span> (initfile.<a class="code" href="structinpfile.html#a23ae60bf94767194de9a09a4eee8cc1d">flag</a>&gt;0)</div><div class="line"><a name="l00267"></a><span class="lineno">  267</span>&#160;                    {</div><div class="line"><a name="l00268"></a><span class="lineno">  268</span>&#160;                        res=strncmp(initfile.<a class="code" href="structinpfile.html#a66e8759aa38ee20350f04cdc58dee702">filename</a>,<span class="stringliteral">&quot;yes&quot;</span>,3);</div><div class="line"><a name="l00269"></a><span class="lineno">  269</span>&#160;                        <span class="keywordflow">if</span> (res==0)</div><div class="line"><a name="l00270"></a><span class="lineno">  270</span>&#160;                        {</div><div class="line"><a name="l00271"></a><span class="lineno">  271</span>&#160;                            printf(<span class="stringliteral">&quot; Initially particles are placed in rock matrix randomly. &quot;</span>);</div><div class="line"><a name="l00272"></a><span class="lineno">  272</span>&#160;                            printf(<span class="stringliteral">&quot; The closest cells to initial particles positions &quot;</span>);</div><div class="line"><a name="l00273"></a><span class="lineno">  273</span>&#160;                            printf(<span class="stringliteral">&quot; will be set as starting point in DFN. &quot;</span>);</div><div class="line"><a name="l00274"></a><span class="lineno">  274</span>&#160;                            flag_in=5;</div><div class="line"><a name="l00275"></a><span class="lineno">  275</span>&#160;                            <a class="code" href="_initial_part_positions_8c.html#a0355cfa365c20964163da103eec2f0e5">InitInMatrix</a>();</div><div class="line"><a name="l00276"></a><span class="lineno">  276</span>&#160;                            k_new=<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>;</div><div class="line"><a name="l00277"></a><span class="lineno">  277</span>&#160;                        } <span class="comment">// end if flag_in=5</span></div><div class="line"><a name="l00278"></a><span class="lineno">  278</span>&#160;                        <span class="keywordflow">else</span></div><div class="line"><a name="l00279"></a><span class="lineno">  279</span>&#160;                        {</div><div class="line"><a name="l00280"></a><span class="lineno">  280</span>&#160;                            initfile=<a class="code" href="_func_def_8h.html#a2afbea1a5d0b14d8a098f551fe055420">Control_File</a>(<span class="stringliteral">&quot;init_fluxw:&quot;</span>,11);</div><div class="line"><a name="l00281"></a><span class="lineno">  281</span>&#160;                            res=strncmp(initfile.<a class="code" href="structinpfile.html#a66e8759aa38ee20350f04cdc58dee702">filename</a>,<span class="stringliteral">&quot;yes&quot;</span>,3);</div><div class="line"><a name="l00282"></a><span class="lineno">  282</span>&#160;                            <span class="keywordflow">if</span> (res==0)</div><div class="line"><a name="l00283"></a><span class="lineno">  283</span>&#160;                            {</div><div class="line"><a name="l00285"></a><span class="lineno">  285</span>&#160;                                flag_in=6;</div><div class="line"><a name="l00286"></a><span class="lineno">  286</span>&#160;                                <a class="code" href="_func_def_8h.html#adef0bd5e77533cdd4cdd0bcf0d5dc96b">flag_w</a>=0;</div><div class="line"><a name="l00287"></a><span class="lineno">  287</span>&#160;                                </div><div class="line"><a name="l00288"></a><span class="lineno">  288</span>&#160;                                initfile = <a class="code" href="_func_def_8h.html#a9bb4185d1b46d02dc42d5889e4d70ebf">Control_Data</a>(<span class="stringliteral">&quot;init_totalnumber:&quot;</span>,17 );</div><div class="line"><a name="l00289"></a><span class="lineno">  289</span>&#160;                                <a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a> = initfile.<a class="code" href="structinpfile.html#a23ae60bf94767194de9a09a4eee8cc1d">flag</a>;</div><div class="line"><a name="l00290"></a><span class="lineno">  290</span>&#160;                                <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>&lt;<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>)</div><div class="line"><a name="l00291"></a><span class="lineno">  291</span>&#160;                                {</div><div class="line"><a name="l00292"></a><span class="lineno">  292</span>&#160;                                    printf(<span class="stringliteral">&quot;\n  The requested number of particles (%d) is less than number of nodes in in-flow boundary (\%d). The initial number of particles is changed to number of nodes in  in-flow boundary. \n&quot;</span>, <a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>, <a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>);</div><div class="line"><a name="l00293"></a><span class="lineno">  293</span>&#160;                                    <a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>=<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>;</div><div class="line"><a name="l00294"></a><span class="lineno">  294</span>&#160;                                }</div><div class="line"><a name="l00295"></a><span class="lineno">  295</span>&#160;                                weight_p=(<a class="code" href="_func_def_8h.html#aa27b8bb0d1e8acff46eaadd07476ed04">totalFluxIn</a>/(float)<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>);</div><div class="line"><a name="l00296"></a><span class="lineno">  296</span>&#160;                </div><div class="line"><a name="l00297"></a><span class="lineno">  297</span>&#160;                                printf(<span class="stringliteral">&quot;\n  Requested number of particles %d with flux weight %5.12e (Total Flux %5.12e) \n&quot;</span>, <a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>, weight_p, <a class="code" href="_func_def_8h.html#aa27b8bb0d1e8acff46eaadd07476ed04">totalFluxIn</a>);</div><div class="line"><a name="l00298"></a><span class="lineno">  298</span>&#160;                                <span class="keywordtype">int</span> npart_alloc=0;</div><div class="line"><a name="l00299"></a><span class="lineno">  299</span>&#160;                                npart_alloc=<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>*3;</div><div class="line"><a name="l00300"></a><span class="lineno">  300</span>&#160;                                <span class="comment">/*  memory allocatation for particle structures */</span></div><div class="line"><a name="l00301"></a><span class="lineno">  301</span>&#160;                                <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>=(<span class="keyword">struct </span><a class="code" href="structcontam.html">contam</a>*) malloc (npart_alloc*<span class="keyword">sizeof</span>(<span class="keyword">struct</span> <a class="code" href="structcontam.html">contam</a>));</div><div class="line"><a name="l00302"></a><span class="lineno">  302</span>&#160;                            }</div><div class="line"><a name="l00303"></a><span class="lineno">  303</span>&#160;                        }<span class="comment">//end if/else flag_in=6</span></div><div class="line"><a name="l00304"></a><span class="lineno">  304</span>&#160;                        </div><div class="line"><a name="l00305"></a><span class="lineno">  305</span>&#160;                    } <span class="comment">//end if/else flag_in=5</span></div><div class="line"><a name="l00306"></a><span class="lineno">  306</span>&#160;                    </div><div class="line"><a name="l00307"></a><span class="lineno">  307</span>&#160;                } <span class="comment">//end if flag_in=4</span></div><div class="line"><a name="l00308"></a><span class="lineno">  308</span>&#160;                </div><div class="line"><a name="l00309"></a><span class="lineno">  309</span>&#160;            }   <span class="comment">//end if/else flag_in=3</span></div><div class="line"><a name="l00310"></a><span class="lineno">  310</span>&#160;            </div><div class="line"><a name="l00311"></a><span class="lineno">  311</span>&#160;        } <span class="comment">//end if /else flag_in=2</span></div><div class="line"><a name="l00312"></a><span class="lineno">  312</span>&#160;        </div><div class="line"><a name="l00313"></a><span class="lineno">  313</span>&#160;    } <span class="comment">//end if/else flag_in=1</span></div><div class="line"><a name="l00314"></a><span class="lineno">  314</span>&#160;    </div><div class="line"><a name="l00315"></a><span class="lineno">  315</span>&#160;    <span class="comment">/* loop on fractures and initial setting of particles */</span></div><div class="line"><a name="l00316"></a><span class="lineno">  316</span>&#160;    </div><div class="line"><a name="l00317"></a><span class="lineno">  317</span>&#160;    k_current=0;</div><div class="line"><a name="l00318"></a><span class="lineno">  318</span>&#160;    <span class="keywordflow">if</span> ((flag_in&lt;=2) || (flag_in==6))</div><div class="line"><a name="l00319"></a><span class="lineno">  319</span>&#160;    {</div><div class="line"><a name="l00320"></a><span class="lineno">  320</span>&#160;        </div><div class="line"><a name="l00321"></a><span class="lineno">  321</span>&#160;        </div><div class="line"><a name="l00322"></a><span class="lineno">  322</span>&#160;        <span class="keywordflow">for</span> (i=0; i&lt;numberf; i++)</div><div class="line"><a name="l00323"></a><span class="lineno">  323</span>&#160;        {</div><div class="line"><a name="l00324"></a><span class="lineno">  324</span>&#160;            <span class="keywordflow">if</span> (lastnode[i]!=firstnode[i])</div><div class="line"><a name="l00325"></a><span class="lineno">  325</span>&#160;            {</div><div class="line"><a name="l00326"></a><span class="lineno">  326</span>&#160;                <span class="keywordflow">if</span> (flag_in==6) <span class="comment">//initial set according to in flow flux</span></div><div class="line"><a name="l00327"></a><span class="lineno">  327</span>&#160;                {</div><div class="line"><a name="l00328"></a><span class="lineno">  328</span>&#160;                    </div><div class="line"><a name="l00329"></a><span class="lineno">  329</span>&#160;                    k_new=<a class="code" href="_initial_part_positions_8c.html#ad3b7970929b5c9502172afc4e0bc65e8">InitParticles_flux</a>(k_current, firstind[i], lastind[i], weight_p);</div><div class="line"><a name="l00330"></a><span class="lineno">  330</span>&#160;                    <span class="comment">//                printf (&quot;IPF %d %d %d %d %12.5e %d %d %d\n&quot;,i+1, k_current, firstind[i], lastind[i], weight_p, firstnode[i], lastnode[i], node[lastnode[i]-1].fracture[0]);</span></div><div class="line"><a name="l00331"></a><span class="lineno">  331</span>&#160;                    </div><div class="line"><a name="l00332"></a><span class="lineno">  332</span>&#160;                    k_current=k_new;</div><div class="line"><a name="l00333"></a><span class="lineno">  333</span>&#160;                }</div><div class="line"><a name="l00334"></a><span class="lineno">  334</span>&#160;                </div><div class="line"><a name="l00335"></a><span class="lineno">  335</span>&#160;                <span class="keywordflow">if</span> (flag_in==2) <span class="comment">//equidistant placing of particles on the entire in-flow boundary</span></div><div class="line"><a name="l00336"></a><span class="lineno">  336</span>&#160;                {</div><div class="line"><a name="l00337"></a><span class="lineno">  337</span>&#160;                    </div><div class="line"><a name="l00338"></a><span class="lineno">  338</span>&#160;                    k_new=<a class="code" href="_initial_part_positions_8c.html#a01e7bd2a5bd9417f1d28849b8577b0b9">InitParticles_eq</a> (k_current, firstnode[i], lastnode[i], parts_dist, firstind[i], lastind[i]);</div><div class="line"><a name="l00339"></a><span class="lineno">  339</span>&#160;                    <span class="comment">//                 printf(&quot;IPE fract first %d fract last %d number parts %d\n&quot;, node[firstn-1].fracture[0], node[lastn-1].fracture[0], k_new);</span></div><div class="line"><a name="l00340"></a><span class="lineno">  340</span>&#160;                }</div><div class="line"><a name="l00341"></a><span class="lineno">  341</span>&#160;                </div><div class="line"><a name="l00342"></a><span class="lineno">  342</span>&#160;                <span class="keywordflow">if</span> (flag_in==1) <span class="comment">//equidistant placing of particles on each fracture</span></div><div class="line"><a name="l00343"></a><span class="lineno">  343</span>&#160;                {</div><div class="line"><a name="l00344"></a><span class="lineno">  344</span>&#160;                    k_new=<a class="code" href="_initial_part_positions_8c.html#a30b78ca10a78d9e5d6de6d5831a5557a">InitParticles_np</a> (k_current, firstnode[i], lastnode[i], parts_fracture, firstind[i], lastind[i]);</div><div class="line"><a name="l00345"></a><span class="lineno">  345</span>&#160;                }</div><div class="line"><a name="l00346"></a><span class="lineno">  346</span>&#160;                </div><div class="line"><a name="l00347"></a><span class="lineno">  347</span>&#160;                k_current=k_new;</div><div class="line"><a name="l00348"></a><span class="lineno">  348</span>&#160;            }</div><div class="line"><a name="l00349"></a><span class="lineno">  349</span>&#160;        }<span class="comment">// end of loop on fractures</span></div><div class="line"><a name="l00350"></a><span class="lineno">  350</span>&#160;    }</div><div class="line"><a name="l00351"></a><span class="lineno">  351</span>&#160;    </div><div class="line"><a name="l00352"></a><span class="lineno">  352</span>&#160;    frc=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[0]-1].<a class="code" href="structvertex.html#aa2cd06573422bc415cb8f38a67232de7">fracture</a>[0];</div><div class="line"><a name="l00353"></a><span class="lineno">  353</span>&#160;    firstn=<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[0];</div><div class="line"><a name="l00354"></a><span class="lineno">  354</span>&#160;    lastn=<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[0];</div><div class="line"><a name="l00355"></a><span class="lineno">  355</span>&#160;    </div><div class="line"><a name="l00356"></a><span class="lineno">  356</span>&#160;    <span class="keywordtype">double</span> thirdcoor=0.0;</div><div class="line"><a name="l00357"></a><span class="lineno">  357</span>&#160;    <span class="keywordtype">int</span> firstcoor=0, secondcoor=0, frc_count=0;</div><div class="line"><a name="l00358"></a><span class="lineno">  358</span>&#160;    <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[1]-1].<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[0]==frc)</div><div class="line"><a name="l00359"></a><span class="lineno">  359</span>&#160;    {</div><div class="line"><a name="l00360"></a><span class="lineno">  360</span>&#160;        <span class="keywordflow">if</span> (fabs(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[0]-1].coord[0])-fabs(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[1]-1].coord[0])&lt;1e-10)</div><div class="line"><a name="l00361"></a><span class="lineno">  361</span>&#160;        {</div><div class="line"><a name="l00362"></a><span class="lineno">  362</span>&#160;            firstcoor=1;</div><div class="line"><a name="l00363"></a><span class="lineno">  363</span>&#160;            secondcoor=2;</div><div class="line"><a name="l00364"></a><span class="lineno">  364</span>&#160;            </div><div class="line"><a name="l00365"></a><span class="lineno">  365</span>&#160;        }</div><div class="line"><a name="l00366"></a><span class="lineno">  366</span>&#160;        <span class="keywordflow">else</span></div><div class="line"><a name="l00367"></a><span class="lineno">  367</span>&#160;        {</div><div class="line"><a name="l00368"></a><span class="lineno">  368</span>&#160;            firstcoor=0;</div><div class="line"><a name="l00369"></a><span class="lineno">  369</span>&#160;            <span class="keywordflow">if</span> (fabs(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[0]-1].coord[1])-fabs(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[1]-1].coord[1])&lt;1e-10)</div><div class="line"><a name="l00370"></a><span class="lineno">  370</span>&#160;                secondcoor=2;</div><div class="line"><a name="l00371"></a><span class="lineno">  371</span>&#160;            <span class="keywordflow">else</span></div><div class="line"><a name="l00372"></a><span class="lineno">  372</span>&#160;                secondcoor=1;</div><div class="line"><a name="l00373"></a><span class="lineno">  373</span>&#160;        }</div><div class="line"><a name="l00374"></a><span class="lineno">  374</span>&#160;    }</div><div class="line"><a name="l00375"></a><span class="lineno">  375</span>&#160;    </div><div class="line"><a name="l00376"></a><span class="lineno">  376</span>&#160;    </div><div class="line"><a name="l00377"></a><span class="lineno">  377</span>&#160;    <span class="comment">/*** loop on all nodes in zone: define the boundary nodes for each fracture ****/</span></div><div class="line"><a name="l00378"></a><span class="lineno">  378</span>&#160;    <span class="comment">//  k_current=0;</span></div><div class="line"><a name="l00379"></a><span class="lineno">  379</span>&#160;    </div><div class="line"><a name="l00380"></a><span class="lineno">  380</span>&#160;    <span class="comment">//  fprintf(inp,&quot;%d\n&quot;, nzone_in);</span></div><div class="line"><a name="l00381"></a><span class="lineno">  381</span>&#160;    <span class="keywordflow">for</span> (i=0; i&lt;<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>-1; i++)</div><div class="line"><a name="l00382"></a><span class="lineno">  382</span>&#160;    {</div><div class="line"><a name="l00383"></a><span class="lineno">  383</span>&#160;        <span class="comment">//        printf(&quot; %d %d %d  %d %d %d\n&quot;, i,nodezonein[i], node[nodezonein[i]-1].fracture[0], node[nodezonein[i]-1].fracture[1], firstn, lastn);</span></div><div class="line"><a name="l00384"></a><span class="lineno">  384</span>&#160;        </div><div class="line"><a name="l00385"></a><span class="lineno">  385</span>&#160;        </div><div class="line"><a name="l00386"></a><span class="lineno">  386</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[0]==frc)</div><div class="line"><a name="l00387"></a><span class="lineno">  387</span>&#160;        {</div><div class="line"><a name="l00388"></a><span class="lineno">  388</span>&#160;            </div><div class="line"><a name="l00389"></a><span class="lineno">  389</span>&#160;            <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].coord[firstcoor]!=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].coord[firstcoor])</div><div class="line"><a name="l00390"></a><span class="lineno">  390</span>&#160;            {</div><div class="line"><a name="l00391"></a><span class="lineno">  391</span>&#160;                <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].coord[firstcoor]&lt;<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].coord[firstcoor])</div><div class="line"><a name="l00392"></a><span class="lineno">  392</span>&#160;                {</div><div class="line"><a name="l00393"></a><span class="lineno">  393</span>&#160;                    firstn=<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i];</div><div class="line"><a name="l00394"></a><span class="lineno">  394</span>&#160;                    first_ind=i;</div><div class="line"><a name="l00395"></a><span class="lineno">  395</span>&#160;                }</div><div class="line"><a name="l00396"></a><span class="lineno">  396</span>&#160;                <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].coord[firstcoor]&gt;<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastn-1].coord[firstcoor])</div><div class="line"><a name="l00397"></a><span class="lineno">  397</span>&#160;                {</div><div class="line"><a name="l00398"></a><span class="lineno">  398</span>&#160;                    lastn=<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i];</div><div class="line"><a name="l00399"></a><span class="lineno">  399</span>&#160;                    last_ind=i;</div><div class="line"><a name="l00400"></a><span class="lineno">  400</span>&#160;                }</div><div class="line"><a name="l00401"></a><span class="lineno">  401</span>&#160;                <span class="comment">//    fprintf(inp,&quot;first coord, %d %d\n&quot;,firstn, lastn);</span></div><div class="line"><a name="l00402"></a><span class="lineno">  402</span>&#160;            }</div><div class="line"><a name="l00403"></a><span class="lineno">  403</span>&#160;            <span class="keywordflow">else</span></div><div class="line"><a name="l00404"></a><span class="lineno">  404</span>&#160;            {</div><div class="line"><a name="l00405"></a><span class="lineno">  405</span>&#160;                <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].coord[secondcoor]&lt;<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].coord[secondcoor])</div><div class="line"><a name="l00406"></a><span class="lineno">  406</span>&#160;                {</div><div class="line"><a name="l00407"></a><span class="lineno">  407</span>&#160;                    firstn=<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i];</div><div class="line"><a name="l00408"></a><span class="lineno">  408</span>&#160;                    first_ind=i;</div><div class="line"><a name="l00409"></a><span class="lineno">  409</span>&#160;                }</div><div class="line"><a name="l00410"></a><span class="lineno">  410</span>&#160;                <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].coord[secondcoor]&gt;<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastn-1].coord[secondcoor])</div><div class="line"><a name="l00411"></a><span class="lineno">  411</span>&#160;                {</div><div class="line"><a name="l00412"></a><span class="lineno">  412</span>&#160;                    lastn=<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i];</div><div class="line"><a name="l00413"></a><span class="lineno">  413</span>&#160;                    last_ind=i;</div><div class="line"><a name="l00414"></a><span class="lineno">  414</span>&#160;                }</div><div class="line"><a name="l00415"></a><span class="lineno">  415</span>&#160;                <span class="comment">//      fprintf(inp,&quot;second coord, %d %d\n&quot;,firstn, lastn);</span></div><div class="line"><a name="l00416"></a><span class="lineno">  416</span>&#160;            }</div><div class="line"><a name="l00417"></a><span class="lineno">  417</span>&#160;            </div><div class="line"><a name="l00418"></a><span class="lineno">  418</span>&#160;            </div><div class="line"><a name="l00419"></a><span class="lineno">  419</span>&#160;        }</div><div class="line"><a name="l00420"></a><span class="lineno">  420</span>&#160;        <span class="comment">//  printf(&quot;fracture in flow-in zone %d first node %d last node %d \n&quot;, frc, firstn, lastn);</span></div><div class="line"><a name="l00421"></a><span class="lineno">  421</span>&#160;        <span class="keywordflow">if</span> ((<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[0]!=frc)||((i==<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>-2)))</div><div class="line"><a name="l00422"></a><span class="lineno">  422</span>&#160;        {</div><div class="line"><a name="l00423"></a><span class="lineno">  423</span>&#160;            <span class="comment">//     printf(&quot;fracture in flow-in zone %d first node %d last node %d \n&quot;, frc, firstn, lastn);</span></div><div class="line"><a name="l00424"></a><span class="lineno">  424</span>&#160;            <span class="keywordflow">if</span> ((i==<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>-2)  &amp;&amp; (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[0]==<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>-1]-1].<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[0]))</div><div class="line"><a name="l00425"></a><span class="lineno">  425</span>&#160;            {</div><div class="line"><a name="l00426"></a><span class="lineno">  426</span>&#160;                lastn=<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>-1];</div><div class="line"><a name="l00427"></a><span class="lineno">  427</span>&#160;                last_ind=<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>-1;</div><div class="line"><a name="l00428"></a><span class="lineno">  428</span>&#160;                </div><div class="line"><a name="l00429"></a><span class="lineno">  429</span>&#160;            }</div><div class="line"><a name="l00430"></a><span class="lineno">  430</span>&#160;            </div><div class="line"><a name="l00431"></a><span class="lineno">  431</span>&#160;            <span class="keywordflow">if</span> (firstn==lastn)</div><div class="line"><a name="l00432"></a><span class="lineno">  432</span>&#160;            {</div><div class="line"><a name="l00433"></a><span class="lineno">  433</span>&#160;                <span class="comment">//printf(&quot;firstn %d lastn %d fracture %d %d \n&quot;, firstn, lastn, node[nodezonein[i]-1].fracture[0],node[nodezonein[i]-1].fracture[1] );</span></div><div class="line"><a name="l00434"></a><span class="lineno">  434</span>&#160;            }</div><div class="line"><a name="l00435"></a><span class="lineno">  435</span>&#160;            <span class="keywordflow">if</span> (firstn!=lastn)</div><div class="line"><a name="l00436"></a><span class="lineno">  436</span>&#160;            {</div><div class="line"><a name="l00437"></a><span class="lineno">  437</span>&#160;                </div><div class="line"><a name="l00438"></a><span class="lineno">  438</span>&#160;                </div><div class="line"><a name="l00439"></a><span class="lineno">  439</span>&#160;                <span class="keywordflow">if</span> (flag_in==3)</div><div class="line"><a name="l00440"></a><span class="lineno">  440</span>&#160;                {</div><div class="line"><a name="l00441"></a><span class="lineno">  441</span>&#160;                    </div><div class="line"><a name="l00442"></a><span class="lineno">  442</span>&#160;                    <span class="keywordtype">double</span> cx1=0, cx2=0, cy1=0, cy2=0;</div><div class="line"><a name="l00443"></a><span class="lineno">  443</span>&#160;                    </div><div class="line"><a name="l00444"></a><span class="lineno">  444</span>&#160;                    <span class="comment">/* define ends points of fracture edge, then calculate an intersection with starting region */</span></div><div class="line"><a name="l00445"></a><span class="lineno">  445</span>&#160;                    inter_p[frc_count][0]=1e-10;</div><div class="line"><a name="l00446"></a><span class="lineno">  446</span>&#160;                    inter_p[frc_count][1]=1e-10;</div><div class="line"><a name="l00447"></a><span class="lineno">  447</span>&#160;                    inter_p[frc_count][2]=1e-10;</div><div class="line"><a name="l00448"></a><span class="lineno">  448</span>&#160;                    inter_p[frc_count][3]=1e-10;</div><div class="line"><a name="l00449"></a><span class="lineno">  449</span>&#160;                    <span class="keywordflow">if</span> ((zonenumb_in==1)||(zonenumb_in==2))</div><div class="line"><a name="l00450"></a><span class="lineno">  450</span>&#160;                    {</div><div class="line"><a name="l00451"></a><span class="lineno">  451</span>&#160;                        cx1=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[0];</div><div class="line"><a name="l00452"></a><span class="lineno">  452</span>&#160;                        cy1=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[1];</div><div class="line"><a name="l00453"></a><span class="lineno">  453</span>&#160;                        <span class="keywordflow">if</span> ((cx1&gt;ixmin) &amp;&amp; (cx1&lt;ixmax) &amp;&amp;(cy1&gt;iymin)&amp;&amp;(cy1&lt;iymax))</div><div class="line"><a name="l00454"></a><span class="lineno">  454</span>&#160;                        {</div><div class="line"><a name="l00455"></a><span class="lineno">  455</span>&#160;                            inter_p[frc_count][0]=cx1;</div><div class="line"><a name="l00456"></a><span class="lineno">  456</span>&#160;                            inter_p[frc_count][1]=cy1;</div><div class="line"><a name="l00457"></a><span class="lineno">  457</span>&#160;                        }</div><div class="line"><a name="l00458"></a><span class="lineno">  458</span>&#160;                        cx2=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastn-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[0];</div><div class="line"><a name="l00459"></a><span class="lineno">  459</span>&#160;                        cy2=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastn-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[1];</div><div class="line"><a name="l00460"></a><span class="lineno">  460</span>&#160;                        <span class="keywordflow">if</span> ((cx2&gt;ixmin) &amp;&amp; (cx2&lt;ixmax) &amp;&amp;(cy2&gt;iymin)&amp;&amp;(cy2&lt;iymax))</div><div class="line"><a name="l00461"></a><span class="lineno">  461</span>&#160;                        {</div><div class="line"><a name="l00462"></a><span class="lineno">  462</span>&#160;                            inter_p[frc_count][0]=cx2;</div><div class="line"><a name="l00463"></a><span class="lineno">  463</span>&#160;                            inter_p[frc_count][1]=cy2;</div><div class="line"><a name="l00464"></a><span class="lineno">  464</span>&#160;                        }</div><div class="line"><a name="l00465"></a><span class="lineno">  465</span>&#160;                        </div><div class="line"><a name="l00466"></a><span class="lineno">  466</span>&#160;                        thirdcoor= <a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[2];</div><div class="line"><a name="l00467"></a><span class="lineno">  467</span>&#160;                    }</div><div class="line"><a name="l00468"></a><span class="lineno">  468</span>&#160;                    </div><div class="line"><a name="l00469"></a><span class="lineno">  469</span>&#160;                    <span class="keywordflow">if</span> ((zonenumb_in==3)||(zonenumb_in==5))</div><div class="line"><a name="l00470"></a><span class="lineno">  470</span>&#160;                    {</div><div class="line"><a name="l00471"></a><span class="lineno">  471</span>&#160;                        cx1=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[2];</div><div class="line"><a name="l00472"></a><span class="lineno">  472</span>&#160;                        cy1=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[1];</div><div class="line"><a name="l00473"></a><span class="lineno">  473</span>&#160;                        <span class="keywordflow">if</span> ((cx1&gt;izmin) &amp;&amp; (cx1&lt;izmax) &amp;&amp;(cy1&gt;iymin)&amp;&amp;(cy1&lt;iymax))</div><div class="line"><a name="l00474"></a><span class="lineno">  474</span>&#160;                        {</div><div class="line"><a name="l00475"></a><span class="lineno">  475</span>&#160;                            inter_p[frc_count][0]=cx1;</div><div class="line"><a name="l00476"></a><span class="lineno">  476</span>&#160;                            inter_p[frc_count][1]=cy1;</div><div class="line"><a name="l00477"></a><span class="lineno">  477</span>&#160;                        }</div><div class="line"><a name="l00478"></a><span class="lineno">  478</span>&#160;                        </div><div class="line"><a name="l00479"></a><span class="lineno">  479</span>&#160;                        cx2=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastn-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[2];</div><div class="line"><a name="l00480"></a><span class="lineno">  480</span>&#160;                        cy2=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastn-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[1];</div><div class="line"><a name="l00481"></a><span class="lineno">  481</span>&#160;                        <span class="keywordflow">if</span> ((cx2&gt;izmin) &amp;&amp; (cx2&lt;izmax) &amp;&amp;(cy2&gt;iymin)&amp;&amp;(cy2&lt;iymax))</div><div class="line"><a name="l00482"></a><span class="lineno">  482</span>&#160;                        {</div><div class="line"><a name="l00483"></a><span class="lineno">  483</span>&#160;                            inter_p[frc_count][0]=cx2;</div><div class="line"><a name="l00484"></a><span class="lineno">  484</span>&#160;                            inter_p[frc_count][1]=cy2;</div><div class="line"><a name="l00485"></a><span class="lineno">  485</span>&#160;                        }</div><div class="line"><a name="l00486"></a><span class="lineno">  486</span>&#160;                        thirdcoor= <a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[0];</div><div class="line"><a name="l00487"></a><span class="lineno">  487</span>&#160;                    }</div><div class="line"><a name="l00488"></a><span class="lineno">  488</span>&#160;                    </div><div class="line"><a name="l00489"></a><span class="lineno">  489</span>&#160;                    <span class="keywordflow">if</span> ((zonenumb_in==4)||(zonenumb_in==6))</div><div class="line"><a name="l00490"></a><span class="lineno">  490</span>&#160;                    {</div><div class="line"><a name="l00491"></a><span class="lineno">  491</span>&#160;                        cx1=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[0];</div><div class="line"><a name="l00492"></a><span class="lineno">  492</span>&#160;                        cy1=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[2];</div><div class="line"><a name="l00493"></a><span class="lineno">  493</span>&#160;                        <span class="keywordflow">if</span> ((cx1&gt;ixmin) &amp;&amp; (cx1&lt;ixmax) &amp;&amp;(cy1&gt;izmin)&amp;&amp;(cy1&lt;izmax))</div><div class="line"><a name="l00494"></a><span class="lineno">  494</span>&#160;                        {</div><div class="line"><a name="l00495"></a><span class="lineno">  495</span>&#160;                            inter_p[frc_count][0]=cx1;</div><div class="line"><a name="l00496"></a><span class="lineno">  496</span>&#160;                            inter_p[frc_count][1]=cy1;</div><div class="line"><a name="l00497"></a><span class="lineno">  497</span>&#160;                        }</div><div class="line"><a name="l00498"></a><span class="lineno">  498</span>&#160;                        cx2=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastn-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[0];</div><div class="line"><a name="l00499"></a><span class="lineno">  499</span>&#160;                        cy2=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastn-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[2];</div><div class="line"><a name="l00500"></a><span class="lineno">  500</span>&#160;                        <span class="keywordflow">if</span> ((cx2&gt;ixmin) &amp;&amp; (cx2&lt;ixmax) &amp;&amp;(cy2&gt;izmin)&amp;&amp;(cy2&lt;izmax))</div><div class="line"><a name="l00501"></a><span class="lineno">  501</span>&#160;                        {</div><div class="line"><a name="l00502"></a><span class="lineno">  502</span>&#160;                            inter_p[frc_count][0]=cx2;</div><div class="line"><a name="l00503"></a><span class="lineno">  503</span>&#160;                            inter_p[frc_count][1]=cy2;</div><div class="line"><a name="l00504"></a><span class="lineno">  504</span>&#160;                        }</div><div class="line"><a name="l00505"></a><span class="lineno">  505</span>&#160;                        thirdcoor= <a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[1];</div><div class="line"><a name="l00506"></a><span class="lineno">  506</span>&#160;                    }</div><div class="line"><a name="l00507"></a><span class="lineno">  507</span>&#160;                    </div><div class="line"><a name="l00508"></a><span class="lineno">  508</span>&#160;                    <span class="keywordtype">double</span> pr1, pr2, pr3, pr4, p_x, p_y;</div><div class="line"><a name="l00509"></a><span class="lineno">  509</span>&#160;                    <span class="keywordtype">int</span> ii;</div><div class="line"><a name="l00510"></a><span class="lineno">  510</span>&#160;                    </div><div class="line"><a name="l00511"></a><span class="lineno">  511</span>&#160;                    <span class="comment">/* define intersection points of boundary fracture edges and starting region sides*/</span></div><div class="line"><a name="l00512"></a><span class="lineno">  512</span>&#160;                    <span class="keywordflow">for</span> (ii=0; ii&lt;4; ii++)</div><div class="line"><a name="l00513"></a><span class="lineno">  513</span>&#160;                    {</div><div class="line"><a name="l00514"></a><span class="lineno">  514</span>&#160;                        <span class="keywordtype">int</span> kk;</div><div class="line"><a name="l00515"></a><span class="lineno">  515</span>&#160;                        kk=ii+1;</div><div class="line"><a name="l00516"></a><span class="lineno">  516</span>&#160;                        <span class="keywordflow">if</span> (ii==3)</div><div class="line"><a name="l00517"></a><span class="lineno">  517</span>&#160;                            kk=0;</div><div class="line"><a name="l00518"></a><span class="lineno">  518</span>&#160;                        pr1=(px[ii]-cx1)*(py[kk]-cy1)-(py[ii]-cy1)*(px[kk]-cx1);</div><div class="line"><a name="l00519"></a><span class="lineno">  519</span>&#160;                        pr2=(cx1-px[ii])*(cy2-py[ii])-(cy1-py[ii])*(cx2-px[ii]);</div><div class="line"><a name="l00520"></a><span class="lineno">  520</span>&#160;                        pr3=(px[ii]-cx2)*(py[kk]-cy2)-(py[ii]-cy2)*(px[kk]-cx2);</div><div class="line"><a name="l00521"></a><span class="lineno">  521</span>&#160;                        pr4=(cx1-px[kk])*(cy2-py[kk])-(cy1-py[kk])*(cx2-px[kk]);</div><div class="line"><a name="l00522"></a><span class="lineno">  522</span>&#160;                        </div><div class="line"><a name="l00523"></a><span class="lineno">  523</span>&#160;                        <span class="keywordflow">if</span> ((pr1*pr3&lt;0)&amp;&amp;(pr2*pr4&lt;0))</div><div class="line"><a name="l00524"></a><span class="lineno">  524</span>&#160;                        {</div><div class="line"><a name="l00525"></a><span class="lineno">  525</span>&#160;                            </div><div class="line"><a name="l00526"></a><span class="lineno">  526</span>&#160;                            pr1=cx1*cy2-cy1*cx2;</div><div class="line"><a name="l00527"></a><span class="lineno">  527</span>&#160;                            pr2=px[ii]*py[kk]-py[ii]*px[kk];</div><div class="line"><a name="l00528"></a><span class="lineno">  528</span>&#160;                            pr3=(cx1-cx2)*(py[ii]-py[kk])-(cy1-cy2)*(px[ii]-px[kk]);</div><div class="line"><a name="l00529"></a><span class="lineno">  529</span>&#160;                            p_x=((px[ii]-px[kk])*pr1-(cx1-cx2)*pr2)/pr3;</div><div class="line"><a name="l00530"></a><span class="lineno">  530</span>&#160;                            p_y=((py[ii]-py[kk])*pr1-(cy1-cy2)*pr2)/pr3;</div><div class="line"><a name="l00531"></a><span class="lineno">  531</span>&#160;                            <span class="keywordflow">if</span> (inter_p[frc_count][0]==1e-10)</div><div class="line"><a name="l00532"></a><span class="lineno">  532</span>&#160;                            {</div><div class="line"><a name="l00533"></a><span class="lineno">  533</span>&#160;                                inter_p[frc_count][0]=p_x;</div><div class="line"><a name="l00534"></a><span class="lineno">  534</span>&#160;                                inter_p[frc_count][1]=p_y;</div><div class="line"><a name="l00535"></a><span class="lineno">  535</span>&#160;                            }</div><div class="line"><a name="l00536"></a><span class="lineno">  536</span>&#160;                            <span class="keywordflow">else</span></div><div class="line"><a name="l00537"></a><span class="lineno">  537</span>&#160;                            {</div><div class="line"><a name="l00538"></a><span class="lineno">  538</span>&#160;                                inter_fr[frc_count]=frc;</div><div class="line"><a name="l00539"></a><span class="lineno">  539</span>&#160;                                inter_p[frc_count][2]=p_x;</div><div class="line"><a name="l00540"></a><span class="lineno">  540</span>&#160;                                inter_p[frc_count][3]=p_y;</div><div class="line"><a name="l00541"></a><span class="lineno">  541</span>&#160;                                frc_count++;</div><div class="line"><a name="l00542"></a><span class="lineno">  542</span>&#160;                                <span class="keywordflow">break</span>;</div><div class="line"><a name="l00543"></a><span class="lineno">  543</span>&#160;                            }</div><div class="line"><a name="l00544"></a><span class="lineno">  544</span>&#160;                            </div><div class="line"><a name="l00545"></a><span class="lineno">  545</span>&#160;                        }</div><div class="line"><a name="l00546"></a><span class="lineno">  546</span>&#160;                        </div><div class="line"><a name="l00547"></a><span class="lineno">  547</span>&#160;                    }</div><div class="line"><a name="l00548"></a><span class="lineno">  548</span>&#160;                    </div><div class="line"><a name="l00549"></a><span class="lineno">  549</span>&#160;                    </div><div class="line"><a name="l00550"></a><span class="lineno">  550</span>&#160;                } <span class="comment">//end of flag</span></div><div class="line"><a name="l00551"></a><span class="lineno">  551</span>&#160;            }</div><div class="line"><a name="l00552"></a><span class="lineno">  552</span>&#160;            </div><div class="line"><a name="l00553"></a><span class="lineno">  553</span>&#160;            <span class="keywordflow">if</span> (i&lt;<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>-2)</div><div class="line"><a name="l00554"></a><span class="lineno">  554</span>&#160;            {</div><div class="line"><a name="l00555"></a><span class="lineno">  555</span>&#160;                frc=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].<a class="code" href="structvertex.html#aa2cd06573422bc415cb8f38a67232de7">fracture</a>[0];</div><div class="line"><a name="l00556"></a><span class="lineno">  556</span>&#160;                </div><div class="line"><a name="l00557"></a><span class="lineno">  557</span>&#160;                firstn=<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i];</div><div class="line"><a name="l00558"></a><span class="lineno">  558</span>&#160;                first_ind=i;</div><div class="line"><a name="l00559"></a><span class="lineno">  559</span>&#160;                lastn=<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i];</div><div class="line"><a name="l00560"></a><span class="lineno">  560</span>&#160;                last_ind=i;</div><div class="line"><a name="l00561"></a><span class="lineno">  561</span>&#160;                <span class="comment">//        fprintf(inp,&quot;p %d %d %d %d\n&quot;, i, frc, firstn, lastn);</span></div><div class="line"><a name="l00562"></a><span class="lineno">  562</span>&#160;                <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i+1]-1].<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[0]==frc)</div><div class="line"><a name="l00563"></a><span class="lineno">  563</span>&#160;                {</div><div class="line"><a name="l00564"></a><span class="lineno">  564</span>&#160;                    <span class="keywordflow">if</span> (abs(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].coord[0])-abs(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i+1]-1].coord[0])&lt;1e-10)</div><div class="line"><a name="l00565"></a><span class="lineno">  565</span>&#160;                    {</div><div class="line"><a name="l00566"></a><span class="lineno">  566</span>&#160;                        firstcoor=1;</div><div class="line"><a name="l00567"></a><span class="lineno">  567</span>&#160;                        secondcoor=2;</div><div class="line"><a name="l00568"></a><span class="lineno">  568</span>&#160;                    }</div><div class="line"><a name="l00569"></a><span class="lineno">  569</span>&#160;                    <span class="keywordflow">else</span></div><div class="line"><a name="l00570"></a><span class="lineno">  570</span>&#160;                    {</div><div class="line"><a name="l00571"></a><span class="lineno">  571</span>&#160;                        firstcoor=0;</div><div class="line"><a name="l00572"></a><span class="lineno">  572</span>&#160;                        </div><div class="line"><a name="l00573"></a><span class="lineno">  573</span>&#160;                        <span class="keywordflow">if</span> (abs(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].coord[1])-abs(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i+1]-1].coord[1])&lt;1e-10)</div><div class="line"><a name="l00574"></a><span class="lineno">  574</span>&#160;                            </div><div class="line"><a name="l00575"></a><span class="lineno">  575</span>&#160;                            secondcoor=2;</div><div class="line"><a name="l00576"></a><span class="lineno">  576</span>&#160;                        <span class="keywordflow">else</span></div><div class="line"><a name="l00577"></a><span class="lineno">  577</span>&#160;                            secondcoor=1;</div><div class="line"><a name="l00578"></a><span class="lineno">  578</span>&#160;                    }</div><div class="line"><a name="l00579"></a><span class="lineno">  579</span>&#160;                }</div><div class="line"><a name="l00580"></a><span class="lineno">  580</span>&#160;            }</div><div class="line"><a name="l00581"></a><span class="lineno">  581</span>&#160;            </div><div class="line"><a name="l00582"></a><span class="lineno">  582</span>&#160;            </div><div class="line"><a name="l00583"></a><span class="lineno">  583</span>&#160;        }</div><div class="line"><a name="l00584"></a><span class="lineno">  584</span>&#160;    }</div><div class="line"><a name="l00585"></a><span class="lineno">  585</span>&#160;    <span class="comment">//   printf(&quot;fracture in flow-in zone %d first node %d last node %d \n&quot;, frc, firstn, lastn);</span></div><div class="line"><a name="l00586"></a><span class="lineno">  586</span>&#160;    <span class="keywordflow">if</span> ((flag_in==3)&amp;&amp;(frc_count==0))</div><div class="line"><a name="l00587"></a><span class="lineno">  587</span>&#160;    {</div><div class="line"><a name="l00588"></a><span class="lineno">  588</span>&#160;        printf(<span class="stringliteral">&quot;\n There is no fracture crosses the given range. Try to increase the range. \n&quot;</span>);</div><div class="line"><a name="l00589"></a><span class="lineno">  589</span>&#160;        printf(<span class="stringliteral">&quot;\n Program is terminated. \n&quot;</span>);</div><div class="line"><a name="l00590"></a><span class="lineno">  590</span>&#160;        exit(1);</div><div class="line"><a name="l00591"></a><span class="lineno">  591</span>&#160;    }</div><div class="line"><a name="l00592"></a><span class="lineno">  592</span>&#160;    </div><div class="line"><a name="l00593"></a><span class="lineno">  593</span>&#160;    <span class="comment">/* place particles in starting region: every fracture (or part of fracture</span></div><div class="line"><a name="l00594"></a><span class="lineno">  594</span>&#160;<span class="comment">     edge will have the same amount of particles */</span></div><div class="line"><a name="l00595"></a><span class="lineno">  595</span>&#160;    <span class="keywordflow">if</span> (flag_in==3)</div><div class="line"><a name="l00596"></a><span class="lineno">  596</span>&#160;    {</div><div class="line"><a name="l00597"></a><span class="lineno">  597</span>&#160;        <span class="keywordflow">for</span> (i=0; i&lt;frc_count; i++)</div><div class="line"><a name="l00598"></a><span class="lineno">  598</span>&#160;        {</div><div class="line"><a name="l00599"></a><span class="lineno">  599</span>&#160;            parts_fracture=(int) <a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>/frc_count;</div><div class="line"><a name="l00600"></a><span class="lineno">  600</span>&#160;            </div><div class="line"><a name="l00601"></a><span class="lineno">  601</span>&#160;            k_new=<a class="code" href="_initial_part_positions_8c.html#a2e978219b990b7dfb64b4f9d567054db">InitParticles_ones</a> (k_current, inter_p, inter_fr[i], parts_fracture, i, thirdcoor, zonenumb_in, first_ind, last_ind);</div><div class="line"><a name="l00602"></a><span class="lineno">  602</span>&#160;            k_current=k_new;</div><div class="line"><a name="l00603"></a><span class="lineno">  603</span>&#160;            <span class="comment">//   printf(&quot; %d intersects at %f %f %f %f\n&quot;,i, inter_p[i][0], inter_p[i][1], inter_p[i][2], inter_p[i][3]);</span></div><div class="line"><a name="l00604"></a><span class="lineno">  604</span>&#160;        }</div><div class="line"><a name="l00605"></a><span class="lineno">  605</span>&#160;    }</div><div class="line"><a name="l00606"></a><span class="lineno">  606</span>&#160;    </div><div class="line"><a name="l00607"></a><span class="lineno">  607</span>&#160;    <span class="comment">// fclose(inp);</span></div><div class="line"><a name="l00608"></a><span class="lineno">  608</span>&#160;    <span class="keywordflow">if</span> (flag_in==0)</div><div class="line"><a name="l00609"></a><span class="lineno">  609</span>&#160;    {</div><div class="line"><a name="l00610"></a><span class="lineno">  610</span>&#160;        printf(<span class="stringliteral">&quot;\n There is no option specified for particles initial positions! \n&quot;</span>);</div><div class="line"><a name="l00611"></a><span class="lineno">  611</span>&#160;        printf(<span class="stringliteral">&quot;\n Program is terminated. \n&quot;</span>);</div><div class="line"><a name="l00612"></a><span class="lineno">  612</span>&#160;        exit(1);</div><div class="line"><a name="l00613"></a><span class="lineno">  613</span>&#160;    }</div><div class="line"><a name="l00614"></a><span class="lineno">  614</span>&#160;    </div><div class="line"><a name="l00615"></a><span class="lineno">  615</span>&#160;    <span class="keywordflow">return</span> k_new;</div><div class="line"><a name="l00616"></a><span class="lineno">  616</span>&#160;}</div><div class="line"><a name="l00617"></a><span class="lineno">  617</span>&#160;</div><div class="line"><a name="l00618"></a><span class="lineno">  618</span>&#160;</div><div class="line"><a name="l00620"></a><span class="lineno">  620</span>&#160;</div><div class="line"><a name="l00621"></a><span class="lineno"><a class="line" href="_initial_part_positions_8c.html#a461479d92dd1f67c7ef94fa953186860">  621</a></span>&#160;<span class="keywordtype">int</span> <a class="code" href="_initial_part_positions_8c.html#a461479d92dd1f67c7ef94fa953186860">InitCell</a> ()</div><div class="line"><a name="l00623"></a><span class="lineno">  623</span>&#160;{</div><div class="line"><a name="l00624"></a><span class="lineno">  624</span>&#160;    </div><div class="line"><a name="l00625"></a><span class="lineno">  625</span>&#160;    <span class="keywordtype">int</span> i, j, k, curcel, insc=0;</div><div class="line"><a name="l00626"></a><span class="lineno">  626</span>&#160;    </div><div class="line"><a name="l00627"></a><span class="lineno">  627</span>&#160;    </div><div class="line"><a name="l00628"></a><span class="lineno">  628</span>&#160;    <span class="keywordflow">for</span> (i=0; i&lt;<a class="code" href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a>; i++)</div><div class="line"><a name="l00629"></a><span class="lineno">  629</span>&#160;    {</div><div class="line"><a name="l00630"></a><span class="lineno">  630</span>&#160;        </div><div class="line"><a name="l00631"></a><span class="lineno">  631</span>&#160;        <span class="keywordflow">for</span> (j=0; j&lt;<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].<a class="code" href="structvertex.html#acf83b561551952403148f73fe032f57b">numneighb</a>; j++)</div><div class="line"><a name="l00632"></a><span class="lineno">  632</span>&#160;        {</div><div class="line"><a name="l00633"></a><span class="lineno">  633</span>&#160;            <span class="keywordflow">for</span> (k=0; k&lt;4; k++)</div><div class="line"><a name="l00634"></a><span class="lineno">  634</span>&#160;                <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].cells[j][k]!=0)</div><div class="line"><a name="l00635"></a><span class="lineno">  635</span>&#160;                {</div><div class="line"><a name="l00636"></a><span class="lineno">  636</span>&#160;                    curcel=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[i]-1].<a class="code" href="structvertex.html#a583032cadaa99e8f5c044d3bd10191cf">cells</a>[j][k];</div><div class="line"><a name="l00637"></a><span class="lineno">  637</span>&#160;                    insc=<a class="code" href="_func_def_8h.html#a5c8fe8978ce1689f5dcf12bc7b56793d">InsideCell</a> (curcel);</div><div class="line"><a name="l00638"></a><span class="lineno">  638</span>&#160;                    <span class="keywordflow">if</span> (insc==1)</div><div class="line"><a name="l00639"></a><span class="lineno">  639</span>&#160;                        <span class="keywordflow">break</span>;</div><div class="line"><a name="l00640"></a><span class="lineno">  640</span>&#160;                }</div><div class="line"><a name="l00641"></a><span class="lineno">  641</span>&#160;            <span class="keywordflow">if</span> (insc==1)</div><div class="line"><a name="l00642"></a><span class="lineno">  642</span>&#160;                <span class="keywordflow">break</span>;</div><div class="line"><a name="l00643"></a><span class="lineno">  643</span>&#160;        }</div><div class="line"><a name="l00644"></a><span class="lineno">  644</span>&#160;        </div><div class="line"><a name="l00645"></a><span class="lineno">  645</span>&#160;        <span class="keywordflow">if</span> (insc==1)</div><div class="line"><a name="l00646"></a><span class="lineno">  646</span>&#160;            <span class="keywordflow">break</span>;</div><div class="line"><a name="l00647"></a><span class="lineno">  647</span>&#160;    }</div><div class="line"><a name="l00648"></a><span class="lineno">  648</span>&#160;    <span class="comment">//   if (insc==0)</span></div><div class="line"><a name="l00649"></a><span class="lineno">  649</span>&#160;    <span class="comment">//  {</span></div><div class="line"><a name="l00650"></a><span class="lineno">  650</span>&#160;    <span class="comment">// printf(&quot;Init Cell: NOT FOUND \n&quot;);</span></div><div class="line"><a name="l00651"></a><span class="lineno">  651</span>&#160;    </div><div class="line"><a name="l00652"></a><span class="lineno">  652</span>&#160;    <span class="comment">//  printf(&quot;%d find cell: particle cell %d   \n&quot;, np, particle[np].cell);</span></div><div class="line"><a name="l00653"></a><span class="lineno">  653</span>&#160;    </div><div class="line"><a name="l00654"></a><span class="lineno">  654</span>&#160;    <span class="comment">//      }</span></div><div class="line"><a name="l00655"></a><span class="lineno">  655</span>&#160;    </div><div class="line"><a name="l00656"></a><span class="lineno">  656</span>&#160;    <span class="keywordflow">return</span> insc;</div><div class="line"><a name="l00657"></a><span class="lineno">  657</span>&#160;}</div><div class="line"><a name="l00659"></a><span class="lineno">  659</span>&#160;</div><div class="line"><a name="l00660"></a><span class="lineno"><a class="line" href="_initial_part_positions_8c.html#a30b78ca10a78d9e5d6de6d5831a5557a">  660</a></span>&#160;<span class="keywordtype">int</span> <a class="code" href="_initial_part_positions_8c.html#a30b78ca10a78d9e5d6de6d5831a5557a">InitParticles_np</a> (<span class="keywordtype">int</span> k_current, <span class="keywordtype">int</span> firstnd, <span class="keywordtype">int</span> lastnd, <span class="keywordtype">int</span> parts_fracture, <span class="keywordtype">int</span> first_ind, <span class="keywordtype">int</span> last_ind)</div><div class="line"><a name="l00662"></a><span class="lineno">  662</span>&#160;{</div><div class="line"><a name="l00663"></a><span class="lineno">  663</span>&#160;    </div><div class="line"><a name="l00664"></a><span class="lineno">  664</span>&#160;    <span class="keywordtype">double</span> deltax, deltay;</div><div class="line"><a name="l00665"></a><span class="lineno">  665</span>&#160;    <span class="keywordtype">int</span> j, pf;</div><div class="line"><a name="l00666"></a><span class="lineno">  666</span>&#160;    pf=parts_fracture;</div><div class="line"><a name="l00667"></a><span class="lineno">  667</span>&#160;    </div><div class="line"><a name="l00668"></a><span class="lineno">  668</span>&#160;    <span class="keywordtype">int</span> ii,jj,  k, curcel, insc=0;</div><div class="line"><a name="l00669"></a><span class="lineno">  669</span>&#160;    </div><div class="line"><a name="l00670"></a><span class="lineno">  670</span>&#160;    deltax=fabs(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastnd-1].coord_xy[0]-<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstnd-1].coord_xy[0]);</div><div class="line"><a name="l00671"></a><span class="lineno">  671</span>&#160;    </div><div class="line"><a name="l00672"></a><span class="lineno">  672</span>&#160;    deltay=fabs(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastnd-1].coord_xy[1]-<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstnd-1].coord_xy[1]);</div><div class="line"><a name="l00673"></a><span class="lineno">  673</span>&#160;    </div><div class="line"><a name="l00674"></a><span class="lineno">  674</span>&#160;    </div><div class="line"><a name="l00675"></a><span class="lineno">  675</span>&#160;    <span class="keywordflow">for</span> (j=0; j&lt;pf; j++)</div><div class="line"><a name="l00676"></a><span class="lineno">  676</span>&#160;    {</div><div class="line"><a name="l00677"></a><span class="lineno">  677</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstnd-1].coord_xy[0]&lt;<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastnd-1].coord_xy[0])</div><div class="line"><a name="l00678"></a><span class="lineno">  678</span>&#160;            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#afdb17185baff7dc4f22cd6c76c6fd6fd">position</a>[0]=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstnd-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[0]+(deltax/pf)*(j)+deltax/(2.0*pf);</div><div class="line"><a name="l00679"></a><span class="lineno">  679</span>&#160;        <span class="keywordflow">else</span></div><div class="line"><a name="l00680"></a><span class="lineno">  680</span>&#160;            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#afdb17185baff7dc4f22cd6c76c6fd6fd">position</a>[0]=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstnd-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[0]-(deltax/pf)*(j)-deltax/(2.0*pf);</div><div class="line"><a name="l00681"></a><span class="lineno">  681</span>&#160;        </div><div class="line"><a name="l00682"></a><span class="lineno">  682</span>&#160;        </div><div class="line"><a name="l00683"></a><span class="lineno">  683</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstnd-1].coord_xy[1]&lt;<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastnd-1].coord_xy[1])</div><div class="line"><a name="l00684"></a><span class="lineno">  684</span>&#160;            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#afdb17185baff7dc4f22cd6c76c6fd6fd">position</a>[1]=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstnd-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[1]+(deltay/pf)*(j)+deltay/(2.0*pf);</div><div class="line"><a name="l00685"></a><span class="lineno">  685</span>&#160;        <span class="keywordflow">else</span></div><div class="line"><a name="l00686"></a><span class="lineno">  686</span>&#160;            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#afdb17185baff7dc4f22cd6c76c6fd6fd">position</a>[1]=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstnd-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[1]-(deltay/pf)*(j)-deltay/(2.0*pf);</div><div class="line"><a name="l00687"></a><span class="lineno">  687</span>&#160;        </div><div class="line"><a name="l00688"></a><span class="lineno">  688</span>&#160;        </div><div class="line"><a name="l00689"></a><span class="lineno">  689</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a70e44829116675e21bda7422f0c2b777">velocity</a>[0]=0.;</div><div class="line"><a name="l00690"></a><span class="lineno">  690</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a70e44829116675e21bda7422f0c2b777">velocity</a>[1]=0.;</div><div class="line"><a name="l00691"></a><span class="lineno">  691</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a0b9253ce268d781958145f9d110df2de">fracture</a>=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstnd-1].<a class="code" href="structvertex.html#aa2cd06573422bc415cb8f38a67232de7">fracture</a>[0];</div><div class="line"><a name="l00692"></a><span class="lineno">  692</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a1747a7b7a63f22656afa1df6462279e8">intcell</a>=0;</div><div class="line"><a name="l00693"></a><span class="lineno">  693</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a11c5e87e6e62d258648d5a02472251d5">time</a>=0.0;</div><div class="line"><a name="l00694"></a><span class="lineno">  694</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#ad50f20b01e37af1fedfd7bbe0cd1d7e3">fl_weight</a>=0.0;</div><div class="line"><a name="l00695"></a><span class="lineno">  695</span>&#160;        </div><div class="line"><a name="l00696"></a><span class="lineno">  696</span>&#160;        </div><div class="line"><a name="l00697"></a><span class="lineno">  697</span>&#160;        k_current++;</div><div class="line"><a name="l00698"></a><span class="lineno">  698</span>&#160;        </div><div class="line"><a name="l00699"></a><span class="lineno">  699</span>&#160;        <span class="keywordflow">if</span> (k_current&gt;<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>)</div><div class="line"><a name="l00700"></a><span class="lineno">  700</span>&#160;        {</div><div class="line"><a name="l00701"></a><span class="lineno">  701</span>&#160;            printf(<span class="stringliteral">&quot; \n Number of particles with allocated memory is less than number of particles set up initially. \n&quot;</span>);</div><div class="line"><a name="l00702"></a><span class="lineno">  702</span>&#160;            printf(<span class="stringliteral">&quot; Increase the number of particles. Program is terminated. \n&quot;</span>);</div><div class="line"><a name="l00703"></a><span class="lineno">  703</span>&#160;            exit(1);</div><div class="line"><a name="l00704"></a><span class="lineno">  704</span>&#160;            </div><div class="line"><a name="l00705"></a><span class="lineno">  705</span>&#160;        }</div><div class="line"><a name="l00706"></a><span class="lineno">  706</span>&#160;    }</div><div class="line"><a name="l00707"></a><span class="lineno">  707</span>&#160;    <span class="keywordflow">return</span> k_current;</div><div class="line"><a name="l00708"></a><span class="lineno">  708</span>&#160;}</div><div class="line"><a name="l00710"></a><span class="lineno">  710</span>&#160;</div><div class="line"><a name="l00711"></a><span class="lineno"><a class="line" href="_initial_part_positions_8c.html#a01e7bd2a5bd9417f1d28849b8577b0b9">  711</a></span>&#160;<span class="keywordtype">int</span> <a class="code" href="_initial_part_positions_8c.html#a01e7bd2a5bd9417f1d28849b8577b0b9">InitParticles_eq</a> (<span class="keywordtype">int</span> k_current, <span class="keywordtype">int</span> firstn, <span class="keywordtype">int</span> lastn, <span class="keywordtype">double</span> parts_dist, <span class="keywordtype">int</span> first_ind, <span class="keywordtype">int</span> last_ind)</div><div class="line"><a name="l00713"></a><span class="lineno">  713</span>&#160;{</div><div class="line"><a name="l00714"></a><span class="lineno">  714</span>&#160;    <span class="keywordtype">double</span> deltax, deltay, edgelength, eqdist_x, eqdist_y;</div><div class="line"><a name="l00715"></a><span class="lineno">  715</span>&#160;    <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> j, pf;</div><div class="line"><a name="l00716"></a><span class="lineno">  716</span>&#160;    </div><div class="line"><a name="l00717"></a><span class="lineno">  717</span>&#160;    deltax=(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastn-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[0]-<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[0]);</div><div class="line"><a name="l00718"></a><span class="lineno">  718</span>&#160;    deltay=(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[lastn-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[1]-<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[1]);</div><div class="line"><a name="l00719"></a><span class="lineno">  719</span>&#160;    edgelength=sqrt(deltax*deltax+deltay*deltay);</div><div class="line"><a name="l00720"></a><span class="lineno">  720</span>&#160;    pf=(int)(edgelength/parts_dist);</div><div class="line"><a name="l00721"></a><span class="lineno">  721</span>&#160;    <span class="keywordflow">if</span> (pf&lt;2)</div><div class="line"><a name="l00722"></a><span class="lineno">  722</span>&#160;    {</div><div class="line"><a name="l00723"></a><span class="lineno">  723</span>&#160;        pf=1;</div><div class="line"><a name="l00724"></a><span class="lineno">  724</span>&#160;        eqdist_x=deltax/2.0;</div><div class="line"><a name="l00725"></a><span class="lineno">  725</span>&#160;        eqdist_y=deltay/2.0;</div><div class="line"><a name="l00726"></a><span class="lineno">  726</span>&#160;    }</div><div class="line"><a name="l00727"></a><span class="lineno">  727</span>&#160;    <span class="keywordflow">else</span></div><div class="line"><a name="l00728"></a><span class="lineno">  728</span>&#160;    {</div><div class="line"><a name="l00729"></a><span class="lineno">  729</span>&#160;        eqdist_x=deltax/pf;</div><div class="line"><a name="l00730"></a><span class="lineno">  730</span>&#160;        eqdist_y=deltay/pf;</div><div class="line"><a name="l00731"></a><span class="lineno">  731</span>&#160;    }</div><div class="line"><a name="l00732"></a><span class="lineno">  732</span>&#160;    </div><div class="line"><a name="l00733"></a><span class="lineno">  733</span>&#160;    <span class="keywordflow">for</span> (j=0; j&lt;pf; j++)</div><div class="line"><a name="l00734"></a><span class="lineno">  734</span>&#160;    {</div><div class="line"><a name="l00735"></a><span class="lineno">  735</span>&#160;        </div><div class="line"><a name="l00736"></a><span class="lineno">  736</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#afdb17185baff7dc4f22cd6c76c6fd6fd">position</a>[0]=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[0]+eqdist_x*(j)+eqdist_x/2.0;</div><div class="line"><a name="l00737"></a><span class="lineno">  737</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#afdb17185baff7dc4f22cd6c76c6fd6fd">position</a>[1]=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[1]+eqdist_y*(j)+eqdist_y/2.0;</div><div class="line"><a name="l00738"></a><span class="lineno">  738</span>&#160;        </div><div class="line"><a name="l00739"></a><span class="lineno">  739</span>&#160;        </div><div class="line"><a name="l00740"></a><span class="lineno">  740</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a70e44829116675e21bda7422f0c2b777">velocity</a>[0]=0.;</div><div class="line"><a name="l00741"></a><span class="lineno">  741</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a70e44829116675e21bda7422f0c2b777">velocity</a>[1]=0.;</div><div class="line"><a name="l00742"></a><span class="lineno">  742</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a0b9253ce268d781958145f9d110df2de">fracture</a>=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[firstn-1].<a class="code" href="structvertex.html#aa2cd06573422bc415cb8f38a67232de7">fracture</a>[0];</div><div class="line"><a name="l00743"></a><span class="lineno">  743</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a1747a7b7a63f22656afa1df6462279e8">intcell</a>=0;</div><div class="line"><a name="l00744"></a><span class="lineno">  744</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a11c5e87e6e62d258648d5a02472251d5">time</a>=0.0;</div><div class="line"><a name="l00745"></a><span class="lineno">  745</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#ad50f20b01e37af1fedfd7bbe0cd1d7e3">fl_weight</a>=0.0;</div><div class="line"><a name="l00746"></a><span class="lineno">  746</span>&#160;        </div><div class="line"><a name="l00747"></a><span class="lineno">  747</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>==13)</div><div class="line"><a name="l00748"></a><span class="lineno">  748</span>&#160;            printf(<span class="stringliteral">&quot;os %lf %lf \n&quot;</span>, <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#afdb17185baff7dc4f22cd6c76c6fd6fd">position</a>[0], <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#afdb17185baff7dc4f22cd6c76c6fd6fd">position</a>[1]);</div><div class="line"><a name="l00749"></a><span class="lineno">  749</span>&#160;        </div><div class="line"><a name="l00750"></a><span class="lineno">  750</span>&#160;        </div><div class="line"><a name="l00751"></a><span class="lineno">  751</span>&#160;        k_current++;</div><div class="line"><a name="l00752"></a><span class="lineno">  752</span>&#160;        </div><div class="line"><a name="l00753"></a><span class="lineno">  753</span>&#160;        <span class="keywordflow">if</span> (k_current&gt;<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>)</div><div class="line"><a name="l00754"></a><span class="lineno">  754</span>&#160;        {</div><div class="line"><a name="l00755"></a><span class="lineno">  755</span>&#160;            printf(<span class="stringliteral">&quot;\n Number of particles with allocated memory is less than number of particles set up initially. \n&quot;</span>);</div><div class="line"><a name="l00756"></a><span class="lineno">  756</span>&#160;            printf(<span class="stringliteral">&quot; Increase the number of particles. Program is terminated. \n&quot;</span>);</div><div class="line"><a name="l00757"></a><span class="lineno">  757</span>&#160;            exit(1);</div><div class="line"><a name="l00758"></a><span class="lineno">  758</span>&#160;        }</div><div class="line"><a name="l00759"></a><span class="lineno">  759</span>&#160;    }</div><div class="line"><a name="l00760"></a><span class="lineno">  760</span>&#160;    <span class="keywordflow">return</span> k_current;</div><div class="line"><a name="l00761"></a><span class="lineno">  761</span>&#160;}</div><div class="line"><a name="l00763"></a><span class="lineno">  763</span>&#160;</div><div class="line"><a name="l00764"></a><span class="lineno"><a class="line" href="_initial_part_positions_8c.html#a2e978219b990b7dfb64b4f9d567054db">  764</a></span>&#160;<span class="keywordtype">int</span> <a class="code" href="_initial_part_positions_8c.html#a2e978219b990b7dfb64b4f9d567054db">InitParticles_ones</a> (<span class="keywordtype">int</span> k_current, <span class="keywordtype">double</span> inter_p[][4], <span class="keywordtype">int</span> fracture_n, <span class="keywordtype">int</span> parts_fracture, <span class="keywordtype">int</span> ii, <span class="keywordtype">double</span> thirdcoor, <span class="keywordtype">int</span> zonenumb_in, <span class="keywordtype">int</span> first_ind, <span class="keywordtype">int</span> last_ind)</div><div class="line"><a name="l00766"></a><span class="lineno">  766</span>&#160;{</div><div class="line"><a name="l00767"></a><span class="lineno">  767</span>&#160;    <span class="keywordtype">int</span> j=fracture_n-1;</div><div class="line"><a name="l00768"></a><span class="lineno">  768</span>&#160;    <span class="keywordtype">double</span> x_1=0, y_1=0, z_1=0, x_2=0, z_2=0, y_2=0;</div><div class="line"><a name="l00769"></a><span class="lineno">  769</span>&#160;    <span class="keywordtype">double</span> x1cor=0.0, y1cor=0.0, z1cor=0, x2cor=0.0, y2cor=0.0, z2cor=0;</div><div class="line"><a name="l00770"></a><span class="lineno">  770</span>&#160;    </div><div class="line"><a name="l00771"></a><span class="lineno">  771</span>&#160;    <span class="keywordflow">if</span> ((zonenumb_in==1)||(zonenumb_in==2))</div><div class="line"><a name="l00772"></a><span class="lineno">  772</span>&#160;    {</div><div class="line"><a name="l00773"></a><span class="lineno">  773</span>&#160;        x1cor=inter_p[ii][0];</div><div class="line"><a name="l00774"></a><span class="lineno">  774</span>&#160;        y1cor=inter_p[ii][1];</div><div class="line"><a name="l00775"></a><span class="lineno">  775</span>&#160;        z1cor=thirdcoor;</div><div class="line"><a name="l00776"></a><span class="lineno">  776</span>&#160;        x2cor=inter_p[ii][2];</div><div class="line"><a name="l00777"></a><span class="lineno">  777</span>&#160;        y2cor=inter_p[ii][3];</div><div class="line"><a name="l00778"></a><span class="lineno">  778</span>&#160;        z2cor=thirdcoor;</div><div class="line"><a name="l00779"></a><span class="lineno">  779</span>&#160;        </div><div class="line"><a name="l00780"></a><span class="lineno">  780</span>&#160;    }</div><div class="line"><a name="l00781"></a><span class="lineno">  781</span>&#160;    </div><div class="line"><a name="l00782"></a><span class="lineno">  782</span>&#160;    <span class="keywordflow">if</span> ((zonenumb_in==3)||(zonenumb_in==5))</div><div class="line"><a name="l00783"></a><span class="lineno">  783</span>&#160;    {</div><div class="line"><a name="l00784"></a><span class="lineno">  784</span>&#160;        x1cor=thirdcoor;</div><div class="line"><a name="l00785"></a><span class="lineno">  785</span>&#160;        x2cor=thirdcoor;</div><div class="line"><a name="l00786"></a><span class="lineno">  786</span>&#160;        z1cor=inter_p[ii][0];</div><div class="line"><a name="l00787"></a><span class="lineno">  787</span>&#160;        z2cor=inter_p[ii][2];</div><div class="line"><a name="l00788"></a><span class="lineno">  788</span>&#160;        y1cor=inter_p[ii][1];</div><div class="line"><a name="l00789"></a><span class="lineno">  789</span>&#160;        y2cor=inter_p[ii][3];</div><div class="line"><a name="l00790"></a><span class="lineno">  790</span>&#160;        </div><div class="line"><a name="l00791"></a><span class="lineno">  791</span>&#160;    }</div><div class="line"><a name="l00792"></a><span class="lineno">  792</span>&#160;    </div><div class="line"><a name="l00793"></a><span class="lineno">  793</span>&#160;    <span class="keywordflow">if</span> ((zonenumb_in==4)||(zonenumb_in==6))</div><div class="line"><a name="l00794"></a><span class="lineno">  794</span>&#160;    {</div><div class="line"><a name="l00795"></a><span class="lineno">  795</span>&#160;        y1cor=thirdcoor;</div><div class="line"><a name="l00796"></a><span class="lineno">  796</span>&#160;        y2cor=thirdcoor;</div><div class="line"><a name="l00797"></a><span class="lineno">  797</span>&#160;        x1cor=inter_p[ii][0];</div><div class="line"><a name="l00798"></a><span class="lineno">  798</span>&#160;        x2cor=inter_p[ii][2];</div><div class="line"><a name="l00799"></a><span class="lineno">  799</span>&#160;        z1cor=inter_p[ii][1];</div><div class="line"><a name="l00800"></a><span class="lineno">  800</span>&#160;        z2cor=inter_p[ii][3];</div><div class="line"><a name="l00801"></a><span class="lineno">  801</span>&#160;        </div><div class="line"><a name="l00802"></a><span class="lineno">  802</span>&#160;    }</div><div class="line"><a name="l00803"></a><span class="lineno">  803</span>&#160;    </div><div class="line"><a name="l00804"></a><span class="lineno">  804</span>&#160;    <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].theta!=0.0)</div><div class="line"><a name="l00805"></a><span class="lineno">  805</span>&#160;    {</div><div class="line"><a name="l00806"></a><span class="lineno">  806</span>&#160;        </div><div class="line"><a name="l00807"></a><span class="lineno">  807</span>&#160;        x_1=<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[0][0]*x1cor+<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[0][1]*y1cor+<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[0][2]*z1cor;</div><div class="line"><a name="l00808"></a><span class="lineno">  808</span>&#160;        y_1=<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[1][0]*x1cor+<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[1][1]*y1cor+<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[1][2]*z1cor;</div><div class="line"><a name="l00809"></a><span class="lineno">  809</span>&#160;        z_1=<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[2][0]*x1cor+<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[2][1]*y1cor+<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[2][2]*z1cor;</div><div class="line"><a name="l00810"></a><span class="lineno">  810</span>&#160;        x_2=<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[0][0]*x2cor+<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[0][1]*y2cor+<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[0][2]*z2cor;</div><div class="line"><a name="l00811"></a><span class="lineno">  811</span>&#160;        y_2=<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[1][0]*x2cor+<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[1][1]*y2cor+<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[1][2]*z2cor;</div><div class="line"><a name="l00812"></a><span class="lineno">  812</span>&#160;        z_2=<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[2][0]*x2cor+<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[2][1]*y2cor+<a class="code" href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a>[j].<a class="code" href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">rot2mat</a>[2][2]*z2cor;</div><div class="line"><a name="l00813"></a><span class="lineno">  813</span>&#160;        </div><div class="line"><a name="l00814"></a><span class="lineno">  814</span>&#160;    }</div><div class="line"><a name="l00815"></a><span class="lineno">  815</span>&#160;    <span class="keywordflow">else</span></div><div class="line"><a name="l00816"></a><span class="lineno">  816</span>&#160;    {</div><div class="line"><a name="l00817"></a><span class="lineno">  817</span>&#160;        <span class="comment">/* if angle =0 and fracture is parallel to xy plane, we use the same x and y coordinates */</span></div><div class="line"><a name="l00818"></a><span class="lineno">  818</span>&#160;        x_1=x1cor;</div><div class="line"><a name="l00819"></a><span class="lineno">  819</span>&#160;        y_1=y1cor;</div><div class="line"><a name="l00820"></a><span class="lineno">  820</span>&#160;        z_1=z1cor;</div><div class="line"><a name="l00821"></a><span class="lineno">  821</span>&#160;        x_2=x2cor;</div><div class="line"><a name="l00822"></a><span class="lineno">  822</span>&#160;        y_2=y2cor;</div><div class="line"><a name="l00823"></a><span class="lineno">  823</span>&#160;        z_2=z2cor;</div><div class="line"><a name="l00824"></a><span class="lineno">  824</span>&#160;        </div><div class="line"><a name="l00825"></a><span class="lineno">  825</span>&#160;    }</div><div class="line"><a name="l00826"></a><span class="lineno">  826</span>&#160;    <span class="keywordtype">double</span> deltax, deltay;</div><div class="line"><a name="l00827"></a><span class="lineno">  827</span>&#160;    <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span>  pf;</div><div class="line"><a name="l00828"></a><span class="lineno">  828</span>&#160;    pf=parts_fracture;</div><div class="line"><a name="l00829"></a><span class="lineno">  829</span>&#160;    deltax=(x_2-x_1);</div><div class="line"><a name="l00830"></a><span class="lineno">  830</span>&#160;    deltay=(y_2-y_1);</div><div class="line"><a name="l00831"></a><span class="lineno">  831</span>&#160;    </div><div class="line"><a name="l00832"></a><span class="lineno">  832</span>&#160;    <span class="keywordflow">for</span> (j=0; j&lt;pf; j++)</div><div class="line"><a name="l00833"></a><span class="lineno">  833</span>&#160;    {</div><div class="line"><a name="l00834"></a><span class="lineno">  834</span>&#160;        </div><div class="line"><a name="l00835"></a><span class="lineno">  835</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#afdb17185baff7dc4f22cd6c76c6fd6fd">position</a>[0]=x_1+(deltax/pf)*(j)+deltax/(2.0*pf);</div><div class="line"><a name="l00836"></a><span class="lineno">  836</span>&#160;        </div><div class="line"><a name="l00837"></a><span class="lineno">  837</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#afdb17185baff7dc4f22cd6c76c6fd6fd">position</a>[1]=y_1+(deltay/pf)*(j)+deltay/(2.0*pf);</div><div class="line"><a name="l00838"></a><span class="lineno">  838</span>&#160;        </div><div class="line"><a name="l00839"></a><span class="lineno">  839</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a70e44829116675e21bda7422f0c2b777">velocity</a>[0]=0.;</div><div class="line"><a name="l00840"></a><span class="lineno">  840</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a70e44829116675e21bda7422f0c2b777">velocity</a>[1]=0.;</div><div class="line"><a name="l00841"></a><span class="lineno">  841</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a0b9253ce268d781958145f9d110df2de">fracture</a>=fracture_n;</div><div class="line"><a name="l00842"></a><span class="lineno">  842</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a1747a7b7a63f22656afa1df6462279e8">intcell</a>=0;</div><div class="line"><a name="l00843"></a><span class="lineno">  843</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a11c5e87e6e62d258648d5a02472251d5">time</a>=0.0;</div><div class="line"><a name="l00844"></a><span class="lineno">  844</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#ad50f20b01e37af1fedfd7bbe0cd1d7e3">fl_weight</a>=0.0;</div><div class="line"><a name="l00845"></a><span class="lineno">  845</span>&#160;        k_current++;</div><div class="line"><a name="l00846"></a><span class="lineno">  846</span>&#160;        </div><div class="line"><a name="l00847"></a><span class="lineno">  847</span>&#160;        <span class="keywordflow">if</span> (k_current&gt;<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>)</div><div class="line"><a name="l00848"></a><span class="lineno">  848</span>&#160;        {</div><div class="line"><a name="l00849"></a><span class="lineno">  849</span>&#160;            printf(<span class="stringliteral">&quot; \n Number of particles with allocated memory is less than number of particles set up initially. \n&quot;</span>);</div><div class="line"><a name="l00850"></a><span class="lineno">  850</span>&#160;            printf(<span class="stringliteral">&quot; Increase the number of particles. Program is terminated. \n&quot;</span>);</div><div class="line"><a name="l00851"></a><span class="lineno">  851</span>&#160;            exit(1);</div><div class="line"><a name="l00852"></a><span class="lineno">  852</span>&#160;            </div><div class="line"><a name="l00853"></a><span class="lineno">  853</span>&#160;        }</div><div class="line"><a name="l00854"></a><span class="lineno">  854</span>&#160;    }</div><div class="line"><a name="l00855"></a><span class="lineno">  855</span>&#160;    </div><div class="line"><a name="l00856"></a><span class="lineno">  856</span>&#160;    </div><div class="line"><a name="l00857"></a><span class="lineno">  857</span>&#160;    </div><div class="line"><a name="l00858"></a><span class="lineno">  858</span>&#160;    <span class="keywordflow">return</span> k_current;</div><div class="line"><a name="l00859"></a><span class="lineno">  859</span>&#160;}</div><div class="line"><a name="l00861"></a><span class="lineno"><a class="line" href="_initial_part_positions_8c.html#a59d54f4d5c97e70b0f25e32de87ee3da">  861</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="_initial_part_positions_8c.html#a59d54f4d5c97e70b0f25e32de87ee3da">FlowInWeight</a>(<span class="keywordtype">int</span> numberpart)</div><div class="line"><a name="l00863"></a><span class="lineno">  863</span>&#160;{</div><div class="line"><a name="l00864"></a><span class="lineno">  864</span>&#160;    </div><div class="line"><a name="l00865"></a><span class="lineno">  865</span>&#160;    <span class="keywordtype">int</span> ind1=0,ind2=0,ver1=0, ver2=0, ver3=0,incell=0, n1in=0, n2in=0, jj;</div><div class="line"><a name="l00866"></a><span class="lineno">  866</span>&#160;    <span class="keywordtype">int</span> ins;</div><div class="line"><a name="l00867"></a><span class="lineno">  867</span>&#160;    <span class="keywordtype">double</span> sumflux1=0, sumflux2=0, particleflux[numberpart], totalflux=0;</div><div class="line"><a name="l00868"></a><span class="lineno">  868</span>&#160;    <span class="keywordflow">for</span> (<a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>=0; <a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>&lt;numberpart; <a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>++)</div><div class="line"><a name="l00869"></a><span class="lineno">  869</span>&#160;    {</div><div class="line"><a name="l00870"></a><span class="lineno">  870</span>&#160;        ins=0;</div><div class="line"><a name="l00871"></a><span class="lineno">  871</span>&#160;        ins=<a class="code" href="_initial_part_positions_8c.html#a461479d92dd1f67c7ef94fa953186860">InitCell</a>();</div><div class="line"><a name="l00872"></a><span class="lineno">  872</span>&#160;        incell=<a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[<a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>].<a class="code" href="structcontam.html#a858c05f3221a70f6e212c1e7e172252c">cell</a>;</div><div class="line"><a name="l00873"></a><span class="lineno">  873</span>&#160;        <span class="keywordflow">if</span> (incell!=0)</div><div class="line"><a name="l00874"></a><span class="lineno">  874</span>&#160;        {</div><div class="line"><a name="l00875"></a><span class="lineno">  875</span>&#160;            incell=<a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[<a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>].<a class="code" href="structcontam.html#a858c05f3221a70f6e212c1e7e172252c">cell</a>;</div><div class="line"><a name="l00876"></a><span class="lineno">  876</span>&#160;            ver1=<a class="code" href="_func_def_8h.html#a26a18460ba7bc9d4e9ace81fda1819ee">cell</a>[incell-1].<a class="code" href="structelement.html#aadc95dc4817009223151c40387f40e06">node_ind</a>[0];</div><div class="line"><a name="l00877"></a><span class="lineno">  877</span>&#160;            ver2=<a class="code" href="_func_def_8h.html#a26a18460ba7bc9d4e9ace81fda1819ee">cell</a>[incell-1].<a class="code" href="structelement.html#aadc95dc4817009223151c40387f40e06">node_ind</a>[1];</div><div class="line"><a name="l00878"></a><span class="lineno">  878</span>&#160;            ver3=<a class="code" href="_func_def_8h.html#a26a18460ba7bc9d4e9ace81fda1819ee">cell</a>[incell-1].<a class="code" href="structelement.html#aadc95dc4817009223151c40387f40e06">node_ind</a>[2];</div><div class="line"><a name="l00879"></a><span class="lineno">  879</span>&#160;            n1in=0;</div><div class="line"><a name="l00880"></a><span class="lineno">  880</span>&#160;            n2in=0;</div><div class="line"><a name="l00881"></a><span class="lineno">  881</span>&#160;            <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[ver1-1].typeN&gt;=300)</div><div class="line"><a name="l00882"></a><span class="lineno">  882</span>&#160;            {</div><div class="line"><a name="l00883"></a><span class="lineno">  883</span>&#160;                n1in=ver1;</div><div class="line"><a name="l00884"></a><span class="lineno">  884</span>&#160;                ind1=0;</div><div class="line"><a name="l00885"></a><span class="lineno">  885</span>&#160;            }</div><div class="line"><a name="l00886"></a><span class="lineno">  886</span>&#160;            <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[ver2-1].typeN&gt;=300)</div><div class="line"><a name="l00887"></a><span class="lineno">  887</span>&#160;            {</div><div class="line"><a name="l00888"></a><span class="lineno">  888</span>&#160;                <span class="keywordflow">if</span> (n1in==0)</div><div class="line"><a name="l00889"></a><span class="lineno">  889</span>&#160;                {</div><div class="line"><a name="l00890"></a><span class="lineno">  890</span>&#160;                    n1in=ver2;</div><div class="line"><a name="l00891"></a><span class="lineno">  891</span>&#160;                    ind1=1;</div><div class="line"><a name="l00892"></a><span class="lineno">  892</span>&#160;                }</div><div class="line"><a name="l00893"></a><span class="lineno">  893</span>&#160;                <span class="keywordflow">else</span></div><div class="line"><a name="l00894"></a><span class="lineno">  894</span>&#160;                {</div><div class="line"><a name="l00895"></a><span class="lineno">  895</span>&#160;                    n2in=ver2;</div><div class="line"><a name="l00896"></a><span class="lineno">  896</span>&#160;                    ind2=1;</div><div class="line"><a name="l00897"></a><span class="lineno">  897</span>&#160;                }</div><div class="line"><a name="l00898"></a><span class="lineno">  898</span>&#160;            }</div><div class="line"><a name="l00899"></a><span class="lineno">  899</span>&#160;            <span class="keywordflow">if</span> (<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[ver3-1].typeN&gt;=300)</div><div class="line"><a name="l00900"></a><span class="lineno">  900</span>&#160;            {</div><div class="line"><a name="l00901"></a><span class="lineno">  901</span>&#160;                <span class="keywordflow">if</span> (n1in==0)</div><div class="line"><a name="l00902"></a><span class="lineno">  902</span>&#160;                {</div><div class="line"><a name="l00903"></a><span class="lineno">  903</span>&#160;                    n1in=ver3;</div><div class="line"><a name="l00904"></a><span class="lineno">  904</span>&#160;                    ind1=2;</div><div class="line"><a name="l00905"></a><span class="lineno">  905</span>&#160;                }</div><div class="line"><a name="l00906"></a><span class="lineno">  906</span>&#160;                <span class="keywordflow">else</span></div><div class="line"><a name="l00907"></a><span class="lineno">  907</span>&#160;                {</div><div class="line"><a name="l00908"></a><span class="lineno">  908</span>&#160;                    n2in=ver3;</div><div class="line"><a name="l00909"></a><span class="lineno">  909</span>&#160;                    ind2=2;</div><div class="line"><a name="l00910"></a><span class="lineno">  910</span>&#160;                }</div><div class="line"><a name="l00911"></a><span class="lineno">  911</span>&#160;                </div><div class="line"><a name="l00912"></a><span class="lineno">  912</span>&#160;            }</div><div class="line"><a name="l00913"></a><span class="lineno">  913</span>&#160;       </div><div class="line"><a name="l00914"></a><span class="lineno">  914</span>&#160;            <span class="keywordflow">if</span> ((n1in!=0) &amp;&amp; (n2in!=0))</div><div class="line"><a name="l00915"></a><span class="lineno">  915</span>&#160;            {</div><div class="line"><a name="l00916"></a><span class="lineno">  916</span>&#160;                sumflux1=0;</div><div class="line"><a name="l00917"></a><span class="lineno">  917</span>&#160;                sumflux2=0;</div><div class="line"><a name="l00918"></a><span class="lineno">  918</span>&#160;                <span class="keywordflow">for</span> (jj=0; jj&lt;<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[n1in-1].<a class="code" href="structvertex.html#acf83b561551952403148f73fe032f57b">numneighb</a>; jj++)</div><div class="line"><a name="l00919"></a><span class="lineno">  919</span>&#160;                    </div><div class="line"><a name="l00920"></a><span class="lineno">  920</span>&#160;                    sumflux1=sumflux1+fabs(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[n1in-1].flux[jj]);</div><div class="line"><a name="l00921"></a><span class="lineno">  921</span>&#160;                </div><div class="line"><a name="l00922"></a><span class="lineno">  922</span>&#160;                <span class="keywordflow">for</span> (jj=0; jj&lt;<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[n2in-1].<a class="code" href="structvertex.html#acf83b561551952403148f73fe032f57b">numneighb</a>; jj++)</div><div class="line"><a name="l00923"></a><span class="lineno">  923</span>&#160;                    </div><div class="line"><a name="l00924"></a><span class="lineno">  924</span>&#160;                    sumflux2=sumflux2+fabs(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[n2in-1].flux[jj]);</div><div class="line"><a name="l00925"></a><span class="lineno">  925</span>&#160;                </div><div class="line"><a name="l00926"></a><span class="lineno">  926</span>&#160;                </div><div class="line"><a name="l00927"></a><span class="lineno">  927</span>&#160;                particleflux[<a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>]=<a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[<a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>].<a class="code" href="structcontam.html#a5b317ab55ab5b48dc82c9a473c76a41e">weight</a>[ind1]*sumflux1+<a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[<a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>].<a class="code" href="structcontam.html#a5b317ab55ab5b48dc82c9a473c76a41e">weight</a>[ind2]*sumflux2;</div><div class="line"><a name="l00928"></a><span class="lineno">  928</span>&#160;                totalflux=totalflux+particleflux[<a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>];</div><div class="line"><a name="l00929"></a><span class="lineno">  929</span>&#160;                </div><div class="line"><a name="l00930"></a><span class="lineno">  930</span>&#160;               </div><div class="line"><a name="l00931"></a><span class="lineno">  931</span>&#160;                </div><div class="line"><a name="l00932"></a><span class="lineno">  932</span>&#160;            }</div><div class="line"><a name="l00933"></a><span class="lineno">  933</span>&#160;            <span class="keywordflow">else</span></div><div class="line"><a name="l00934"></a><span class="lineno">  934</span>&#160;            {</div><div class="line"><a name="l00935"></a><span class="lineno">  935</span>&#160;                <span class="keywordtype">int</span> ncent=0;</div><div class="line"><a name="l00936"></a><span class="lineno">  936</span>&#160;                ncent=n1in+n2in;</div><div class="line"><a name="l00937"></a><span class="lineno">  937</span>&#160;                <span class="keywordflow">if</span> (ncent!=0)</div><div class="line"><a name="l00938"></a><span class="lineno">  938</span>&#160;                {</div><div class="line"><a name="l00939"></a><span class="lineno">  939</span>&#160;                    sumflux1=0;</div><div class="line"><a name="l00940"></a><span class="lineno">  940</span>&#160;                    </div><div class="line"><a name="l00941"></a><span class="lineno">  941</span>&#160;                    <span class="keywordflow">for</span> (jj=0; jj&lt;<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[ncent-1].<a class="code" href="structvertex.html#acf83b561551952403148f73fe032f57b">numneighb</a>; jj++)</div><div class="line"><a name="l00942"></a><span class="lineno">  942</span>&#160;                        </div><div class="line"><a name="l00943"></a><span class="lineno">  943</span>&#160;                        sumflux1=sumflux1+fabs(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[ncent-1].flux[jj]);</div><div class="line"><a name="l00944"></a><span class="lineno">  944</span>&#160;                    </div><div class="line"><a name="l00945"></a><span class="lineno">  945</span>&#160;                    </div><div class="line"><a name="l00946"></a><span class="lineno">  946</span>&#160;                    particleflux[<a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>]=sumflux1;</div><div class="line"><a name="l00947"></a><span class="lineno">  947</span>&#160;                    totalflux=totalflux+particleflux[<a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>];</div><div class="line"><a name="l00948"></a><span class="lineno">  948</span>&#160;                    </div><div class="line"><a name="l00949"></a><span class="lineno">  949</span>&#160;                    </div><div class="line"><a name="l00950"></a><span class="lineno">  950</span>&#160;                   </div><div class="line"><a name="l00951"></a><span class="lineno">  951</span>&#160;                }</div><div class="line"><a name="l00952"></a><span class="lineno">  952</span>&#160;            }</div><div class="line"><a name="l00953"></a><span class="lineno">  953</span>&#160;            </div><div class="line"><a name="l00954"></a><span class="lineno">  954</span>&#160;        }</div><div class="line"><a name="l00955"></a><span class="lineno">  955</span>&#160;    }</div><div class="line"><a name="l00956"></a><span class="lineno">  956</span>&#160;    </div><div class="line"><a name="l00957"></a><span class="lineno">  957</span>&#160;    <span class="keywordflow">for</span> (<a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>=0; <a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>&lt;numberpart; <a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>++)</div><div class="line"><a name="l00958"></a><span class="lineno">  958</span>&#160;    {</div><div class="line"><a name="l00959"></a><span class="lineno">  959</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[<a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>].<a class="code" href="structcontam.html#ad50f20b01e37af1fedfd7bbe0cd1d7e3">fl_weight</a>=particleflux[<a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>]/totalflux;</div><div class="line"><a name="l00960"></a><span class="lineno">  960</span>&#160;        <span class="comment">// printf(&quot;%d   %5.12e %5.12e  %5.12e  \n&quot;, np+1,particleflux[np], totalflux, particle[np].fl_weight);</span></div><div class="line"><a name="l00961"></a><span class="lineno">  961</span>&#160;    }</div><div class="line"><a name="l00962"></a><span class="lineno">  962</span>&#160;    </div><div class="line"><a name="l00963"></a><span class="lineno">  963</span>&#160;    </div><div class="line"><a name="l00964"></a><span class="lineno">  964</span>&#160;    </div><div class="line"><a name="l00965"></a><span class="lineno">  965</span>&#160;    <span class="keywordflow">return</span>;</div><div class="line"><a name="l00966"></a><span class="lineno">  966</span>&#160;}</div><div class="line"><a name="l00968"></a><span class="lineno"><a class="line" href="_initial_part_positions_8c.html#a0355cfa365c20964163da103eec2f0e5">  968</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="_initial_part_positions_8c.html#a0355cfa365c20964163da103eec2f0e5">InitInMatrix</a>()</div><div class="line"><a name="l00972"></a><span class="lineno">  972</span>&#160;{</div><div class="line"><a name="l00973"></a><span class="lineno">  973</span>&#160;  </div><div class="line"><a name="l00974"></a><span class="lineno">  974</span>&#160;    <span class="keyword">struct </span><a class="code" href="structinpfile.html">inpfile</a> inputfile;</div><div class="line"><a name="l00975"></a><span class="lineno">  975</span>&#160;    </div><div class="line"><a name="l00976"></a><span class="lineno">  976</span>&#160;    inputfile = <a class="code" href="_func_def_8h.html#a2afbea1a5d0b14d8a098f551fe055420">Control_File</a>(<span class="stringliteral">&quot;inm_coord:&quot;</span>,10 );</div><div class="line"><a name="l00977"></a><span class="lineno">  977</span>&#160;    </div><div class="line"><a name="l00978"></a><span class="lineno">  978</span>&#160;    FILE *mc= <a class="code" href="_func_def_8h.html#a205add90a27292ac5f18188fa35bcdf5">OpenFile</a> (inputfile.filename,<span class="stringliteral">&quot;r&quot;</span>);</div><div class="line"><a name="l00979"></a><span class="lineno">  979</span>&#160;    </div><div class="line"><a name="l00980"></a><span class="lineno">  980</span>&#160;    printf(<span class="stringliteral">&quot;\n OPEN AND READ FILE: %s \n \n&quot;</span>, inputfile.filename);</div><div class="line"><a name="l00981"></a><span class="lineno">  981</span>&#160;    </div><div class="line"><a name="l00982"></a><span class="lineno">  982</span>&#160;    inputfile = <a class="code" href="_func_def_8h.html#a2afbea1a5d0b14d8a098f551fe055420">Control_File</a>(<span class="stringliteral">&quot;inm_nodeID:&quot;</span>,10 );</div><div class="line"><a name="l00983"></a><span class="lineno">  983</span>&#160;    </div><div class="line"><a name="l00984"></a><span class="lineno">  984</span>&#160;    FILE *mn= <a class="code" href="_func_def_8h.html#a205add90a27292ac5f18188fa35bcdf5">OpenFile</a> (inputfile.filename,<span class="stringliteral">&quot;r&quot;</span>);</div><div class="line"><a name="l00985"></a><span class="lineno">  985</span>&#160;    </div><div class="line"><a name="l00986"></a><span class="lineno">  986</span>&#160;    printf(<span class="stringliteral">&quot;\n OPEN AND READ FILE: %s \n \n&quot;</span>, inputfile.filename);</div><div class="line"><a name="l00987"></a><span class="lineno">  987</span>&#160;    </div><div class="line"><a name="l00988"></a><span class="lineno">  988</span>&#160;    <span class="comment">//FILE *mdf=OpenFile(&quot;distance_time.dat&quot;,&quot;w&quot;);</span></div><div class="line"><a name="l00989"></a><span class="lineno">  989</span>&#160;    <span class="keywordtype">int</span> i, ii;</div><div class="line"><a name="l00990"></a><span class="lineno">  990</span>&#160;    <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> number, no;</div><div class="line"><a name="l00991"></a><span class="lineno">  991</span>&#160;    <span class="keywordtype">char</span> cs;</div><div class="line"><a name="l00992"></a><span class="lineno">  992</span>&#160;    <span class="keywordflow">if</span> (fscanf(mn,<span class="stringliteral">&quot;%d  %d %d %d %d \n&quot;</span>, &amp;no, &amp;no, &amp;number, &amp;no, &amp;no)!=5)</div><div class="line"><a name="l00993"></a><span class="lineno">  993</span>&#160;        printf(<span class="stringliteral">&quot;error&quot;</span>);</div><div class="line"><a name="l00994"></a><span class="lineno">  994</span>&#160;    </div><div class="line"><a name="l00995"></a><span class="lineno">  995</span>&#160;    <span class="keywordflow">for</span> (i=0; i&lt;(number+1); i++)</div><div class="line"><a name="l00996"></a><span class="lineno">  996</span>&#160;    {</div><div class="line"><a name="l00997"></a><span class="lineno">  997</span>&#160;        <span class="keywordflow">do</span></div><div class="line"><a name="l00998"></a><span class="lineno">  998</span>&#160;            cs=fgetc(mn);</div><div class="line"><a name="l00999"></a><span class="lineno">  999</span>&#160;        <span class="keywordflow">while</span> (cs!=<span class="charliteral">&#39;\n&#39;</span>);</div><div class="line"><a name="l01000"></a><span class="lineno"> 1000</span>&#160;    }</div><div class="line"><a name="l01001"></a><span class="lineno"> 1001</span>&#160;    </div><div class="line"><a name="l01002"></a><span class="lineno"> 1002</span>&#160;    </div><div class="line"><a name="l01003"></a><span class="lineno"> 1003</span>&#160;    </div><div class="line"><a name="l01004"></a><span class="lineno"> 1004</span>&#160;    <span class="keywordflow">if</span> (fscanf(mc,<span class="stringliteral">&quot; %d \n&quot;</span>, &amp;<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>)!=1)</div><div class="line"><a name="l01005"></a><span class="lineno"> 1005</span>&#160;        printf(<span class="stringliteral">&quot;error&quot;</span>);</div><div class="line"><a name="l01006"></a><span class="lineno"> 1006</span>&#160;    </div><div class="line"><a name="l01007"></a><span class="lineno"> 1007</span>&#160;    <span class="comment">/*  memory allocatation for particle structures */</span></div><div class="line"><a name="l01008"></a><span class="lineno"> 1008</span>&#160;    </div><div class="line"><a name="l01009"></a><span class="lineno"> 1009</span>&#160;    <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>=(<span class="keyword">struct </span><a class="code" href="structcontam.html">contam</a>*) malloc ((<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>+1)*<span class="keyword">sizeof</span>(<span class="keyword">struct </span><a class="code" href="structcontam.html">contam</a>));</div><div class="line"><a name="l01010"></a><span class="lineno"> 1010</span>&#160;    </div><div class="line"><a name="l01011"></a><span class="lineno"> 1011</span>&#160;    <span class="keywordtype">double</span>* distance=malloc ((<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>+1)*<span class="keyword">sizeof</span>(<span class="keywordtype">double</span>));</div><div class="line"><a name="l01012"></a><span class="lineno"> 1012</span>&#160;    <span class="keywordtype">double</span> xp, yp, zp,  sum_distance=0.0, xp2,yp2,zp2;</div><div class="line"><a name="l01013"></a><span class="lineno"> 1013</span>&#160;    </div><div class="line"><a name="l01014"></a><span class="lineno"> 1014</span>&#160;    <span class="keywordflow">for</span> (ii=0; ii&lt;<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>; ii++)</div><div class="line"><a name="l01015"></a><span class="lineno"> 1015</span>&#160;    {</div><div class="line"><a name="l01016"></a><span class="lineno"> 1016</span>&#160;        </div><div class="line"><a name="l01017"></a><span class="lineno"> 1017</span>&#160;        <span class="keywordflow">if</span> (fscanf(mn,<span class="stringliteral">&quot; %d  %d   %d  %d  %d  %d \n&quot;</span>, &amp;no, &amp;no, &amp;no, &amp;no, &amp;no, &amp;number)!=6)</div><div class="line"><a name="l01018"></a><span class="lineno"> 1018</span>&#160;            printf(<span class="stringliteral">&quot;error&quot;</span>);</div><div class="line"><a name="l01019"></a><span class="lineno"> 1019</span>&#160;        </div><div class="line"><a name="l01020"></a><span class="lineno"> 1020</span>&#160;        <span class="keywordflow">if</span> (fscanf(mc,<span class="stringliteral">&quot; %lf %lf %lf \n &quot;</span>, &amp;xp, &amp;yp, &amp;zp)!=1)</div><div class="line"><a name="l01021"></a><span class="lineno"> 1021</span>&#160;            printf (<span class="stringliteral">&quot;error&quot;</span>);</div><div class="line"><a name="l01022"></a><span class="lineno"> 1022</span>&#160;        </div><div class="line"><a name="l01023"></a><span class="lineno"> 1023</span>&#160;        xp2=(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[number-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[0]-xp)*(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[number-1].coord[0]-xp);</div><div class="line"><a name="l01024"></a><span class="lineno"> 1024</span>&#160;        yp2=(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[number-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[1]-yp)*(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[number-1].coord[1]-yp);</div><div class="line"><a name="l01025"></a><span class="lineno"> 1025</span>&#160;        zp2=(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[number-1].<a class="code" href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">coord</a>[2]-zp)*(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[number-1].coord[2]-zp);</div><div class="line"><a name="l01026"></a><span class="lineno"> 1026</span>&#160;        distance[ii]=sqrt(xp2+yp2+zp2);</div><div class="line"><a name="l01027"></a><span class="lineno"> 1027</span>&#160;        </div><div class="line"><a name="l01028"></a><span class="lineno"> 1028</span>&#160;        </div><div class="line"><a name="l01029"></a><span class="lineno"> 1029</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[ii].<a class="code" href="structcontam.html#a70e44829116675e21bda7422f0c2b777">velocity</a>[0]=0.;</div><div class="line"><a name="l01030"></a><span class="lineno"> 1030</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[ii].<a class="code" href="structcontam.html#a70e44829116675e21bda7422f0c2b777">velocity</a>[1]=0.;</div><div class="line"><a name="l01031"></a><span class="lineno"> 1031</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[ii].<a class="code" href="structcontam.html#a0b9253ce268d781958145f9d110df2de">fracture</a>=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[number-1].<a class="code" href="structvertex.html#aa2cd06573422bc415cb8f38a67232de7">fracture</a>[0];</div><div class="line"><a name="l01032"></a><span class="lineno"> 1032</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[ii].<a class="code" href="structcontam.html#a858c05f3221a70f6e212c1e7e172252c">cell</a>=0;</div><div class="line"><a name="l01033"></a><span class="lineno"> 1033</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[ii].<a class="code" href="structcontam.html#afdb17185baff7dc4f22cd6c76c6fd6fd">position</a>[0]=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[number-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[0];</div><div class="line"><a name="l01034"></a><span class="lineno"> 1034</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[ii].<a class="code" href="structcontam.html#afdb17185baff7dc4f22cd6c76c6fd6fd">position</a>[1]=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[number-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[1];</div><div class="line"><a name="l01035"></a><span class="lineno"> 1035</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[ii].<a class="code" href="structcontam.html#a11c5e87e6e62d258648d5a02472251d5">time</a>=0.0;</div><div class="line"><a name="l01036"></a><span class="lineno"> 1036</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[ii].<a class="code" href="structcontam.html#a5b317ab55ab5b48dc82c9a473c76a41e">weight</a>[0]=0.0;</div><div class="line"><a name="l01037"></a><span class="lineno"> 1037</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[ii].<a class="code" href="structcontam.html#a5b317ab55ab5b48dc82c9a473c76a41e">weight</a>[1]=0.0;</div><div class="line"><a name="l01038"></a><span class="lineno"> 1038</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[ii].<a class="code" href="structcontam.html#a5b317ab55ab5b48dc82c9a473c76a41e">weight</a>[2]=0.0;</div><div class="line"><a name="l01039"></a><span class="lineno"> 1039</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[ii].<a class="code" href="structcontam.html#aacd8d2e19b4a3c689e96cf9693218777">prev_pos</a>[0]=0.0;</div><div class="line"><a name="l01040"></a><span class="lineno"> 1040</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[ii].<a class="code" href="structcontam.html#aacd8d2e19b4a3c689e96cf9693218777">prev_pos</a>[1]=0.0;</div><div class="line"><a name="l01041"></a><span class="lineno"> 1041</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[ii].<a class="code" href="structcontam.html#a1747a7b7a63f22656afa1df6462279e8">intcell</a>=0;</div><div class="line"><a name="l01042"></a><span class="lineno"> 1042</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[ii].<a class="code" href="structcontam.html#ad50f20b01e37af1fedfd7bbe0cd1d7e3">fl_weight</a>=0.0;</div><div class="line"><a name="l01043"></a><span class="lineno"> 1043</span>&#160;        <span class="comment">// define a time that took for particle to reach fracture from a rock matrix</span></div><div class="line"><a name="l01044"></a><span class="lineno"> 1044</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[ii].<a class="code" href="structcontam.html#a11c5e87e6e62d258648d5a02472251d5">time</a>=<a class="code" href="_initial_part_positions_8c.html#a0610d22ad530758b4de42acd23a5df92">TimeFromMatrix</a>(distance[ii]);</div><div class="line"><a name="l01045"></a><span class="lineno"> 1045</span>&#160;        sum_distance=sum_distance+distance[ii];</div><div class="line"><a name="l01046"></a><span class="lineno"> 1046</span>&#160;      </div><div class="line"><a name="l01047"></a><span class="lineno"> 1047</span>&#160;        </div><div class="line"><a name="l01048"></a><span class="lineno"> 1048</span>&#160;    }</div><div class="line"><a name="l01049"></a><span class="lineno"> 1049</span>&#160;    fclose(mc);</div><div class="line"><a name="l01050"></a><span class="lineno"> 1050</span>&#160;    fclose(mn);</div><div class="line"><a name="l01051"></a><span class="lineno"> 1051</span>&#160;    </div><div class="line"><a name="l01052"></a><span class="lineno"> 1052</span>&#160;    <span class="comment">//define weights according to distance from fracture</span></div><div class="line"><a name="l01053"></a><span class="lineno"> 1053</span>&#160;    <span class="keywordflow">for</span> (i=0; i&lt;<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>; i++)</div><div class="line"><a name="l01054"></a><span class="lineno"> 1054</span>&#160;    {</div><div class="line"><a name="l01055"></a><span class="lineno"> 1055</span>&#160;        <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[i].<a class="code" href="structcontam.html#ad50f20b01e37af1fedfd7bbe0cd1d7e3">fl_weight</a>=distance[i]/sum_distance;</div><div class="line"><a name="l01056"></a><span class="lineno"> 1056</span>&#160;    }</div><div class="line"><a name="l01057"></a><span class="lineno"> 1057</span>&#160;    </div><div class="line"><a name="l01058"></a><span class="lineno"> 1058</span>&#160;    free(distance);</div><div class="line"><a name="l01059"></a><span class="lineno"> 1059</span>&#160;    <span class="keywordflow">return</span>;</div><div class="line"><a name="l01060"></a><span class="lineno"> 1060</span>&#160;}</div><div class="line"><a name="l01062"></a><span class="lineno"><a class="line" href="_initial_part_positions_8c.html#a0610d22ad530758b4de42acd23a5df92"> 1062</a></span>&#160;<span class="keywordtype">double</span> <a class="code" href="_initial_part_positions_8c.html#a0610d22ad530758b4de42acd23a5df92">TimeFromMatrix</a>(<span class="keywordtype">double</span> pdist)</div><div class="line"><a name="l01064"></a><span class="lineno"> 1064</span>&#160;{</div><div class="line"><a name="l01065"></a><span class="lineno"> 1065</span>&#160;    </div><div class="line"><a name="l01066"></a><span class="lineno"> 1066</span>&#160;    <span class="keyword">struct </span><a class="code" href="structinpfile.html">inpfile</a> inputfile;</div><div class="line"><a name="l01067"></a><span class="lineno"> 1067</span>&#160;    <span class="keywordtype">double</span> ptime=0.0, ptime1=0.0, ptime2=0.0;</div><div class="line"><a name="l01068"></a><span class="lineno"> 1068</span>&#160;    <span class="keywordtype">double</span> randomnumber=0.0;</div><div class="line"><a name="l01069"></a><span class="lineno"> 1069</span>&#160;    randomnumber=drand48();</div><div class="line"><a name="l01070"></a><span class="lineno"> 1070</span>&#160;   </div><div class="line"><a name="l01071"></a><span class="lineno"> 1071</span>&#160;    <span class="keywordtype">double</span> mporosity=0.0;</div><div class="line"><a name="l01072"></a><span class="lineno"> 1072</span>&#160;    <span class="keywordtype">double</span> mdiffcoeff=0.0;</div><div class="line"><a name="l01073"></a><span class="lineno"> 1073</span>&#160;    </div><div class="line"><a name="l01074"></a><span class="lineno"> 1074</span>&#160;    inputfile=<a class="code" href="_func_def_8h.html#a922a07b425353141b29d8d57e830a7c1">Control_Param</a>(<span class="stringliteral">&quot;inm_porosity:&quot;</span>,13);</div><div class="line"><a name="l01075"></a><span class="lineno"> 1075</span>&#160;    mporosity=inputfile.<a class="code" href="structinpfile.html#a57ff2cce78c35a8a128927da9eeccfbc">param</a>;</div><div class="line"><a name="l01076"></a><span class="lineno"> 1076</span>&#160;    inputfile=<a class="code" href="_func_def_8h.html#a922a07b425353141b29d8d57e830a7c1">Control_Param</a>(<span class="stringliteral">&quot;inm_diffcoeff:&quot;</span>,14);</div><div class="line"><a name="l01077"></a><span class="lineno"> 1077</span>&#160;    mdiffcoeff=inputfile.<a class="code" href="structinpfile.html#a57ff2cce78c35a8a128927da9eeccfbc">param</a>;</div><div class="line"><a name="l01078"></a><span class="lineno"> 1078</span>&#160;    <span class="comment">//  printf(&quot;diff coeff %5.9e porosity %lf \n&quot;,mdiffcoeff, mporosity);</span></div><div class="line"><a name="l01079"></a><span class="lineno"> 1079</span>&#160;    </div><div class="line"><a name="l01080"></a><span class="lineno"> 1080</span>&#160;    <span class="keywordtype">double</span> retardation_factor=1.0;</div><div class="line"><a name="l01081"></a><span class="lineno"> 1081</span>&#160;    <span class="keywordtype">double</span> inverse_erfc=0.0;</div><div class="line"><a name="l01082"></a><span class="lineno"> 1082</span>&#160;    <span class="keywordtype">double</span> z;</div><div class="line"><a name="l01083"></a><span class="lineno"> 1083</span>&#160;    z=1-randomnumber;</div><div class="line"><a name="l01084"></a><span class="lineno"> 1084</span>&#160;    inverse_erfc=0.5*sqrt(<a class="code" href="_func_def_8h.html#a1daf785e3f68d293c7caa1c756d5cb74">pi</a>)*(z+(<a class="code" href="_func_def_8h.html#a1daf785e3f68d293c7caa1c756d5cb74">pi</a>/12)*pow(z,3)+((7*pow(<a class="code" href="_func_def_8h.html#a1daf785e3f68d293c7caa1c756d5cb74">pi</a>,2))/480)*pow(z,5)+((127*pow(<a class="code" href="_func_def_8h.html#a1daf785e3f68d293c7caa1c756d5cb74">pi</a>,3))/40320)*pow(z,7)+((4369*pow(<a class="code" href="_func_def_8h.html#a1daf785e3f68d293c7caa1c756d5cb74">pi</a>,4))/5806080)*pow(z,9)+((34807*pow(<a class="code" href="_func_def_8h.html#a1daf785e3f68d293c7caa1c756d5cb74">pi</a>,5))/182476800)*pow(z,11));</div><div class="line"><a name="l01085"></a><span class="lineno"> 1085</span>&#160;    ptime2=(1.0/inverse_erfc)*(1.0/inverse_erfc);</div><div class="line"><a name="l01086"></a><span class="lineno"> 1086</span>&#160;    ptime1=(pdist*pdist)/(4.0*(mporosity*mdiffcoeff/retardation_factor));</div><div class="line"><a name="l01087"></a><span class="lineno"> 1087</span>&#160;    ptime=(ptime1*ptime2)/<a class="code" href="_func_def_8h.html#a5799b26fb87c6d98cd44352eb84a2e7e">timeunit</a>;</div><div class="line"><a name="l01088"></a><span class="lineno"> 1088</span>&#160;    </div><div class="line"><a name="l01089"></a><span class="lineno"> 1089</span>&#160;    </div><div class="line"><a name="l01090"></a><span class="lineno"> 1090</span>&#160;    </div><div class="line"><a name="l01091"></a><span class="lineno"> 1091</span>&#160;    <span class="keywordflow">return</span> ptime;</div><div class="line"><a name="l01092"></a><span class="lineno"> 1092</span>&#160;}</div><div class="line"><a name="l01094"></a><span class="lineno"> 1094</span>&#160;</div><div class="line"><a name="l01095"></a><span class="lineno"><a class="line" href="_initial_part_positions_8c.html#ad3b7970929b5c9502172afc4e0bc65e8"> 1095</a></span>&#160;<span class="keywordtype">int</span> <a class="code" href="_initial_part_positions_8c.html#ad3b7970929b5c9502172afc4e0bc65e8">InitParticles_flux</a> (<span class="keywordtype">int</span> k_current, <span class="keywordtype">int</span> first_ind, <span class="keywordtype">int</span> last_ind, <span class="keywordtype">double</span> weight_p)</div><div class="line"><a name="l01098"></a><span class="lineno"> 1098</span>&#160;{</div><div class="line"><a name="l01099"></a><span class="lineno"> 1099</span>&#160;    </div><div class="line"><a name="l01100"></a><span class="lineno"> 1100</span>&#160;    <span class="keywordtype">double</span> deltax, deltay, edgelength, eqdist_x, eqdist_y, sumflux=0, startpos_x, startpos_y, sumdelta=0.0;</div><div class="line"><a name="l01101"></a><span class="lineno"> 1101</span>&#160;    <span class="keywordtype">double</span> sumfluxfract=0.0;</div><div class="line"><a name="l01102"></a><span class="lineno"> 1102</span>&#160;    <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i, k, j, jj, pf, curcel=0, insc=0, <a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>, np_init;</div><div class="line"><a name="l01103"></a><span class="lineno"> 1103</span>&#160;    </div><div class="line"><a name="l01104"></a><span class="lineno"> 1104</span>&#160;    </div><div class="line"><a name="l01105"></a><span class="lineno"> 1105</span>&#160;    <span class="keyword">struct </span><a class="code" href="structposit3d.html">posit3d</a> particle3dposition;</div><div class="line"><a name="l01106"></a><span class="lineno"> 1106</span>&#160;    </div><div class="line"><a name="l01107"></a><span class="lineno"> 1107</span>&#160;    np_init=k_current;</div><div class="line"><a name="l01108"></a><span class="lineno"> 1108</span>&#160;    </div><div class="line"><a name="l01109"></a><span class="lineno"> 1109</span>&#160;    <span class="comment">//  printf(&quot;firstn %d lastn %d \n&quot;, first_ind, last_ind);</span></div><div class="line"><a name="l01110"></a><span class="lineno"> 1110</span>&#160;    <span class="keywordflow">for</span> (j=first_ind; j&lt;=last_ind; j++)</div><div class="line"><a name="l01111"></a><span class="lineno"> 1111</span>&#160;    {</div><div class="line"><a name="l01112"></a><span class="lineno"> 1112</span>&#160;        sumflux=0.0;</div><div class="line"><a name="l01113"></a><span class="lineno"> 1113</span>&#160;        </div><div class="line"><a name="l01114"></a><span class="lineno"> 1114</span>&#160;        <span class="comment">// calculate the flux of the current cell</span></div><div class="line"><a name="l01115"></a><span class="lineno"> 1115</span>&#160;        <span class="keywordflow">for</span> (jj=0; jj&lt;<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j]-1].<a class="code" href="structvertex.html#acf83b561551952403148f73fe032f57b">numneighb</a>; jj++)</div><div class="line"><a name="l01116"></a><span class="lineno"> 1116</span>&#160;            sumflux=sumflux+(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j]-1].flux[jj]);</div><div class="line"><a name="l01117"></a><span class="lineno"> 1117</span>&#160;</div><div class="line"><a name="l01118"></a><span class="lineno"> 1118</span>&#160;        sumflux=sumflux/<a class="code" href="_func_def_8h.html#a6f8c052f8417728038991f7f2826d38d">density</a>;         </div><div class="line"><a name="l01119"></a><span class="lineno"> 1119</span>&#160;        sumfluxfract=sumfluxfract+sumflux;</div><div class="line"><a name="l01120"></a><span class="lineno"> 1120</span>&#160;        </div><div class="line"><a name="l01121"></a><span class="lineno"> 1121</span>&#160;        <span class="comment">// ceiling number of particles in cell nodezonein[j]</span></div><div class="line"><a name="l01122"></a><span class="lineno"> 1122</span>&#160;        pf=ceilf(fabs(sumflux)/weight_p);</div><div class="line"><a name="l01123"></a><span class="lineno"> 1123</span>&#160;      </div><div class="line"><a name="l01124"></a><span class="lineno"> 1124</span>&#160;   <span class="comment">//     printf(&quot;%5.12e  %d  %5.12e  %5.12e\n&quot;, sumflux, pf, fabs(sumflux)/weight_p, density);    </span></div><div class="line"><a name="l01125"></a><span class="lineno"> 1125</span>&#160;   </div><div class="line"><a name="l01126"></a><span class="lineno"> 1126</span>&#160;        <span class="keywordflow">if</span> ((j&gt;first_ind) &amp;&amp; (j&lt;last_ind))</div><div class="line"><a name="l01127"></a><span class="lineno"> 1127</span>&#160;        {</div><div class="line"><a name="l01128"></a><span class="lineno"> 1128</span>&#160;            <span class="comment">//calculate distance between nodes</span></div><div class="line"><a name="l01129"></a><span class="lineno"> 1129</span>&#160;            deltax=(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j+1]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[0]-<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j-1]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[0])/2.0;</div><div class="line"><a name="l01130"></a><span class="lineno"> 1130</span>&#160;            deltay=(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j+1]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[1]-<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j-1]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[1])/2.0;</div><div class="line"><a name="l01131"></a><span class="lineno"> 1131</span>&#160;            </div><div class="line"><a name="l01132"></a><span class="lineno"> 1132</span>&#160;            </div><div class="line"><a name="l01133"></a><span class="lineno"> 1133</span>&#160;            <span class="comment">//define starting position of first particle</span></div><div class="line"><a name="l01134"></a><span class="lineno"> 1134</span>&#160;            startpos_x=((<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[0])-(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j-1]-1].coord_xy[0]))/2.0+<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j-1]-1].coord_xy[0];</div><div class="line"><a name="l01135"></a><span class="lineno"> 1135</span>&#160;            </div><div class="line"><a name="l01136"></a><span class="lineno"> 1136</span>&#160;            startpos_y=((<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[1])-(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j-1]-1].coord_xy[1]))/2.0+<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j-1]-1].coord_xy[1];</div><div class="line"><a name="l01137"></a><span class="lineno"> 1137</span>&#160;        }</div><div class="line"><a name="l01138"></a><span class="lineno"> 1138</span>&#160;        </div><div class="line"><a name="l01139"></a><span class="lineno"> 1139</span>&#160;        <span class="keywordflow">if</span> (j==first_ind)</div><div class="line"><a name="l01140"></a><span class="lineno"> 1140</span>&#160;        {</div><div class="line"><a name="l01141"></a><span class="lineno"> 1141</span>&#160;            <span class="comment">//calculate distance between nodes</span></div><div class="line"><a name="l01142"></a><span class="lineno"> 1142</span>&#160;            deltax=(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j+1]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[0]-<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[0])/2.0;</div><div class="line"><a name="l01143"></a><span class="lineno"> 1143</span>&#160;            deltay=(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j+1]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[1]-<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[1])/2.0;</div><div class="line"><a name="l01144"></a><span class="lineno"> 1144</span>&#160;            </div><div class="line"><a name="l01145"></a><span class="lineno"> 1145</span>&#160;            </div><div class="line"><a name="l01146"></a><span class="lineno"> 1146</span>&#160;            </div><div class="line"><a name="l01147"></a><span class="lineno"> 1147</span>&#160;            <span class="comment">//define starting position of first particle</span></div><div class="line"><a name="l01148"></a><span class="lineno"> 1148</span>&#160;            startpos_x=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[0];</div><div class="line"><a name="l01149"></a><span class="lineno"> 1149</span>&#160;            startpos_y=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[1];</div><div class="line"><a name="l01150"></a><span class="lineno"> 1150</span>&#160;            </div><div class="line"><a name="l01151"></a><span class="lineno"> 1151</span>&#160;        }</div><div class="line"><a name="l01152"></a><span class="lineno"> 1152</span>&#160;        </div><div class="line"><a name="l01153"></a><span class="lineno"> 1153</span>&#160;        <span class="keywordflow">if</span> (j==last_ind)</div><div class="line"><a name="l01154"></a><span class="lineno"> 1154</span>&#160;        {</div><div class="line"><a name="l01155"></a><span class="lineno"> 1155</span>&#160;            <span class="comment">//calculate distance between nodes</span></div><div class="line"><a name="l01156"></a><span class="lineno"> 1156</span>&#160;            deltax=(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[0]-<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j-1]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[0])/2.0;</div><div class="line"><a name="l01157"></a><span class="lineno"> 1157</span>&#160;            deltay=(<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[1]-<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j-1]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[1])/2.0;</div><div class="line"><a name="l01158"></a><span class="lineno"> 1158</span>&#160;            </div><div class="line"><a name="l01159"></a><span class="lineno"> 1159</span>&#160;            </div><div class="line"><a name="l01160"></a><span class="lineno"> 1160</span>&#160;            </div><div class="line"><a name="l01161"></a><span class="lineno"> 1161</span>&#160;            <span class="comment">//define starting position of first particle</span></div><div class="line"><a name="l01162"></a><span class="lineno"> 1162</span>&#160;            </div><div class="line"><a name="l01163"></a><span class="lineno"> 1163</span>&#160;            startpos_x=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j-1]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[0];</div><div class="line"><a name="l01164"></a><span class="lineno"> 1164</span>&#160;            startpos_y=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[j-1]-1].<a class="code" href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">coord_xy</a>[1];</div><div class="line"><a name="l01165"></a><span class="lineno"> 1165</span>&#160;            </div><div class="line"><a name="l01166"></a><span class="lineno"> 1166</span>&#160;        }</div><div class="line"><a name="l01167"></a><span class="lineno"> 1167</span>&#160;        </div><div class="line"><a name="l01168"></a><span class="lineno"> 1168</span>&#160;        </div><div class="line"><a name="l01169"></a><span class="lineno"> 1169</span>&#160;        <span class="keywordflow">if</span> (pf&lt;2)</div><div class="line"><a name="l01170"></a><span class="lineno"> 1170</span>&#160;        {</div><div class="line"><a name="l01171"></a><span class="lineno"> 1171</span>&#160;            pf=1;</div><div class="line"><a name="l01172"></a><span class="lineno"> 1172</span>&#160;            eqdist_x=deltax/2.0;</div><div class="line"><a name="l01173"></a><span class="lineno"> 1173</span>&#160;            eqdist_y=deltay/2.0;</div><div class="line"><a name="l01174"></a><span class="lineno"> 1174</span>&#160;        }</div><div class="line"><a name="l01175"></a><span class="lineno"> 1175</span>&#160;        <span class="keywordflow">else</span></div><div class="line"><a name="l01176"></a><span class="lineno"> 1176</span>&#160;        {</div><div class="line"><a name="l01177"></a><span class="lineno"> 1177</span>&#160;            eqdist_x=deltax/pf;</div><div class="line"><a name="l01178"></a><span class="lineno"> 1178</span>&#160;            eqdist_y=deltay/pf;</div><div class="line"><a name="l01179"></a><span class="lineno"> 1179</span>&#160;        }</div><div class="line"><a name="l01180"></a><span class="lineno"> 1180</span>&#160;        </div><div class="line"><a name="l01181"></a><span class="lineno"> 1181</span>&#160;        </div><div class="line"><a name="l01182"></a><span class="lineno"> 1182</span>&#160;        <span class="keywordflow">for</span> (jj=0; jj&lt;pf; jj++)</div><div class="line"><a name="l01183"></a><span class="lineno"> 1183</span>&#160;        {</div><div class="line"><a name="l01184"></a><span class="lineno"> 1184</span>&#160;            </div><div class="line"><a name="l01185"></a><span class="lineno"> 1185</span>&#160;            </div><div class="line"><a name="l01186"></a><span class="lineno"> 1186</span>&#160;            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#afdb17185baff7dc4f22cd6c76c6fd6fd">position</a>[0]=startpos_x+eqdist_x*(jj);</div><div class="line"><a name="l01187"></a><span class="lineno"> 1187</span>&#160;            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#afdb17185baff7dc4f22cd6c76c6fd6fd">position</a>[1]=startpos_y+eqdist_y*(jj);</div><div class="line"><a name="l01188"></a><span class="lineno"> 1188</span>&#160;            </div><div class="line"><a name="l01189"></a><span class="lineno"> 1189</span>&#160;            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a70e44829116675e21bda7422f0c2b777">velocity</a>[0]=0.;</div><div class="line"><a name="l01190"></a><span class="lineno"> 1190</span>&#160;            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a70e44829116675e21bda7422f0c2b777">velocity</a>[1]=0.;</div><div class="line"><a name="l01191"></a><span class="lineno"> 1191</span>&#160;            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a0b9253ce268d781958145f9d110df2de">fracture</a>=<a class="code" href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a>[<a class="code" href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a>[first_ind]-1].<a class="code" href="structvertex.html#aa2cd06573422bc415cb8f38a67232de7">fracture</a>[0];</div><div class="line"><a name="l01192"></a><span class="lineno"> 1192</span>&#160;            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a1747a7b7a63f22656afa1df6462279e8">intcell</a>=0;</div><div class="line"><a name="l01193"></a><span class="lineno"> 1193</span>&#160;            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a11c5e87e6e62d258648d5a02472251d5">time</a>=0.0;</div><div class="line"><a name="l01194"></a><span class="lineno"> 1194</span>&#160;            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#ad50f20b01e37af1fedfd7bbe0cd1d7e3">fl_weight</a>=weight_p;</div><div class="line"><a name="l01195"></a><span class="lineno"> 1195</span>&#160;            <a class="code" href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a>[k_current].<a class="code" href="structcontam.html#a858c05f3221a70f6e212c1e7e172252c">cell</a>=0;</div><div class="line"><a name="l01196"></a><span class="lineno"> 1196</span>&#160;            </div><div class="line"><a name="l01197"></a><span class="lineno"> 1197</span>&#160;            </div><div class="line"><a name="l01198"></a><span class="lineno"> 1198</span>&#160;            <a class="code" href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a>=k_current;</div><div class="line"><a name="l01199"></a><span class="lineno"> 1199</span>&#160;            </div><div class="line"><a name="l01200"></a><span class="lineno"> 1200</span>&#160;            </div><div class="line"><a name="l01201"></a><span class="lineno"> 1201</span>&#160;            k_current++;</div><div class="line"><a name="l01202"></a><span class="lineno"> 1202</span>&#160;            </div><div class="line"><a name="l01203"></a><span class="lineno"> 1203</span>&#160;            <span class="keywordflow">if</span> (k_current&gt;<a class="code" href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a>*3)</div><div class="line"><a name="l01204"></a><span class="lineno"> 1204</span>&#160;            {</div><div class="line"><a name="l01205"></a><span class="lineno"> 1205</span>&#160;<span class="comment">//                   printf(&quot;\n Number of particles with allocated memory (%d) is less than number of particles set up initially (%d).  \n&quot;, k_current, npart);</span></div><div class="line"><a name="l01206"></a><span class="lineno"> 1206</span>&#160;                printf(<span class="stringliteral">&quot; Increase the number of particles. Program is terminated. \n&quot;</span>);</div><div class="line"><a name="l01207"></a><span class="lineno"> 1207</span>&#160;                exit(1);</div><div class="line"><a name="l01208"></a><span class="lineno"> 1208</span>&#160;              </div><div class="line"><a name="l01209"></a><span class="lineno"> 1209</span>&#160;            }</div><div class="line"><a name="l01210"></a><span class="lineno"> 1210</span>&#160;            </div><div class="line"><a name="l01211"></a><span class="lineno"> 1211</span>&#160;        }</div><div class="line"><a name="l01212"></a><span class="lineno"> 1212</span>&#160;        </div><div class="line"><a name="l01213"></a><span class="lineno"> 1213</span>&#160;    }</div><div class="line"><a name="l01214"></a><span class="lineno"> 1214</span>&#160;    </div><div class="line"><a name="l01215"></a><span class="lineno"> 1215</span>&#160;   </div><div class="line"><a name="l01216"></a><span class="lineno"> 1216</span>&#160;    <span class="keywordflow">return</span> k_current;</div><div class="line"><a name="l01217"></a><span class="lineno"> 1217</span>&#160;}</div><div class="line"><a name="l01219"></a><span class="lineno"> 1219</span>&#160;</div><div class="ttc" id="structcontam_html_afdb17185baff7dc4f22cd6c76c6fd6fd"><div class="ttname"><a href="structcontam.html#afdb17185baff7dc4f22cd6c76c6fd6fd">contam::position</a></div><div class="ttdeci">double position[2]</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00183">FuncDef.h:183</a></div></div>
<div class="ttc" id="_func_def_8h_html_a6f8c052f8417728038991f7f2826d38d"><div class="ttname"><a href="_func_def_8h.html#a6f8c052f8417728038991f7f2826d38d">density</a></div><div class="ttdeci">double density</div><div class="ttdef"><b>Definition:</b> <a href="main_8c_source.html#l00047">main.c:47</a></div></div>
<div class="ttc" id="structvertex_html_a7588fad77cf759d8db0957c9e957e23a"><div class="ttname"><a href="structvertex.html#a7588fad77cf759d8db0957c9e957e23a">vertex::aperture</a></div><div class="ttdeci">double aperture</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00157">FuncDef.h:157</a></div></div>
<div class="ttc" id="structposit3d_html_aefc0a716be214a76cae4e2176134c28e"><div class="ttname"><a href="structposit3d.html#aefc0a716be214a76cae4e2176134c28e">posit3d::cord3</a></div><div class="ttdeci">double cord3[3]</div><div class="ttdef"><b>Definition:</b> <a href="_initial_part_positions_8c_source.html#l00017">InitialPartPositions.c:17</a></div></div>
<div class="ttc" id="_func_def_8h_html_a922a07b425353141b29d8d57e830a7c1"><div class="ttname"><a href="_func_def_8h.html#a922a07b425353141b29d8d57e830a7c1">Control_Param</a></div><div class="ttdeci">struct inpfile Control_Param(char fileobject[], int ctr)</div><div class="ttdef"><b>Definition:</b> <a href="_read_grid_init_8c_source.html#l01245">ReadGridInit.c:1245</a></div></div>
<div class="ttc" id="_func_def_8h_html_a1daf785e3f68d293c7caa1c756d5cb74"><div class="ttname"><a href="_func_def_8h.html#a1daf785e3f68d293c7caa1c756d5cb74">pi</a></div><div class="ttdeci">#define pi</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00002">FuncDef.h:2</a></div></div>
<div class="ttc" id="structvertex_html_a1eff0deab63f96c97ba5098608d65ad3"><div class="ttname"><a href="structvertex.html#a1eff0deab63f96c97ba5098608d65ad3">vertex::coord</a></div><div class="ttdeci">double coord[3]</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00115">FuncDef.h:115</a></div></div>
<div class="ttc" id="_func_def_8h_html_a5c8fe8978ce1689f5dcf12bc7b56793d"><div class="ttname"><a href="_func_def_8h.html#a5c8fe8978ce1689f5dcf12bc7b56793d">InsideCell</a></div><div class="ttdeci">int InsideCell(unsigned int numc)</div><div class="ttdef"><b>Definition:</b> <a href="_tracking_part_8c_source.html#l01487">TrackingPart.c:1487</a></div></div>
<div class="ttc" id="_initial_part_positions_8c_html_a01e7bd2a5bd9417f1d28849b8577b0b9"><div class="ttname"><a href="_initial_part_positions_8c.html#a01e7bd2a5bd9417f1d28849b8577b0b9">InitParticles_eq</a></div><div class="ttdeci">int InitParticles_eq(int k_current, int firstn, int lastn, double parts_dist, int first_ind, int last_ind)</div><div class="ttdef"><b>Definition:</b> <a href="_initial_part_positions_8c_source.html#l00711">InitialPartPositions.c:711</a></div></div>
<div class="ttc" id="structvertex_html_aa2cd06573422bc415cb8f38a67232de7"><div class="ttname"><a href="structvertex.html#aa2cd06573422bc415cb8f38a67232de7">vertex::fracture</a></div><div class="ttdeci">unsigned int fracture[2]</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00112">FuncDef.h:112</a></div></div>
<div class="ttc" id="structinpfile_html_a66e8759aa38ee20350f04cdc58dee702"><div class="ttname"><a href="structinpfile.html#a66e8759aa38ee20350f04cdc58dee702">inpfile::filename</a></div><div class="ttdeci">char filename[120]</div><div class="ttdef"><b>Definition:</b> <a href="_initial_part_positions_8c_source.html#l00011">InitialPartPositions.c:11</a></div></div>
<div class="ttc" id="_initial_part_positions_8c_html_a461479d92dd1f67c7ef94fa953186860"><div class="ttname"><a href="_initial_part_positions_8c.html#a461479d92dd1f67c7ef94fa953186860">InitCell</a></div><div class="ttdeci">int InitCell()</div><div class="ttdef"><b>Definition:</b> <a href="_initial_part_positions_8c_source.html#l00621">InitialPartPositions.c:621</a></div></div>
<div class="ttc" id="_func_def_8h_html_a26a18460ba7bc9d4e9ace81fda1819ee"><div class="ttname"><a href="_func_def_8h.html#a26a18460ba7bc9d4e9ace81fda1819ee">cell</a></div><div class="ttdeci">struct element * cell</div><div class="ttdef"><b>Definition:</b> <a href="main_8c_source.html#l00042">main.c:42</a></div></div>
<div class="ttc" id="_initial_part_positions_8c_html_a59d54f4d5c97e70b0f25e32de87ee3da"><div class="ttname"><a href="_initial_part_positions_8c.html#a59d54f4d5c97e70b0f25e32de87ee3da">FlowInWeight</a></div><div class="ttdeci">void FlowInWeight(int numberpart)</div><div class="ttdef"><b>Definition:</b> <a href="_initial_part_positions_8c_source.html#l00861">InitialPartPositions.c:861</a></div></div>
<div class="ttc" id="_initial_part_positions_8c_html_ae2df2bd90c8d4aa96e508002a36327c9"><div class="ttname"><a href="_initial_part_positions_8c.html#ae2df2bd90c8d4aa96e508002a36327c9">InitPos</a></div><div class="ttdeci">int InitPos()</div><div class="ttdef"><b>Definition:</b> <a href="_initial_part_positions_8c_source.html#l00021">InitialPartPositions.c:21</a></div></div>
<div class="ttc" id="_func_def_8h_html_a9bb4185d1b46d02dc42d5889e4d70ebf"><div class="ttname"><a href="_func_def_8h.html#a9bb4185d1b46d02dc42d5889e4d70ebf">Control_Data</a></div><div class="ttdeci">struct inpfile Control_Data(char fileobject[], int ctr)</div><div class="ttdef"><b>Definition:</b> <a href="_read_grid_init_8c_source.html#l01201">ReadGridInit.c:1201</a></div></div>
<div class="ttc" id="_func_def_8h_html_a5225ceefbfdde0f9cb32cb84174a4863"><div class="ttname"><a href="_func_def_8h.html#a5225ceefbfdde0f9cb32cb84174a4863">ncells</a></div><div class="ttdeci">unsigned int ncells</div><div class="ttdef"><b>Definition:</b> <a href="main_8c_source.html#l00030">main.c:30</a></div></div>
<div class="ttc" id="_func_def_8h_html_a69f9dfd9f5eba61f40439fc5dc92f3b6"><div class="ttname"><a href="_func_def_8h.html#a69f9dfd9f5eba61f40439fc5dc92f3b6">Moving2Center</a></div><div class="ttdeci">void Moving2Center(int nnp, int cellnumber)</div><div class="ttdef"><b>Definition:</b> <a href="_tracking_part_8c_source.html#l02440">TrackingPart.c:2440</a></div></div>
<div class="ttc" id="_func_def_8h_html_adef0bd5e77533cdd4cdd0bcf0d5dc96b"><div class="ttname"><a href="_func_def_8h.html#adef0bd5e77533cdd4cdd0bcf0d5dc96b">flag_w</a></div><div class="ttdeci">unsigned int flag_w</div><div class="ttdef"><b>Definition:</b> <a href="main_8c_source.html#l00038">main.c:38</a></div></div>
<div class="ttc" id="_initial_part_positions_8c_html_ad3b7970929b5c9502172afc4e0bc65e8"><div class="ttname"><a href="_initial_part_positions_8c.html#ad3b7970929b5c9502172afc4e0bc65e8">InitParticles_flux</a></div><div class="ttdeci">int InitParticles_flux(int k_current, int first_ind, int last_ind, double weight_p)</div><div class="ttdef"><b>Definition:</b> <a href="_initial_part_positions_8c_source.html#l01095">InitialPartPositions.c:1095</a></div></div>
<div class="ttc" id="structvertex_html_acf83b561551952403148f73fe032f57b"><div class="ttname"><a href="structvertex.html#acf83b561551952403148f73fe032f57b">vertex::numneighb</a></div><div class="ttdeci">unsigned int numneighb</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00127">FuncDef.h:127</a></div></div>
<div class="ttc" id="_initial_part_positions_8c_html_a0610d22ad530758b4de42acd23a5df92"><div class="ttname"><a href="_initial_part_positions_8c.html#a0610d22ad530758b4de42acd23a5df92">TimeFromMatrix</a></div><div class="ttdeci">double TimeFromMatrix(double pdist)</div><div class="ttdef"><b>Definition:</b> <a href="_initial_part_positions_8c_source.html#l01062">InitialPartPositions.c:1062</a></div></div>
<div class="ttc" id="structcontam_html_a70e44829116675e21bda7422f0c2b777"><div class="ttname"><a href="structcontam.html#a70e44829116675e21bda7422f0c2b777">contam::velocity</a></div><div class="ttdeci">double velocity[2]</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00180">FuncDef.h:180</a></div></div>
<div class="ttc" id="_func_def_8h_html_a821d6a33fadb1aacb1a231c2c38cf7ca"><div class="ttname"><a href="_func_def_8h.html#a821d6a33fadb1aacb1a231c2c38cf7ca">particle</a></div><div class="ttdeci">struct contam * particle</div><div class="ttdef"><b>Definition:</b> <a href="main_8c_source.html#l00041">main.c:41</a></div></div>
<div class="ttc" id="structcontam_html_aacd8d2e19b4a3c689e96cf9693218777"><div class="ttname"><a href="structcontam.html#aacd8d2e19b4a3c689e96cf9693218777">contam::prev_pos</a></div><div class="ttdeci">double prev_pos[2]</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00186">FuncDef.h:186</a></div></div>
<div class="ttc" id="structcontam_html_a5b317ab55ab5b48dc82c9a473c76a41e"><div class="ttname"><a href="structcontam.html#a5b317ab55ab5b48dc82c9a473c76a41e">contam::weight</a></div><div class="ttdeci">double weight[3]</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00198">FuncDef.h:198</a></div></div>
<div class="ttc" id="_initial_part_positions_8c_html_a30b78ca10a78d9e5d6de6d5831a5557a"><div class="ttname"><a href="_initial_part_positions_8c.html#a30b78ca10a78d9e5d6de6d5831a5557a">InitParticles_np</a></div><div class="ttdeci">int InitParticles_np(int k_current, int firstnd, int lastnd, int parts_fracture, int first_ind, int last_ind)</div><div class="ttdef"><b>Definition:</b> <a href="_initial_part_positions_8c_source.html#l00660">InitialPartPositions.c:660</a></div></div>
<div class="ttc" id="structcontam_html_a1747a7b7a63f22656afa1df6462279e8"><div class="ttname"><a href="structcontam.html#a1747a7b7a63f22656afa1df6462279e8">contam::intcell</a></div><div class="ttdeci">unsigned int intcell</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00195">FuncDef.h:195</a></div></div>
<div class="ttc" id="_func_def_8h_html_ab71244130eef93109128572e8ef83c2b"><div class="ttname"><a href="_func_def_8h.html#ab71244130eef93109128572e8ef83c2b">fracture</a></div><div class="ttdeci">struct material * fracture</div><div class="ttdef"><b>Definition:</b> <a href="main_8c_source.html#l00039">main.c:39</a></div></div>
<div class="ttc" id="structinpfile_html_a57ff2cce78c35a8a128927da9eeccfbc"><div class="ttname"><a href="structinpfile.html#a57ff2cce78c35a8a128927da9eeccfbc">inpfile::param</a></div><div class="ttdeci">double param</div><div class="ttdef"><b>Definition:</b> <a href="_initial_part_positions_8c_source.html#l00013">InitialPartPositions.c:13</a></div></div>
<div class="ttc" id="_func_def_8h_html_a6d971fbeb47db134f84ff99053b7c47d"><div class="ttname"><a href="_func_def_8h.html#a6d971fbeb47db134f84ff99053b7c47d">nzone_in</a></div><div class="ttdeci">unsigned int nzone_in</div><div class="ttdef"><b>Definition:</b> <a href="main_8c_source.html#l00035">main.c:35</a></div></div>
<div class="ttc" id="_func_def_8h_html_aa27b8bb0d1e8acff46eaadd07476ed04"><div class="ttname"><a href="_func_def_8h.html#aa27b8bb0d1e8acff46eaadd07476ed04">totalFluxIn</a></div><div class="ttdeci">double totalFluxIn</div><div class="ttdef"><b>Definition:</b> <a href="main_8c_source.html#l00052">main.c:52</a></div></div>
<div class="ttc" id="_func_def_8h_html_ae04e94bac8fc762dd97706043f25eaa4"><div class="ttname"><a href="_func_def_8h.html#ae04e94bac8fc762dd97706043f25eaa4">nodezonein</a></div><div class="ttdeci">unsigned int * nodezonein</div><div class="ttdef"><b>Definition:</b> <a href="main_8c_source.html#l00036">main.c:36</a></div></div>
<div class="ttc" id="structcontam_html_a11c5e87e6e62d258648d5a02472251d5"><div class="ttname"><a href="structcontam.html#a11c5e87e6e62d258648d5a02472251d5">contam::time</a></div><div class="ttdeci">double time</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00201">FuncDef.h:201</a></div></div>
<div class="ttc" id="_func_def_8h_html_a67b3bde25c933874d93fd30f06b55f65"><div class="ttname"><a href="_func_def_8h.html#a67b3bde25c933874d93fd30f06b55f65">np</a></div><div class="ttdeci">unsigned int np</div><div class="ttdef"><b>Definition:</b> <a href="main_8c_source.html#l00034">main.c:34</a></div></div>
<div class="ttc" id="_initial_part_positions_8c_html_a0355cfa365c20964163da103eec2f0e5"><div class="ttname"><a href="_initial_part_positions_8c.html#a0355cfa365c20964163da103eec2f0e5">InitInMatrix</a></div><div class="ttdeci">void InitInMatrix()</div><div class="ttdef"><b>Definition:</b> <a href="_initial_part_positions_8c_source.html#l00968">InitialPartPositions.c:968</a></div></div>
<div class="ttc" id="structinpfile_html_a23ae60bf94767194de9a09a4eee8cc1d"><div class="ttname"><a href="structinpfile.html#a23ae60bf94767194de9a09a4eee8cc1d">inpfile::flag</a></div><div class="ttdeci">long int flag</div><div class="ttdef"><b>Definition:</b> <a href="_initial_part_positions_8c_source.html#l00012">InitialPartPositions.c:12</a></div></div>
<div class="ttc" id="structmaterial_html_a572b56f5488c560ee0fe118592b660a7"><div class="ttname"><a href="structmaterial.html#a572b56f5488c560ee0fe118592b660a7">material::rot2mat</a></div><div class="ttdeci">double rot2mat[3][3]</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00098">FuncDef.h:98</a></div></div>
<div class="ttc" id="_initial_part_positions_8c_html_a2e978219b990b7dfb64b4f9d567054db"><div class="ttname"><a href="_initial_part_positions_8c.html#a2e978219b990b7dfb64b4f9d567054db">InitParticles_ones</a></div><div class="ttdeci">int InitParticles_ones(int k_current, double inter_p[][4], int fracture_n, int parts_fracture, int ii, double thirdcoor, int zonenumb_in, int first_ind, int last_ind)</div><div class="ttdef"><b>Definition:</b> <a href="_initial_part_positions_8c_source.html#l00764">InitialPartPositions.c:764</a></div></div>
<div class="ttc" id="_func_def_8h_html_a9d0d7da8d3c9e75857a16e168969f89d"><div class="ttname"><a href="_func_def_8h.html#a9d0d7da8d3c9e75857a16e168969f89d">node</a></div><div class="ttdeci">struct vertex * node</div><div class="ttdef"><b>Definition:</b> <a href="main_8c_source.html#l00040">main.c:40</a></div></div>
<div class="ttc" id="_func_def_8h_html_a5799b26fb87c6d98cd44352eb84a2e7e"><div class="ttname"><a href="_func_def_8h.html#a5799b26fb87c6d98cd44352eb84a2e7e">timeunit</a></div><div class="ttdeci">double timeunit</div><div class="ttdef"><b>Definition:</b> <a href="main_8c_source.html#l00051">main.c:51</a></div></div>
<div class="ttc" id="_func_def_8h_html_a484100bd5902d4efcdd15c21d1cb9048"><div class="ttname"><a href="_func_def_8h.html#a484100bd5902d4efcdd15c21d1cb9048">Control_File_Optional</a></div><div class="ttdeci">struct inpfile Control_File_Optional(char fileobject[], int ctr)</div><div class="ttdef"><b>Definition:</b> <a href="_read_grid_init_8c_source.html#l01135">ReadGridInit.c:1135</a></div></div>
<div class="ttc" id="_func_def_8h_html_a9c5371905330e52a346417559120efec"><div class="ttname"><a href="_func_def_8h.html#a9c5371905330e52a346417559120efec">npart</a></div><div class="ttdeci">unsigned int npart</div><div class="ttdef"><b>Definition:</b> <a href="main_8c_source.html#l00033">main.c:33</a></div></div>
<div class="ttc" id="structcontam_html_a0b9253ce268d781958145f9d110df2de"><div class="ttname"><a href="structcontam.html#a0b9253ce268d781958145f9d110df2de">contam::fracture</a></div><div class="ttdeci">unsigned int fracture</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00189">FuncDef.h:189</a></div></div>
<div class="ttc" id="structcontam_html_ad50f20b01e37af1fedfd7bbe0cd1d7e3"><div class="ttname"><a href="structcontam.html#ad50f20b01e37af1fedfd7bbe0cd1d7e3">contam::fl_weight</a></div><div class="ttdeci">double fl_weight</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00204">FuncDef.h:204</a></div></div>
<div class="ttc" id="structinpfile_html"><div class="ttname"><a href="structinpfile.html">inpfile</a></div><div class="ttdef"><b>Definition:</b> <a href="_initial_part_positions_8c_source.html#l00010">InitialPartPositions.c:10</a></div></div>
<div class="ttc" id="_func_def_8h_html_a2afbea1a5d0b14d8a098f551fe055420"><div class="ttname"><a href="_func_def_8h.html#a2afbea1a5d0b14d8a098f551fe055420">Control_File</a></div><div class="ttdeci">struct inpfile Control_File(char fileobject[], int ctr)</div><div class="ttdef"><b>Definition:</b> <a href="_read_grid_init_8c_source.html#l01090">ReadGridInit.c:1090</a></div></div>
<div class="ttc" id="structcontam_html_a858c05f3221a70f6e212c1e7e172252c"><div class="ttname"><a href="structcontam.html#a858c05f3221a70f6e212c1e7e172252c">contam::cell</a></div><div class="ttdeci">unsigned int cell</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00192">FuncDef.h:192</a></div></div>
<div class="ttc" id="structposit3d_html"><div class="ttname"><a href="structposit3d.html">posit3d</a></div><div class="ttdef"><b>Definition:</b> <a href="_initial_part_positions_8c_source.html#l00016">InitialPartPositions.c:16</a></div></div>
<div class="ttc" id="_func_def_8h_html"><div class="ttname"><a href="_func_def_8h.html">FuncDef.h</a></div></div>
<div class="ttc" id="structvertex_html_a3473aebf2ce957df8c02095cc98ebb39"><div class="ttname"><a href="structvertex.html#a3473aebf2ce957df8c02095cc98ebb39">vertex::coord_xy</a></div><div class="ttdeci">double coord_xy[6]</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00118">FuncDef.h:118</a></div></div>
<div class="ttc" id="structelement_html_aadc95dc4817009223151c40387f40e06"><div class="ttname"><a href="structelement.html#aadc95dc4817009223151c40387f40e06">element::node_ind</a></div><div class="ttdeci">unsigned int node_ind[3]</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00167">FuncDef.h:167</a></div></div>
<div class="ttc" id="structvertex_html_a583032cadaa99e8f5c044d3bd10191cf"><div class="ttname"><a href="structvertex.html#a583032cadaa99e8f5c044d3bd10191cf">vertex::cells</a></div><div class="ttdeci">unsigned int ** cells</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00139">FuncDef.h:139</a></div></div>
<div class="ttc" id="structcontam_html"><div class="ttname"><a href="structcontam.html">contam</a></div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00177">FuncDef.h:177</a></div></div>
<div class="ttc" id="_func_def_8h_html_a205add90a27292ac5f18188fa35bcdf5"><div class="ttname"><a href="_func_def_8h.html#a205add90a27292ac5f18188fa35bcdf5">OpenFile</a></div><div class="ttdeci">FILE * OpenFile(char filen[120], char fileopt[2])</div><div class="ttdef"><b>Definition:</b> <a href="_read_grid_init_8c_source.html#l00728">ReadGridInit.c:728</a></div></div>
<div class="ttc" id="structelement_html_af51b08146a1c495aadf45e7f085e1267"><div class="ttname"><a href="structelement.html#af51b08146a1c495aadf45e7f085e1267">element::fracture</a></div><div class="ttdeci">unsigned int fracture</div><div class="ttdef"><b>Definition:</b> <a href="_func_def_8h_source.html#l00164">FuncDef.h:164</a></div></div>
</div><!-- fragment --></div><!-- contents -->
</div><!-- doc-content -->
<!-- start footer part -->
<div id="nav-path" class="navpath"><!-- id is needed for treeview function! -->
  <ul>
    <li class="navelem"><a class="el" href="dir_f80ddb2df819b297fd9c7beccba09984.html">Docs_PT</a></li><li class="navelem"><a class="el" href="dir_966eb28a5fc4cd2ac9e3086fcb212f97.html">ParticleTracking_source</a></li><li class="navelem"><a class="el" href="_initial_part_positions_8c.html">InitialPartPositions.c</a></li>
    <li class="footer">Generated by
    <a href="http://www.doxygen.org/index.html">
    <img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.8.15 </li>
  </ul>
</div>
</body>
</html>
