init_neoclassics Subroutine

    subroutine init_neoclassics(r_eff, eps_eff, iota)
    !! Constructor of the neoclassics object from the effective radius,
    !! epsilon effective and iota only.
    !
    ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

    real(dp), intent(in) :: r_eff
    real(dp), dimension(4,no_roots) :: mynu
    real(dp), intent(in)::eps_eff, iota

    !! This should be called as the standard constructor
    call init_profile_values_from_PROCESS(r_eff, densities, temperatures, dr_densities, dr_temperatures)
    call gauss_laguerre_30_roots(roots)
    call gauss_laguerre_30_weights(weights)


    KT = neoclassics_calc_KT()
    nu = neoclassics_calc_nu()
    nu_star = neoclassics_calc_nu_star()
    nu_star_averaged = neoclassics_calc_nu_star_fromT(iota)
    vd = neoclassics_calc_vd()

    D11_plateau = neoclassics_calc_D11_plateau()

    D11_mono = neoclassics_calc_D11_mono(eps_eff) !for using epseff

    !alternatively use:  = myneo%interpolate_D11_mono() !

    D111 = neoclassics_calc_D111()

    D112 = neoclassics_calc_D112()
    D113 = neoclassics_calc_D113()

    Gamma_flux = neoclassics_calc_Gamma_flux(densities, temperatures, dr_densities, dr_temperatures)
    q_flux = neoclassics_calc_q_flux()

    ! Return:

    end subroutine init_neoclassics