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_expt
    pdf_params%min_pitch_angle = min_pitch_angle_expt
    pdf_params%min_energy = min_energy_expt*C_E ! In Joules
    pdf_params%max_energy = max_energy_expt*C_E ! In Joules
    pdf_params%Zeff = Zeff_expt
    pdf_params%E = E_expt
    pdf_params%k = k_expt
    pdf_params%t = t_expt
    pdf_params%Bo = Bo_expt
    pdf_params%lambda = lambda_expt

    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_expt
  END SUBROUTINE initialize_params