span class="c">!          Vden(ii)=Vpart(ii)/(C_PI*dr**2/3)</span>
<span class="c">!       else</span>
<span class="c">!          Vden(ii)=Vpart(ii)/(2*C_PI*dr**2*(ii-1))</span>
<span class="c">!       end if</span>

       <span class="c">! Gaussian weighting</span>
<span class="c">!       arg=MIN(r_1D(ii)**2._rp/(2._rp*sigr**2._rp),100._rp)</span>
<span class="c">!       arg1=MIN((ar-r_1D(ii))**2._rp/(2._rp*sigr**2._rp),100._rp)</span>
<span class="c">!       arg2=MIN((ar-r_1D(ii))/(sqrt(2._rp)*sigr),10._