neoclassics_calc_q_flux Function

    function neoclassics_calc_q_flux()
        !! Calculates the Energy flux by neoclassicsal energy transport
        !
        ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

        real(dp),dimension(4) ::  z, neoclassics_calc_q_flux


        ! densities = densities
        ! temps = temperatures
        ! dr_densities = dr_densities
        ! dr_temps = dr_temperatures

        z = (/-1.0,1.0,1.0,2.0/)

        q_flux = - densities * temperatures * D112 * ((dr_densities/densities - z * Er/temperatures) + &
                        (D113/D112-3.0/2.0) * dr_temperatures/temperatures )

        neoclassics_calc_q_flux = q_flux
    end function neoclassics_calc_q_flux