Re: [ferret_users] Vertical regridding without interpolation

[Ansley Manke <ansley.b.manke@xxxxxxxx>]

Hi Matthieu,
Do you mean that the levels in your grid correspond to known depths, so that you would like to simply assign a depth to each level?  If that is the case, then use the @ASN transformation to put the data onto the depth axis. Make sure there are the same number of depths as levels, for instance:

yes? define axis/z=/units=meters/depth zdepths = {2,5,10,15,20,30,50,100}
yes? let/units="`temp,return=units`" temp_on_depth = temp[gz=zdepths@ASN]

and note that on the DEFINE AXIS command, you could use /EDGES so that you define the edges of the grid cells instead of the grid cell centers. For unequally spaced grids, this is often the better choice.

-Ansley

On 8/9/2013 1:57 AM, mcaillau wrote:

Hi Ferret users,

I have data on a vertical level axis (level-grid) and i would like to regrid it to a regular z axis.
I have the thickness of each level (which vary in time).

I know that ZAXREPLACE would do it but i dont want my data to be interpolated on the z-axis. I am looking for the following transformation :
for each grid cell of my z-axis (so at a given depth) it will affect the value at the same depth from my level-grid.

Thank in advance for your help

Matthieu