...
# Retrieve land-sea mask and interpolate to LL grid
lsm = retrieve(
type : "an",
date : -1,
param : "lsm",
grid : [1.5,1.5],
levtype : "sfc"
)
# The following line forces the values to 0 or 1.
lsm = lsm > 0.5
# Now compute the average value on land and on sea
land = integrate(field, lsm)
sea = integrate(field, not lsm)
number or list Anchor FieldsetFunctions-interpolate FieldsetFunctions-interpolate interpolate
( fieldset,list )
number or list interpolate
( fieldset,number,number )
...