This component provides a mechanism to apply the constrained two-stream radiative transfer model described in Müller 2014. This correction handles both the loading and scattering effects on the filter.

Usage

da.corr.mueller2014 [switches…​] [[station] variables times [archive]|[file]]

Switches

--asymmetry-end=INPUT

This is the input asymmetry parameter at the end of each sampling interval. That is, the asymmetry parameter calculated from the scattering and back scattering integrated optical depths. When this option is set, this is used for all absorption channels. Therefor, normally only a single absorption can be corrected while using this option if there is any wavelength dependence.

Default: Resolved automatically by instrument suffix

--asymmetry-start=INPUT

This is the input asymmetry parameter at the start of each sampling interval. That is, the asymmetry parameter calculated from the scattering and back scattering integrated optical depths. When this option is set, this is used for all absorption channels. Therefor, normally only a single absorption can be corrected while using this option if there is any wavelength dependence.

Default: Resolved automatically by instrument suffix

--correct=SELECTION

These are the variables to that the correction is applied to. This option is mutually exclusive with instrument specification.

Default: All absorptions

--delta-af=INPUT

This is the absorption optical of an unloaded filter used in the model calculation.

Default: Wavelength dependent

--delta-sf=INPUT

This is the scattering optical of an unloaded filter used in the model calculation.

Default: Wavelength dependent

--deltaf-end=INPUT

This is the input end filter optical depth used in the correction. When this option is set, this value is used for all absorption channels. Therefor, normally only a single absorption can be corrected while using this option.

Default: Resolved automatically by instrument suffix

--deltaf-start=INPUT

This is the input start filter optical depth used in the correction. When this option is set, this value is used for all absorption channels. Therefor, normally only a single absorption can be corrected while using this option.

Default: Resolved automatically by instrument suffix

--deltas-end=INPUT

This is the input end scattering optical depth used in the correction. That is, this is the scattering integrated with the sample length at the start of each calculation interval.When this option is set, this is used for all absorption channels. Therefor, normally only a single absorption can be corrected while using this option if there is any wavelength dependence.

Default: Resolved automatically by instrument suffix

--deltas-start=INPUT

This is the input start scattering optical depth used in the correction. That is, this is the scattering integrated with the sample length at the start of each calculation interval.When this option is set, this is used for all absorption channels. Therefor, normally only a single absorption can be corrected while using this option if there is any wavelength dependence.

Default: Resolved automatically by instrument suffix

--dl=INPUT

This is the change in sample length used for each calculation interval. When omitted it is calculated from the start and end sample lengths.

--gf=INPUT

This is the average asymmetry parameter of the filter fibers used in the model calculation.

Default: 0.75

--hf=INPUT

This is the fractional thickness of the particle layer used in the two stream model calculation.

Default: 0.2

--instruments=SUFFIX…​

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

Default: All instrument suffixes

--length-end=INPUT

This is the input end sample length used in the correction. When this option is set, this length is used for all absorption channels.

Default: Resolved automatically by instrument suffix

--length-start=INPUT

This is the input start sample length used in the correction. When this option is set, this length is used for all absorption channels.

Default: Resolved automatically by instrument suffix

--transmittance-end=INPUT

This is the input end transmittance used in the correction. When this option is set, this transmittance is used for all absorption channels. Therefor, normally only a single absorption can be corrected while using this option.

Default: Resolved automatically by instrument suffix

--transmittance-start=INPUT

This is the input start transmittance used in the correction. When this option is set, this transmittance is used for all absorption channels. Therefor, normally only a single absorption can be corrected while using this option.

Default: Resolved automatically by instrument suffix

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 correct all absorptions present using pre-calculated scattering optical depths associated with each wavelength.

da.corr.mueller2014 bnd S11a 2015-05-01 2015-05-03

Or:

da.corr.mueller2014 input_file.c3d

Or:

da.get bnd S11a 2015-05-01 2015-05-03 | da.corr.mueller2014
Single instrument with a manual constant

This will correct all absorptions from A11 using an alternate setting for the absorption optical depth of an unloaded filter.

da.corr.mueller2014 --delta-af=0.02 --instruments=A11 bnd S11a 2015-05-01 2015-05-03