Skip to main content
U.S. FlagAn official website of the United States government
icon dot gov
Official websites use .gov

A .gov website belongs to an official government organization in the United States.

icon https
Secure .gov websites use HTTPS

A lock () or https:// means you’ve safely connected to the .gov website. Share sensitive information only on official, secure websites.

Select data for a particular time period from an input file

da.select --include=.*,2012:1,2012:2 inputfile.cpd3

Select a particular variable from a data stream

... | da.select --include=X1c_A81 |...

Remove size-cut distinction from variables in data stream

... | da.select --include=':::cover:!' --set-flavor-remove=cover | ...
--alternatively--
... |da.select --set-flavor-remove=pm1 | da.select --set-flavor-remove=pm10 | ...

Remove neph zero data to enable calculating the average for an entire period. Normally, da.avg will break the data into hourly averages because the metadata for the zeros specifies that they shouldn't be averaged (normally correct, you don't want zeros being folded into hourly averages if you put the two together, since the zero specify “segments” rather than values). This is, unfortunately, hard to do because it was designed to “do the right thing” by default, so it ends up being complicated when you want to override the intelligence (also applies to the cut sizes above). However, da.select can do this:

... | da.select --path="^Smoothing" --path-apply=:raw_meta: --explicit-meta=true --path-remove=true | da.avg --interval=undef | ...

What is the percentage of hours that the nephelometer RH was above 40% at EGB and is there any any seasonal dependence? (Change the –interval=1mo to –interval=undef to get a single coverage number for the whole time.)

data.get egb S11a 2010 2014 avgh | da.select --include=U_S11 | da.where --include='U_S11 > 40' | da.avg --interval=1mo | da.select --include=':::cover:!' --set-flavor-remove=cover | da.export

Is there any annual cycle in the % of hours for which Bsp(PM1) < 1.0 Mm-1 at EGB?

data.get egb S11a 2010 2014 avgh | da.select --include=::BsG_S11:pm1 | da.where --include='BsG_S11 < 1' | da.avg --interval=1mo | da.select --include=':::cover:!' --set-flavor-remove=cover | da.export