Revision 66ff08f00acc06131fe610be0f9878a6c78bfe44 authored by Software Heritage on 10 January 2019, 12:27:59 UTC, committed by Software Heritage on 10 January 2019, 12:27:59 UTC
0 parent
plotPreHarmo.sce
/*
This file is part of libszdist.
libszdist is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
libszdist is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with libszdist. If not, see <http://www.gnu.org/licenses/>.
*/
// Compares results of inversion algorithms : libszdist vs. IfT, with data
// circulated before the "hamozition paper"
//
// testLib -i ../../files/ift/Le_ambient_080913_dat_STD.raw -o ../../files/out/preharmo.csv -t ift
info="Before harmonization paper, 2008-09-13 ambient, eq len = 0 m, flow ratio = 1:10";
exec('plotResults.sci');
scf();
//-----------------------------------------------------------------------------
subplot(2,1,1);
scRef=avgScan("../files/ift/dNdLogDp_NG2.csv","ift");
plotScan(scRef,color("green"),"IfT");
errL=[scRef(:,1), scRef(:,2)*0.90];
plotScan(errL,color("grey"),"IfT -10%");
errH=[scRef(:,1), scRef(:,2)*1.10];
plotScan(errH,color("grey"),"IfT +10%");
scTst=avgScan("../files/out/preharmo.csv","lamp");
plotScan(scTst,color("red"),"LaMP (libszdist 2013-10)");
xgrid(2);
a=gca();
a.grid=[color("grey"),color("grey"),-1];
a.data_bounds(1,1)=10;
a.data_bounds(2,1)=1000;
xtitle(info,"Dp (nm)","dN/dLog");
//-----------------------------------------------------------------------------
subplot(2,1,2);
plotRatio(scTst,scRef,color("red"),"LaMP (libszdist 2013-10) / IfT");
xgrid(2);
a=gca();
a.grid=[color("grey"),color("grey"),-1];
a.data_bounds(1,1)=10; // xmin
a.data_bounds(2,1)=1000; // xmax
a.data_bounds(1,2)=0.8; // ymin
a.data_bounds(2,2)=1.2; // ymax
xtitle("","Dp (nm)","ratio");
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...