Could somebody please explain the difference between rho, theta and pressure levels within UM diagnostics?
On the one hand, MetOffice documentation seems to suggest that rho = upper boundary of the gridcell and theta = lower boundary of the gridcell. On the other hand, the UM documentation and the vertlevs file ($UMDIR/vn13.2/ctldata/vert/vertlevs_L85_50t_35s_85km) suggest that rho = centre of gridcell and theta = bounds of the gridcell.
And I cannot find any definition of pressure levels.
"I believe Rho is the middle and Theta is the edges of the cell, but some variables are only on particular grids, i.e. density is on the rho grid and potential temperature on the theta grid. Pressure level fields are interpolated at runtime and need to have the relevant Heaviside function considered to properly process them."
If you want further information or have further scientific questions that you can’t find answers to in the UM documentation, I suggest posting on the UM github discussions board: MetOffice/simulation-systems Um · Discussions · GitHub