differential equations – Solution of coupled ODE eigenvalues

I want a value solution of the following coupled ODEs: But the code indicating errors: the BCs are incomplete.
enter the description of the image here

                gamma = 1.67; ep = 0.01; d = 0.05; H = 1.1; gpara = 0.38;
Eqns = {Rho[x]/ Exp[-x/H]* gpara -
1 / gamma / Exp[-x/H]*(RE[Rho[Rho[Rho[Rho[x], X]+ D[T[T[T[T[x], X]+
1.3 * ep * D[V[V[V[V[x], {x, 2}

+ I * lambda * V[x] == 0,
re[Exp[Exp[Exp[Exp[-x/H]* V[x], X]- I * lambda * Rho[x] == 0,
re[Exp[Exp[Exp[Exp[-x/H]* V[x], X]- V[x]*RE[Exp[Exp[Exp[Exp[-x/H]X]-
gamma * Exp[-x/H]*RE[V[V[V[V[x], X]+ gamma * d * D[T[T[T[T[x], {x, 2}

+
I * lambda * t[x] == 0};
BC = {V[0] == 0, V[1] == 0, V & # 39;[0] == 0, V & # 39;[1] == 0, T[0] == 0,
T[1] == 0, T & # 39;[0] == 0, T & # 39;[1] == 0, Rho[0] == 0, Rho[1] == 0};
ground =
ParametricNDSolve[{Eqns, BCs}, {V, Rho, T}, {x, 0, 1}, {lambda}];
FindRoot[sol[lambda][1], {lambda, 1}]

Q1: I have to find the numerical values ​​of lambda (eigenvalues),
Q2: In the equations above, gpara is considered constant, but it is a variable with 50 rows and 1 column in the file gpara.txt, which is in another directory. I import Gpara by Import[“E\ ……\ gpara.txt”] to use in the equations above, but does not work too. Please help to get the solution.