   vpll(cc)*(rm(cc)-r_1D(rind))/dr</span>

       <span class="c">! Gaussian weighting</span>
       
<span class="c">!       do ii=1_idef,F%dim_1D</span>
<span class="c">!          arg=MIN((r_1D(ii)-rm(cc))**2._rp/(2._rp*sigr**2._rp),100._rp)</span>
<span class="c">!          Vpart(ii)=Vpart(ii)+1/sqrt(2._rp*C_PI*sigr**2._rp)* &amp;</span>
<span class="c">!               exp(-arg)*vpll(cc)           </span>
<span class="c">!       end do</span>

    <span class="k">end do</span>

<span class="k">    </span