n class="n">Ct</span> <span class="o">+</span> <span class="n">Cr</span><span class="p">)</span><span class="o">/</span><span class="mi">2</span>
    <span class="n">G</span> <span class="o">=</span> <span class="mf">0.5</span><span class="o">*</span><span class="p">(</span><span class="mi">1</span> <span class="o">-</span> <span class="n">np</span><span class="o">.</span><span class="n">sqrt</span><span class="p">((</span><span class="n">Cavg</span><span class="o">**</span><span class="mf">7.0</span><span 