Re: [ferret_users] POP vertical section

["William S. Kessler" <william.s.kessler@xxxxxxxx>]

! Better to define an intermediate variable with your axis points,  
then use that to define the axis:

let depth_in_meters={2.5,.....,6100.}

! Then define the axis using the FROM_DATA qualifier, and specify  
that it is depth (positive down), and in meters (better labels):

def axis/z/from_data/name=zmeters/depth/units=m depth_in_meters

! Then do not redefine the entire grid (unnecessary) unless you need  
to change the x and y specifications also:

let temp_n=temp[gz=zmeters]

Billy K

On 6 Mar 10, at 10:12 PM, Osinski, Robert FORNATL, PL wrote:

> Hi Ferreters,
> I have a problem with creating vertical sections with POP (Parallel  
> Ocean Program) model data. The ocean model is z-coordinates so  
> levels are numbered 1, 2.... and my vertical section is on Z axis  
> is in number of levels not in meters. I've tried to define the new  
> grid:
>
> use hist.Rgrid12.19790320.nc ! my model output file with "temp"  
> variable
>
> define axis/z my_zaxis= 
> {2.5,7.5,12.5,17.5,23.,29.6,37.7,47.4,59.,73.1,90.1,110.6,135.3,165.2, 
> 201.1,244.5,296.9,360.1,463.2,528.1,639.2,775.,950.,1150.,1350.,1550., 
> 1750.,1950.,2150.,2350.,2550.,2750.,2950.,3150.,3375.,3620.,3875.,4120 
> .5,4375.,4625.,4900.,5200.,5500.,5800.,6100.}
> define axis/x=1:1280:1 my_xaxis
> define axis/y=1:720:1 my_yaxis
> define grid/x=my_xaxis/y=my_yaxis/z=my_zaxis my_grid
>
> let temp_n=temp[d=1,g=my_grid]
> fill/x=1:100/y=10/vlimits=1:100:1 temp_n
>
> but the results is still in levels not meters (figure is attached).
> Any hints, pleas. I struggle with this problem for weeks and have  
> no answer.
> Thanks in advance,
> Dr. Robert Osinski
> NPS, Monterey, CA
>
> <vertical_POP.tiff>