Description: Update pv-vector discarding velocity.
The source code for the Fortran routine PVUP is available: