<span class="n">cc</span><span class="p">))</span><span class="o">*</span><span class="n">vpll</span><span class="p">(</span><span class="n">cc</span><span class="p">)</span>

       <span class="c">! First-order weighting</span>
<span class="c">!       rind=FLOOR(rm(cc)/dr)+1_ip</span>
<span class="c">!       Vpart(rind)=Vpart(rind)+ &amp;</span>
<span class="c">!            vpll(cc)*(r_1D(rind+1)-rm(cc))/dr</span>
<span class="c">!       Vpart(rind+1)=Vpart(rind+1)+ &amp;</span>
<span class="c">!         