Page History

...

# 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)

...