function [q_1e] = q1e_func_GenI_alt(x,d,p,u)

h_2 = x(1);

q_1e_nom = p(21);
h_2_nom = p(30);

Cv_1 = q_1e_nom/h_2_nom^(0.5);

q_1e = Cv_1 * h_2^0.5;

end