Switched to NetCDF4 storage of predictive output

Changes:

1. Where xarrays were pickled to disc, we now use xarray.to_netcdf
2. Where xarrays were read from Python pickles, we now use xarray.load_dataset

(cherry picked from commit 63519eb3)