Interpolating 3-d Cube on model levels onto single altitude level (horizontal slice)

[Andrew Elvidge]

Hi,

I have a 3-d cube on model vertical levels. In concert with orography data, IRIS has calculated a 3-d 'altitude' field, which is the altitude in m of every data point. I would like to interpolate this 3d cube onto a single vertical altitude in metres, Note that my domain is over orography so using the standard model-altitude is not appropriate.

Can this be done is IRIS?

If not, could anyone recommend a good 3-d interpolation tool that can handle the irregular altitude coordinate of the raw data?

Thanks, Andy

[Andrew Elvidge]

p.s. as mentioned here http://scitools.org.uk/iris/docs/latest/userguide/interpolation_and_regridding.html, the iris.analysis.Linear scheme requires 1D numeric, monotonic, coordinates.

p.p.s. 3-d interpolation tool for python, that is

[Andrew Elvidge]

Given the lack of reply I take it this cannot be done in IRIS...? I have spoken to another IRIS and Unified Model user who agrees that this would be very useful functionality to have in IRIS.

I've also found it difficult to do in Python without very slow for loops.  For information, my solution has been to use a call to Cython to execute a simple linear interpolation script. For large datasets, this is considerably quicker than anything I could muster just in Python.

On Thursday, February 9, 2017 at 4:52:30 PM UTC, Andrew Elvidge wrote:

[Andrew Dawson]

You may be interested in this: https://github.com/SciTools-incubator/python-stratify

It is still under development but I'm sure user feedback would be welcomed.

[Andrew Elvidge]

Hi Andrew, Yes that looks precisely what I'm after, thanks for drawing my attention to this.