Background on the surface interpolation:
The original grid (at 0.2° resolution) has invalid (NaN) values in non-land oceanic Arctic regions.
Surface will interpolate values for every cell of a grid. It’s trivial to mask the land values.
After interpolating with surface, how can the newly interpolated cells be re-set to match (almost) those NaN-valued non-land cells found in the original grid?