hyu :: (real ~ "1" field (x(half), y(grid)) on
(i = 0 .. nlon + 1 by j = 1 .. nlat + 1)
)