Re: list output

[Arne Biastoch <abiastoch@xxxxxxxx>]

Jonathan, Muyin

It was clear to me that I could format the output to display bad values
by another variable.

But what I want instead is to get ONLY the defined values in the
listing. In you example that would be:

> yes? let myvar = "do something fancy" 
> yes? list myvar

>  26.00 26.00 26.00 27.00 27.00 27.00 28.00 28.00 27.00 28.00 
>  28.10 27.53 27.00 27.71 27.75 29.00 27.60 27.50 27.00 27.80 
>  27.21 27.00 27.00 27.00 27.30 28.00 

Muyin Wang gave another hint:

> Try to use the external function
> compressi/j/k/l(dat) which returns data, compressed along the I/J/K/L 
> axis, and
> missing points moved to the end. But the problem is that your new data set
> i/j/k/l index may not be true (compared with the original).

That works well. I get the missing values at the end - but then it is
possible to get rid of these in an editor.

Thanks to both of you.

-Arne

Jonathan Callahan wrote:
> 
> Arne,
> 
> Here's an example using the COADS SST time series from CDC.  In it I
> define a new variable which resets the value of the bad flag and use the
> /FORMAT qualifier to the LIST command.
> 
> yes? use
> "http://www.cdc.noaa.gov/cgi-bin/nph-nc/Datasets/coads/1degree/global/enh/sst.mean.nc";
> yes? set region/x=180:190/y=0:10/t="15-jan-1961"
> yes? shade sst
> yes? let mysst = if sst then sst else 99.99
> yes? list/format=(10f6.2) mysst
>              IF SST THEN SST ELSE 99.99
>              LONGITUDE: 180E to 170W
>              LATITUDE: 0 to 10N
>              TIME: 01-JAN-1961 00:00
>              DATA SET:
> http://www.cdc.noaa.gov/cgi-bin/nph-nc/Datasets/coads/1degree/global/enh/sst.mean.nc
>  99.99 99.99 99.99 99.99 99.99 99.99 99.99 99.99 99.99 99.99
>  99.99 99.99 99.99 99.99 99.99 99.99 99.99 99.99 99.99 26.00
>  99.99 99.99 99.99 99.99 99.99 99.99 99.99 26.00 26.00 99.99
>  99.99 99.99 99.99 99.99 99.99 99.99 27.00 99.99 99.99 99.99
>  99.99 99.99 99.99 27.00 99.99 27.00 99.99 99.99 28.00 28.00
>  99.99 99.99 27.00 99.99 99.99 99.99 99.99 28.00 99.99 99.99
>  28.10 27.53 27.00 27.71 27.75 29.00 27.60 99.99 99.99 27.50
>  99.99 99.99 99.99 27.00 27.80 99.99 99.99 27.21 99.99 27.00
>  27.00 27.00 27.30 99.99 99.99 99.99 99.99 99.99 99.99 99.99
>  28.00 99.99 99.99 99.99 99.99 99.99 99.99 99.99 99.99 99.99
> yes?
> 
> Note: low latitudes appear at the bottom of the image but at the top of
> the listing.
> 
> I could have used the original SST variable and
> 
> yes? list/format=(10g) sst
> 
> to get the bad flags to appear as their floating point values instead of
> "***".
> 
> -- Jonathan Callahan
> 
> Arne Biastoch wrote:
> >
> > Dear Ferreters,
> >
> > is there any possibility to suppress output that contains missing
> > values? I want to subsample a field and do a
> >
> > let newvar = if oldvar gt xx then oldvar
> >
> > Now newvar conatins a lot of missing values that show up as stars when
> > listed to an ASCII file. It is possible to prevent FERRET from listing
> > these values?
> >
> > OR: How should the let command be defined to ignore those values and
> > enable the newvar smaller in dimension that oldvar?
> >
> > -Arne
> >
> > --
> >
> > Dr. Arne Biastoch
> > Scripps Institution of Oceanography        phone: +1-858-822-3787
> > Physical Oceanography Research Division    fax  : +1-858-534-9820
> > MS: 0230                                   email: abiastoch@ucsd.edu
> > 8605 La Jolla Shores Dr.
> > La Jolla, CA 92093-0230, U.S.A.   http://www.ecco.ucsd.edu/~biastoch

-- 

Dr. Arne Biastoch
Scripps Institution of Oceanography        phone: +1-858-822-3787
Physical Oceanography Research Division    fax  : +1-858-534-9820
MS: 0230                                   email: abiastoch@ucsd.edu
8605 La Jolla Shores Dr.            
La Jolla, CA 92093-0230, U.S.A.   http://www.ecco.ucsd.edu/~biastoch