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.

• 介质单元的速度用双线性函数表示,推导了计算Jacobi矩阵元素的解析公式,提高了精度和速度;