...
# 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 interpolate interpolate interpolate
( fieldset,list )
number or list interpolate
( fieldset,number,number )
...
With n fields in the input fieldsets, if xik, yik are the ith value of the kth input fieldsets and zi is the ith value of the resulting field:
number or list Anchor nearest_gridpoint nearest_gridpoint nearest_gridpoint
( fieldset,list )
number or list nearest_gridpoint
( fieldset,number,number )vector or list nearest_gridpoint ( fieldset,vector,vector )
...