Integrate obsevations (netcdf)

### What problem does your feature request solve?

Load in observational data (netcdf input) to produce plots (eg. box-whisker profile plots) in order to compare with model output

### Describe the solution you'd like

Basic recipe to handle obs data:

- [ ] fetch data
- [ ] read data
- [ ] manipulate/generate statistics and plots (eg. below)


### Describe alternatives you've considered

No files, tests, or entry points are named. Start by clarifying the observational data source and how NetCDF input should connect to existing model-output workflows. Done would cover fetching and reading the data, then generating the requested statistics and comparison plots.

Written by the indexing model from the issue text.