data.process.psap_int

data.process.psap_int takes high frequency PSAP/CLAP intensities or transmittance and uses that to calculate a longer time base average.

Command Line Usage

data.process.psap_int [--cpd2] [--cpd2out] [--instrument=A11] 
                      [--wavelengths=...] intensities flow[;offset;slope] 
                      time area average [inputFile] 

Arguments

intensities

A series of comma separated lists of the form 'Sample;SampleDark;Reference;ReferenceDark'. Each list represents an absorption that will be calculated. The dark specifiers are optional and if omitted default to zero (no offset). Each individual component is either a field index (starting from one) in CPD1 mode (the default) or a variable name base in CPD2 mode. If the reference is omitted it is assumed to be one; this allows just a transmittance to be used.

flow[;offset;slope]

A flow rate and an optional slope and offset calibration. If the flow rate begins with a ”+” it is assumed be to a constant value. Otherwise if running in CPD1 mode (the default) it is a field index (starting from one), if in CPD2 mode it is a variable name base.

time

An instrument time field specifier. In CPD1 mode this is a field index (starting from one). In CPD2 mode it is a variable name base. If the argument is blank (””) or invalid the record time is used.

area

The filter spot area in m^2.

inputFile

A file to read from. If omitted or '-' then standard input is used.

--cpd2

Use CPD2 input mode. A field specifiers are treated as variable name bases.

--cpd2out

Output CPD2 style data instead of CPD1.

--instrument=A11

Set the default instrument for CPD2 output, defaulting to null.

--wavelengths=...

Comma separated list of wavelength strings to use in the CPD2 output.

Example Usage

Average PSAP3W Raw Data to 60 Seconds

cnvt_psap3w --record=m lr_X.sfa | data.process.psap_int --cpd2 'IpB;;IfB,IpG;;IfG,IpR;;IfR' 'Q;0.0;1.0' 'Sp' '2.248E-05' 60

Average CPD1 Logged CLAP3W Data Raw to Four Seconds

data.process.psap_int '33;32;13;12,34;32;14;12,35;32;15;12' '9;0.0;1.0' 6 2.248E-05 4

Note the output is red, green, then blue in this example.

Three ways to calculate hourly averages for PSAP-3W data

Calculating absorptions using data.avg and backing out the Weiss correction:

data.get cpr A12a 2009:40 2009:50 | data.edit.corr.psap_weiss_undo | \
  data.avg --interval=1h --contam --nocut --nostddev > instrument.cpd2

Calculating absorptions using instrument transmittances:

data.get cpr A12a,A12m 2009:40 2009:50 | data.process.psap_int --cpd2 --cpd2out \
  --instrument=A13 IrB,IrG,IrR Q '' 1.93E-5 3600 > transmittance.cpd2

Calculating absorptions using instrument intensities:

data.get cpr A12m 2009:40 2009:50 | data.process.psap_int --cpd2 --cpd2out \
  --instrument=A14 'IpB;;IfB,IpG;;IfG,IpR;;IfR' Q '' \
  1.93E-5 3600 > intensities.cpd2

Combining all three files:

data.multiplex instrument.cpd2 transmittance.cpd2 intensities.cpd2 | \
  data.consolidate --source=- 'Ba*'