bdata.getByName has still problems with doble, need to set in main omega_nu by hand. Allocate IPPL vector with neutrinos and deduce from the size of the vector the array lenghts in the initializer