VelocityPenalty2 to constrain velocities to 1D

@nicholasmalaya, this should be straightforward to use: just add VelocityPenalty2 (and its AdjointStabilization) to the physics line and add a VelocityPenalty2 section for the vertical penalty_val.