generate histograms for 2D fields on vertical levels

[Sylviabohnenstengel]

• [x] develop a script to look into data range for any 2D field for the whole model simulation.
• [x] as option enable selection of level
• [x] as an option generate plots for selected model or pressure levels using the slider option for time.
• [x] either a pdf or box whisker
• [ ] this script can get extended to display the stats against forecast range in ticket #493
• [ ] develop test for plotting function
• [x] beautify the histogram
• [x] ensure documentation is up to date.
• [ ] develop LFRic recipe
• [ ] ensure it works for ensembles using the postage stamp option
• [ ] Should the postage stamp option not only realisation coordinate, but maybe also vertical level?
• [ ] add min, max, man,median to PDF curve

[Sylviabohnenstengel]

Recipe:

• read cube (lat, lon, plevel, time)
• use sequence and serial coordinate to calculate pdf for each level and time
• plot histogram against sequence coordinate

[Sylviabohnenstengel]

The recipe currently plots a simple histogram of the data range for a chosen variable for every pressure level and every validity time.

using tox the test for def test_plot_vertical_line_series(tmp_working_dir is failing as it first requires ticket #567 to be reviewed and included.

[Sylviabohnenstengel]

@jfrost-mo Do you have ensemble data that work with CSET at the moment that I can point to? i would like to try the new internal histogram plotting function.

[Sylviabohnenstengel]

@jfrost-mo how can we ensure that we use the same bins and limits for the xaxis in the histogram? is there a way to read out infromation from all cubes across all realisation coordinates, time coordinates and pressure level coordinates to set the min and max limits for the xrange and also the bin width? In the example here different bins are used for the plots

I assume this can be done in the collate step if we write out all cubes into a single netcdf. Would this be written into the intermediate directory or collate from there?

Imanaged to solve this in the collate step

[Sylviabohnenstengel]

Needed to use the following to link the plot ino the figure.

Put the shared colorbar in its own axes.

colorbar_axes = fig.add_axes([0.15, 0.07, 0.7, 0.03])
colorbar = fig.colorbar(plot, colorbar_axes, orientation="horizontal")
colorbar.set_label(f"{cube.name()} / {cube.units}")

How can I do this properly if I don't need a common colourbar?

all dealt with in the collate step now as it has access to the whole cube

[Sylviabohnenstengel]

PS: https://github.com/MetOffice/CSET/pull/594