vn_avg (t_fwd*cd.vbtot+t_rev*cd2.vbtot) "Time-averaged pulse reverse plating velocity"
vn_bnd -i_avg/(2*F_const)*M/rho "Boundary velocity for upper boundary"