CL_FFT3D_spectra_2D_peak_line.m

%% FFT 3D spectra
%% 2D peak line fft
%%
setlatex

%cd G:\orange_backup % any dir
%cd /Volumes/GROUP_BLACK/orange_backup/
clear all

%Gi=[0.05 0.1 0.2 0.3 0.4 0.5 1 1.25 3]; %All the G

Gi=[ 1.0];
mi=[1.5 1.7 2.0 2.2 2.5];
ri='*';

Gi=[0.4]
mi=[1.5]
ri=141
for i3=Gi
    for i4=mi

        close all
        clearvars -except i3 i4 Gi mi ri
        %sys=computer;
        sysdir
        gstr=num2str(i3,['%10.1','f']);
        mstr=num2str(i4,['%10.1','f']);
        rstr=num2str(ri);
        list0=dir(['g',gstr,'*m',mstr,'*r',rstr]);
        if isempty(list0); continue; end
        f1pmax=2.5;
        %A1din = input('input: digits after decimal point for A_1 legend? (default is 3 digits)','s');
        %if size(A1din)==0;
            A1din= '3';
        %end

        maximumA1=inf;

        A2switch=false; % A2=true or A2dA1=false

        mrswitch=0; %mass ratio is a variable or not

        n1=1;%starting number of files
        disp(['starting number= ',num2str(n1)])

        n0=length(list0);
        s=0;sp=0;
        countl=0;
        clear gamall aryaxismax aryaxismin arxaxismax arxaxismin
        for ai=n1:n0
            close all
            

            

            str=list0(ai).name;
            str(str=='r')=[' '];
            str(str=='u')=[];
            str(str=='n')=[];
            str(str=='_')=[];
            str(str=='z')=[];
            str(str=='l')=[' '];
            str1=str;
            str(str=='g')=' ';
            str(str=='a')=' ';
            str(str=='m')=' ';
            gamr=str2num(str);
            gamall(1:length(gamr),ai)=gamr;
            A1=gamr(2);
%             if A1>maximumA1
%                 continue
%             end
            countl=countl+1;
            A2m=[];f2m=[];f1temp=[]; %initialise A2m and f2m
            
            cd (list0(ai).name)
            listf=dir('0*');
            if isempty(listf); continue; end
            n=length(listf);
            
            clear A2dA1 A2p A2peak f1 f2p f2peak h
            for f1i=1:n
                dirnamef=[list0(ai).name,'/',listf(f1i).name,'/'];
                disp(dirnamef);
                f1(f1i)=str2double(listf(f1i).name);
                T1=1/f1(f1i);

    %             if A1==0.1 && f(i)==0.4
    %                 stop
    %             else
    %                 continue
    %             end

                cd(listf(f1i).name)
                %[t,CL]=textscanty;
                [t,CL]=textscan_ft(A1,f1(f1i));
                [tt,y]=textscanty;
                
                cd ../
                [t,CL]=autocut(f1(f1i),t,CL,tt,y);
                
                [A2p(f1i),f2p(f1i),A2, f2]=Contfft(t,CL,f1(f1i));
                
                f2(f2>f1pmax)=[];
                A2=A2(1:numel(f2));
                A2m(f1i,1:numel(A2))=A2; f2m(f1i,1:numel(f2))=f2;
                
                if A2switch==false
                    A2dA1(f1i)=A2p(f1i)/A1;
                end
                clear ymax ymin yamp yampds t v02 y ;
            end
            %% 3D spectra
            
          
            Lf1=length(f1);
            figure(1)
            %waterfall(f1,[],A2m)
            hold on
            for itemp=1:Lf1
                f1temp(1:length(f2m(itemp,:)))=f1(itemp);
                h(itemp)=plot3(f2m(itemp,:),f1temp,A2m(itemp,:),'color','k','LineWidth',1.5);
                %h(itemp)=plot3(f1temp,f2m(itemp,:),A2m(itemp,:),'color','k');
                clear f1temp
            end
            grid on; %grid minor
            
            xlim([0 2.5]);ylim([0 2.5])
            xlabel('$f_{CY2}/f_n$','Interpreter','LaTex','FontSize',22);
            ylabel('$f_1/f_n$','Interpreter','LaTex','FontSize',22);
            zlabel('$C_{YA2}$','Interpreter','LaTex','FontSize',22);
            
            campos([12.5839  -11.2161    0.5939*300])
            [A2peak, IA2peak]=max(A2m,[],2);
            for itemp=1:Lf1
                f2peak(itemp)=f2m(itemp,IA2peak(itemp));
            end

            hp3=plot3(f2peak,f1,zeros(1,Lf1),'color','r','LineStyle','--','LineWidth',1.5);
            %zlim([0 2.5])
            cd ../
            %
            
            temp=['CL_FFT_spectra_fig/',list0(ai).name];
            mkdir(temp)
            cd (temp)
            temp=[list0(ai).name,'_CL_3DFFT'];
            printpdf(gcf,[temp,'.eps']);
            saveas(gcf,[temp,'.fig'])
            %% 2D f1-f2
            hold off
            close all
            figure(1)
            hp3=plot(f1,f2peak,'color','r','LineStyle','-','LineWidth',1.5);
            grid on; grid minor;
            %set(gca,'XTick',0:0.2:2.5);
            xlabel('$f_1/f_n$','Interpreter','LaTex','FontSize',20);
            ylabel('$f_{CY2}/f_n$','Interpreter','LaTex','FontSize',20);
            setfigure;
            
            printpdf(gcf,[list0(ai).name,'_2Df1_f2peak.eps']);
            
            cd ../;cd ../;
            
        end


    end
end