For cylindrical unit vectors: now check for scalar product and make sure the radial component is never negative.
For cylindrical unit vectors: now check for scalar product and make sure the radial component is never negative.

mzemp/iof