This calculates intensive optical intensive properties from its input. This include name normalization, single scattering albedos, Ångström exponents, asymmetry parameters, and radiative forcing efficiency.

Usage

da.calc.intensives [switches…​] [[station] variables times [archive]|[file]]

Switches

--absorption=SELECTION

This is the input light absorption to use in intensives generation. This option is mutually exclusive with absorption instrument specification.

Default: All light absorptions

--absorption-instruments=SUFFIX…​

These are the instrument suffixes to use absorption data from. For example A11 would usually specifies the primary absorption instrument. This option is mutually exclusive with manual absorption variable specification.

Default: All instrument suffixes

--angstrom=INPUT

This is the Ångström exponent used to perform wavelength adjustment when interpolation is not possible. When not set, data are simply left as undefined when they cannot be adjusted to the target wavelengths.

--angstrom-distance=NUMBERS

This is the maximum distance that the fallback Ågström exponent is valid for. This is used in conjunction with angstrom-wavelength to limit how far data are extrapolated.

This option only accepts numbers greater than or equal to 0. Undefined values are permitted, to indicate that the option is set but not to a specific value.

--angstrom-wavelength=NUMBERS

This is the wavelength of the Ångström exponent used to perform wavelength adjustment when interpolation is not possible. This is used to limit how far data are extrapolated when used in conjunction with angstrom-distance.

This option only accepts numbers strictly greater than 0. Undefined values are permitted, to indicate that the option is set but not to a specific value.

--backscattering=SELECTION

This is the input light backscattering to use in intensives generation. This option is mutually exclusive with scattering instrument specification.

Default: All light backscatterings

--bfr=INPUT

When set, this specifies a backscatter fraction to use. This allows other parameters to be calculated when there is no backscattering data available by assuming a fixed backscatter ratio.

Default: Calculated from the scattering data

--counts=INPUT

This is the input concentration to forward and rename in the output. This option is mutually exclusive with concentration instrument specification.

Default: All concentrations

--counts-instruments=SUFFIX…​

These are the instrument suffixes to use particle concentrations from. For example N61 or N71 usually specifies the system CPC. This option is mutually exclusive with manual counts variable specification.

Default: All instrument suffixes

--extinction=SELECTION

This is the input light extinction to use in intensives generation. This option is mutually exclusive with extinction instrument specification.

Default: All light extinctions

--extinction-instruments=SUFFIX…​

These are the instrument suffixes to use extinction data from. For example E11 would usually specifies the first extinction instrument. This option is mutually exclusive with manual extinction variable specification.

Default: All instrument suffixes

--mac=INPUT

This is the calibration factor used to calculate equivalent black carbon concentrations. The efficiency for each wavelength is this value divided by the wavelength in nm.

Default: 14625

--output=SUFFIX…​

These are the instrument suffixes that output data will be generated with first. If the number of outputs exceeds the number of suffixes defined then the defaults are used.

Default: XI, XJ, etc

--scattering=SELECTION

This is the input light scattering to use in intensives generation. This option is mutually exclusive with scattering instrument specification.

Default: All light scatterings

--scattering-instruments=SUFFIX…​

These are the instrument suffixes to use scattering data from. For example S11 would usually specifies the reference nephelometer. This option is mutually exclusive with manual scattering variable specification.

Default: All instrument suffixes

--ssa=INPUT

When set, this specifies a single scattering albedo to use. This allows other parameters to be calculated assuming a fixed value.

Default: Calculated from data

--wavelengths=NUMBER…​

These are the output wavelengths to calculate the intensive parameters at. Input data will be adjusted to these wavelengths before use.

Default: 450,550,700

This option only accepts numbers strictly greater than 0.

Arguments

If no bare word input specification is supplied then data are read from standard input.

station

This argument is used to specify the station used to look up variables that do not include a station as part of an archive read specification. The station is the three letter GAW station code of the location, such as BND. The argument accepts multiple stations specified as regular expressions and separated by : or ; or ,. For example BND,MLO and (BND|MLO) are two ways of selecting both the BND and MLO stations.

The special value allstations may also be specified to select all stations.

variables

This argument may be split into multiple actual program arguments. Each part consists of one or more variable specifications separated by , (commas). The station and archive used if none are explicitly given are defined in the other arguments or inferred from system defaults. For example simply BsB_S11 selects the blue scattering from the S11 instrument and returns all cut sizes and metadata. If instead ::BsB_S11:pm10 is used then only PM10 scattering is returned. That is, the selection specifies "any" station and archive (the defaults are set in the other arguments) then requires that it has the "PM10" flavor. A specification such as bnd:BsB_S11,thd:BsB_S11 allows for data selection from multiple stations.

This may also be used to select a "record" type alias such as "S11a". In this mode all the variables defined in the alias are returned.

Finally the special value everything can be used to select all available data for the given station(s) and archive(s).

times

This argument may be split into multiple actual program arguments. This full list of time arguments defines a time bounds list that sets the range of data queried. The final resulting data are clipped to the given time bounds but any values that intersect them are returned.

archive

This argument is used to specify the archive used to look up variables that do not include an archive as part of an archive read specification. The archive is the internal name, such as raw or clean_meta. The argument accepts multiple archives specified as regular expressions and separated by : or ; or ,. For example raw(_meta)? selects both the raw and raw_meta archives.

The special value allarchives may also be used to select all available archives.

file

This argument is used to specify the the file to read data from. If it is present and exists then data is read from the given file name instead of from standard input. Alternatively - (a single dash) may be used to explicitly specify standard input.

Examples

Defaults

This will generate intensives for all input data. Note that the output variable assignment is un-ordered in the case of multiple sources of extensive parameters,

da.calc.intensives bnd S11a 2015-05-01 2015-05-03

Or:

da.calc.intensives input_file.c3d

Or:

da.get bnd S11a 2015-05-01 2015-05-03 | da.calc.intensives
Restricted instruments

This will calculate intensives from only S12, A12, and N61 regardless of what data is present in the input. Output variables will have the XJ suffix (for example BsG_XJ).

da.calc.intensives --absorption-instruments=A12 --counts-instruments=N61 --extinction-instruments= --output=XJ --scattering-instruments=S12 bnd S11a 2015-05-01 2015-05-03
Non-standard wavelengths and manual backscatter fraction

This will calculate intensives at 467, 528, and 652nm using a fixed intensives fraction of 0.2 instead of calculating it from the data.

da.calc.intensives --bfr=0.2 --wavelengths=467,528,652 bnd S11a 2015-05-01 2015-05-03
Explicit scattering

This will calculate intensives using only green scattering from the S11 instrument.

da.calc.intensives --backscattering=bnd:raw:BbsG_S11 --scattering=bnd:raw:BsG_S11 bnd S11a 2015-05-01 2015-05-03