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