initialize_params Subroutine

  SUBROUTINE initialize_params(params)
    TYPE(KORC_PARAMS), INTENT(IN) :: params
    REAL(rp) :: max_pitch_angle
    REAL(rp) :: min_pitch_angle
    REAL(rp) :: max_energy
    REAL(rp) :: min_energy
    REAL(rp) :: Zeff
    REAL(rp) :: E
    REAL(rp) :: k
    REAL(rp) :: t
    REAL(rp) :: Bo
    REAL(rp) :: lambda
    REAL(rp) :: A_fact

    NAMELIST /ExperimentalPDF/ max_pitch_angle,min_pitch_angle,max_energy,min_energy,Zeff,E,k,t,Bo,lambda, &
         A_fact

    open(unit=default_unit_open,file=TRIM(params%path_to_inputs),status='OLD',form='formatted')
    read(default_unit_open,nml=ExperimentalPDF)
    close(default_unit_open)

    pdf_params%max_pitch_angle = max_pitch_angle
    pdf_params%min_pitch_angle = min_pitch_angle
    pdf_params%min_energy = min_energy*C_E ! In Joules
    pdf_params%max_energy = max_energy*C_E ! In Joules
    pdf_params%Zeff = Zeff
    pdf_params%E = E
    pdf_params%k = k
    pdf_params%t = t
    pdf_params%Bo = Bo
    pdf_params%lambda = lambda

    pdf_params%max_p = SQRT((pdf_params%max_energy/(C_ME*C_C**2))**2 - 1.0_rp) ! In units of mc
    pdf_params%min_p = SQRT((pdf_params%min_energy/(C_ME*C_C**2))**2 - 1.0_rp) ! In units of mc

    pdf_params%fGo = &
         IntGamma(SQRT(pdf_params%min_p**2.0_rp + 1.0_rp),SQRT(pdf_params%max_p**2.0_rp + 1.0_rp),pdf_params%k,pdf_params%t/xo)

    pdf_params%A_fact = A_fact
  END SUBROUTINE initialize_params