In Joules
function VTe_wu(Te) REAL(rp), INTENT(IN) :: Te !! In Joules REAL(rp) :: VTe_wu VTe_wu = SQRT(2.0_rp*Te/C_ME) end function VTe_wu