{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "##
Notebook S06. Size Distribution Inversion of Ambient Data
\n", "\n", " Notebook Version 2: This notebook has been updated to reflect changes made to the package DifferentialMobilityAnalyzers.jl to work with the Julia v1 series (tested with Julia 1.1.0). To read the original supplement published with the paper please switch to v1.0.0 of the package DifferentialMobilityAnalyzers.jl and/or download the virtual machine on zenodo.org which contains a complete installation that works with Julia 0.6.4 \n", "\n", "This notebook demonstrates how to use JDL to invert size distribution data obtained from the TSI Aerosol Instrument Manager (AIM) software. The notebook also demonstrates that the software produces comparable results to the TSI software. This notebook is a supplement to the manuscript
\n", "\n", "Petters, M. D. (2018) A language to simplify computation of differential mobility analyzer response functions, Aerosol Science & Technology.
\n", "\n", "## Block 1. Dependencies" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "data": { "text/html": [ " \n" ], "text/plain": [ "HTML{String}(\" \\n\")" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "" ], "text/plain": [ "HTML{String}(\"\")" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "" ], "text/plain": [ "HTML{String}(\"\")" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ " \n" ], "text/plain": [ "HTML{String}(\" \\n\")" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "using Plots, Plots.PlotMeasures, LambertW, ProgressMeter, DifferentialMobilityAnalyzers\n", "plotlyjs(); " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Block 2. Setup the DMA \n", "Sheath flow = 4 L min-1 and sample aerosol flow = 1 L min-1 and specified in the data file loaded below." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "t,p = 295.15, 1e5 # Temperature [K], Pressure [Pa]\n", "qsa,β = 1.66e-5, 1/4 # Qsample [m3 s-1], Sample-to-sheath ratio,\n", "r₁,r₂,l = 9.37e-3,1.961e-2,0.44369 # DMA geometry [m]\n", "leff = 0.0 # DMA effective diffusion length [m]\n", "m = 6 # Upper number of charges\n", "Λ = DMAconfig(t,p,qsa,qsa/β,r₁,r₂,l,leff,:-,m,:cylindrical) # Specify DMA with negative polarity\n", "tscan, tc = 120,1 # SMPS scan time, bin average time\n", "v1,v2 = 10,10000 # Voltage range\n", "z1,z2 = vtoz(Λ,v2), vtoz(Λ,v1)\n", "δ = setupSMPS(Λ, v1, v2, tscan, tc);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Block 3. Response function\n", " loadtsidata Loads txt file exported from TSI Aerosol Instrument Manager (20130618_1410_SMPS3_BP.txt). The file was edited to remove a unicode character in line 154 that results in a load error. The load function will need to be adapted for files with different scan length and export parameters. The data corresponds to ambient data collected by Dr. Grieshop and Saha during the SOAS campaign (alse see Acknowledgements). It contains a total of 160 spectra. The instrument was operated at 4 L min-1 sheath flow and 1 L min-1 sample flow with a TSI 3081 long column DMA and TSI 3010 CPC. The data were exported with multiple charge correction turned on and diffusion correction turned off.
\n", "𝕟: Array of TSI inverted size distribution
\n", "rawdp: raw diameter grid obtained at 10Hz acquisition
\n", "rawc: 2D array containing raw counts vs diameter for each spectrum
\n", "Nt_TSI: array of integrated number concentration computed by TSI software" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "𝕟, rawdp, rawc, Nt_TSI = loadtsidata();" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The response function is computed from the 10 Hz raw count data. The raw counts falling within a bin are summed together. Concentration is computed from the number of milliseconds in each bin. The coincidence corrected concentration is obtained using Eq. (11) in Collins et al. (2003)

\n", "\n", "
$N_a = \\frac{-W(-N_m Q_{cpc}\\tau)}{Q_{cpc}\\tau}$
\n", "\n", "where $N_a$ is the coincidence corrected concentration, $N_m$ is the measured raw concentration, $Q_{cpc}$ is the volumetric CPC flow rate, $\\tau$ is the beam transit time, and $W$ is the zero branch of the LambertW function. Here $\\tau = 0.4\\mu s$, where the value is for the TSI 3010 CPC and obtained from the manufacturer manual." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "application/vnd.plotly.v1+json": { "data": [ { "colorbar": { "title": "" }, "hoverinfo": "text", "legendgroup": "y1", "line": { "color": "rgba(0, 0, 0, 1.000)", "dash": "solid", "shape": "linear", "width": 1 }, "mode": "lines", "name": "y1", "showlegend": true, "text": [ null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ], "type": "scatter", "x": [ 641.2116606343327, 641.2116606343327, 611.4599387459494, 611.4599387459494, 583.306474390341, 583.306474390341, 556.6613542884901, 556.6613542884901, 531.4394703614233, 531.4394703614233, 507.56026030320544, 507.56026030320544, 484.94747088813676, 484.94747088813676, 463.52894137710354, 463.52894137710354, 443.2364043391983, 443.2364043391983, 424.00530127516816, 424.00530127516816, 405.7746106002292, 405.7746106002292, 388.48668579141435, 388.48668579141435, 372.08710180294423, 372.08710180294423, 356.5245081767597, 356.5245081767597, 341.75048760168943, 341.75048760168943, 327.7194189851799, 327.7194189851799, 314.38834438227576, 314.38834438227576, 301.71683936864554, 301.71683936864554, 289.6668866433689, 289.6668866433689, 278.20275280215174, 278.20275280215174, 267.2908683347188, 267.2908683347188, 256.899710975487, 256.899710975487, 246.99969257957602, 246.99969257957602, 237.56304971266977, 237.56304971266977, 228.56373813904784, 228.56373813904784, 219.97733137279684, 219.97733137279684, 211.78092342763756, 211.78092342763756, 203.95303586511605, 203.95303586511605, 196.4735292024492, 196.4735292024492, 189.32351870267976, 189.32351870267976, 182.48529453289518, 182.48529453289518, 175.94224624241147, 175.94224624241147, 169.67879148283004, 169.67879148283004, 163.68030886617655, 163.68030886617655, 157.93307483604227, 157.93307483604227, 152.4242044096698, 152.4242044096698, 147.14159563600944, 147.14159563600944, 142.07387760555852, 142.07387760555852, 137.21036184188898, 137.21036184188898, 132.5409969017305, 132.5409969017305, 128.0563260098899, 128.0563260098899, 123.74744755674116, 123.74744755674116, 119.60597828913852, 119.60597828913852, 115.62401903004877, 115.62401903004877, 111.79412276766202, 111.79412276766202, 108.10926496096282, 108.10926496096282, 104.56281591549961, 104.56281591549961, 101.14851509018342, 101.14851509018342, 97.86044720323027, 97.86044720323027, 94.69302001269416, 94.69302001269416, 91.64094365432143, 91.64094365432143, 88.69921142660979, 88.69921142660979, 85.86308191991205, 85.86308191991205, 83.12806239313848, 83.12806239313848, 80.48989330804814, 80.48989330804814, 77.94453393726194, 77.94453393726194, 75.48814896795007, 75.48814896795007, 73.11709602865756, 73.11709602865756, 70.82791407191412, 70.82791407191412, 68.61731255014551, 68.61731255014551, 66.48216132696216, 66.48216132696216, 64.41948127016396, 64.41948127016396, 62.426435476770884, 62.426435476770884, 60.50032108408952, 60.50032108408952, 58.638561624262884, 58.638561624262884, 56.838699882942535, 56.838699882942535, 55.098391225681986, 55.098391225681986, 53.415397358389946, 53.415397358389946, 51.78758049072226, 51.78758049072226, 50.212897873635384, 50.212897873635384, 48.68939668449611, 48.68939668449611, 47.21520923514558, 47.21520923514558, 45.7885484801673, 45.7885484801673, 44.40770380431811, 44.40770380431811, 43.07103706965834, 43.07103706965834, 41.77697890437404, 41.77697890437404, 40.524025216626946, 40.524025216626946, 39.31073391800763, 39.31073391800763, 38.1357218423103, 38.1357218423103, 36.99766184640428, 36.99766184640428, 35.89528008094817, 35.89528008094817, 34.82735341959329, 34.82735341959329, 33.792707036151306, 33.792707036151306, 32.79021211996428, 32.79021211996428, 31.81878372042438, 31.81878372042438, 30.87737871224071, 30.87737871224071, 29.96499387365312, 29.96499387365312, 29.080664070349634, 29.080664070349634, 28.223460538357696, 28.223460538357696, 27.392489259653978, 27.392489259653978, 26.586889424678173, 26.586889424678173, 25.805831976338997, 25.805831976338997, 25.04851823047987, 25.04851823047987, 24.314178568116613, 24.314178568116613, 23.60207119508215, 23.60207119508215, 22.911480965010405, 22.911480965010405, 22.2417182618682, 22.2417182618682, 21.592117938498276, 21.592117938498276, 20.962038307875016, 20.962038307875016, 20.350860183991877, 20.350860183991877, 19.757985969505228, 19.757985969505228, 19.1828387874474, 19.1828387874474, 18.624861654496534, 18.624861654496534, 18.083516693454932, 18.083516693454932, 17.558284382738123, 17.558284382738123, 17.04866284081756, 17.04866284081756, 16.554167143690737, 16.554167143690737, 16.074328673574186, 16.074328673574186, 15.608694497127134, 15.608694497127134, 15.156826771620556, 15.156826771620556, 14.718302177562125, 14.718302177562125, 14.292711376380993, 14.292711376380993, 13.879658491860328, 13.879658491860328, 13.478760614084793, 13.478760614084793, 13.089647324745192, 13.089647324745192, 12.711960242710555, 12.711960242710555, 12.345352588843555, 12.345352588843555, 11.989488769095177, 11.989488769095177, 11.644043974971833, 11.644043974971833, 11.308703800521819 ], "xaxis": "x1", "y": [ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.7200322574396036, 0.7200322574396036, 1.5000750037651924, 1.5000750037651924, 1.2600609869538875, 1.2600609869538875, 1.8000936049106635, 1.8000936049106635, 1.740089788669434, 1.740089788669434, 3.1801946283562716, 3.1801946283562716, 3.840251921391259, 3.840251921391259, 4.800345626728641, 4.800345626728641, 5.940472865357314, 5.940472865357314, 4.920358204078665, 4.920358204078665, 7.2006336626757195, 7.2006336626757195, 7.560683493867237, 7.560683493867237, 8.400806488715006, 8.400806488715006, 9.060909729677427, 9.060909729677427, 11.581357366023152, 11.581357366023152, 11.461334129244213, 11.461334129244213, 13.861835361062166, 13.861835361062166, 12.361513087099313, 12.361513087099313, 14.221917172958557, 14.221917172958557, 15.135465823229493, 15.135465823229493, 17.94286380083212, 17.94286380083212, 19.803406347998827, 19.803406347998827, 19.503315718367265, 19.503315718367265, 19.323261916969948, 19.323261916969948, 22.38423533316738, 22.38423533316738, 21.483936065107628, 21.483936065107628, 23.464608721525924, 23.464608721525924, 22.924420081861328, 22.924420081861328, 23.764715202646418, 23.764715202646418, 22.4442556686553, 22.4442556686553, 23.34456646534051, 23.34456646534051, 19.743388126008508, 19.743388126008508, 22.26419480629913, 22.26419480629913, 21.904074378548984, 21.904074378548984, 24.665041851733996, 24.665041851733996, 22.144154471573778, 22.144154471573778, 21.904074378548984, 21.904074378548984, 24.877847911754355, 24.877847911754355, 21.603975343059435, 21.603975343059435, 22.4442556686553, 22.4442556686553, 21.84405447538101, 21.84405447538101, 25.445333690730784, 25.445333690730784, 21.30387750843918, 21.30387750843918, 24.60501973885711, 24.60501973885711, 25.085197986974226, 25.085197986974226, 24.90513078363427, 24.90513078363427, 25.085197986974226, 25.085197986974226, 25.085197986974226, 25.085197986974226, 25.02517553782124, 25.02517553782124, 26.4057040225599, 26.4057040225599, 24.424953688460665, 24.424953688460665, 25.265265622672285, 25.265265622672285, 24.168496726708252, 24.168496726708252, 24.544997674019314, 24.544997674019314, 22.864399362083123, 22.864399362083123, 20.223535246836864, 20.223535246836864, 20.40359120968603, 20.40359120968603, 18.483016559491443, 18.483016559491443, 16.382444351236078, 16.382444351236078, 17.402714932561945, 17.402714932561945, 15.362187649784662, 15.362187649784662, 13.441742098922902, 13.441742098922902, 13.474476554293103, 13.474476554293103, 12.54155017524025, 12.54155017524025, 11.52134572362356, 11.52134572362356, 9.78102898249031, 9.78102898249031, 9.000900104053622, 9.000900104053622, 9.600998521053487, 9.600998521053487, 9.000900104053622, 9.000900104053622, 8.340797391269184, 8.340797391269184, 9.3009487123405, 9.3009487123405, 7.146080808738024, 7.146080808738024, 8.10076148164165, 8.10076148164165, 6.000480042244097, 6.000480042244097, 5.34040373714643, 5.34040373714643, 5.34040373714643, 5.34040373714643, 4.680333241420863, 4.680333241420863, 4.920358204078665, 4.920358204078665, 4.680333241420863, 4.680333241420863, 5.127653108198008, 5.127653108198008, 4.56032104815464, 4.56032104815464, 3.42021479006593, 3.42021479006593, 4.320297237744097, 4.320297237744097, 4.86035189139833, 4.86035189139833, 5.280397088389138, 5.280397088389138, 6.3005166468641, 6.3005166468641, 5.1603839349082214, 5.1603839349082214, 4.691243450331413, 4.691243450331413, 3.1201897079496437, 3.1201897079496437, 4.680333241420863, 4.680333241420863, 5.34040373714643, 5.34040373714643, 3.840251921391259, 3.840251921391259, 2.9401752347782644, 2.9401752347782644, 2.5201431443169056, 2.5201431443169056, 2.8201658260369027, 2.8201658260369027, 2.5637827325224656, 2.5637827325224656, 2.7001566093263794, 2.7001566093263794, 1.920101381411932, 1.920101381411932, 2.2801258630724237, 2.2801258630724237, 1.3200644191483035, 1.3200644191483035, 1.2000576027649343, 1.2000576027649343, 0.3000126005274222, 0.3000126005274222, 0.06000242409784236, 0.06000242409784236, 0, 0, 0, 0, 0.06000242409784236, 0.06000242409784236, 0, 0, 0.06000242409784236, 0.06000242409784236, 0.3000126005274222, 0.3000126005274222, 0.3000126005274222, 0.3000126005274222, 0.1800074163048173, 0.1800074163048173, 0.8400383057410972, 0.8400383057410972, 0.49093033479978115, 0.49093033479978115, 1.560078627982324, 1.560078627982324, 1.5000750037651924, 1.5000750037651924 ], "yaxis": "y1", "zmax": 0.1, "zmin": -0.1 } ], "layout": { "annotations": [], "height": 288, "margin": { "b": 20, "l": 0, "r": 0, "t": 20 }, "paper_bgcolor": "rgba(255, 255, 255, 1.000)", "plot_bgcolor": "rgba(255, 255, 255, 1.000)", "showlegend": false, "width": 432, "xaxis": { "anchor": "y1", "domain": [ 0.10942208612812285, 0.9421296296296298 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ 0.9030899869919435, 3 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 17 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "10^1", "10^2", "10^3" ], "tickvals": [ 10, 100, 1000 ], "title": "Apparent +1 Mobility Diam. (nm)", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 17 }, "type": "log", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" }, "yaxis": { "anchor": "x1", "domain": [ 0.1294109069699621, 0.9863298337707787 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ -0.792171120676797, 27.197875143236697 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 17 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "0", "5", "10", "15", "20", "25" ], "tickvals": [ 0, 5, 10, 15, 20, 25 ], "title": "Concentration (cm⁻³)", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 17 }, "type": "-", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" } } }, "image/svg+xml": [ "10251002510000510152025Apparent +1 Mobility Diam. (nm)Concentration (cm⁻³)" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Qcpc,τ = 16.6666, 0.4e-6 # 1 L min-1 = 16.6666 cm3 s-1, 0.4 μs \n", "correct = @. x -> -lambertw(-x*Qcpc*τ,0)/(Qcpc*τ)\n", "\n", "a,b = size(rawc)\n", "𝕣 = SizeDistribution[]\n", "for j = 1:b\n", " R = Float64[]\n", " for i = 1:length(δ.De)-1\n", " ii = (rawdp .<= δ.De[i]) .& (rawdp .> δ.De[i+1])\n", " un = rawdp[ii]\n", " c = rawc[ii,j]\n", " Nm = length(c) > 0 ? sum(c)./(length(c)*1.6666) : 0\n", " Na = correct(Nm)\n", " push!(R,Na)\n", " end\n", " push!(𝕣, SizeDistribution([],δ.De,δ.Dp,δ.ΔlnD,R./δ.ΔlnD,R,:raw))\n", "end\n", "\n", "j = 110 # plot example distribtution\n", "figure(\"Nimbus Sans L\", 2, 3, 2, 8)\n", "p = plot(𝕣[j].Dp, 𝕣[j].N,xaxis = :log10,xticks = [10, 100, 1000],lt = :steppre,legend = :none,color = :black,\n", " ylabel = \"Concentration (cm⁻³)\", xlim = (8,1000), xlabel = \"Apparent +1 Mobility Diam. (nm)\", \n", " right_margin = 25px, fmt = :svg)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " Figure 1. Coincidence corrected size distribution obtained from raw counts." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Block 4. Inversion using regularization\n", "Theory of the inversion is presented in detail in Notebook S3. Here, the same routines are applied to actual data. The block makes use of the wrapper function rinv(), which computes the inverse of a response function and returns a size distribution\n", "```julia\n", "function rinv(R, δ;λ₁= 1e-2, λ₂=1e1)\n", " eyeM = Matrix{Float64}(I, length(R), length(R))\n", " setupRegularization(δ.𝐀,eyeM,R,inv(δ.𝐒)*R) # setup the system\n", " λopt = lcorner(λ₁,λ₂;n=10,r=3) # compute the optimal λ\n", " N = clean((reginv(λopt, r = :Nλ))[1]) # find the inverted size\n", " return SizeDistribution([],δ.De,δ.Dp,δ.ΔlnD,N./δ.ΔlnD,N,:regularized)\n", "end\n", "```" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "\u001b[32mInverting TSI SMPS ...100%|█████████████████████████████| Time: 0:03:52\u001b[39m\n" ] } ], "source": [ "λ₁,λ₂ = 1e-2, 1e1 # bounds [λ₁,λ₂] within which the optimal distribution lies\n", "\n", "𝕟ᵢₙᵥ = SizeDistribution[]\n", "@showprogress 1 \"Inverting TSI SMPS ...\" for i = 1:b\n", " push!(𝕟ᵢₙᵥ, rinv(𝕣[i].N, δ,λ₁=λ₁,λ₂=λ₂))\n", "end\n", "\n", "# Compute size distribution statistics.\n", "Nt, Ninv = Float64[],Float64[]\n", "At, Ainv = Float64[],Float64[]\n", "Vt, Vinv = Float64[],Float64[]\n", "for i = 1:b\n", " push!(Nt, sum(𝕟[i].N))\n", " push!(Ninv,sum(𝕟ᵢₙᵥ[i].N))\n", " push!(At, sum(π/4.0.*(𝕟[i].Dp./1000.0).^2 .* 𝕟[i].N))\n", " push!(Ainv,sum(π/4.0.*(𝕟ᵢₙᵥ[i].Dp./1000.0).^2 .* 𝕟ᵢₙᵥ[i].N))\n", " push!(Vt, sum(π/6.0.*(𝕟[i].Dp./1000.0).^3 .* 𝕟[i].N))\n", " push!(Vinv,sum(π/6.0.*(𝕟ᵢₙᵥ[i].Dp./1000.0).^3 .* 𝕟ᵢₙᵥ[i].N))\n", "end" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Block 5. Visualization" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### (a) Example size distribution comparison" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "application/vnd.plotly.v1+json": { "data": [ { "colorbar": { "title": "" }, "hoverinfo": "text", "legendgroup": "TSI AIM Inversion", "line": { "color": "rgba(0, 0, 0, 1.000)", "dash": "solid", "shape": "linear", "width": 1.5 }, "mode": "lines", "name": "TSI AIM Inversion", "showlegend": true, "text": [ null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ], "type": "scatter", "x": [ 552.3, 532.8, 514, 495.8, 478.3, 461.4, 445.1, 429.4, 414.2, 399.5, 385.4, 371.8, 358.7, 346, 333.8, 322, 310.6, 299.6, 289, 278.8, 269, 259.5, 250.3, 241.4, 232.9, 224.7, 216.7, 209.1, 201.7, 194.6, 187.7, 181.1, 174.7, 168.5, 162.5, 156.8, 151.2, 145.9, 140.7, 135.8, 131, 126.3, 121.9, 117.6, 113.4, 109.4, 105.5, 101.8, 98.2, 94.7, 91.4, 88.2, 85.1, 82, 79.1, 76.4, 73.7, 71, 68.5, 66.1, 63.8, 61.5, 59.4, 57.3, 55.2, 53.3, 51.4, 49.6, 47.8, 46.1, 44.5, 42.9, 41.4, 40, 38.5, 37.2, 35.9, 34.6, 33.4, 32.2, 31.1, 30, 28.9, 27.9, 26.9, 25.9, 25, 24.1, 23.3, 22.5, 21.7, 20.9, 20.2, 19.5, 18.8, 18.1, 17.5, 16.8, 16.3, 15.7, 15.1, 14.6, 14.1, 13.6, 13.1, 12.6, 12.2 ], "xaxis": "x1", "y": [ 8.00626220333464, 41.66960009075682, 17.100779519422446, 42.83989343114151, 57.95312425413374, 39.925734028122505, 60.60318918270739, 67.87762175458684, 57.48321762471443, 135.8881376206361, 129.8288610091219, 151.39939933629262, 204.7659395670461, 170.98912053150258, 173.76122220949105, 247.37674265898366, 222.09124933361258, 262.6817144957361, 263.8725499651149, 344.0220308948419, 302.9490645633238, 362.53209600804047, 376.07687230963904, 376.8008412109717, 349.04290940012544, 397.34991891670603, 466.35410055734985, 512.5326328181227, 462.623510957801, 461.84612383519413, 581.6462566682061, 497.78833515750716, 530.7947157821543, 534.885769801683, 502.43528611387205, 494.153290343977, 510.10492666428337, 436.7308739467291, 463.1794078946371, 474.2539171831701, 510.41327574643486, 461.6333195390615, 502.42225727941496, 562.1247197066549, 451.50991516589687, 506.29182111317294, 607.5823231274683, 563.0671387323849, 624.7065545489137, 675.8925021860308, 689.1905992219083, 715.2222104671893, 705.7849913754317, 821.893621112266, 778.8029226178255, 829.5849763867727, 747.642293541267, 819.4659149584269, 782.220820190404, 700.4040827446504, 671.6885316012073, 580.2174278227445, 622.6566845943303, 586.7752744994835, 478.4578877679935, 499.48208363692993, 471.4657466093512, 416.49839691830226, 341.6381884836749, 401.3345707881683, 375.80587255293136, 371.31657049349747, 426.4116027622259, 364.1576602538043, 333.4830067024957, 267.8810880330818, 273.8773919447201, 265.9358830486372, 256.2024751202215, 284.35257484822654, 314.1008782696355, 217.5949985624682, 283.17954545260585, 345.6406464288953, 427.03134098790173, 387.4792739320089, 364.60976080946557, 283.11700704721176, 411.70856307739126, 366.941922177286, 291.62440165321453, 228.95223355872008, 271.92523824856494, 290.14432605888834, 210.310142923023, 247.9018046876047, 135.9150638785141, 127.64435976514854, 13.622601872755679, 0, 0, 7.438205021005183, 0, 36.128263078360476, 61.71845741223492, 41.98906711164486, 196.53388766256998 ], "yaxis": "y1", "zmax": 0.1, "zmin": -0.1 }, { "colorbar": { "title": "" }, "hoverinfo": "text", "legendgroup": "N=(𝐀ᵀ𝐀+λ²𝐈)⁻¹(𝐀ᵀR-λ²𝐒⁻¹R)", "line": { "color": "rgba(163, 0, 0, 1.000)", "dash": "solid", "shape": "linear", "width": 1 }, "mode": "lines", "name": "N=(𝐀ᵀ𝐀+λ²𝐈)⁻¹(𝐀ᵀR-λ²𝐒⁻¹R)", "showlegend": true, "text": [ null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ], "type": "scatter", "x": [ 641.2116606343327, 641.2116606343327, 611.4599387459494, 611.4599387459494, 583.306474390341, 583.306474390341, 556.6613542884901, 556.6613542884901, 531.4394703614233, 531.4394703614233, 507.56026030320544, 507.56026030320544, 484.94747088813676, 484.94747088813676, 463.52894137710354, 463.52894137710354, 443.2364043391983, 443.2364043391983, 424.00530127516816, 424.00530127516816, 405.7746106002292, 405.7746106002292, 388.48668579141435, 388.48668579141435, 372.08710180294423, 372.08710180294423, 356.5245081767597, 356.5245081767597, 341.75048760168943, 341.75048760168943, 327.7194189851799, 327.7194189851799, 314.38834438227576, 314.38834438227576, 301.71683936864554, 301.71683936864554, 289.6668866433689, 289.6668866433689, 278.20275280215174, 278.20275280215174, 267.2908683347188, 267.2908683347188, 256.899710975487, 256.899710975487, 246.99969257957602, 246.99969257957602, 237.56304971266977, 237.56304971266977, 228.56373813904784, 228.56373813904784, 219.97733137279684, 219.97733137279684, 211.78092342763756, 211.78092342763756, 203.95303586511605, 203.95303586511605, 196.4735292024492, 196.4735292024492, 189.32351870267976, 189.32351870267976, 182.48529453289518, 182.48529453289518, 175.94224624241147, 175.94224624241147, 169.67879148283004, 169.67879148283004, 163.68030886617655, 163.68030886617655, 157.93307483604227, 157.93307483604227, 152.4242044096698, 152.4242044096698, 147.14159563600944, 147.14159563600944, 142.07387760555852, 142.07387760555852, 137.21036184188898, 137.21036184188898, 132.5409969017305, 132.5409969017305, 128.0563260098899, 128.0563260098899, 123.74744755674116, 123.74744755674116, 119.60597828913852, 119.60597828913852, 115.62401903004877, 115.62401903004877, 111.79412276766202, 111.79412276766202, 108.10926496096282, 108.10926496096282, 104.56281591549961, 104.56281591549961, 101.14851509018342, 101.14851509018342, 97.86044720323027, 97.86044720323027, 94.69302001269416, 94.69302001269416, 91.64094365432143, 91.64094365432143, 88.69921142660979, 88.69921142660979, 85.86308191991205, 85.86308191991205, 83.12806239313848, 83.12806239313848, 80.48989330804814, 80.48989330804814, 77.94453393726194, 77.94453393726194, 75.48814896795007, 75.48814896795007, 73.11709602865756, 73.11709602865756, 70.82791407191412, 70.82791407191412, 68.61731255014551, 68.61731255014551, 66.48216132696216, 66.48216132696216, 64.41948127016396, 64.41948127016396, 62.426435476770884, 62.426435476770884, 60.50032108408952, 60.50032108408952, 58.638561624262884, 58.638561624262884, 56.838699882942535, 56.838699882942535, 55.098391225681986, 55.098391225681986, 53.415397358389946, 53.415397358389946, 51.78758049072226, 51.78758049072226, 50.212897873635384, 50.212897873635384, 48.68939668449611, 48.68939668449611, 47.21520923514558, 47.21520923514558, 45.7885484801673, 45.7885484801673, 44.40770380431811, 44.40770380431811, 43.07103706965834, 43.07103706965834, 41.77697890437404, 41.77697890437404, 40.524025216626946, 40.524025216626946, 39.31073391800763, 39.31073391800763, 38.1357218423103, 38.1357218423103, 36.99766184640428, 36.99766184640428, 35.89528008094817, 35.89528008094817, 34.82735341959329, 34.82735341959329, 33.792707036151306, 33.792707036151306, 32.79021211996428, 32.79021211996428, 31.81878372042438, 31.81878372042438, 30.87737871224071, 30.87737871224071, 29.96499387365312, 29.96499387365312, 29.080664070349634, 29.080664070349634, 28.223460538357696, 28.223460538357696, 27.392489259653978, 27.392489259653978, 26.586889424678173, 26.586889424678173, 25.805831976338997, 25.805831976338997, 25.04851823047987, 25.04851823047987, 24.314178568116613, 24.314178568116613, 23.60207119508215, 23.60207119508215, 22.911480965010405, 22.911480965010405, 22.2417182618682, 22.2417182618682, 21.592117938498276, 21.592117938498276, 20.962038307875016, 20.962038307875016, 20.350860183991877, 20.350860183991877, 19.757985969505228, 19.757985969505228, 19.1828387874474, 19.1828387874474, 18.624861654496534, 18.624861654496534, 18.083516693454932, 18.083516693454932, 17.558284382738123, 17.558284382738123, 17.04866284081756, 17.04866284081756, 16.554167143690737, 16.554167143690737, 16.074328673574186, 16.074328673574186, 15.608694497127134, 15.608694497127134, 15.156826771620556, 15.156826771620556, 14.718302177562125, 14.718302177562125, 14.292711376380993, 14.292711376380993, 13.879658491860328, 13.879658491860328, 13.478760614084793, 13.478760614084793, 13.089647324745192, 13.089647324745192, 12.711960242710555, 12.711960242710555, 12.345352588843555, 12.345352588843555, 11.989488769095177, 11.989488769095177, 11.644043974971833, 11.644043974971833, 11.308703800521819 ], "xaxis": "x1", "y": [ 9.794843599707866, 9.539089078501105, 9.539089078501105, 1.8703637563411324, 1.8703637563411324, 0, 0, 0, 0, 28.93810444595935, 28.93810444595935, 70.62723651631723, 70.62723651631723, 66.91796660849091, 66.91796660849091, 75.73498792631153, 75.73498792631153, 66.40331798773849, 66.40331798773849, 107.92415485454899, 107.92415485454899, 133.2263626014492, 133.2263626014492, 173.7243758509, 173.7243758509, 215.5914235295337, 215.5914235295337, 183.1121360327545, 183.1121360327545, 229.35317888649533, 229.35317888649533, 223.51603047437567, 223.51603047437567, 244.93292285170898, 244.93292285170898, 271.58845845130514, 271.58845845130514, 346.571532209115, 346.571532209115, 343.35061874787397, 343.35061874787397, 382.3862762094065, 382.3862762094065, 326.4931690888932, 326.4931690888932, 348.6906035217081, 348.6906035217081, 354.83320107701917, 354.83320107701917, 441.0444990535429, 441.0444990535429, 516.4755986602988, 516.4755986602988, 505.08960198759644, 505.08960198759644, 469.102776191186, 469.102776191186, 510.7130394238147, 510.7130394238147, 492.8891317781741, 492.8891317781741, 549.447766786939, 549.447766786939, 543.0065106923687, 543.0065106923687, 543.0887601887559, 543.0887601887559, 499.14618780079127, 499.14618780079127, 493.8978546283724, 493.8978546283724, 383.1162491364551, 383.1162491364551, 431.36311884618226, 431.36311884618226, 457.7541300012219, 457.7541300012219, 557.3882622612988, 557.3882622612988, 499.0316826456282, 499.0316826456282, 472.3969005705528, 472.3969005705528, 565.339154082934, 565.339154082934, 453.0476181887128, 453.0476181887128, 488.53666002872853, 488.53666002872853, 502.09373018230343, 502.09373018230343, 623.5009289767337, 623.5009289767337, 536.7624557466448, 536.7624557466448, 658.1068646734118, 658.1068646734118, 695.3428626944735, 695.3428626944735, 692.5019344451716, 692.5019344451716, 722.4735700953453, 722.4735700953453, 719.4956104937633, 719.4956104937633, 728.2231912202011, 728.2231912202011, 803.228662557507, 803.228662557507, 764.2506954248317, 764.2506954248317, 808.4665948524668, 808.4665948524668, 789.37888606167, 789.37888606167, 810.6224368820817, 810.6224368820817, 750.9587879336658, 750.9587879336658, 647.4597272870468, 647.4597272870468, 658.4458311286688, 658.4458311286688, 586.254247237369, 586.254247237369, 526.4491098130969, 526.4491098130969, 599.5274597477704, 599.5274597477704, 532.547289306425, 532.547289306425, 464.0460887416672, 464.0460887416672, 476.5624496351967, 476.5624496351967, 443.59650291071483, 443.59650291071483, 411.2548361402344, 411.2548361402344, 342.54130257285885, 342.54130257285885, 317.0534581555238, 317.0534581555238, 368.8322729208582, 368.8322729208582, 375.05404351662753, 375.05404351662753, 377.88400552111176, 377.88400552111176, 442.43650599040035, 442.43650599040035, 339.5705477661618, 339.5705477661618, 386.51223059997466, 386.51223059997466, 273.3075432440192, 273.3075432440192, 235.0776020856017, 235.0776020856017, 250.2970507090045, 250.2970507090045, 232.50291704812645, 232.50291704812645, 269.4378201958652, 269.4378201958652, 271.16869038818527, 271.16869038818527, 301.5337669625609, 301.5337669625609, 263.85813289797227, 263.85813289797227, 191.50737758254954, 191.50737758254954, 262.43005950011514, 262.43005950011514, 323.7538166661209, 323.7538166661209, 381.72399039140015, 381.72399039140015, 472.6565893683877, 472.6565893683877, 392.45009219032767, 392.45009219032767, 348.0513218575341, 348.0513218575341, 225.71260471933954, 225.71260471933954, 357.29918581743334, 357.29918581743334, 431.6208981474245, 431.6208981474245, 322.3073766002152, 322.3073766002152, 246.72351349392247, 246.72351349392247, 205.76620887596673, 205.76620887596673, 237.25022747523127, 237.25022747523127, 230.19666968543098, 230.19666968543098, 260.4962424399532, 260.4962424399532, 195.02834274159397, 195.02834274159397, 233.68035425565927, 233.68035425565927, 132.64051540797703, 132.64051540797703, 115.23911199591511, 115.23911199591511, 12.57244644073366, 12.57244644073366, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4.772627799839306, 4.772627799839306, 39.5376914077545, 39.5376914077545, 42.64246526621957, 42.64246526621957, 29.527056747312802, 29.527056747312802, 138.29848397352842, 138.29848397352842, 96.5087414558226, 96.5087414558226, 328.1992718457263, 328.1992718457263, 391.0841480648158, 391.0841480648158 ], "yaxis": "y1", "zmax": 0.1, "zmin": -0.1 }, { "colorbar": { "title": "" }, "hoverinfo": "text", "legendgroup": "R measured", "line": { "color": "rgba(58, 129, 252, 1.000)", "dash": "solid", "shape": "linear", "width": 1 }, "mode": "lines", "name": "R measured", "showlegend": true, "text": [ null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ], "type": "scatter", "x": [ 641.2116606343327, 641.2116606343327, 611.4599387459494, 611.4599387459494, 583.306474390341, 583.306474390341, 556.6613542884901, 556.6613542884901, 531.4394703614233, 531.4394703614233, 507.56026030320544, 507.56026030320544, 484.94747088813676, 484.94747088813676, 463.52894137710354, 463.52894137710354, 443.2364043391983, 443.2364043391983, 424.00530127516816, 424.00530127516816, 405.7746106002292, 405.7746106002292, 388.48668579141435, 388.48668579141435, 372.08710180294423, 372.08710180294423, 356.5245081767597, 356.5245081767597, 341.75048760168943, 341.75048760168943, 327.7194189851799, 327.7194189851799, 314.38834438227576, 314.38834438227576, 301.71683936864554, 301.71683936864554, 289.6668866433689, 289.6668866433689, 278.20275280215174, 278.20275280215174, 267.2908683347188, 267.2908683347188, 256.899710975487, 256.899710975487, 246.99969257957602, 246.99969257957602, 237.56304971266977, 237.56304971266977, 228.56373813904784, 228.56373813904784, 219.97733137279684, 219.97733137279684, 211.78092342763756, 211.78092342763756, 203.95303586511605, 203.95303586511605, 196.4735292024492, 196.4735292024492, 189.32351870267976, 189.32351870267976, 182.48529453289518, 182.48529453289518, 175.94224624241147, 175.94224624241147, 169.67879148283004, 169.67879148283004, 163.68030886617655, 163.68030886617655, 157.93307483604227, 157.93307483604227, 152.4242044096698, 152.4242044096698, 147.14159563600944, 147.14159563600944, 142.07387760555852, 142.07387760555852, 137.21036184188898, 137.21036184188898, 132.5409969017305, 132.5409969017305, 128.0563260098899, 128.0563260098899, 123.74744755674116, 123.74744755674116, 119.60597828913852, 119.60597828913852, 115.62401903004877, 115.62401903004877, 111.79412276766202, 111.79412276766202, 108.10926496096282, 108.10926496096282, 104.56281591549961, 104.56281591549961, 101.14851509018342, 101.14851509018342, 97.86044720323027, 97.86044720323027, 94.69302001269416, 94.69302001269416, 91.64094365432143, 91.64094365432143, 88.69921142660979, 88.69921142660979, 85.86308191991205, 85.86308191991205, 83.12806239313848, 83.12806239313848, 80.48989330804814, 80.48989330804814, 77.94453393726194, 77.94453393726194, 75.48814896795007, 75.48814896795007, 73.11709602865756, 73.11709602865756, 70.82791407191412, 70.82791407191412, 68.61731255014551, 68.61731255014551, 66.48216132696216, 66.48216132696216, 64.41948127016396, 64.41948127016396, 62.426435476770884, 62.426435476770884, 60.50032108408952, 60.50032108408952, 58.638561624262884, 58.638561624262884, 56.838699882942535, 56.838699882942535, 55.098391225681986, 55.098391225681986, 53.415397358389946, 53.415397358389946, 51.78758049072226, 51.78758049072226, 50.212897873635384, 50.212897873635384, 48.68939668449611, 48.68939668449611, 47.21520923514558, 47.21520923514558, 45.7885484801673, 45.7885484801673, 44.40770380431811, 44.40770380431811, 43.07103706965834, 43.07103706965834, 41.77697890437404, 41.77697890437404, 40.524025216626946, 40.524025216626946, 39.31073391800763, 39.31073391800763, 38.1357218423103, 38.1357218423103, 36.99766184640428, 36.99766184640428, 35.89528008094817, 35.89528008094817, 34.82735341959329, 34.82735341959329, 33.792707036151306, 33.792707036151306, 32.79021211996428, 32.79021211996428, 31.81878372042438, 31.81878372042438, 30.87737871224071, 30.87737871224071, 29.96499387365312, 29.96499387365312, 29.080664070349634, 29.080664070349634, 28.223460538357696, 28.223460538357696, 27.392489259653978, 27.392489259653978, 26.586889424678173, 26.586889424678173, 25.805831976338997, 25.805831976338997, 25.04851823047987, 25.04851823047987, 24.314178568116613, 24.314178568116613, 23.60207119508215, 23.60207119508215, 22.911480965010405, 22.911480965010405, 22.2417182618682, 22.2417182618682, 21.592117938498276, 21.592117938498276, 20.962038307875016, 20.962038307875016, 20.350860183991877, 20.350860183991877, 19.757985969505228, 19.757985969505228, 19.1828387874474, 19.1828387874474, 18.624861654496534, 18.624861654496534, 18.083516693454932, 18.083516693454932, 17.558284382738123, 17.558284382738123, 17.04866284081756, 17.04866284081756, 16.554167143690737, 16.554167143690737, 16.074328673574186, 16.074328673574186, 15.608694497127134, 15.608694497127134, 15.156826771620556, 15.156826771620556, 14.718302177562125, 14.718302177562125, 14.292711376380993, 14.292711376380993, 13.879658491860328, 13.879658491860328, 13.478760614084793, 13.478760614084793, 13.089647324745192, 13.089647324745192, 12.711960242710555, 12.711960242710555, 12.345352588843555, 12.345352588843555, 11.989488769095177, 11.989488769095177, 11.644043974971833, 11.644043974971833, 11.308703800521819 ], "xaxis": "x2", "y": [ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.7200322574396036, 0.7200322574396036, 1.5000750037651924, 1.5000750037651924, 1.2600609869538875, 1.2600609869538875, 1.8000936049106635, 1.8000936049106635, 1.740089788669434, 1.740089788669434, 3.1801946283562716, 3.1801946283562716, 3.840251921391259, 3.840251921391259, 4.800345626728641, 4.800345626728641, 5.940472865357314, 5.940472865357314, 4.920358204078665, 4.920358204078665, 7.2006336626757195, 7.2006336626757195, 7.560683493867237, 7.560683493867237, 8.400806488715006, 8.400806488715006, 9.060909729677427, 9.060909729677427, 11.581357366023152, 11.581357366023152, 11.461334129244213, 11.461334129244213, 13.861835361062166, 13.861835361062166, 12.361513087099313, 12.361513087099313, 14.221917172958557, 14.221917172958557, 15.135465823229493, 15.135465823229493, 17.94286380083212, 17.94286380083212, 19.803406347998827, 19.803406347998827, 19.503315718367265, 19.503315718367265, 19.323261916969948, 19.323261916969948, 22.38423533316738, 22.38423533316738, 21.483936065107628, 21.483936065107628, 23.464608721525924, 23.464608721525924, 22.924420081861328, 22.924420081861328, 23.764715202646418, 23.764715202646418, 22.4442556686553, 22.4442556686553, 23.34456646534051, 23.34456646534051, 19.743388126008508, 19.743388126008508, 22.26419480629913, 22.26419480629913, 21.904074378548984, 21.904074378548984, 24.665041851733996, 24.665041851733996, 22.144154471573778, 22.144154471573778, 21.904074378548984, 21.904074378548984, 24.877847911754355, 24.877847911754355, 21.603975343059435, 21.603975343059435, 22.4442556686553, 22.4442556686553, 21.84405447538101, 21.84405447538101, 25.445333690730784, 25.445333690730784, 21.30387750843918, 21.30387750843918, 24.60501973885711, 24.60501973885711, 25.085197986974226, 25.085197986974226, 24.90513078363427, 24.90513078363427, 25.085197986974226, 25.085197986974226, 25.085197986974226, 25.085197986974226, 25.02517553782124, 25.02517553782124, 26.4057040225599, 26.4057040225599, 24.424953688460665, 24.424953688460665, 25.265265622672285, 25.265265622672285, 24.168496726708252, 24.168496726708252, 24.544997674019314, 24.544997674019314, 22.864399362083123, 22.864399362083123, 20.223535246836864, 20.223535246836864, 20.40359120968603, 20.40359120968603, 18.483016559491443, 18.483016559491443, 16.382444351236078, 16.382444351236078, 17.402714932561945, 17.402714932561945, 15.362187649784662, 15.362187649784662, 13.441742098922902, 13.441742098922902, 13.474476554293103, 13.474476554293103, 12.54155017524025, 12.54155017524025, 11.52134572362356, 11.52134572362356, 9.78102898249031, 9.78102898249031, 9.000900104053622, 9.000900104053622, 9.600998521053487, 9.600998521053487, 9.000900104053622, 9.000900104053622, 8.340797391269184, 8.340797391269184, 9.3009487123405, 9.3009487123405, 7.146080808738024, 7.146080808738024, 8.10076148164165, 8.10076148164165, 6.000480042244097, 6.000480042244097, 5.34040373714643, 5.34040373714643, 5.34040373714643, 5.34040373714643, 4.680333241420863, 4.680333241420863, 4.920358204078665, 4.920358204078665, 4.680333241420863, 4.680333241420863, 5.127653108198008, 5.127653108198008, 4.56032104815464, 4.56032104815464, 3.42021479006593, 3.42021479006593, 4.320297237744097, 4.320297237744097, 4.86035189139833, 4.86035189139833, 5.280397088389138, 5.280397088389138, 6.3005166468641, 6.3005166468641, 5.1603839349082214, 5.1603839349082214, 4.691243450331413, 4.691243450331413, 3.1201897079496437, 3.1201897079496437, 4.680333241420863, 4.680333241420863, 5.34040373714643, 5.34040373714643, 3.840251921391259, 3.840251921391259, 2.9401752347782644, 2.9401752347782644, 2.5201431443169056, 2.5201431443169056, 2.8201658260369027, 2.8201658260369027, 2.5637827325224656, 2.5637827325224656, 2.7001566093263794, 2.7001566093263794, 1.920101381411932, 1.920101381411932, 2.2801258630724237, 2.2801258630724237, 1.3200644191483035, 1.3200644191483035, 1.2000576027649343, 1.2000576027649343, 0.3000126005274222, 0.3000126005274222, 0.06000242409784236, 0.06000242409784236, 0, 0, 0, 0, 0.06000242409784236, 0.06000242409784236, 0, 0, 0.06000242409784236, 0.06000242409784236, 0.3000126005274222, 0.3000126005274222, 0.3000126005274222, 0.3000126005274222, 0.1800074163048173, 0.1800074163048173, 0.8400383057410972, 0.8400383057410972, 0.49093033479978115, 0.49093033479978115, 1.560078627982324, 1.560078627982324, 1.5000750037651924, 1.5000750037651924 ], "yaxis": "y2", "zmax": 0.1, "zmin": -0.1 } ], "layout": { "annotations": [], "height": 252, "legend": { "bgcolor": "rgba(255, 255, 255, 1.000)", "bordercolor": "rgba(0, 0, 0, 1.000)", "font": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 14 }, "tracegroupgap": 0, "x": 1, "y": 1 }, "margin": { "b": 20, "l": 0, "r": 0, "t": 20 }, "paper_bgcolor": "rgba(255, 255, 255, 1.000)", "plot_bgcolor": "rgba(255, 255, 255, 1.000)", "showlegend": true, "width": 828, "xaxis1": { "anchor": "y1", "domain": [ 0.10466988727858294, 0.47181964573268925 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ 0.9030899869919435, 3 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 17 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "10^1", "10^2", "10^3" ], "tickvals": [ 10, 100, 1000 ], "title": "Mobility Diameter (nm)", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 17 }, "type": "log", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" }, "xaxis2": { "anchor": "y2", "domain": [ 0.5966183574879227, 0.963768115942029 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ 0.9030899869919435, 3 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 17 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "10^1", "10^2", "10^3" ], "tickvals": [ 10, 100, 1000 ], "title": "Apparent +1 Mobility Diam. (nm)", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 17 }, "type": "log", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" }, "yaxis1": { "anchor": "x1", "domain": [ 0.14789817939424243, 1 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ -24.887549291603182, 854.4725256783759 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 17 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "0", "200", "400", "600", "800" ], "tickvals": [ 0, 200, 400, 600, 800 ], "title": "dN/dlnD (cm⁻³)", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 17 }, "type": "-", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" }, "yaxis2": { "anchor": "x2", "domain": [ 0.14789817939424243, 1 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ -0.792171120676797, 27.197875143236697 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 17 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "0", "5", "10", "15", "20", "25" ], "tickvals": [ 0, 5, 10, 15, 20, 25 ], "title": "Concentration (cm⁻³)", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 17 }, "type": "-", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" } } }, "image/svg+xml": [ "1025100251000020040060080010251002510000510152025TSI AIM InversionN=(𝐀ᵀ𝐀+λ²𝐈)⁻¹(𝐀ᵀR-λ²𝐒⁻¹R)R measuredMobility Diameter (nm)Apparent +1 Mobility Diam. (nm)dN/dlnD (cm⁻³)Concentration (cm⁻³)" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "figure(\"Nimbus Sans L\", 2, 5.75, 1.75, 8)\n", "\n", "p2 = plot(𝕟[j].Dp, 𝕟[j].S, xaxis = :log10, xticks = [10, 100, 1000],\n", " ylabel = \"dN/dlnD (cm⁻³)\", label = \"TSI AIM Inversion\", xlim = (8,1000), \n", " color = :black, xlabel = \"Mobility Diameter (nm)\", left_margin = 30px, lw= 1.5)\n", "\n", "p2 = plot!(𝕟ᵢₙᵥ[j].Dp, 𝕟ᵢₙᵥ[j].S, color = RGBA(163/255,0,0,1),\n", " lt = :steppre, label = \"N=(𝐀ᵀ𝐀+λ²𝐈)⁻¹(𝐀ᵀR-λ²𝐒⁻¹R)\", xlim = (8,1000))\n", "\n", "p3 = plot(𝕣[j].Dp, 𝕣[j].N, xaxis = :log10, xticks = [10, 100, 1000], lt = :steppre,\n", " label = \"R measured\", ylabel = \"Concentration (cm⁻³)\", xlim = (8,1000),left_margin = 30px,\n", " color = RGBA(58/255,129/255,252/255,1), xlabel = \"Apparent +1 Mobility Diam. (nm)\")\n", "\n", "plot(p2,p3, layout = grid(1,2), right_margin = 30px, top_margin = 0px, fmt = :svg)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " Figure 2. Left: Size distribution expressed as dN/dlnD output by TSI AIM software and the inversion performed in Block 3. Right: Measured response function as shown in Figure 1." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### (b) Comparison of integrated number, area, and volume concentration" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "application/vnd.plotly.v1+json": { "data": [ { "colorbar": { "title": "" }, "hoverinfo": "text", "legendgroup": "TSI AIM", "line": { "color": "rgba(153, 0, 153, 1.000)", "dash": "solid", "shape": "linear", "width": 1 }, "mode": "lines", "name": "TSI AIM", "showlegend": true, "text": [ null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ], "type": "scatter", "x": [ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160 ], "xaxis": "x1", "y": [ 196.50052644636776, 191.72960798667748, 191.27371555365312, 216.21038329834005, 244.8780176538852, 219.98531666487438, 230.8096129947479, 216.22432836727307, 247.17433122520742, 4977.57192991733, 20328.379275196057, 155068.20847996254, 210033.07859874598, 2975.1149058986543, 1203.132417796781, 1218.3448314618354, 1262.3266938720164, 1579.7410516298105, 1246.1188000164245, 1273.6312662050811, 1225.3348712932955, 1265.4128502239535, 1437.0671062926574, 1880.601341102498, 2008.3620294063937, 1879.9282004951701, 2031.3549202519575, 1815.9845991448021, 1930.538585710011, 1524.3595981860365, 1451.9087679442098, 1463.6594627444445, 1438.8494285447953, 1365.5473616894903, 1326.728315554913, 1365.0678451299525, 1370.2786794037313, 1369.3029820072127, 1426.0508659980305, 1440.5628812044802, 1448.6175121744654, 1340.3435735348396, 1587.887857147671, 1445.6175564199275, 1382.0026554727967, 1089.3571906628435, 1075.8053869990292, 1357.1852763483128, 1262.7397284749832, 1493.1182323290893, 1219.4746175224618, 1302.1658362216529, 1457.21405231707, 1370.3641310618816, 1305.8549007153097, 1293.4387594094424, 1521.341892833851, 1427.621199078114, 1546.5375540194625, 1297.1135447666607, 1360.9468784197768, 1325.5850573282487, 1374.9111320319216, 1264.3178119390816, 1161.68769806865, 1180.1259277617862, 1233.3373403721291, 1198.5358600361353, 1166.2520575111307, 1141.8030325532711, 1117.432622707541, 1135.3876198087517, 1205.2717224860678, 1151.2620590102288, 1144.0449505303484, 1211.6652759944836, 1132.8812172993344, 1130.3651094461263, 1185.8544559175994, 1386.9693012587068, 1420.4914003434426, 1423.4393953647602, 1323.7994589537752, 1259.811251370198, 2352.6341457741178, 1350.973057789779, 1437.5299119411195, 1871.7293866507478, 1569.958458577866, 1360.3047092711736, 1784.9944588601895, 1347.5429996017072, 1248.1181186685078, 1244.3600326227756, 1223.0429120486222, 1384.31806747665, 1306.7688588110211, 1303.1314033750693, 1280.112424630013, 1497.8874190053955, 1131.73310533136, 1357.9135084805998, 1403.9180571279953, 1260.6890339499516, 1390.0397577378035, 1213.0859012930118, 1355.572114329724, 1312.2222162798248, 1294.9932353555737, 1392.0529740324966, 1359.2111046055863, 1433.5800051479018, 1334.4322829078849, 1391.616048322139, 1366.033185547155, 1417.334497065981, 1434.4063315002832, 1450.653505218748, 1510.556651281714, 1463.129484228177, 1493.5752061062378, 1381.611125611063, 1396.0818762309514, 1420.690837702979, 1415.5574139766338, 1921.6273179025015, 1330.0195728049175, 1382.9082659462642, 1382.517009665803, 1499.451744754616, 1380.8300737958616, 1371.2407997705016, 1374.8316804289568, 1440.701119933659, 1377.6053721871067, 1379.2501588417279, 1638.1950038515527, 1751.7781323473632, 1506.7427743446221, 1526.6923817992617, 2116.58042021365, 1551.612943006435, 1759.7195948941514, 2178.170342018214, 1422.6153240893318, 1560.665336366599, 1620.0200797475406, 1570.5199816370261, 1506.9354860780065, 1544.3344894162638, 1555.3296515902102, 1523.1022797316393, 1670.9202108183422, 1491.1948442884557, 1494.3242950197239, 1590.4054148443515, 1464.1097607598049, 1445.9014571197654, 1516.2505271241062, 1580.0862547821894 ], "yaxis": "y1", "zmax": 0.1, "zmin": -0.1 }, { "colorbar": { "title": "" }, "hoverinfo": "text", "legendgroup": "This work", "line": { "color": "rgba(0, 0, 153, 1.000)", "dash": "dash", "shape": "linear", "width": 1 }, "mode": "lines", "name": "This work", "showlegend": true, "text": [ null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ], "type": "scatter", "x": [ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160 ], "xaxis": "x1", "y": [ 193.23799229021432, 188.16471951573416, 188.3585222074517, 212.5572970235677, 239.38944599038788, 217.85872671746142, 226.22865336239124, 212.04543783223502, 243.47621874145858, 5044.599430749679, 20259.01587377313, 149427.11669844866, 212201.87905364487, 3496.9815189246747, 1190.0005712037917, 1212.707453357946, 1263.872119515561, 1567.5791223722272, 1245.1883024752688, 1258.734407608199, 1214.4825339125116, 1253.8514176033552, 1428.894327268066, 1856.1083666090458, 1982.205730451339, 1857.0287364477688, 1997.9921036236876, 1815.4505198192342, 1904.6919425310991, 1515.4193860829664, 1431.8990336363725, 1450.2817467401592, 1426.2492230239727, 1350.908234220685, 1318.0011238461348, 1351.389970146644, 1369.2934664878373, 1354.2288389939038, 1426.7384843709453, 1430.8170743872104, 1433.004230584908, 1329.3016052552312, 1567.4086494501348, 1436.851389978293, 1372.1805764215178, 1083.1192238254241, 1103.6059969610694, 1373.7776425787367, 1302.3158801497843, 1550.9038228321376, 1293.273065352025, 1318.8626151497833, 1510.7866806741968, 2151.7900698498156, 1351.6097792692176, 1332.8384494776562, 1550.0474879249455, 1486.8021149747708, 1573.6771265749082, 1359.9819766254097, 1383.220311530664, 1371.4011977481514, 1400.4282539327648, 1284.4936663713088, 1160.929190466488, 1182.7214645567487, 1253.532548964498, 1206.9254003758258, 1153.021861247311, 1145.206091441279, 1112.1884489592294, 1147.522444594584, 1214.2428566503027, 1149.814549687478, 1141.8471391031921, 1199.993365052746, 1126.6068416067956, 1118.834530322003, 1195.835940381806, 1394.3288764741123, 1460.8704318551092, 1425.5686047745153, 1343.4708741785796, 1271.6989711303797, 2361.334697803266, 1359.2049958044104, 1443.527341544435, 1910.3375180019762, 1579.5095661952691, 1404.5538517087718, 1775.456788198148, 1341.644278712517, 1234.0460893172244, 1244.2246845063792, 1235.7324766373893, 1378.2165968397592, 1316.6475095140986, 1299.1864452929774, 1294.3288185336837, 1473.9699898220065, 1120.2596234688235, 1355.0050375273088, 1396.901866807664, 1262.0376693437486, 1385.5327630558688, 1219.8744722961806, 1344.2755295517557, 1317.8584052285214, 1286.7974690113726, 1385.1691713993864, 1355.3117202592052, 1421.8339506815034, 1329.8975106957598, 1378.3227740774207, 1372.9783873897427, 1405.7152687986072, 1427.0777011319512, 1444.3709370147874, 1500.961429153926, 1459.552426529107, 1485.2249472763608, 1370.8291975087154, 1387.1541371129904, 1428.4990417448212, 1402.4185414078124, 1915.3141794799185, 1315.95321192799, 1370.602122286995, 1377.955493317029, 1485.4932891609064, 1375.5161879412194, 1357.8676525230067, 1358.5965780723984, 1430.4857194258313, 1364.9249019256952, 1363.3405207859762, 1624.9108976039129, 1729.965192598937, 1485.8288348319509, 1519.6468405266723, 2089.5662699155478, 1560.9130049952923, 1744.4549647394476, 2140.975011442188, 1403.1557354561826, 1546.0448068063254, 1620.896586306945, 1558.559014213663, 1510.9115342775076, 1535.6877198037168, 1550.5378277013185, 1508.7882427446991, 1653.2286572518603, 1467.9694256245327, 1480.3058284933127, 1583.338894010587, 1448.8214673144876, 1442.766895306091, 1504.592179004853, 1570.3769627486622 ], "yaxis": "y1", "zmax": 0.1, "zmin": -0.1 }, { "colorbar": { "title": "" }, "hoverinfo": "text", "legendgroup": "TSI AIM", "line": { "color": "rgba(204, 0, 0, 1.000)", "dash": "solid", "shape": "linear", "width": 1 }, "mode": "lines", "name": "TSI AIM", "showlegend": true, "text": [ null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ], "type": "scatter", "x": [ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160 ], "xaxis": "x2", "y": [ 5.111134561301739, 4.900138170363288, 4.944048991751053, 5.329773453889977, 5.625970762338133, 5.676650513546044, 5.763567558653393, 5.393581742542189, 6.060791094971256, 361.7369180880948, 359.03340292744144, 3777.939594017695, 1151.6952149833442, 24.437646370291272, 22.328151159905754, 23.078432501520293, 22.31525493399014, 25.831221657470117, 22.18568453973351, 22.943126712842528, 22.797604664838722, 23.074078328268197, 25.274743290279314, 26.524376695285035, 26.516535416685795, 25.675510502439018, 27.13877192720197, 26.159182212879045, 26.269259551521785, 25.395276764281522, 24.936582842052754, 24.62406226466108, 25.16722217556518, 24.35947058533022, 24.045717576889903, 24.94902546995615, 25.03578701300408, 24.83567496662721, 24.780892494538932, 25.19852909512568, 24.57217335958428, 25.678626283362444, 25.662273876579096, 25.638633519394777, 24.365593364209143, 12.673728423486022, 11.303458868804144, 9.719332292969884, 9.65221431600658, 9.728734781592921, 9.958699551774183, 9.963986329600626, 10.310948965340797, 10.690206813359142, 10.926314239247514, 11.858584699262526, 13.380929390454929, 13.792747818199295, 13.223569642471059, 12.59966705646303, 12.914684903386854, 12.51948136528437, 12.73429705851746, 12.915327203397556, 13.38713593221909, 13.654144812595664, 14.40451603997972, 14.859740559742963, 14.770866364468338, 15.23144461677763, 16.783634816807844, 16.312751379380018, 16.699200483334018, 17.004471759602623, 17.420229350764583, 17.60560093511016, 18.532632159149443, 18.49376450771169, 18.24466549499038, 18.15676486878093, 17.765269906509108, 17.951610634215022, 18.320816125164402, 18.583610183305854, 19.945987563354077, 18.71322338198299, 18.303762399365674, 20.723767279342564, 20.257904861785146, 19.624188313601437, 20.55204871397247, 19.740819199813377, 19.99219734073468, 20.053995379931703, 20.156040754158337, 20.53812815693389, 19.953733360549826, 20.815342195373272, 20.470125940193107, 20.983923291760117, 21.35385531705423, 21.766661771144765, 21.010420636071636, 21.921132451591948, 21.749544497751476, 21.34017748718946, 22.23995131334944, 22.08354176456735, 22.584835822613265, 23.036365780441933, 23.14712604960178, 22.721092144731646, 22.760595523745838, 23.863372042355177, 23.02416402300976, 23.22548669713661, 23.952318204792025, 23.559322018311878, 23.962904121323117, 24.078203197137018, 24.219978724610055, 24.19313720138893, 24.139519106099275, 24.40129391591202, 24.47730945459479, 26.035617010709906, 24.55573932793573, 24.37148085326492, 24.66731842137711, 25.272070175608, 24.941643546495857, 24.510742427872945, 25.055668150103998, 25.230342596278117, 24.79529474314666, 25.631297253682707, 27.066927806070723, 25.657837904204122, 26.155535813894364, 26.33075333766881, 27.094391841737174, 26.488646922898333, 27.473089738386356, 26.78062375517149, 25.571911588063635, 25.637671637776503, 26.84166301046554, 25.399943180034956, 26.077859764429675, 25.461919231788453, 26.084040212187382, 26.35968399224973, 25.64619667859604, 27.287178760005244, 26.643265765065646, 27.915452197498734, 26.919477153679, 26.900767244976528, 28.498736707640642, 27.106371935027614 ], "yaxis": "y2", "zmax": 0.1, "zmin": -0.1 }, { "colorbar": { "title": "" }, "hoverinfo": "text", "legendgroup": "This work", "line": { "color": "rgba(0, 0, 0, 1.000)", "dash": "dash", "shape": "linear", "width": 1 }, "mode": "lines", "name": "This work", "showlegend": true, "text": [ null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ], "type": "scatter", "x": [ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160 ], "xaxis": "x2", "y": [ 5.18877907940097, 4.946675479293642, 4.995785067803004, 5.405763951889549, 5.641788734452565, 5.68805993812453, 5.666594682182609, 5.416120732262276, 6.189771025047425, 324.82011244087613, 318.39550400331, 3339.8374245774194, 1314.8122912321653, 24.671279029463733, 22.611732725324316, 23.211826822157207, 22.509654625035125, 26.158370894989808, 22.59720160815715, 22.954948720804378, 23.022036024200304, 23.274343907970408, 25.739822738025687, 26.69153437260644, 26.832813564629806, 25.993205104791986, 27.4410629106996, 26.373737210910104, 26.672611444821527, 25.488175827111313, 25.328153565605223, 25.031260579435784, 25.60688973454421, 24.576820559264238, 24.323776235918174, 25.33432742651706, 25.1372827744119, 25.03738215781513, 24.976488107804393, 25.455863685537974, 24.784655143859702, 25.923829999238915, 26.037680263613435, 26.05927759673925, 24.70556928748869, 12.869942616351775, 11.43745089515352, 9.85063639715437, 9.79881413701325, 10.055859768360007, 9.966563557785584, 10.203899718041809, 10.451917025471138, 23.197093296128262, 11.147086946567342, 12.105087698968921, 13.66238971670205, 14.102857105512653, 13.540064670261366, 12.897159708854334, 13.183790325012835, 12.764964195344735, 13.003219406553638, 13.215191863246142, 13.67516598203967, 13.905065270353587, 14.624260327748067, 15.059056575813571, 15.126376344692998, 15.477665086678373, 17.18026668191269, 16.58522394249865, 16.968615188920644, 17.417185995020294, 17.831714482067202, 17.958720839463098, 18.78740418312768, 18.892392040155286, 18.709756361575785, 18.715324478058463, 18.298447329730347, 18.404471392379644, 18.720304850944462, 19.03119108084694, 20.257031923852146, 19.26574415383203, 18.706162951131937, 21.0662981276273, 20.78368269328872, 20.08285674521012, 21.160755147917705, 20.182612106580937, 20.44595492542503, 20.425279987491397, 20.652734707372982, 20.966284339874683, 20.36454819053668, 21.292586545740097, 20.82427824836683, 21.61616979502778, 21.62675352835975, 22.38385815600228, 21.5723786456251, 22.53825977388569, 22.157736202493773, 21.789822091559103, 22.648674307379046, 22.595482340894126, 23.06062280372447, 23.6237868782938, 23.81258102449379, 23.188920663259857, 23.43486005629615, 24.164653555926012, 23.67420214207086, 23.81386913460375, 24.454038128274426, 24.07549637027948, 24.665949375255245, 24.54333344105244, 24.901549477748883, 24.651728416145023, 24.418971843665563, 24.82679698401061, 25.165716650970133, 26.9404760437114, 24.95175439981794, 24.78827086573102, 25.104390745850587, 25.59780327235449, 25.179224435360037, 25.058619821256208, 25.708746524304726, 25.755448207087348, 25.19215685837868, 26.28201601790959, 27.642541280732527, 26.177543409830278, 26.675065766367975, 26.76667578370924, 27.468124680325523, 26.85822187438905, 27.82577936299636, 27.198803137275068, 26.167155868347734, 26.38290064720897, 27.44588754282129, 26.011095183523047, 26.604113830407687, 26.05636505905741, 26.487127248143945, 27.014153636085325, 25.99638417676017, 27.84910359903737, 27.081659183601623, 28.110538639909784, 27.484159172294536, 27.4332626770125, 29.060057697033688, 27.807911125937757 ], "yaxis": "y2", "zmax": 0.1, "zmin": -0.1 }, { "colorbar": { "title": "" }, "hoverinfo": "text", "legendgroup": "TSI AIM", "line": { "color": "rgba(0, 0, 204, 1.000)", "dash": "solid", "shape": "linear", "width": 1 }, "mode": "lines", "name": "TSI AIM", "showlegend": true, "text": [ null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ], "type": "scatter", "x": [ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160 ], "xaxis": "x3", "y": [ 1.0378039008278503, 0.9757623181961994, 0.9898452959182209, 1.080321854763355, 1.1631731213419203, 1.178366527892624, 1.1994466566798578, 1.1035351090608532, 1.2157152958493498, 100.29798458848424, 97.40242352175201, 1026.5614807808624, 151.29862440411648, 3.901771013880543, 3.7301247955449974, 3.8759621011700065, 3.659557126210126, 4.086835266824942, 3.6340888739813773, 3.858324504680583, 3.7551828212045004, 3.815760524425657, 4.2259033713728975, 4.371737983114335, 4.307380361442847, 4.187602191775825, 4.437147867992696, 4.356108312396685, 4.298134322136188, 4.318494549060967, 4.144414878295289, 4.064397655100121, 4.201834039385025, 4.090374027710338, 3.9812634489272587, 4.165582641164475, 4.2508450891419445, 4.1992154570223095, 4.118323968726959, 4.254132727086158, 4.0768109644365, 4.334334559710399, 4.162406395217695, 4.293655162083104, 4.047267856390056, 1.9404908700418344, 1.7322784562657245, 1.4262722646573889, 1.4002468768220941, 1.4188457776352552, 1.4942248325184317, 1.4610398206721549, 1.4946762633477086, 1.579364101613276, 1.628202105324122, 1.8138519188905389, 2.0332905300668207, 2.119034438365694, 1.9823157691427975, 1.902569630066241, 2.0262156931632433, 1.894360004246963, 1.9326606197284226, 1.985607969208476, 2.055836884050246, 2.069751067566363, 2.2360626692816243, 2.3146263492456645, 2.3107996803407076, 2.42964276340863, 2.700027303300703, 2.6020378567303237, 2.6780875422740755, 2.7304187730744607, 2.7806438474301753, 2.8686015472805813, 3.072208621791253, 3.0058409222007043, 2.941624606650066, 2.934913046036969, 2.813442674066407, 2.8614137527887276, 2.980381973943654, 3.0365828033385696, 3.0562505983088197, 3.00128116217715, 2.919862179921546, 3.32580923458729, 3.3279063002519633, 3.2173439434520508, 3.1882897757128092, 3.2360860950335164, 3.286528071204994, 3.349785053135752, 3.3537105645306524, 3.460444914991341, 3.289072289917036, 3.503778532597925, 3.442636619561169, 3.4800351828883778, 3.6782473459680296, 3.5297646841783856, 3.5130960805312323, 3.7128028254228576, 3.6902510915178013, 3.5840906974377296, 3.813603765807182, 3.737248002543765, 3.8420491813898714, 3.8867667888319763, 3.9041698223314447, 3.86073056841989, 3.8209979056671757, 4.081339568377051, 3.894585159244175, 3.8889170848112498, 4.0731343393503705, 3.9617419922521324, 4.001120438466195, 4.067291872336939, 4.100552386247567, 4.130323647660298, 4.145861854944665, 4.183502711005124, 4.204352540879397, 4.311600564251246, 4.224798341878524, 4.203110275563693, 4.235702653397426, 4.306961930433635, 4.3160658380236105, 4.23257363292066, 4.287429097210338, 4.307907571422263, 4.256150607357722, 4.415548127206195, 4.532412710834881, 4.30850977965567, 4.515983885502413, 4.531079288338962, 4.507612869099098, 4.603585119686571, 4.610085159300437, 4.409948478389974, 4.346258230693706, 4.40341434713326, 4.646708086233395, 4.293545427600795, 4.488340299468738, 4.359718588290123, 4.492473163733969, 4.557691540875151, 4.318873893158243, 4.734451002540511, 4.5897528897842435, 4.879455374451865, 4.654380060994856, 4.704406091735462, 4.973844823291304, 4.603959370273221 ], "yaxis": "y3", "zmax": 0.1, "zmin": -0.1 }, { "colorbar": { "title": "" }, "hoverinfo": "text", "legendgroup": "This work", "line": { "color": "rgba(0, 0, 0, 1.000)", "dash": "dash", "shape": "linear", "width": 1 }, "mode": "lines", "name": "This work", "showlegend": true, "text": [ null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ], "type": "scatter", "x": [ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160 ], "xaxis": "x3", "y": [ 1.0857790679367962, 1.002480959875828, 1.0175758856743218, 1.117979201944521, 1.1860396502812893, 1.1931514036238238, 1.1704642716048816, 1.1207859454857054, 1.2723664436364805, 84.49191304279364, 80.31750108389734, 842.8986886296851, 206.41371313703783, 4.1230410294249635, 3.7978766186593984, 3.8766949117675593, 3.6872019479634157, 4.14782768042253, 3.760214925849103, 3.8232645414402566, 3.7906398901878218, 3.847668877273603, 4.374300551915247, 4.4197699356296285, 4.405286530242607, 4.27631611943179, 4.524543175030498, 4.404658280542666, 4.427456015697982, 4.299595128721327, 4.245655827104088, 4.166581523856618, 4.317817995227425, 4.137981544346955, 4.027347631245824, 4.263145642963551, 4.253198620692602, 4.23870920846987, 4.14595560706737, 4.30987462720298, 4.10950743757508, 4.3849533609624665, 4.298607851503398, 4.411629731372064, 4.130104752146885, 2.014013673108549, 1.7899167339767463, 1.4985958635818057, 1.4644903601248604, 1.5114576718625015, 1.5079959196873594, 1.5678906277340756, 1.5439091053363423, 4.605777709038979, 1.7144956422946347, 1.9199837445389945, 2.13138883789023, 2.2395944066174978, 2.1066600560736926, 2.019937597534639, 2.158333387812875, 1.9794869608610792, 2.033774846206659, 2.088768110455931, 2.1593152640300355, 2.1504908250245975, 2.30145874074324, 2.3704693805466173, 2.4364360262938187, 2.5006308712520235, 2.8217434524802307, 2.677185797618641, 2.752029904582889, 2.854177758636442, 2.9097410843454368, 2.971997154721216, 3.130711867585389, 3.1192124235421947, 3.094158611329475, 3.1323889894674144, 2.9980865477136103, 3.017380147538852, 3.114126235840075, 3.1847493708438277, 3.1405786909221645, 3.179574128537354, 3.043685629590943, 3.4227395890782564, 3.521603661669303, 3.3714117402482704, 3.398161371277994, 3.3661263968308788, 3.4335689633006203, 3.44173991301157, 3.518830168177832, 3.6088062634567333, 3.4114979261961675, 3.679133353679267, 3.54108997987351, 3.7306290544108807, 3.7447579006386014, 3.770190146012226, 3.6999024091455044, 3.9095729971261086, 3.8182391601391474, 3.7217248939880387, 3.9544250813661908, 3.8979232432295974, 3.9783179793580854, 4.080221723237815, 4.140085681327101, 4.003383050244376, 4.035889220723589, 4.153505830012789, 4.104076691736123, 4.073618978698611, 4.246879989248927, 4.124424756521174, 4.2401967117291886, 4.20930769078522, 4.345464934643725, 4.262678572167572, 4.213058074907197, 4.298605649711017, 4.452242917787, 4.592950973359897, 4.328941650649829, 4.316890741721304, 4.3729451426185735, 4.402521373706199, 4.351492047767091, 4.429301711906779, 4.508639054871161, 4.4860204127156456, 4.368184026885298, 4.63239595780465, 4.679366253248214, 4.478440778244224, 4.686761127951312, 4.653635715428732, 4.597897612425877, 4.710732745178608, 4.682686059854525, 4.510711363333386, 4.540099667166455, 4.652505245649606, 4.854929239221948, 4.496427172907768, 4.663530671727717, 4.554855879437058, 4.609797429033649, 4.7812987697637865, 4.395722554679775, 4.894705247142196, 4.7077447157422885, 4.885332610946976, 4.835121580894191, 4.870532808946831, 5.144704221001939, 4.80466985043328 ], "yaxis": "y3", "zmax": 0.1, "zmin": -0.1 } ], "layout": { "annotations": [], "height": 432, "legend": { "bgcolor": "rgba(255, 255, 255, 1.000)", "bordercolor": "rgba(0, 0, 0, 1.000)", "font": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 17 }, "tracegroupgap": 0, "x": 1, "y": 1 }, "margin": { "b": 20, "l": 0, "r": 0, "t": 20 }, "paper_bgcolor": "rgba(255, 255, 255, 1.000)", "plot_bgcolor": "rgba(255, 255, 255, 1.000)", "showlegend": true, "width": 864, "xaxis1": { "anchor": "y1", "domain": [ 0.1292438271604938, 0.9954432779235928 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ -2, 165 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "", "", "", "" ], "tickvals": [ 0, 50, 100, 150 ], "title": "", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "type": "-", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" }, "xaxis2": { "anchor": "y2", "domain": [ 0.1292438271604938, 0.9954432779235928 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ -2, 165 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "", "", "", "" ], "tickvals": [ 0, 50, 100, 150 ], "title": "", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "type": "-", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" }, "xaxis3": { "anchor": "y3", "domain": [ 0.1292438271604938, 0.9954432779235928 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ -2, 165 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "0", "50", "100", "150" ], "tickvals": [ 0, 50, 100, 150 ], "title": "Scan number", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "type": "-", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" }, "yaxis1": { "anchor": "x1", "domain": [ 0.7498285322359396, 0.9908865558471858 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ 0, 2500 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "0", "1000", "2000" ], "tickvals": [ 0, 1000, 2000 ], "title": "N (cm⁻³)", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "type": "-", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" }, "yaxis2": { "anchor": "x2", "domain": [ 0.4314986282578875, 0.6725566518691336 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ 0, 30 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "0", "10", "20" ], "tickvals": [ 0, 10, 20 ], "title": "A (μm² cm⁻³)", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "type": "-", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" }, "yaxis3": { "anchor": "x3", "domain": [ 0.11316872427983533, 0.35422674789108144 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ 0, 5 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "0", "2", "4" ], "tickvals": [ 0, 2, 4 ], "title": "V (μm³ cm⁻³)", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "type": "-", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" } } }, "image/svg+xml": [ "01000200001020050100150024TSI AIMTSI AIMTSI AIMThis workThis workThis workScan numberN (cm⁻³)A (μm² cm⁻³)V (μm³ cm⁻³)" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "figure(\"Nimbus Sans L\", 2, 6, 3, 10)\n", "\n", "p1 = plot(Nt, ylabel = \"N (cm⁻³)\", label = \"TSI AIM\",ylim = (0,2500), color = RGBA(0.6,0,0.6,1))\n", "p1 = plot!(Ninv, ylabel = \"N (cm⁻³)\", label = \"This work\", ytick = [0,1000, 2000],\n", " color = RGBA(0,0,0.6,1), ls = :dash, xticks = ([0,50,100,150],[\"\",\"\",\"\",\"\"]))\n", "\n", "p2 = plot(At, ylabel = \"A(μm² cm⁻³)\", label = \"TSI AIM\", ylim = (0,30), color = RGBA(0.8,0,0,1))\n", "p2 = plot!(Ainv, ylabel = \"A (μm² cm⁻³)\", label = \"This work\", ytick = [0, 10, 20],\n", " color = RGBA(0,0,0,1), ls = :dash, xticks = ([0,50,100,150],[\"\",\"\",\"\",\"\"]))\n", "\n", "p3 = plot(Vt, ylabel = \"V(μm³ cm⁻³)\", label = \"TSI AIM\", ylim = (0,5), color = RGBA(0.0,0,0.8,1))\n", "p3 = plot!(Vinv, ylabel = \"V (μm³ cm⁻³)\", label = \"This work\", xlabel = \"Scan number\", \n", " color = RGBA(0,0,0,1), ls = :dash, ytick = [0,2,4],xticks = [0,50,100,150])\n", "\n", "\n", "plot(p1,p2,p3, layout = grid(3,1), left_margin = 30px, bottom_margin = 10px, xlim = (-2,165), fmt = :svg)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " Figure 3. Top to bottom: integrated number concentration, surface area, and volume concentration from TSI AIM software provided dn/dlnD distribution and from the inversion of this work. Note that for scan 54, the agreement is poor. The spectrum can be visualized using the block above if desired. The disagreement is due to failed convergence of the L-curve algorithm. This can be fixed by changing the bounds $\\lambda_1$ and $\\lambda_2$ in the lcorner() call or perhaps a more robust L-curve algorithm. It is left in place here to alert potential users of this issue. It is apparent from graphing the inverted spectra from scan 54 that inversion failed." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### (c) Comparison of relative error" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "application/vnd.plotly.v1+json": { "data": [ { "colorbar": { "title": "" }, "hoverinfo": "text", "legendgroup": "N: Rel. Error %", "line": { "color": "rgba(153, 0, 153, 1.000)", "dash": "solid", "shape": "linear", "width": 0 }, "marker": { "color": "rgba(153, 0, 153, 1.000)", "line": { "color": "rgba(0, 0, 0, 1.000)", "width": 1 }, "size": 4, "symbol": "square" }, "mode": "lines+markers", "name": "N: Rel. Error %", "showlegend": true, "text": [ null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ], "type": "scatter", "x": [ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160 ], "xaxis": "x1", "y": [ -1.660318277591944, -1.8593312260832588, -1.5240951103831595, -1.6895979827812373, -2.241349270988853, -0.966696313941989, -1.984735199248788, -1.9326643614033545, -1.4961555536199227, 1.3465903009755786, -0.3412146166889112, -3.637813215751962, 1.0325994692684717, 17.541057388786363, -1.0914714289751173, -0.46270792622195506, 0.12242675775192711, -0.7698685328861921, -0.07467165579586063, -1.1696366909450076, -0.885662983648676, -0.9136490607434622, -0.568712413554262, -1.302401203175433, -1.3023697208010727, -1.2181031191175105, -1.6423922917483953, -0.029409903906645773, -1.3388306957566427, -0.5864897045099308, -1.3781674682059786, -0.9139910166809772, -0.8757139747114461, -1.0720336679273739, -0.6577979535416832, -1.001992320902297, -0.07189872620091141, -1.1008624980289068, 0.04821836228356262, -0.6765276923643314, -1.0778056635613105, -0.8238162585797215, -1.2897137291749947, -0.6063959587862235, -0.7107134716697562, -0.5726282335019798, 2.5841672014294614, 1.2225571938908697, 3.1341495624436777, 3.8701282491815503, 6.051659195621092, 1.282231376656112, 3.6763733009534683, 57.02323353884159, 3.503825618669021, 3.04611948432742, 1.8868602269029502, 4.1454214839954, 1.7548602350398663, 4.846794801611495, 1.6366129688139917, 3.4562957817468405, 1.855910633520915, 1.5957897802043564, -0.06529358995734796, 0.21993727397254914, 1.6374440253528175, 0.6999824218390548, -1.134419971961622, 0.29804255120939344, -0.4693055886989369, 1.0687825526824317, 0.7443246196575856, -0.12573239180620077, -0.19210883507133114, -0.9632949935086542, -0.5538423266912507, -1.0200756399649675, 0.8417124390263507, 0.5306227909101232, 2.842610064510346, 0.1495820206106796, 1.4859815126643943, 0.9436111756624106, 0.36982171855222007, 0.6093339883549583, 0.4172038128387305, 2.0626983594200743, 0.6083669007430297, 3.252884602693624, -0.5343249450831283, -0.4377389731484359, -1.127459744458756, -0.010876925716682995, 1.0375404218247597, -0.4407564114230343, 0.7559600641283791, -0.30272910865892927, 1.1105582314599907, -1.5967441130702795, -1.013797494169551, -0.2141867604325844, -0.4997578230943499, 0.10697605495714017, -0.3242349477305231, 0.5596117303756379, -0.8333444350582558, 0.4295148244536892, -0.6328810159344723, -0.49450723223335735, -0.2868858511505929, -0.8193511645125489, -0.3398278256760441, -0.9552400793843895, 0.5084211654642833, -0.8197943598654129, -0.5109173187116867, -0.43308537713237527, -0.635211007786198, -0.2444799136118197, -0.5590785650256079, -0.7803880485964478, -0.6394853532561784, 0.5496061377059849, -0.9281765924217241, -0.3285308427793367, -1.0576055544252319, -0.8898741848830262, -0.32994287353300483, -0.9309039548981117, -0.38483271442911704, -0.9752588494838472, -1.180879273271688, -0.7090575808185635, -0.9204718940141636, -1.153499091790014, -0.8108989599167122, -1.2451884942299851, -1.3880232159578847, -0.46149056329775917, -1.27631107422677, 0.5993802791331048, -0.8674467340702647, -1.707641035161736, -1.3678742456683446, -0.9368138844111037, 0.05410467255078257, -0.761592820417076, 0.26384992829714576, -0.5599026423229893, -0.30809056356589265, -0.9397948632485961, -1.0587910453137261, -1.5575039541533064, -0.9381140742429119, -0.44432198028300546, -1.0442040518453621, -0.21678944980927858, -0.7688932607572219, -0.6144786086292209 ], "yaxis": "y1", "zmax": 0.1, "zmin": -0.1 }, { "colorbar": { "title": "" }, "hoverinfo": "text", "legendgroup": "A: Rel. Error %", "line": { "color": "rgba(204, 0, 0, 1.000)", "dash": "solid", "shape": "linear", "width": 0 }, "marker": { "color": "rgba(204, 0, 0, 1.000)", "line": { "color": "rgba(0, 0, 0, 1.000)", "width": 1 }, "size": 4, "symbol": "square" }, "mode": "lines+markers", "name": "A: Rel. Error %", "showlegend": true, "text": [ null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ], "type": "scatter", "x": [ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160 ], "xaxis": "x2", "y": [ 1.5191249059867329, 0.9497142185054716, 1.0464312982794077, 1.425773509080578, 0.28115987058307995, 0.2009886737127764, -1.682514787654209, 0.417885382960076, 2.1281038738191462, -10.205429360745594, -11.318695863054312, -11.596325418595978, 14.163215590956513, 0.9560358458107765, 1.2700629057357142, 0.5780042497605767, 0.8711515580710656, 1.26648767084186, 1.8548765880386995, 0.05152744919999445, 0.9844514924312515, 0.8679245032156491, 1.8400956338307997, 0.6302039789350433, 1.1927581902158686, 1.2373448322392253, 1.1138712698883628, 0.8201900055018758, 1.5354520842456565, 0.3658123661816237, 1.5702661669109252, 1.6536601897692933, 1.7469848516134792, 0.8922606637638062, 1.1563749683873483, 1.544356740606534, 0.4054027195354519, 0.8121671404500348, 0.7893001162430482, 1.0212286179119605, 0.8647252368197385, 0.954894211125858, 1.4628726543868613, 1.6406649637792512, 1.395311487792208, 1.54819628691225, 1.1854072979304873, 1.350958072289174, 1.518820616773496, 3.3624617600432107, 0.0789661940348406, 2.4078052749676844, 1.3671686340819884, 116.99386832385454, 2.0205597467333343, 2.0786881905201153, 2.103443774599813, 2.248350302643639, 2.3934159712351617, 2.361115187077139, 2.08371651061659, 1.9608067051489115, 2.11179578111307, 2.3217736192521885, 2.15154347635608, 1.8376870994253265, 1.5255235730131187, 1.341315585351345, 2.406832283581197, 1.616527362280058, 2.363206000571721, 1.6703041490784292, 1.6133389491042132, 2.427092362834535, 2.3621108713161636, 2.0057248011837174, 1.374721204146157, 2.1554699275930265, 2.549188236490884, 3.0763168070647406, 3.0012345775049893, 2.522674802791724, 2.180518176978729, 2.408471191152948, 1.559433241949565, 2.9525686760144563, 2.198458125637647, 1.652840642667173, 2.5954205782425763, 2.3372606513910252, 2.961779832350239, 2.2379664303482234, 2.269673397860107, 1.851424619012547, 2.4642436442393776, 2.0846894111732506, 2.0588369232149275, 2.292752844932348, 1.7300934503697682, 3.0130042627249334, 1.2779809886956126, 2.8355123598957577, 2.6746632982143614, 2.815216429427335, 1.8767827748507415, 2.1070331052286884, 1.8377872697242457, 2.3181995976214975, 2.106665662075878, 2.549972957759646, 2.874892431423293, 2.05900541905371, 2.962420433361943, 1.2625269934026488, 2.8232865193779433, 2.533348149555371, 2.094661231504617, 2.1909558839019128, 2.933890025902709, 1.9317481462683521, 2.8140848548569557, 1.895542570352415, 1.1576566059084392, 1.7437725620817162, 2.812430008503729, 3.475466061085758, 1.6127189924666037, 1.71015464745657, 1.7718680117847834, 1.2889054774027944, 0.9525470461530934, 2.2352541747582166, 2.606509514287358, 2.081246454761547, 1.6005541347384595, 2.5387664065008986, 2.1266302507102486, 2.0255233802884094, 1.9863097287329374, 1.6555638969007316, 1.3793734170945218, 1.3952201959067734, 1.2837639594545283, 1.5614997840475047, 2.32772696024002, 2.906773360551159, 2.2510696603268, 2.406115632449416, 2.0180109515575553, 2.334646582834211, 1.545339727578808, 2.4828432845705604, 1.3654558707193813, 2.0592998784313266, 1.645419230516379, 0.6988475093680551, 2.0976708254467784, 1.9794804630912932, 1.9696346373225493, 2.5880969706742443 ], "yaxis": "y2", "zmax": 0.1, "zmin": -0.1 }, { "colorbar": { "title": "" }, "hoverinfo": "text", "legendgroup": "V: Rel. Error %", "line": { "color": "rgba(0, 0, 204, 1.000)", "dash": "solid", "shape": "linear", "width": 0 }, "marker": { "color": "rgba(0, 0, 204, 1.000)", "line": { "color": "rgba(0, 0, 0, 1.000)", "width": 1 }, "size": 4, "symbol": "square" }, "mode": "lines+markers", "name": "V: Rel. Error %", "showlegend": true, "text": [ null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null, null ], "type": "scatter", "x": [ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160 ], "xaxis": "x3", "y": [ 4.622758410396836, 2.7382325778905687, 2.8015074547964436, 3.4857526037381446, 1.9658749432748712, 1.2546924391717722, -2.4163129651135313, 1.5632340360727914, 4.659902526565801, -15.759111821181481, -17.540551682514604, -17.891066009166124, 36.4280170754955, 5.671014899573882, 1.8163419946517356, 0.01890654703077735, 0.7554144067131728, 1.4924118447517662, 3.4706375171707533, -0.9086836319183332, 0.944216850990713, 0.8362252464139851, 3.511608465721695, 1.0986923896357639, 2.2729863765029625, 2.11848985632384, 1.969628005147789, 1.1145261932036072, 3.008786693700174, -0.43763909216347946, 2.442828524214776, 2.5141208471143046, 2.760317393672576, 1.1638915246894037, 1.1575265718972418, 2.3421213838120423, 0.055366203691326894, 0.9405030975849397, 0.6709437759203909, 1.3102999763479948, 0.8020110185094076, 1.1678563469140586, 3.272180641519985, 2.7476489106713413, 2.0467361858949236, 3.788876526130201, 3.3273101967262666, 5.070812965839313, 4.588011183325684, 6.527269960347595, 0.921620820992328, 7.313339824835423, 3.293879965576234, 191.62228673770076, 5.299927858362181, 5.851184682891435, 4.824608503939936, 5.68938220488694, 6.2726780902653365, 6.168918373006508, 6.520416118353866, 4.493705337067409, 5.231866652948364, 5.195393191767747, 5.033394467362824, 3.900940491030616, 2.9246081677408977, 2.412615380411271, 5.436920691220934, 2.921750839773742, 4.507959939172955, 2.888041797468195, 2.7610136390846214, 4.532600888274305, 4.642710249806747, 3.6043907017568566, 1.9042732117594858, 3.771706629720322, 5.185366084257601, 6.728510873502661, 6.562915795271155, 5.450676072208009, 4.4874872773254, 4.879385055541923, 2.7592008541448725, 5.94056193758499, 4.24072925499258, 2.9144893063295227, 5.820397088784458, 4.788664174676721, 6.582575936601109, 4.018443823139864, 4.4740494804815185, 2.7450973246697883, 4.923489981321267, 4.287348942405092, 3.7221935393283636, 5.00473472994661, 2.859824349538546, 7.200900518325038, 1.8082132171855811, 6.811373656479409, 5.317427258807683, 5.2997743471831305, 3.4682753408171103, 3.840142679667789, 3.692604796062302, 4.299292971097144, 3.5467739098258604, 4.977271467938385, 6.042663862781854, 3.694960818850162, 5.623957938780692, 1.7682003770256032, 5.379046135239829, 4.749442836123764, 4.26564987606737, 4.106344244203577, 5.975233111319234, 3.4916554529607318, 5.9726721018745135, 3.204468603380503, 1.6208021953840197, 2.751353271580367, 5.896041649631817, 6.525428432341578, 2.465048041204202, 2.7070540313708764, 3.240135119283373, 2.2187204070072224, 0.8207986410073559, 4.6479540829716735, 5.15950124527436, 4.134555775405781, 2.632271032277392, 4.911005935194242, 3.242280696593151, 3.944078284119263, 3.781617622621344, 2.7047954646297483, 2.0029391597868718, 2.327482227575976, 1.5748277536180166, 2.2848993687155237, 4.459961331883642, 5.65676720108139, 4.481046562951552, 4.725273057617144, 3.903232833742917, 4.4759148370506985, 2.6115740934591276, 4.906150995152674, 1.7793680348776062, 3.3848537985859863, 2.570766417962681, 0.12044861657886843, 3.883256578335861, 3.531300529161681, 3.4351573838922804, 4.359518927469338 ], "yaxis": "y3", "zmax": 0.1, "zmin": -0.1 } ], "layout": { "annotations": [], "height": 432, "legend": { "bgcolor": "rgba(255, 255, 255, 1.000)", "bordercolor": "rgba(0, 0, 0, 1.000)", "font": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 17 }, "tracegroupgap": 0, "x": 1, "y": 1 }, "margin": { "b": 20, "l": 0, "r": 0, "t": 20 }, "paper_bgcolor": "rgba(255, 255, 255, 1.000)", "plot_bgcolor": "rgba(255, 255, 255, 1.000)", "showlegend": true, "width": 864, "xaxis1": { "anchor": "y1", "domain": [ 0.0932355967078189, 0.9954432779235929 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ -2, 165 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "", "", "", "", "" ], "tickvals": [ 0, 50, 100, 150 ], "title": "", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "type": "-", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" }, "xaxis2": { "anchor": "y2", "domain": [ 0.0932355967078189, 0.9954432779235929 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ -2, 165 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "", "", "", "", "" ], "tickvals": [ 0, 50, 100, 150 ], "title": "", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "type": "-", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" }, "xaxis3": { "anchor": "y3", "domain": [ 0.0932355967078189, 0.9954432779235929 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ -2, 165 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "0", "50", "100", "150" ], "tickvals": [ 0, 50, 100, 150 ], "title": "Scan number", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "type": "-", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" }, "yaxis1": { "anchor": "x1", "domain": [ 0.7357938282406058, 0.9908865558471858 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ -5, 5 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "-4", "-2", "0", "2", "4" ], "tickvals": [ -4, -2, 0, 2, 4 ], "title": "N (cm⁻³)", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "type": "-", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" }, "yaxis2": { "anchor": "x2", "domain": [ 0.4174639242625537, 0.6725566518691337 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ -5, 5 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "-4", "-2", "0", "2", "4" ], "tickvals": [ -4, -2, 0, 2, 4 ], "title": "A(μm² cm⁻³)", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "type": "-", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" }, "yaxis3": { "anchor": "x3", "domain": [ 0.09913402028450162, 0.35422674789108166 ], "gridcolor": "rgba(0, 0, 0, 0.100)", "gridwidth": 0.5, "linecolor": "rgba(0, 0, 0, 1.000)", "mirror": false, "range": [ -5, 5 ], "showgrid": true, "showline": true, "showticklabels": true, "tickangle": 0, "tickcolor": "rgb(0, 0, 0)", "tickfont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "tickmode": "array", "ticks": "inside", "ticktext": [ "-4", "-2", "0", "2", "4" ], "tickvals": [ -4, -2, 0, 2, 4 ], "title": "V(μm³ cm⁻³)", "titlefont": { "color": "rgba(0, 0, 0, 1.000)", "family": "Nimbus Sans L", "size": 20 }, "type": "-", "visible": true, "zeroline": false, "zerolinecolor": "rgba(0, 0, 0, 1.000)" } } }, "image/svg+xml": [ "-4-2024-4-2024050100150-4-2024N: Rel. Error %A: Rel. Error %V: Rel. Error %Scan numberN (cm⁻³)A(μm² cm⁻³)V(μm³ cm⁻³)" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "figure(\"Nimbus Sans L\", 2, 6, 3, 10)\n", "\n", "p1 = plot(100*(Ninv - Nt)./Nt, ylabel = \"N (cm⁻³)\", \n", " label = \"N: Rel. Error %\", color = RGBA(0.6,0,0.6,1), marker = :square,\n", " xticks = ([0,50,100,150],[\"\",\"\",\"\",\"\",\"\"]))\n", "p2 = plot(100*(Ainv - At)./At, ylabel = \"A(μm² cm⁻³)\", \n", " label = \"A: Rel. Error %\", color = RGBA(0.8,0,0,1), marker = :square,\n", " xticks = ([0,50,100,150],[\"\",\"\",\"\",\"\",\"\"]))\n", "p3 = plot(100*(Vinv - Vt)./Vt, ylabel = \"V(μm³ cm⁻³)\", xlabel = \"Scan number\",\n", " label = \"V: Rel. Error %\", color = RGBA(0,0,0.8,1), marker = :square)\n", "\n", "plot(p1,p2,p3, layout = grid(3,1), left_margin = 30px, ylim = (-5,5), \n", " ms = 2, lw = 0, ytick = [-4, -2, 0, 2, 4], xlim = (-2,165), fmt = :svg)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " Figure 4. Top to bottom: relative error in integrated number concentration, surface area, and volume concentration between TSI AIM software provided dn/dlnD distribution and from the inversion of this work." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Summary and Conclusions\n", "Figures 2-4 demonstrate that the code can be used to invert size distribution data. Figure 2 shows that the inverted spectra is in good qualitative agreement with the inverted spectra from the TSI software. The error of the regularization routine here appears to be similar to that of the TSI algorithm. Number surface area, and volume concentration agree within ~2%, 2%, and 4% for the entire dataset. This agreement is comparable to the range given by Wiedensohler et al. (2012), who compared the performance of several commercial and non-commercial inversion routines. " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Acknowlegements and References\n", "\n", "The data were collected by Drs. Andrew Grieshop and Provat Saha during the 2013 SOAS campaign (Saha et al., 2017). I thank Andrew Grieshop for sharing the datafile.
\n", "\n", "Collins, A. M., W. D. Dick & F. J. Romay (2013) A New Coincidence Correction Method for Condensation Particle Counters, Aerosol Science and Technology, 47:2, 177-182, DOI:10.1080/02786826.2012.737049\n", "\n", "Reineking A. & J. Porstendörfer (1986) Measurements of Particle Loss Functions in a Differential Mobility Analyzer (TSI, Model 3071) for Different Flow Rates, Aerosol Science and Technology, 5:4, 483-486, DOI:10.1080/02786828608959112\n", "\n", "Saha, P. K., Khlystov, A., Yahya, K., Zhang, Y., Xu, L., Ng, N. L., and Grieshop, A. P.: Quantifying the volatility of organic aerosol in the southeastern US, Atmos. Chem. Phys., 17, 501-520, https://doi.org/10.5194/acp-17-501-2017, 2017. \n", "\n", "Wiedensohler, A., Birmili, W., Nowak, A., Sonntag, A., Weinhold, K., Merkel, M., Wehner, B., Tuch, T., Pfeifer, S., Fiebig, M., Fjäraa, A. M., Asmi, E., Sellegri, K., Depuy, R., Venzac, H., Villani, P., Laj, P., Aalto, P., Ogren, J. A., Swietlicki, E., Williams, P., Roldin, P., Quincey, P., Hüglin, C., Fierz-Schmidhauser, R., Gysel, M., Weingartner, E., Riccobono, F., Santos, S., Grüning, C., Faloon, K., Beddows, D., Harrison, R., Monahan, C., Jennings, S. G., O'Dowd, C. D., Marinoni, A., Horn, H.-G., Keck, L., Jiang, J., Scheckman, J., McMurry, P. H., Deng, Z., Zhao, C. S., Moerman, M., Henzing, B., de Leeuw, G., Löschau, G., and Bastian, S.: Mobility particle size spectrometers: harmonization of technical standards and data structure to facilitate high quality long-term observations of atmospheric particle number size distributions, Atmos. Meas. Tech., 5, 657-685, https://doi.org/10.5194/amt-5-657-2012, 2012. " ] } ], "metadata": { "kernelspec": { "display_name": "Julia 1.1.0", "language": "julia", "name": "julia-1.1" }, "language_info": { "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", "version": "1.1.0" } }, "nbformat": 4, "nbformat_minor": 2 }