[Thread Prev][Thread Next][Index]
Re: [ferret_users] POP vertical section
! 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>
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Privacy Policy | Disclaimer | Accessibility Statement