Re: Listing data subsets. (visualizing box boundaries)

[Steve Hankin <hankin@xxxxxxxxxxxxx>]

    *** Visualizing the box boundaries on a time (or other) axis ***
=======

Hi Paulo,

Your expression appears to compute month by month averages from 1985 to
'95 using a regularly spaced destination time axis. Thus, the range of
time associated with each destination time point, extends from the
midpoint to midpoint forward and backward from each time point.

As stated in Chapter 4., 2.3.1 Regridding transformations

     @AVE--averaging
     Computes the length-weighted average of all points on the
     source grid that lie partly or completely within each grid
     cell of the destination grid.

To visually see the bounds of the time points from your source
(irregular) axis, the SHADE command is handy. A plot of

    LET my_time_series = time-series-at-a-single-point
    SHADE 0*my_time_series+0*i[i=1:2] + L

should reveal the boundaries that you are after quite well. (Use
/LEVELS=n if necessary to specify enough of color levels to see distinct
time cells)

    - steve

===================================

"Paulo B. Oliveira" wrote:

> Hi,
>
> I'm computing monthly climatologies from a
> set of unevenly spaced series of temperature
> records using:
> def ax/t=15-jan-1985:1-dec-1995:`365.25/12`/units=days t85_95m
> def grid/x=x/y=y/t=t85_95m mgrid
> def var mclim=SST[gt=mgrid@ave]
>
> Is there a way of listing the actual SST values used
> to compute the average for a given month ?
>
> Thanks!
>
> Paulo.

--
Steve Hankin
NOAA/PMEL, 7600 Sand Point Way NE, Seattle, WA 98115-0070
ph. (206) 526-6080 -- FAX (206) 526-6744