Skip to content

Commit 2469821

Browse files
HolgerPollyNetHolgerPollyNet
committed
raman bsc error calculation in smart mode
1 parent cbd5270 commit 2469821

2 files changed

Lines changed: 131 additions & 50 deletions

File tree

‎lib/interface/picassoProcV3.m‎

Lines changed: 19 additions & 50 deletions
Original file line numberDiff line numberDiff line change
@@ -439,6 +439,10 @@
439439
flag532RR = data.flag532nmChannel & data.flagRotRamanChannel;
440440
flag1064RR = data.flag1064nmChannel & data.flagRotRamanChannel;
441441

442+
%%%% to be put in config file:
443+
sigma_angstroem=0.2
444+
MC_count=3;
445+
442446
if isempty(PollyConfig.H)
443447
print_msg('Using transmission ratios instead of GHK.\n', 'flagTimestamp', true);
444448
flagGHK = 0;
@@ -1569,10 +1573,6 @@
15691573
bg387 = squeeze(sum(data.bg(flag387FR, :, flagClFre), 3));
15701574

15711575
[thisAerExt355_raman, thisAerExtStd355_raman] = pollyRamanExt_smart_MC(data.distance0, sig387, 355, 387, mExt355(iGrp,:), mExt387(iGrp,:), number_density(iGrp, :), PollyConfig.angstrexp, PollyConfig.smoothWin_raman_355, 'moving',15,bg387);
1572-
% thisAerExt355_raman = pollyRamanExt_smart(data.distance0, sig387, 355, 387, mExt355(iGrp,:), mExt387(iGrp,:), number_density(iGrp, :), PollyConfig.angstrexp, PollyConfig.smoothWin_raman_355, 'moving');
1573-
% thisAerExt355_raman = pollyRamanExt(data.distance0, sig387, 355, 387, PollyConfig.angstrexp, data.pressure(iGrp, :), data.temperature(iGrp, :) + 273.17, PollyConfig.smoothWin_raman_355, 380, 70, 'moving');
1574-
1575-
%thisAerExtStd355_raman = pollyRamanExtStd(data.distance0, sig387, bg387, 355, 387, PollyConfig.angstrexp, data.pressure(iGrp, :), data.temperature(iGrp, :) + 273.17, PollyConfig.smoothWin_raman_355, 380, 70, 15);
15761576
data.aerExt355_raman(iGrp, :) = thisAerExt355_raman;
15771577
data.aerExtStd355_raman(iGrp, :) = thisAerExtStd355_raman;
15781578

@@ -1595,14 +1595,12 @@
15951595
continue;
15961596
end
15971597
%here the lower end of the exitncion profiles is set to contant values
1598-
%according to the value at fullüverlap + smoothing window/2 --> should be
1598+
%according to the value at full overlap + smoothing window/2 --> should be
15991599
%the mean value in future?
16001600
thisAerExt355_raman_tmp = thisAerExt355_raman;
16011601
thisAerExt355_raman(1:hBaseInd355) = thisAerExt355_raman(hBaseInd355);
1602-
[thisAerBsc355_raman, ~] = pollyRamanBsc_smart(data.distance0, sig355, sig387, thisAerExt355_raman, PollyConfig.angstrexp, mExt355(iGrp,:), mBsc355(iGrp,:),mExt387(iGrp,:), mBsc387(iGrp,:), refH355, 355, PollyConfig.refBeta355, PollyConfig.smoothWin_raman_355, true, 355, 387);
1603-
%[thisAerBsc355_raman, ~] = pollyRamanBsc(data.distance0, sig355, sig387, thisAerExt355_raman, PollyConfig.angstrexp, mExt355(iGrp,:), mBsc355(iGrp,:), refH355, 355, PollyConfig.refBeta355, PollyConfig.smoothWin_raman_355, true);
1604-
thisAerBscStd355_raman = pollyRamanBscStd(data.distance0, sig355, bg355, sig387, bg387, thisAerExt355_raman, thisAerExtStd355_raman, PollyConfig.angstrexp, 0.2, mExt355(iGrp,:), mBsc355(iGrp,:), refH355, 355, PollyConfig.refBeta355, PollyConfig.smoothWin_raman_355, true);
1605-
1602+
[thisAerBsc355_raman, thisAerBscStd355_raman, ~] = pollyRamanBsc_smart_MC(data.distance0, sig355, sig387, thisAerExt355_raman, PollyConfig.angstrexp, mExt355(iGrp,:), mBsc355(iGrp,:),mExt387(iGrp,:), mBsc387(iGrp,:), refH355, PollyConfig.refBeta355, PollyConfig.smoothWin_raman_355, true, 355, 387, bg355, bg387, thisAerExtStd355_raman, sigma_angstroem, MC_count, 'monte-carlo');
1603+
16061604
% lidar ratio
16071605
[thisLR355_raman, thisLRStd355_raman] = pollyLR(thisAerExt355_raman_tmp, thisAerBsc355_raman, ...
16081606
'hRes', data.hRes, ...
@@ -1615,7 +1613,6 @@
16151613
data.LRStd355_raman(iGrp, :) = thisLRStd355_raman;
16161614

16171615
end
1618-
% here delete temp variables
16191616
%clear vars here
16201617
clearvars thisAerBsc355_raman thisAerBscStd355_raman thisAerExt355_raman thisAerExt355_raman_tmp thisAerExtStd355_raman
16211618
%% Raman method (532 nm)
@@ -1642,10 +1639,7 @@
16421639
sig607 = squeeze(sum(data.signal(flag607FR, :, flagClFre), 3));
16431640
bg607 = squeeze(sum(data.bg(flag607FR, :, flagClFre), 3));
16441641

1645-
%[thisAerExt355_raman, thisAerExtStd355_raman]=pollyRamanExt_smart_MC(data.distance0, sig387, 355, 387, mExt355(iGrp,:), mExt387(iGrp,:), number_density(iGrp, :), PollyConfig.angstrexp, PollyConfig.smoothWin_raman_355, 'moving',15,bg387);
16461642
[thisAerExt532_raman,thisAerExtStd532_raman] = pollyRamanExt_smart_MC(data.distance0, sig607, 532, 607, mExt532(iGrp,:), mExt607(iGrp,:), number_density(iGrp, :), PollyConfig.angstrexp, PollyConfig.smoothWin_raman_532, 'moving',15,bg607);
1647-
%thisAerExt532_raman = pollyRamanExt(data.distance0, sig607, 532, 607, PollyConfig.angstrexp, data.pressure(iGrp, :), data.temperature(iGrp, :) + 273.17, PollyConfig.smoothWin_raman_532, 380, 70, 'moving');
1648-
%thisAerExtStd532_raman = pollyRamanExtStd(data.distance0, sig607, bg607, 532, 607, PollyConfig.angstrexp, data.pressure(iGrp, :), data.temperature(iGrp, :) + 273.17, PollyConfig.smoothWin_raman_532, 380, 70, 15);
16491643
data.aerExt532_raman(iGrp, :) = thisAerExt532_raman;
16501644
data.aerExtStd532_raman(iGrp, :) = thisAerExtStd532_raman;
16511645

@@ -1670,9 +1664,7 @@
16701664

16711665
thisAerExt532_raman_tmp = thisAerExt532_raman;
16721666
thisAerExt532_raman(1:hBaseInd532) = thisAerExt532_raman(hBaseInd532);
1673-
[thisAerBsc532_raman, ~] = pollyRamanBsc_smart(data.distance0, sig532, sig607, thisAerExt532_raman, PollyConfig.angstrexp, mExt532(iGrp,:), mBsc532(iGrp,:),mExt607(iGrp,:), mBsc607(iGrp,:), refH532, 532, PollyConfig.refBeta532, PollyConfig.smoothWin_raman_532, true, 532, 607);
1674-
%[thisAerBsc532_raman, ~] = pollyRamanBsc(data.distance0, sig532, sig607, thisAerExt532_raman, PollyConfig.angstrexp, mExt532(iGrp,:), mBsc532(iGrp,:), refH532, 532, PollyConfig.refBeta532, PollyConfig.smoothWin_raman_532, true);
1675-
thisAerBscStd532_raman = pollyRamanBscStd(data.distance0, sig532, bg532, sig607, bg607, thisAerExt532_raman, thisAerExtStd532_raman, PollyConfig.angstrexp, 0.2, mExt532(iGrp,:), mBsc532(iGrp,:), refH532, 532, PollyConfig.refBeta532, PollyConfig.smoothWin_raman_532, true);
1667+
[thisAerBsc532_raman, thisAerBscStd532_raman, ~] = pollyRamanBsc_smart_MC(data.distance0, sig532, sig607, thisAerExt532_raman, PollyConfig.angstrexp, mExt532(iGrp,:), mBsc532(iGrp,:),mExt607(iGrp,:), mBsc607(iGrp,:), refH532, PollyConfig.refBeta532, PollyConfig.smoothWin_raman_532, true, 532, 607, bg532, bg607, thisAerExtStd532_raman, sigma_angstroem, MC_count, 'monte-carlo');
16761668

16771669
% lidar ratio
16781670
[thisLR532_raman, thisLRStd532_raman] = pollyLR(thisAerExt532_raman_tmp, thisAerBsc532_raman, ...
@@ -1711,13 +1703,11 @@
17111703
sig607 = squeeze(sum(data.signal(flag607FR, :, flagClFre), 3));
17121704
bg607 = squeeze(sum(data.bg(flag607FR, :, flagClFre), 3));
17131705

1714-
%[thisAerExt355_raman, thisAerExtStd355_raman]= pollyRamanExt_smart_MC(data.distance0, sig387, 355, 387, mExt355(iGrp,:), mExt387(iGrp,:), number_density(iGrp, :), PollyConfig.angstrexp, PollyConfig.smoothWin_raman_355, 'moving',15,bg387);
1715-
[thisAerExt532_raman, thisAerExtStd532_raman] = pollyRamanExt_smart_MC(data.distance0, sig607, 532, 607, mExt532(iGrp,:), mExt607(iGrp,:), number_density(iGrp, :), PollyConfig.angstrexp, PollyConfig.smoothWin_raman_1064, 'moving',15,bg607);
1716-
% thisAerExt532_raman = pollyRamanExt(data.distance0, sig607, 532, 607, PollyConfig.angstrexp, data.pressure(iGrp, :), data.temperature(iGrp, :) + 273.17, PollyConfig.smoothWin_raman_1064, 380, 70, 'moving');
1717-
%thisAerExtStd532_raman = pollyRamanExtStd(data.distance0, sig607, bg607, 532, 607, PollyConfig.angstrexp, data.pressure(iGrp, :), data.temperature(iGrp, :) + 273.17, PollyConfig.smoothWin_raman_1064, 380, 70, 15);
1706+
[thisAerExt532_raman, thisAerExtStd532_raman] = pollyRamanExt_smart_MC(data.distance0, sig607, 532, 607, mExt532(iGrp,:), mExt607(iGrp,:), number_density(iGrp, :), PollyConfig.angstrexp, PollyConfig.smoothWin_raman_1064, 'moving',15,bg607);
17181707
thisAerExt1064_raman = thisAerExt532_raman / (1064/532).^PollyConfig.angstrexp;
17191708
data.aerExt1064_raman(iGrp, :) = thisAerExt1064_raman;
1720-
data.aerExtStd1064_raman(iGrp, :) = thisAerExtStd532_raman / (1064/532).^PollyConfig.angstrexp;
1709+
thisAerExtStd1064_raman = thisAerExtStd532_raman / (1064/532).^PollyConfig.angstrexp;
1710+
data.aerExtStd1064_raman(iGrp, :) = thisAerExtStd1064_raman;
17211711

17221712
if isnan(data.refHInd1064(iGrp, 1))
17231713
continue;
@@ -1740,9 +1730,7 @@
17401730

17411731
thisAerExt1064_raman_tmp = thisAerExt1064_raman;
17421732
thisAerExt1064_raman(1:hBaseInd1064) = thisAerExt1064_raman(hBaseInd1064);
1743-
[thisAerBsc1064_raman, ~] = pollyRamanBsc_smart(data.distance0, sig1064, sig607, thisAerExt1064_raman, PollyConfig.angstrexp, mExt1064(iGrp,:), mBsc1064(iGrp,:),mExt607(iGrp,:), mBsc607(iGrp,:), refH1064, 1064, PollyConfig.refBeta1064, PollyConfig.smoothWin_raman_1064, true, 1064, 607);
1744-
% [thisAerBsc1064_raman, ~] = pollyRamanBsc(data.distance0, sig1064, sig607, thisAerExt1064_raman, PollyConfig.angstrexp, mExt1064(iGrp,:), mBsc1064(iGrp,:), refH1064, 1064, PollyConfig.refBeta1064, PollyConfig.smoothWin_raman_1064, true);
1745-
thisAerBscStd1064_raman = pollyRamanBscStd(data.distance0, sig1064, bg1064, sig607, bg607, thisAerExt1064_raman, data.aerExtStd1064_raman(iGrp, :), PollyConfig.angstrexp, 0.2, mExt1064(iGrp,:), mBsc1064(iGrp,:), refH1064, 1064, PollyConfig.refBeta1064, PollyConfig.smoothWin_raman_1064, true);
1733+
[thisAerBsc1064_raman, thisAerBscStd1064_raman, ~] = pollyRamanBsc_smart_MC(data.distance0, sig1064, sig607, thisAerExt1064_raman, PollyConfig.angstrexp, mExt1064(iGrp,:), mBsc1064(iGrp,:),mExt607(iGrp,:), mBsc607(iGrp,:), refH1064, PollyConfig.refBeta1064, PollyConfig.smoothWin_raman_1064, true, 1064, 607, bg1064, bg607, thisAerExtStd1064_raman, sigma_angstroem, MC_count, 'monte-carlo');
17461734

17471735
% lidar ratio
17481736
[thisLR1064_raman, thisLRStd1064_raman] = pollyLR(thisAerExt1064_raman_tmp, thisAerBsc1064_raman, ...
@@ -1811,11 +1799,8 @@
18111799

18121800
thisAerExt355_RR_tmp = thisAerExt355_RR;
18131801
thisAerExt355_RR(1:hBaseInd355) = thisAerExt355_RR(hBaseInd355);
1814-
%[thisAerBsc355_raman, ~] = pollyRamanBsc_smart(data.distance0, sig355, sig387, thisAerExt355_raman, PollyConfig.angstrexp, mExt355(iGrp,:), mBsc355(iGrp,:),mExt387(iGrp,:), mBsc387(iGrp,:), refH355, 355, PollyConfig.refBeta355, PollyConfig.smoothWin_raman_355, true, 355, 387);
1815-
[thisAerBsc355_RR, ~] = pollyRamanBsc_smart(data.distance0, sig355, sig355RR, thisAerExt355_RR, PollyConfig.angstrexp, mExt355(iGrp,:), mBsc355(iGrp,:), mExt355(iGrp,:), mBsc355(iGrp,:), refH355, 355, PollyConfig.refBeta355, PollyConfig.smoothWin_raman_355, true,355, 355);
1816-
%[thisAerBsc355_RR, ~] = pollyRamanBsc(data.distance0, sig355, sig355RR, thisAerExt355_RR, PollyConfig.angstrexp, mExt355(iGrp,:), mBsc355(iGrp,:), refH355, 355, PollyConfig.refBeta355, PollyConfig.smoothWin_raman_355, true);
1817-
thisAerBscStd355_RR = pollyRamanBscStd(data.distance0, sig355, bg355, sig355RR, bg355RR, thisAerExt355_RR, data.aerExtStd355_RR(iGrp, :), PollyConfig.angstrexp, 0.2, mExt355(iGrp,:), mBsc355(iGrp,:), refH355, 355, PollyConfig.refBeta355, PollyConfig.smoothWin_raman_355, true);
1818-
1802+
[thisAerBsc355_RR, thisAerBscStd355_RR, ~] = pollyRamanBsc_smart_MC(data.distance0, sig355, sig355RR, thisAerExt355_RR, PollyConfig.angstrexp, mExt355(iGrp,:), mBsc355(iGrp,:), mExt355(iGrp,:), mBsc355(iGrp,:), refH355, PollyConfig.refBeta355, PollyConfig.smoothWin_raman_355, true, 355, 355, bg355, bg355RR, thisAerExtStd355_RR, sigma_angstroem, MC_count, 'monte-carlo');
1803+
18191804
% lidar ratio
18201805
[thisLR355_RR, thisLRStd355_RR] = pollyLR(thisAerExt355_RR_tmp, thisAerBsc355_RR, ...
18211806
'hRes', data.hRes, ...
@@ -1853,10 +1838,7 @@
18531838
sig532RR = squeeze(sum(data.signal(flag532RR, :, flagClFre), 3));
18541839
bg532RR = squeeze(sum(data.bg(flag532RR, :, flagClFre), 3));
18551840

1856-
%[thisAerExt355_raman, thisAerExtStd355_raman] = pollyRamanExt_smart_MC(data.distance0, sig387, 355, 387, mExt355(iGrp,:), mExt387(iGrp,:), number_density(iGrp, :), PollyConfig.angstrexp, PollyConfig.smoothWin_raman_355, 'moving',15,bg387);
18571841
[thisAerExt532_RR,thisAerExtStd532_RR] = pollyRamanExt_smart_MC(data.distance0, sig532RR, 532, 532, mExt532(iGrp,:), mExt532(iGrp,:),number_density(iGrp, :),PollyConfig.angstrexp, PollyConfig.smoothWin_raman_532, 'moving',15,bg532RR);
1858-
%thisAerExt532_RR = pollyRamanExt(data.distance0, sig532RR, 532, 532, PollyConfig.angstrexp, data.pressure(iGrp, :), data.temperature(iGrp, :) + 273.17, PollyConfig.smoothWin_raman_532, 380, 70, 'moving');
1859-
%thisAerExtStd532_RR = pollyRamanExtStd(data.distance0, sig532RR, bg532RR, 532, 532, PollyConfig.angstrexp, data.pressure(iGrp, :), data.temperature(iGrp, :) + 273.17, PollyConfig.smoothWin_raman_532, 380, 70, 15);
18601842
data.aerExt532_RR(iGrp, :) = thisAerExt532_RR;
18611843
data.aerExtStd532_RR(iGrp, :) = thisAerExtStd532_RR;
18621844

@@ -1881,11 +1863,8 @@
18811863

18821864
thisAerExt532_RR_tmp = thisAerExt532_RR;
18831865
thisAerExt532_RR(1:hBaseInd532) = thisAerExt532_RR(hBaseInd532);
1884-
[thisAerBsc532_RR, ~] = pollyRamanBsc_smart(data.distance0, sig532, sig532RR, thisAerExt532_RR, PollyConfig.angstrexp, mExt532(iGrp,:), mBsc532(iGrp,:), mExt532(iGrp,:), mBsc532(iGrp,:), refH532, 532, PollyConfig.refBeta532, PollyConfig.smoothWin_raman_532, true, 532, 532);
1885-
%[thisAerBsc355_raman, ~] = pollyRamanBsc_smart(data.distance0, sig355, sig387, thisAerExt355_raman, PollyConfig.angstrexp, mExt355(iGrp,:), mBsc355(iGrp,:),mExt387(iGrp,:), mBsc387(iGrp,:), refH355, 355, PollyConfig.refBeta355, PollyConfig.smoothWin_raman_355, true, 355, 387);
1886-
%[thisAerBsc532_RR, ~] = pollyRamanBsc(data.distance0, sig532, sig532RR, thisAerExt532_RR, PollyConfig.angstrexp, mExt532(iGrp,:), mBsc532(iGrp,:), refH532, 532, PollyConfig.refBeta532, PollyConfig.smoothWin_raman_532, true);
1887-
thisAerBscStd532_RR = pollyRamanBscStd(data.distance0, sig532, bg532, sig532RR, bg532RR, thisAerExt532_RR, data.aerExtStd532_RR(iGrp, :), PollyConfig.angstrexp, 0.2, mExt532(iGrp,:), mBsc532(iGrp,:), refH532, 532, PollyConfig.refBeta532, PollyConfig.smoothWin_raman_532, true);
1888-
1866+
[thisAerBsc532_RR, thisAerBscStd532_RR, ~] = pollyRamanBsc_smart_MC(data.distance0, sig532, sig532RR, thisAerExt532_RR, PollyConfig.angstrexp, mExt532(iGrp,:), mBsc532(iGrp,:), mExt532(iGrp,:), mBsc532(iGrp,:), refH532, PollyConfig.refBeta532, PollyConfig.smoothWin_raman_532, true, 532, 532, bg532, bg532RR, thisAerExtStd532_RR, sigma_angstroem, MC_count, 'monte-carlo');
1867+
18891868
% lidar ratio
18901869
[thisLR532_RR, thisLRStd532_RR] = pollyLR(thisAerExt532_RR_tmp, thisAerBsc532_RR, ...
18911870
'hRes', data.hRes, ...
@@ -1950,10 +1929,7 @@
19501929

19511930
thisAerExt1064_RR_tmp = thisAerExt1064_RR;
19521931
thisAerExt1064_RR(1:hBaseInd1064) = thisAerExt1064_RR(hBaseInd1064);
1953-
%[thisAerBsc532_RR, ~] = pollyRamanBsc_smart(data.distance0, sig532, sig532RR, thisAerExt532_RR, PollyConfig.angstrexp, mExt532(iGrp,:), mBsc532(iGrp,:), mExt532(iGrp,:), mBsc532(iGrp,:), refH532, 532, PollyConfig.refBeta532, PollyConfig.smoothWin_raman_532, true, 532, 532);
1954-
[thisAerBsc1064_RR, ~] = pollyRamanBsc_smart(data.distance0, sig1064, sig1064RR, thisAerExt1064_RR, PollyConfig.angstrexp, mExt1064(iGrp,:), mBsc1064(iGrp,:), mExt1058(iGrp,:), mBsc1058(iGrp,:), refH1064, 1064, PollyConfig.refBeta1064, PollyConfig.smoothWin_raman_1064, true, 1064,1064);
1955-
%[thisAerBsc1064_RR, ~] = pollyRamanBsc(data.distance0, sig1064, sig1064RR, thisAerExt1064_RR, PollyConfig.angstrexp, mExt1064(iGrp,:), mBsc1064(iGrp,:), refH1064, 1064, PollyConfig.refBeta1064, PollyConfig.smoothWin_raman_1064, true);
1956-
thisAerBscStd1064_RR = pollyRamanBscStd(data.distance0, sig1064, bg1064, sig1064RR, bg1064RR, thisAerExt1064_RR, data.aerExtStd1064_RR(iGrp, :), PollyConfig.angstrexp, 0.2, mExt1064(iGrp,:), mBsc1064(iGrp,:), refH1064, 1064, PollyConfig.refBeta1064, PollyConfig.smoothWin_raman_1064, true);
1932+
[thisAerBsc1064_RR, thisAerBscStd1064_RR, ~] = pollyRamanBsc_smart_MC(data.distance0, sig1064, sig1064RR, thisAerExt1064_RR, PollyConfig.angstrexp, mExt1064(iGrp,:), mBsc1064(iGrp,:), mExt1058(iGrp,:), mBsc1058(iGrp,:), refH1064, PollyConfig.refBeta1064, PollyConfig.smoothWin_raman_1064, true, 1064, 1058, bg1064, bg1064RR, thisAerExtStd1064_RR, sigma_angstroem, MC_count, 'monte-carlo');
19571933

19581934
% lidar ratio
19591935
[thisLR1064_RR, thisLRStd1064_RR] = pollyLR(thisAerExt1064_RR_tmp, thisAerBsc1064_RR, ...
@@ -2042,10 +2018,7 @@
20422018

20432019
thisAerExt355_NR_raman_tmp = thisAerExt355_NR_raman;
20442020
thisAerExt355_NR_raman(1:hBaseInd355) = thisAerExt355_NR_raman(hBaseInd355);
2045-
%[thisAerBsc355_raman, ~] = pollyRamanBsc_smart(data.distance0, sig355, sig387, thisAerExt355_raman, PollyConfig.angstrexp, mExt355(iGrp,:), mBsc355(iGrp,:),mExt387(iGrp,:), mBsc387(iGrp,:), refH355, 355, PollyConfig.refBeta355, PollyConfig.smoothWin_raman_355, true, 355, 387);
2046-
[thisAerBsc355_NR_raman, ~] = pollyRamanBsc_smart(data.distance0, sig355, sig387, thisAerExt355_NR_raman, PollyConfig.angstrexp, mExt355(iGrp,:), mBsc355(iGrp,:),mExt387(iGrp,:), mBsc387(iGrp,:), refH355_NR, 355, refBeta355_NF, PollyConfig.smoothWin_raman_NR_355, true,355, 387);
2047-
%[thisAerBsc355_NR_raman, ~] = pollyRamanBsc(data.distance0, sig355, sig387, thisAerExt355_NR_raman, PollyConfig.angstrexp, mExt355(iGrp,:), mBsc355(iGrp,:), refH355, 355, refBeta355, PollyConfig.smoothWin_raman_NR_355, true);
2048-
thisAerBscStd355_NR_raman = pollyRamanBscStd(data.distance0, sig355, bg355, sig387, bg387, thisAerExt355_NR_raman, thisAerExtStd355_NR_raman, PollyConfig.angstrexp, 0.2, mExt355(iGrp,:), mBsc355(iGrp,:), refH355_NR, 355, refBeta355_NF, PollyConfig.smoothWin_raman_NR_355, true);
2021+
[thisAerBsc355_NR_raman, thisAerBscStd355_NR_raman, ~] = pollyRamanBsc_smart_MC(data.distance0, sig355, sig387, thisAerExt355_NR_raman, PollyConfig.angstrexp, mExt355(iGrp,:), mBsc355(iGrp,:),mExt387(iGrp,:), mBsc387(iGrp,:), refH355_NR, refBeta355_NF, PollyConfig.smoothWin_raman_NR_355, true, 355, 387, bg355, bg387, thisAerExtStd355_NR_raman, sigma_angstroem, MC_count, 'monte-carlo');
20492022

20502023
% lidar ratio
20512024
[thisLR355_NR_raman, thisLRStd355_NR_raman] = pollyLR(thisAerExt355_NR_raman_tmp, thisAerBsc355_NR_raman, ...
@@ -2136,11 +2109,7 @@
21362109
thisAerExt532_NR_raman_tmp = thisAerExt532_NR_raman;
21372110
thisAerExt532_NR_raman(1:hBaseInd532) = thisAerExt532_NR_raman(hBaseInd532);
21382111

2139-
%thisAerBsc355_raman, ~] = pollyRamanBsc_smart(data.distance0, sig355, sig387, thisAerExt355_raman, PollyConfig.angstrexp, mExt355(iGrp,:), mBsc355(iGrp,:),mExt387(iGrp,:), mBsc387(iGrp,:), refH355, 355, PollyConfig.refBeta355, PollyConfig.smoothWin_raman_355, true, 355, 387);
2140-
[thisAerBsc532_NR_raman, ~] = pollyRamanBsc_smart(data.distance0, sig532, sig607, thisAerExt532_NR_raman, PollyConfig.angstrexp, mExt532(iGrp,:), mBsc532(iGrp,:),mExt607(iGrp,:), mBsc607(iGrp,:), refH532_NR, 532, refBeta532_NR, PollyConfig.smoothWin_raman_NR_532, true, 532, 607);
2141-
%[thisAerBsc532_NR_raman, ~] = pollyRamanBsc(data.distance0, sig532, sig607, thisAerExt532_NR_raman, PollyConfig.angstrexp, mExt532(iGrp,:), mBsc532(iGrp,:), refH532, 532, refBeta532, PollyConfig.smoothWin_raman_NR_532, true);
2142-
thisAerBscStd532_NR_raman = pollyRamanBscStd(data.distance0, sig532, bg532, sig607, bg607, thisAerExt532_NR_raman, thisAerExtStd532_NR_raman, PollyConfig.angstrexp, 0.2, mExt532(iGrp,:), mBsc532(iGrp,:), refH532_NR, 532, refBeta532_NR, PollyConfig.smoothWin_raman_NR_532, true);
2143-
2112+
[thisAerBsc532_NR_raman, thisAerBscStd532_NR_raman, ~] = pollyRamanBsc_smart_MC(data.distance0, sig532, sig607, thisAerExt532_NR_raman, PollyConfig.angstrexp, mExt532(iGrp,:), mBsc532(iGrp,:),mExt607(iGrp,:), mBsc607(iGrp,:), refH532_NR, refBeta532_NR, PollyConfig.smoothWin_raman_NR_532, true, 532, 607, bg532, bg607, thisAerExtStd532_NR_raman, sigma_angstroem, MC_count, 'monte-carlo');
21442113
% lidar ratio
21452114
[thisLR532_NR_raman, thisLRStd532_NR_raman] = pollyLR(thisAerExt532_NR_raman_tmp, thisAerBsc532_NR_raman, ...
21462115
'hRes', data.hRes, ...

0 commit comments

Comments
 (0)