    linetypes={'-','--^','--o','--d'};
    Expheads={'EQ1M40_','EQ1M40lon60180_','EQ1M40lon180300_','EQ1M40lon30060_'};
    appends={'_2','_2','',''};
    labels={'circumglobal','60-180E','180-300E','300-60E'};
    Amps=[0,1,2,3,4,5,6,7,10];
    
    Uk1s=zeros(length(Expheads),length(Amps));
    Cov_exit=zeros(length(Expheads),length(Amps));
    for iexp=1:length(Expheads)
    Exphead=char(Expheads(iexp));
    linetype=char(linetypes(iexp));
    ExpName=cell(length(Amps),1);
    for ia=1:length(Amps)
    ExpName(ia)={strcat(Exphead,sprintf('%d',Amps(ia)))};
    end
    file2=strcat('./All_WA1Jan35_',ExpName,appends(iexp),'.nc');
    for ifile=1:length(file2)
    Utropo__=double(ncread(char(file2(ifile)),'U'));
    lat=double(squeeze(ncread(char(file2(ifile)),'lat')));
    lev=double(squeeze(ncread(char(file2(ifile)),'lev')));
    lon=double(squeeze(ncread(char(file2(ifile)),'lon')));
    lat0=45.;
    lev0=400.;
    [~,ilat]=min(abs(lat-lat0));
    [~,ilev]=min(abs(lev-lev0));
    r1=-5;
    r2=5;
    Utropolon=smoothn(squeeze(mean(Utropo__(:,ilat+r1:ilat+r2,ilev-1:ilev+1),3)),'robust');
    Uk1s(iexp,ifile)=max(abs(sum(Utropolon.*repmat(exp(1i*lon/180*pi),[1,r2-r1+1]),1))/length(lon));
    [Cov_exit(iexp,ifile),imax_conv]=max(max(Utropolon(1:end-1,:)-Utropolon(2:end,:)));
    end
    
    if length(ExpName)>5
    figure(9)
    if iexp==1
        clf
    end
    hold on
    plot(Amps,Uk1s(iexp,:)-Uk1s(iexp,1),strcat('k',linetype),'LineWidth',2)
    xlabel('MJO Amplitude (K/day)')
    title('k=1 zonal asymmetry amplitude of U wind')
    set(gca,'FontSize',18)
    end
    end
    hold on
    plot(Amps,sum(Uk1s(2:4,:)-Uk1s(2:4,1),1),'b','LineWidth',2)
    labels(end+1)={'sum'};
    legend(labels,'Location','southwest')
    