u :: (real ~ " m / s" field (x(half), y(grid), z(half)) on
      (i = 0 .. nlon + 1 by
       j = 0 .. nlat + 1 by
       k = 1 .. nlev
      )
     )