theta atan2(y,x)[rad] "Angle along circumference"
h_f C*(1+ep*cos(theta)) "Lubricant film thickness"
u_b -U*sin(theta) "X-component of journal velocity"
v_b U*cos(theta) "Y-component of journal velocity"