epsl "if(z<-H_active/2,epsl_lower,if(z>H_active/2,epsl_upper, (z+H_active/2)/H_active*(epsl_upper-epsl_lower)+epsl_lower))" "Electrolyte volume fraction"
Eeq_cell intop_point_o2(we.Eeq_o2er1)-intop_point_h2(we.Eeq_h2er1) "Cell equilibrium voltage"
Etherm_cell intop_point_o2(we.Etherm_o2er1)-intop_point_h2(we.Etherm_h2er1) "Cell thermoneutral voltage"
I_stack intop_cc(we.nIs) "Stack current"
I_h2 -intop_h2_electrodes(we.iloc_h2er1) "Total hydrogen production current"
Eff_coulombic I_h2/(I_stack*N_cells) "Current-to-hydrogen coloumbic efficiency"
Eff_energy I_h2*Eeq_cell/(E_stack*I_stack) "Electric energy-to-hydrogen energy efficiency"