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])
)
)