I cooked up the following NCL script to convert full 3D model output from sigma-pressure hybrid coordinates to user defined pure-pressure coordinates. The script uses the inbuilt vinth2p() function for interpolation. While it is also possible to interpolate the whole file (with several time slices) at once, to ensure efficient memory usage, in this script interpolatation is implemented one time-slice at a time. The "efficient" method to write the interpolated variables to the netCDF files as described on the NCAR website has been used, thus ensuring speed. A 'sequence' of well numbered files can be converted at once. I have extensively used this script to interpolate my model outputs onto pressure coordinates.

Click to Download