Skip to content

POLIPHON 2 - first part #323

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
AthenaAugousta wants to merge 4 commits into
base: dev
Choose a base branch
from
+1,599 −3
Open

POLIPHON 2 - first part #323

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
3544606
POLIPHON 2 - first part
Feb 3, 2026
2e7befa
POLIPHON 2 - first part
Feb 9, 2026
36ea95a
POLIPHON 2 - first part
Feb 9, 2026
7565518
final changes for poliphon2
HolgerPollyNet Feb 11, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion lib/config/polly_global_config.json
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -314,5 +314,5 @@
"logbookFileName": "",
"radiosondeFolder": "",
"imgFormat": "png",
"prodSaveList": ["overlap", "aerProfFR", "aerProfNR", "aerProfOC", "aerAttBetaFR", "aerAttBetaOC", "aerAttBetaNR", "WVMR_RH", "volDepol", "quasiV1", "quasiV2", "TC", "TCV2", "cloudinfo", "poliphon_one","RCS"]
"prodSaveList": ["overlap", "aerProfFR", "aerProfNR", "aerProfOC", "aerAttBetaFR", "aerAttBetaOC", "aerAttBetaNR", "WVMR_RH", "volDepol", "quasiV1", "quasiV2", "TC", "TCV2", "cloudinfo", "poliphon_one","RCS", "poliphon_two"]
}
25 changes: 23 additions & 2 deletions lib/interface/picassoProcV3.m
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -56,6 +56,7 @@
% - 2024-08-28: GHK formalism for depol calculation implemented by Moritz Haarig
% - 2025-02-19: Smoothing added into meteo profiles, read all meteo data for HR products and interpolate, recalculation of signals using mean shot_number by Cristofer Jimenez
% - 2025-03-14: Compute and save attenuated backscatter co and cross polarized by Cristofer Jimenez
% - 2026-02-09: POLIPHON step 2 added
%
% .. Authors: - zhenping@tropos.de, jimenez@tropos.de, floutsi@tropos.de, haarig@tropos.de

Expand Down Expand Up @@ -3151,7 +3152,7 @@
for iGrp = 1:size(clFreGrps, 1)

% Angstroem exponent 355-532 (based on parameters by Klett method)
if (~ isnan(data.aerExt355_NR_klett(iGrp, 60))) && (~ isnan(data.aerExt355_NR_klett(iGrp, 60))) %check what the 60 mean HB
if (~ isnan(data.aerExt355_NR_klett(iGrp, 60))) && (~ isnan(data.aerExt532_NR_klett(iGrp, 60))) %check what the 60 mean HB
[thisAE_Bsc_355_532_NR_klett, thisAEStd_Bsc_355_532_NR_klett] = pollyAE(data.aerBsc355_NR_klett(iGrp, :), zeros(size(data.height)), data.aerBsc532_NR_klett(iGrp, :), zeros(size(data.height)), 355, 532, PollyConfig.smoothWin_klett_NR_532);
data.AE_Bsc_355_532_NR_klett(iGrp, :) = thisAE_Bsc_355_532_NR_klett;
data.AEStd_Bsc_355_532_NR_klett(iGrp, :) = thisAEStd_Bsc_355_532_NR_klett;
Expand Down Expand Up @@ -3287,6 +3288,16 @@

print_msg('Finish.\n', 'flagTimestamp', true);

%% POLIPHON (2-step)
[data.POLIPHON2] = poliphon_two ...
(data.aerBsc355_klett, data.pdr355_klett, ...
data.aerBsc532_klett, data.pdr532_klett, data.aerBsc1064_klett, data.pdr1064_klett,...
data.aerBsc355_raman, data.pdr355_raman, data.aerBsc532_raman, data.pdr532_raman,...
data.aerBsc1064_raman, data.pdr1064_raman);

print_msg('Finish. \n', 'flagTimestamp', true);


%% Signal status
data.SNR = NaN(size(data.signal));
for iCh = 1:size(data.signal, 1)
Expand Down Expand Up @@ -5096,6 +5107,16 @@
%end
end

case 'poliphon_two'
if PicassoConfig.flagSaveProfiles
print_msg('--> start saving POLIPHON 2 products.\n', 'flagSimpleMsg', true, 'flagTimestamp', true);
% try
pollySavePOLIPHON2(data, data.POLIPHON2);
print_msg('--> finish!\n', 'flagSimpleMsg', true, 'flagTimestamp', true);
% catch
% print_msg('--> WARNING, could not save POLIPHON 2 products.\n', 'flagSimpleMsg', true, 'flagTimestamp', true);
end

otherwise
warning('Unknow product %s', PollyConfig.prodSaveList{iProd});
end
Expand Down Expand Up @@ -5266,4 +5287,4 @@
end
%% Clean
% clear;
end
end
Loading