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