FIXME `isna.sum()``: count only cells in cfactual which are nan but not in factual something like: `isnan = np.isnan(v) & (~np.isnan(v_orig))`. The problem is, currently this doesn't work because v and v_orig are masked arrays. Should read the v and v_orig from the netcdf file differently such that it is a normal array. #94