function plotGouyWFS(phi,wfs,probeNames,iPrb,portNames,cDrvNames,saveDir);
% plot Gouy phase dependance of WFS signal and save
% phi: Gouy phase [rad]
%
% Author: Yuta Michimura
colororder = [0.0   0.0   1.0;
              1.0   0.0   0.0;
              0.0   0.7   1.0;
              1.0   0.7   0.0;
              0.3   0.3   0.3
              0.0   1.0   0.0
              0.0   0.5   0.0
              0.5   0.5   0.0
              1.0   0.0   1.0
              0.7   0.0   0.5
              0.5   0.1   0.1];

ph=phi/pi*180;
Nprb=length(probeNames);
Nport=length(portNames);

for kk=1:Nprb
    for jj=1:Nport
        if strfind(probeNames{kk},portNames{jj}) == 1
            wfssignal=squeeze(wfs(iPrb(kk),:,:));
            figure();
            set(gcf,'DefaultAxesColorOrder',colororder)
            plot(ph,real(wfssignal),'.')
            legend(cDrvNames);
            title(['WFS signal at ',strrep(probeNames{kk},'_A',' ')]);
            xlabel('Gouy phase [deg]');
            ylabel('WFS signal intensity [W/rad/sqrt(2/pi)]');
            xlim([-180,180]);
            set(gca,'XTick',linspace(-180,180,9));
            grid on;        
            saveas(gcf,[saveDir,'WFS_',strrep(probeNames{kk},'_A','_'),'.png']);
        end
    end
end