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