g_j Function

  1. korc_collisions.f90
  2. korc_collisions
  3. g_j

public function g_j(i, v)

Arguments

Type IntentOptional AttributesName
integer, intent(in) :: i
real(kind=rp), intent(in) :: v

Return Value real(kind=rp)


Contents

Source Code

g_j

Source Code

  function g_j(i,v)
    INTEGER, INTENT(IN) 	:: i
    REAL(rp), INTENT(IN) 	:: v   
    REAL(rp)  :: gam
    REAL(rp)  :: p
    REAL(rp)  :: g_j

    gam=1/sqrt(1-v**2)
    p=v*gam
    
    g_j=2._rp/3._rp*((cparams_ms%Zo(i)**2-cparams_ms%Zj(i)**2)* &
         log((p*cparams_ms%aZj(i))**(3._rp/2._rp)+1)- &
         (cparams_ms%Zo(i)-cparams_ms%Zj(i))**2* &
         (p*cparams_ms%aZj(i))**(3._rp/2._rp)/ &
         ((p*cparams_ms%aZj(i))**(3._rp/2._rp)+1))

!    write(6,'("g_j: ",E17.10)') g_j
    
  end function g_j