rp)</span>
<span class="c">!       arg3=MIN((r_1D(ii))/(sqrt(2._rp)*sigr),10._rp)</span>

<span class="c">!       Ai(ii)=sqrt(C_PI*sigr)*(sqrt(2._rp)*sigr*(exp(-arg)- &amp;</span>
<span class="c">!            exp(-arg1))+r_1D(ii)*sqrt(C_PI)* &amp;</span>
<span class="c">!            (erf(arg2)-erf(-arg3)))</span>
<span class="c">!       Vden(ii)=Vpart(ii)/Ai(ii)</span>

    <span class="k">end do</span>
<span class="k">    </span>
<span class="k">  end subroutine </span><span class="n">calculate_SC_p</span>