load lcgt_optmat.mat
%opt = lcgt; % opt can be changed in extractTF so we must not reload lcgt.

%%%%%%%%%%%%%%
% SYM

freq = 100;
idx = find(f>freq, 1 );

nCARM_phi = real( atan(CARM_nCARM_Q(idx) / CARM_nCARM_P(idx)) );
disp(['nCARM_PHASE = ', num2str(180/pi * nCARM_phi), ';']);

%%%%%%%%%%%%%%
% POBS
%
% SB1: PRCL
% SB1_Q: MICH
% SB2: SRCL
%

freq = 100;
idx = find(f>freq, 1 );

nPRCL_phi = real( atan(PRCL_nPRCL_Q(idx) / PRCL_nPRCL_P(idx)) );
disp(['nPRCL_PHASE = ', num2str(180/pi * nPRCL_phi), ';']);

nMICH_phi = real( atan(MICH_nMICH_Q(idx) / MICH_nMICH_P(idx)) );
disp(['nMICH_PHASE = ', num2str(180/pi * nMICH_phi), ';']);

nSRCL_phi = real( atan(SRCL_nSRCL_Q(idx) / SRCL_nSRCL_P(idx)) );
disp(['nSRCL_PHASE = ', num2str(180/pi * nSRCL_phi), ';']);

% save results
%save lcgt_optmat_tuned.mat f *_PP *_QQ *_phi *ASY_DC
save lcgt_optmat_tuned.mat