```Z :: (real(- infinity .. infinity) ~ " m s / kg" field ( x(half)                                                       , y(half)                                                       , z(half)                                                       ) on      (i = 0 .. nlon by j = 0 .. nlat by k = 1 .. nlev)     )

Scalar Form

Z := ( ( (df(h_y * v, x) + f * ave(ave(h_x * h_y, x), y))       - df(h_x * u, y)       )     / ave(ave(h_x * h_y * df(p, z ~ "1"), x), y)     )
Dimensioned

Z := ( ( (df(h_y * v, x) + f * ave(ave(h_x * h_y, x), y))       - df(h_x * u, y)       )     / ave(ave(h_x * h_y * df(p, z ~ "1"), x), y)     )
Discretized

Z := ( ( (df_g(hyv * v, x) + f * ave_g(ave_g(hxt * hyt, x), y))       - df_g(hxu * u, y)       )     / ave_g(ave_g(hxt * hyt * df_g(p, z), x), y)     )
CSE Eliminated

Z := 4 * ( ( varidx(t_36, [i, j])           - ( 1.56961230576048E-7             * ( ( rdlam                 * ( varidx(hyv, [i, j]) * varidx(v, [i, j, k])                   - varidx(hyv, [i + 1, j]) * varidx(v, [i + 1, j, k])                   )                 )               - rdth * ( varidx(hxu, [i, j]) * varidx(u, [i, j, k])                        - varidx(hxu, [i, j + 1]) * varidx(u, [i, j + 1, k])                        )               )             )           )         / ( varidx(t_31, [i, j, k])           + varidx(t_31, [i, j + 1, k])           + varidx(t_31, [i + 1, j, k])           + varidx(t_31, [i + 1, j + 1, k])           )         )```