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