phi_s :: (real ~ " m^2 s^-2" field (x(grid), y(grid)) on
(i = 1 .. nlon + 1 by j = 1 .. nlat + 1)
)