mean_v_shear = 0.0
mean_v_shear = mean_v_shear +
& v_shear(1) + vmid - mean_v_shear
_RL ub, vb, fb, mean_u_shear, mean_v_shear, umid, vmid