eta_p :: (real ~ " kg s^-3 / m" field (x(grid), y(grid), z(grid)) on
(i = 1 .. nlon + 1 by
j = 1 .. nlat + 1 by
k = 1 .. nlev
)
)
Scalar Form
eta_p := ( int(( 1 / (h_x * h_y)
* ( df(h_y * u_aux, x)
+ df(h_x * v_aux, y)
)
), z ~ "1" = point(z) ~ "1" .. nlev + 1)
+ p_s_t * (1 - B)
)
Dimensioned
eta_p := ( int(( 1 / (h_x * h_y)
* ( df(h_y * u_aux, x)
+ df(h_x * v_aux, y)
)
), z ~ "1" = point(z) ~ "1" .. nlev + 1 ~ "1")
+ p_s_t * (1 ~ "1" - B)
)
Discretized
eta_p := ( p_s_t * (1 - B)
+ ( 1 / (hxt * hyt)
* ( df_h(hyu * mid_h(u_aux, z = point(z) .. nlev + 1), x)
+ df_h(hxv * mid_h(v_aux, z = point(z) .. nlev + 1), y)
)
)
)
CSE Eliminated
(eta_p :=
varidx(t_70, [i, j, k]) + varidx(p_s_t, [i, j]) * (1 - varidx(B, [k])))