[Thread Prev][Thread Next][Index]

Re: [ferret_users] bound variable detected as variable

Hi -
Thanks for the report. I don't see any reason that what Patrick wants to do couldn't be implemented, so we'll look into doing that. As Marco says an axis with more than one point does save with the correct bounds as expected.


Marco Steinacher wrote:
Hi Patrick,

It seems that ferret does not use the bounds variable because your depth
axis has only one point. When I add another one (eg. at z=15 with bounds
at z=10 and z=20) the bounds are used correctly.

With only one grid cell Ferret apparently assumes a regular grid with
boundaries at z=4.5 and z=5.5 (probably in order to have weight=1). When
you save your data Ferret is using these bounds values. You can check
which boundaries Ferret is using with

 list zboxlo[g=SOSALINE_GLOBAL],zboxhi[g=SOSALINE_GLOBAL]

As a workaround you could extend the z-axis with

define axis/z/bounds/unit=meters `SOSALINE_GLOBAL,r=zaxis` = \

Then, the boundaries are treated as you are expecting.


Brockmann Patrick wrote:
Hi all,

I have noticed that when a boundarie variable is present in a netcdf file,
ferret detects this "associate" coordinate variable as a "classic"

Consequence in my scripts, the use of
yes? shade `..varnames[i=1]`
to produce a shade with the first variable found in a file becomes not
since ferret takes the boundaries variable.

I have also saved a file with /bounds qualifier
and noticed that boundaries have been changed

Initially I have:
DEPTH1_1 = 5 ;
DEPTH1_1_bnds =
 -0.000244140625, 10.000244140625 ;

After a
yes? save/clobber/bounds/file=new.nc SOSALINE_GLOBAL

I get a
DEPTH1_1 = 5 ;
DEPTH1_1_bnds =
 4.5, 5.5 ;

I do not understand why ferret has changed those values.
Note that when I use /bounds qualifier, I also save bounds
even when axis is regular. Is there a way to control the write
of boundarie variables for each axis (ie qualifiers XBOUNDS, YBOUNDS,

and try

Tests made with ferret v620 (Linux 32bits)

Thanks gor any help on this topic.



[Thread Prev][Thread Next][Index]

Contact Us
Dept of Commerce / NOAA / OAR / PMEL / TMAP

Privacy Policy | Disclaimer | Accessibility Statement