We very soon have mask capability in CSET which will replace data with np.nan where data is masked out. For subsequent operators (as masking is unlikely to be the only operator used in a recipe), we should check that other operators can handle these (e.g. domain averages do the appropriate nanmean).