><span class="n">ar</span><span class="o">=</span><span class="n">F</span><span class="p">%</span><span class="n">AB</span><span class="p">%</span><span class="n">a</span>
    <span class="c">! Calculating density of minor radial annulus    </span>
    <span class="k">do </span><span class="n">ii</span><span class="o">=</span><span class="mi">1_idef</span><span class="p">,</span><span class="n">F</span><span class="p">%</span><span class="n">dim_1D</span>
       <span class="c">! NGP weighting</span>
      