Initialise module variables
subroutine init_neoclassics_module
!! Initialise module variables
implicit none
species = (/"e","D","T","a"/)
densities = 0
temperatures = 0
dr_densities = 0
dr_temperatures =0
roots = 0
weights = 0
nu = 0
nu_star = 0
nu_star_averaged = 0
vd = 0
KT = 0
Er = 0.0
iota = 1.0d0
D11_mono = 0
D11_plateau = 0
D111 = 0
D112 = 0
D113 = 0
q_flux = 0
Gamma_flux = 0
D31_mono = 0
eps_eff = 1d-5
end subroutine init_neoclassics_module