The velocity in a rectangular cell is represented as a bilinear function with depth and lateral position, and the analytic formulas for computing Jacobin matrix elements are deduced to improve the computation precision.