where</span><span class="p">(</span><span class="n">Cpprod</span> <span class="o">==</span> <span class="mi">0</span><span class="p">)]</span> <span class="o">=</span> <span class="mi">0</span>

    <span class="n">T</span> <span class="o">=</span> <span class="mi">1</span> <span class="o">-</span> <span class="mf">0.17</span><span class="o">*</span><span class="n">np</span><span class="o">.</span><span class="n">cos</span><span class="p">((</span><span class="n">hp</span> <span class="o">-</span> <span cla