Non-monotonic axes?

[David Oxilia <oxilia@xxxxxxxxxxxx>]

I'd like to convert a binary dataset into netCDF. The data consists
of global 1x1 degree gridded fields with 33 vertical levels at a
single time point.

With exception of the vertical axis the conversion is straight
forward. The problem with the vertical axis is that it corresponds
to 33 levels of non-monotonic increasing depths. 

Is it somehow possible to define a non-monotonic depth axis? At
the moment the code to read in the binary data is as follows:

  ! Define temperature grid
  def axis/x=0.5E:0.5W/npoints=360/unit=deg xax
  def axis/y=89.5S:89.5N/npoints=180/unit=deg yax
  def axis/z=1:33/npoints=33/unit="levels" zax
  def grid/x=xax/y=yax/z=zax gg

  ! Read temperature data
  file/form=unf/order=xyz/grid=gg/var=temp temp_u.bin
  set var/tit="Temperature"/unit="deg C" temp

Is it possible to modify this in order to have a real depth axis
at the 33 vertical levels without having to re-grid onto a
monotonic depth axis?

thanks, David.


*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*
David Oxilia			Tel: (541) 737-2368
COAS:Oregon State University	Fax: (541) 737-2064
Ocean Admin. Bldg. 104		WWW: http://www.oce.orst.edu
Corvallis, OR 97331-5503	Email: oxilia@hejira.oce.orst.edu
*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*