Add support for writing HEATING_CALCULATED forcing files (in which the net surface heat flux is calculated within the model rather than from an external source).
Update the search for the surface pressure level based on the Met Office documentation.
Remove some redundancy and order things more logically.
Ignore the third dimension when finding out if the input data is the wrong way around.