Added w_velocity function, and some very basic value checking.  omegaHz and w calculations appear to be correct.