Save the number of observations in the Netcdf

Presently we are only saving the obs, it would be good to update the netcdf to also contain the num of obs.

Guillaime suggested a list of statistics that the user can select like [num_obs, mean etc.]