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