span><span class="p">%</span><span class="n">Ro</span><span class="o">*</span><span class="n">h_params</span><span class="p">%</span><span class="n">Eo</span><span class="o">/</span><span class="n">R</span>
    <span class="n">endif</span>
       
<span class="c">!    E_G=h_params%E*exp(-(rm/h_params%sigma_E)**2/2)</span>
    <span class="n">Z_G</span><span class="o">=</span><span class="n">h_params</span><span class="p">%</span><span class="n">Zeff</span><span class="o">*</span><span class="nb">exp</span><