function CA(v) REAL(rp), INTENT(IN) :: v REAL(rp) :: CA REAL(rp) :: x x = v/cparams_ss%VTe CA = cparams_ss%Gammac*psi(x)/v end function CA