Comparison with Regal Black of A11 (Pall) and A13 (Azumi) measurements on black runs

Creating data files
[aero:/aer/sfa/work/GMD_lab_absorption_study_2015/runs]$ da.get sfa '.*_A1[13]' `query 11` > run11_A11_A13.cpd3
[aero:/aer/sfa/work/GMD_lab_absorption_study_2015/runs]$ da.get sfa '.*_A1[13]' `query 12` > run12_A11_A13.cpd3
[aero:/aer/sfa/work/GMD_lab_absorption_study_2015/runs]$ da.get sfa '.*_A1[13]' `query 22` > run22_A11_A13.cpd3

# create pseudo-instruments with different corrections applied
# A11 = no corrections
# A81 = Weiss correction (same for A11 and A13)
# A91 = Irwin corrections (different for A11 and A13)
cat run11_A11_A13.cpd3 \
|da.corr.weiss --instruments=A11,A13 --constants=0.866,1.317 \
|da.select --capture-variable='(.*)_A1(.)' --set-apply='.*_A(.+)' --set-variable='${1}_A8${2}' \
> run11_A81_A83.cpd3

cat run11_A11_A13.cpd3 \
|da.corr.weiss --instruments=A11 --constants=1.2107,1.4965,-0.52417 \
|da.corr.weiss --instruments=A13 --constants=1.2661,1.5254,-0.60845 \
|da.select --capture-variable='(.*)_A1(.)' --set-apply='.*_A(.+)' --set-variable='${1}_A9${2}' \
> run11_A91_A93.cpd3

cat run12_A11_A13.cpd3 \
|da.corr.weiss --instruments=A11,A13 --constants=0.866,1.317 \
|da.select --capture-variable='(.*)_A1(.)' --set-apply='.*_A(.+)' --set-variable='${1}_A8${2}' \
> run12_A81_A83.cpd3

cat run12_A11_A13.cpd3 \
|da.corr.weiss --instruments=A11 --constants=1.2107,1.4965,-0.52417 \
|da.corr.weiss --instruments=A13 --constants=1.2661,1.5254,-0.60845 \
|da.select --capture-variable='(.*)_A1(.)' --set-apply='.*_A(.+)' --set-variable='${1}_A9${2}' \
> run12_A91_A93.cpd3

cat run22_A11_A13.cpd3 \
|da.corr.weiss --instruments=A11,A13 --constants=0.866,1.317 \
|da.select --capture-variable='(.*)_A1(.)' --set-apply='.*_A(.+)' --set-variable='${1}_A8${2}' \
> run22_A81_A83.cpd3

cat run22_A11_A13.cpd3 \
|da.corr.weiss --instruments=A11 --constants=1.2107,1.4965,-0.52417 \
|da.corr.weiss --instruments=A13 --constants=1.2661,1.5254,-0.60845 \
|da.select --capture-variable='(.*)_A1(.)' --set-apply='.*_A(.+)' --set-variable='${1}_A9${2}' \
> run22_A91_A93.cpd3

da.multiplex run11_A*.cpd3 > run11_all.cpd3
da.multiplex run12_A*.cpd3 > run12_all.cpd3
da.multiplex run22_A*.cpd3 > run22_all.cpd3

#Plotting raw data
da.show.timeseries --y='BaG.*' run11_A11_A13.cpd3 
da.show.timeseries --y='BaG.*' run12_A11_A13.cpd3 
da.show.timeseries --y='BaG.*' run22_A11_A13.cpd3 
for I in 11 12 22; do da.plot.timeseries --y='(Ba|Ir)G.*' --legend=topright --title="CLAP Filter Comparison with Regal Black - no transmittance correction, run $I" --file=run${I}_green.png run${I}_A11_A13.cpd3 ;done

#Compare raw data
# generate density plot with no transmittance correction applied
da.multiplex run*_A11_A13.cpd3 |da.show.density --x='BaG_A11' --y='BaG_A13' --quantile=95 --fit=linearzero --grid=both --x-title="Pallflex E70-2075W (1/Mm)" --y-title="Azumi 371M (1/Mm)" --title="CLAP Filter Comparison with Regal Black - no transmittance correction"
da.multiplex run*_A11_A13.cpd3 |da.show.density --x='BaG_A11' --y='BaG_A13' --quantile=95 --fit=linear --grid=both --x-title="Pallflex E70-2075W (1/Mm)" --y-title="Azumi 371M (1/Mm)" --title="CLAP Filter Comparison with Regal Black - no transmittance correction"



# generate density plot with no transmittance correction applied, green
da.multiplex run*_A11_A13.cpd3 \
|da.plot.density --x='BaG_A11' --y='BaG_A13' --quantile=95 --fit=linearzero --grid=both --x-title="Pallflex E70-2075W (1/Mm)" --y-title="Azumi 371M (1/Mm)" \
--title="CLAP Filter Comparison with Regal Black - no transmittance correction, green wavelength" \
--file=black_filter_compare_no_Tr_corr_green.png

# generate density plot with Bond transmittance correction applied (PSAP/COSMOS default correction), green
da.multiplex run*_A11_A13.cpd3 \
|da.corr.weiss --instruments=A11,A13 --constants=0.866,1.317 \
|da.plot.density --x='BaG_A11' --y='BaG_A13' --quantile=95 --fit=linearzero --grid=both --x-title="Pallflex E70-2075W (1/Mm)" --y-title="Azumi 371M (1/Mm)" \
--title="CLAP Filter Comparison with Regal Black - Bond transmittance correction, green wavelength" \
--file=black_filter_compare_Bond_Tr_corr_green.png


# generate density plot with Irwin transmittance corrections applied, using different correction constants for Pallflex and Azumi filters, green
da.multiplex run*_A11_A13.cpd3 \
|da.corr.weiss --instruments=A11 --constants=1.2107,1.4965,-0.52417 \
|da.corr.weiss --instruments=A13 --constants=1.2661,1.5254,-0.60845 \
|da.plot.density --x='BaG_A11' --y='BaG_A13' --quantile=95 --fit=linearzero --grid=both --x-title="Pallflex E70-2075W (1/Mm)" --y-title="Azumi 371M (1/Mm)" \
--title="CLAP Filter Comparison with Regal Black - Irwin transmittance corrections, green wavelength" \
--file=black_filter_compare_Irwin_Tr_corr_green.png

# generate density plot with no transmittance correction applied, blue
da.multiplex run*_A11_A13.cpd3 \
|da.plot.density --x='BaB_A11' --y='BaB_A13' --quantile=95 --fit=linearzero --grid=both --x-title="Pallflex E70-2075W (1/Mm)" --y-title="Azumi 371M (1/Mm)" \
--title="CLAP Filter Comparison with Regal Black - no transmittance correction, blue wavelength" \
--file=black_filter_compare_no_Tr_corr_blue.png

# generate density plot with Bond transmittance correction applied (PSAP/COSMOS default correction), blue
da.multiplex run*_A11_A13.cpd3 \
|da.corr.weiss --instruments=A11,A13 --constants=0.866,1.317 \
|da.plot.density --x='BaB_A11' --y='BaB_A13' --quantile=95 --fit=linearzero --grid=both --x-title="Pallflex E70-2075W (1/Mm)" --y-title="Azumi 371M (1/Mm)" \
--title="CLAP Filter Comparison with Regal Black - Bond transmittance correction, blue wavelength" \
--file=black_filter_compare_Bond_Tr_corr_blue.png

# generate density plot with Irwin transmittance corrections applied, using different correction constants for Pallflex and Azumi filters, blue
da.multiplex run*_A11_A13.cpd3 \
|da.corr.weiss --instruments=A11 --constants=1.2107,1.4965,-0.52417 \
|da.corr.weiss --instruments=A13 --constants=1.2661,1.5254,-0.60845 \
|da.plot.density --x='BaB_A11' --y='BaB_A13' --quantile=95 --fit=linearzero --grid=both --x-title="Pallflex E70-2075W (1/Mm)" --y-title="Azumi 371M (1/Mm)" \
--title="CLAP Filter Comparison with Regal Black - Irwin transmittance corrections, blue wavelength" \
--file=black_filter_compare_Irwin_Tr_corr_blue.png

# generate density plot with no transmittance correction applied, red
da.multiplex run*_A11_A13.cpd3 \
|da.plot.density --x='BaR_A11' --y='BaR_A13' --quantile=95 --fit=linearzero --grid=both --x-title="Pallflex E70-2075W (1/Mm)" --y-title="Azumi 371M (1/Mm)" \
--title="CLAP Filter Comparison with Regal Black - no transmittance correction, red wavelength" \
--file=black_filter_compare_no_Tr_corr_red.png

# generate density plot with Bond transmittance correction applied (PSAP/COSMOS default correction), red
da.multiplex run*_A11_A13.cpd3 \
|da.corr.weiss --instruments=A11,A13 --constants=0.866,1.317 \
|da.plot.density --x='BaR_A11' --y='BaR_A13' --quantile=95 --fit=linearzero --grid=both --x-title="Pallflex E70-2075W (1/Mm)" --y-title="Azumi 371M (1/Mm)" \
--title="CLAP Filter Comparison with Regal Black - Bond transmittance correction, red wavelength" \
--file=black_filter_compare_Bond_Tr_corr_red.png

# generate density plot with Irwin transmittance corrections applied, using different correction constants for Pallflex and Azumi filters, red
da.multiplex run*_A11_A13.cpd3 \
|da.corr.weiss --instruments=A11 --constants=1.2107,1.4965,-0.52417 \
|da.corr.weiss --instruments=A13 --constants=1.2661,1.5254,-0.60845 \
|da.plot.density --x='BaR_A11' --y='BaR_A13' --quantile=95 --fit=linearzero --grid=both --x-title="Pallflex E70-2075W (1/Mm)" --y-title="Azumi 371M (1/Mm)" \
--title="CLAP Filter Comparison with Regal Black - Irwin transmittance corrections, red wavelength" \
--file=black_filter_compare_Irwin_Tr_corr_red.png