h_y :: real ~ "1" field (x(grid), y(grid));
h_y :: real ~ "1" field (x(grid), y(half));
h_y :: real ~ "1" field (x(half), y(grid))