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