u_profile u_bulk*(y-h1)/h2 Inlet/outlet velocity profile i_loc -tds.ntflux_c*F_const*2 Local current density i_loc_norm -i_loc/maxop1(-i_loc) Normalized local current density