u_aux :: (real ~ " kg s^-3" field (x(half), y(grid), z(half)) on
          (i = 0 .. nlon + 1 by
           j = 1 .. nlat + 1 by
           k = 1 .. nlev
          )
         )


Scalar Form

u_aux := u * df(p, z ~ "1")

Dimensioned

u_aux := u * df(p, z ~ "1")

Discretized

u_aux := u * df_g(ave_g(p, x), z)