clear all;close all;
x=['UVQ_1980.mat';'UVQ_1981.mat';'UVQ_1982.mat';'UVQ_1983.mat';'UVQ_1984.mat';...,
    'UVQ_1985.mat';'UVQ_1986.mat';'UVQ_1987.mat';'UVQ_1988.mat';'UVQ_1989.mat';...,
    'UVQ_1990.mat';'UVQ_1991.mat';'UVQ_1992.mat';'UVQ_1993.mat';'UVQ_1994.mat';...,
    'UVQ_1995.mat';'UVQ_1996.mat';'UVQ_1997.mat';'UVQ_1998.mat';'UVQ_1999.mat';...,
    'UVQ_2000.mat';'UVQ_2001.mat';'UVQ_2002.mat'; 'UVQ_2003.mat'; 'UVQ_2004.mat';...,
    'UVQ_2005.mat';'UVQ_2006.mat';'UVQ_2007.mat';'UVQ_2008.mat';'UVQ_2009.mat'];
x1=NaN(length(x),151,241,92,'single'); %x1=years,151,241,days
y1=NaN(length(x),151,241,92,'single'); %y1=years,151,241,days
for i=1:length(x)
    fprintf('STEP 1.%d\n',i);
    cd narr_data\NARR_UVQ_BIG\;
    load(x(i,:),'newVQ','lat','lon','newUQ');
    cd ..\..
    x1(i,:,:,:)=squeeze(nmean(newUQ(:,:,:,:),4))/(60*60*3);
    y1(i,:,:,:)=squeeze(nmean(newVQ(:,:,:,:),4))/(60*60*3);
    clear newVQ newUQ
end;

vspeed=NaN(30,151,241,92,'single');
vspeed(:,:,:,:)=sqrt(y1(:,:,:,:).^2+x1(:,:,:,:).^2); %vspeed=yrs,151,241,days
vsmax=max(vspeed(:));
vscale=double(vsmax/500);

avg=NaN(151,241,92,'single'); %avg=151,241,days
avgx=NaN(151,241,92,'single'); %avg=151,241,days
avgy=NaN(151,241,92,'single'); %avg=151,241,days
avg(:,:,:)=squeeze(nmean(vspeed(:,:,:,:),1));
avgx(:,:,:)=squeeze(nmean(x1(:,:,:,:),1));
avgy(:,:,:)=squeeze(nmean(y1(:,:,:,:),1));

anom=NaN(30,151,241,92,'single'); %year,151,241,days
anomx=NaN(30,151,241,92,'single'); %year,151,241,days
anomy=NaN(30,151,241,92,'single'); %year,151,241,days

for i=1:30
    fprintf('STEP 2.%d\n',i);
    for j=1:92
        anom(i,:,:,j)=squeeze(vspeed(i,:,:,j))-squeeze(avg(:,:,j));
        anomx(i,:,:,j)=squeeze(x1(i,:,:,j))-squeeze(avgx(:,:,j));
        anomy(i,:,:,j)=squeeze(y1(i,:,:,j))-squeeze(avgy(:,:,j));
    end
end

load e90events.mat
type=e90r; %choose r for ridge... or t for trough
% c=NaN(length(type),6,length(P),29,57);
for i=1:length(type)
    fprintf('run 3.%d\n',i);
    yr=type(i,1)-1979;
    a=type(i,2)-4;
    b=type(i,2)+1;
    c(i,:,:,:)=squeeze(anom(yr,:,:,a:b));%c=events,151,241,day
    clear a b yr
end;
c_avg(:,:,:)=squeeze(nmean(c(:,:,:,:),1)); %c_avg=151,241,day

tc=NaN(151,241,6,'single');
for i=1:6
    fprintf('run 4.%d\n',i);
    for j=1:151
        for k=1:241
            tc(j,k,i)=ttest(c(:,j,k,i)); %tc=151,241,day
        end
    end
end

ff=find(tc(:,:,:)==0);
tc_avg=c_avg;
tc_avg(ff)=NaN;

cx1=NaN(length(type),151,241,6,'single');
cy1=NaN(length(type),151,241,6,'single');
for i=1:length(type)
    fprintf('run 5.%d\n',i);
    yr=type(i,1)-1979;fprintf('yr \n');
    a=type(i,2)-4;fprintf('a \n');
    b=type(i,2)+1;fprintf('b \n');
    cx1(i,:,:,:)=squeeze(x1(yr,:,:,a:b));fprintf('cx1 \n');%c=events,151,241,day
    cy1(i,:,:,:)=squeeze(y1(yr,:,:,a:b));fprintf('cy1 \n');%c=events,151,241,day
    clear a b yr
end;

cx1_avg(:,:,:)=squeeze(nmean(cx1(:,:,:,:),1)); %c_avg=151,241,day
cy1_avg(:,:,:)=squeeze(nmean(cy1(:,:,:,:),1)); %c_avg=151,241,day


a=squeeze(min(min(min(c_avg(:,:,:)))));
b=squeeze(max(max(max(c_avg(:,:,:)))));
if abs(a)>abs(b)
    lim=abs(a);
else
    lim=abs(b);
end

%START HERE AGAIN
for j = 1:6
    fprintf('PLOTTING\n');
    subplot(2,3,j);
    plotANOM(squeeze(tc_avg(:,:,j)),lon,lat);
    hold on;
    clim=[-1*lim lim];
    caxis(clim);
    m_quiver(lon(3:4:151,3:4:241),lat(3:4:151,3:4:241),...,
        squeeze(cx1_avg(3:4:151,3:4:241,j)),...,
        squeeze(cy1_avg(3:4:151,3:4:241,j)),vscale,'k');
    if j==5
        filename2=char(cellstr([' Day **' num2str(j) '**']));
    else
        filename2=char(cellstr([' Day ' num2str(j)]));
    end;
    title(filename2,'fontsize',16);
end;
subplotspace('vertical',-40);
subplotspace('horizontal',-15);
cd temp\
screen_size = get(0, 'ScreenSize');
set(gcf,'paperpositionmode','auto')
f1 = figure(1);
set(f1, 'Position', [0 0 screen_size(3) screen_size(4) ] );
h=colorbar;
ylabel(h,'Moisture Flux Anomaly (kg/m/s)','fontsize',14)
set(h,'position',[0.91 0.275 0.02 0.50])
% text(-5,4,['Series of Meridional Moisture Flux Anomalies for Trough Events'],'fontsize',24);
suptitle('Series of Moisture Flux Anomalies for RIDGE Events (with VECTORS) [STATSIG]  ');
saveas(gcf,['RIDGE Events UVQ Anom with Vectors [STATSIG]  '],'png');
clear a b h;
clf;close;
cd ..