class="o">+</span> <span class="mi">360</span>
    <span class="n">dhp</span><span class="p">[</span><span class="n">np</span><span class="o">.</span><span class="n">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="c1">#dH = 2*np.sqrt(Cpprod)*np.sin(dhp/2*np.pi/180)</span>
    <span class="n">dH</span> <span class="o">=</span> <span class="n">delt