[Thread Prev][Thread Next][Index]
Re: DEFINE AXIS/EDGES, @AVE
Hi Andrew,
This looks to be a subtle bug. It does not have to do with the
/EDGES definition of the axis, per se. If in your example you do
not plot sst, but only plot sst[gt=t2@ave], or just the two averages
together, the average is computed correctly.
When Ferret plots the series sst and sst[gt=t2@ave] together, it uses
an incorrect range of indices for computing the average. It happens
only for some cases where the last point on the axis you are regridding
to is just above a point on the original sst axis. In any event, we have
a fix for this which will be in the next release of Ferret.
Ansley Manke
Andrew Wittenberg wrote:
> Hello Ferreters,
>
> I'm running into a problem with axis definitions and @AVE regridding.
> Can anyone tell me why the following gives two different results?
>
> use coads_climatology
> def axis/t=731:6575:1461/t0="`sst,r=t0`"/unit="`sst,r=tunit`" t1
> def axis/edges/t=1:7306:1461/t0="`sst,r=t0`"/unit="`sst,r=tunit`" t2
> plot/x=140w/y=0/sym/line sst,sst[gt=t1@ave],sst[gt=t2@ave]
>
> Axes t1 & t2 should be nearly identical, but when I regrid using @AVE the
> last point seems to be wrong for axis t2.
>
> Thanks for any insights!
>
> Andrew
>
> ----------------------------------------------------------------------
> Andrew Wittenberg Program in Atmospheric and Oceanic Sciences
> andyw@splash.princeton.edu Princeton University
> ----------------------------------------------------------------------
--
Ansley Manke Pacific Marine Environmental Laboratory Seattle WA (206)526-6246
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement