p_s :: (real(10000 .. 107000) ~ "Pa" field (x(grid), y(grid)) on
        (i = 0 .. nlon + 2 by j = 0 .. nlat + 2)
       )