ev_ext_MJO_EN_bs743.m

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%          Calculating the frequency of extreme events over SA            %
%                      in each MJO phase for DJF                          %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%for El Nino:
%
clc
clear
%
cd G:/Lais_jasmin/u-bs743/data/ENSO
%
load MJO; clear amp_LN amp_NT phase_LN phase_NT;
%
cd G:/Lais_jasmin/u-bs743/data/ev_ext
%
load ev_extr_MJO; 
%
%Putting NaN in ev ext when the MJO is inactive (amplitude < 1)
%
for j = 1:I*J
for i = 1:T
if amp_EN(i) < 1
ev_extr(i,j) = NaN;
end
end
end
%
for j = 1:I*J
for i = 1:T
if isnan(amp_EN(i))
ev_extr(i,j) = NaN;
end
end
end
%
%Separating preca according to the MJO phases:
%
%MJO phase 1-------------------
%
ev_extr_MJO_1 = ev_extr;
%
for j = 1:I*J
for i = 1:T
    if phase_EN(i) ~= 1
        ev_extr_MJO_1(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(phase_EN(:) ~= 1); b = find(amp_EN(:) < 1); 
%
c = [a ; b]; c = unique(c);
%
for i = c(:)
ev_extr_MJO_1(c,:) = [];
end
%
%MJO phase 2-------------------
%
ev_extr_MJO_2 = ev_extr;
%
for j = 1:I*J
for i = 1:T
    if phase_EN(i) ~= 2
        ev_extr_MJO_2(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(phase_EN(:) ~= 2); b = find(amp_EN(:) < 1); 
%
c = [a ; b]; c = unique(c);
%
for i = c(:)
ev_extr_MJO_2(c,:) = [];
end
%
%MJO phase 3-------------------
%
ev_extr_MJO_3 = ev_extr;
%
for j = 1:I*J
for i = 1:T
    if phase_EN(i) ~= 3
        ev_extr_MJO_3(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(phase_EN(:) ~= 3); b = find(amp_EN(:) < 1); 
%
c = [a ; b]; c = unique(c);
%
for i = c(:)
ev_extr_MJO_3(c,:) = [];
end
%
%MJO phase 4-------------------
%
ev_extr_MJO_4 = ev_extr;
%
for j = 1:I*J
for i = 1:T
    if phase_EN(i) ~= 4
        ev_extr_MJO_4(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(phase_EN(:) ~= 4); b = find(amp_EN(:) < 1); 
%
c = [a ; b]; c = unique(c);
%
for i = c(:)
ev_extr_MJO_4(c,:) = [];
end
%
%MJO phase 5-------------------
%
ev_extr_MJO_5 = ev_extr;
%
for j = 1:I*J
for i = 1:T
    if phase_EN(i) ~= 5
        ev_extr_MJO_5(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(phase_EN(:) ~= 5); b = find(amp_EN(:) < 1); 
%
c = [a ; b]; c = unique(c);
%
for i = c(:)
ev_extr_MJO_5(c,:) = [];
end
%
%MJO phase 6-------------------
%
ev_extr_MJO_6 = ev_extr;
%
for j = 1:I*J
for i = 1:T
    if phase_EN(i) ~= 6
        ev_extr_MJO_6(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(phase_EN(:) ~= 6); b = find(amp_EN(:) < 1); 
%
c = [a ; b]; c = unique(c);
%
for i = c(:)
ev_extr_MJO_6(c,:) = [];
end
%
%MJO phase 7-------------------
%
ev_extr_MJO_7 = ev_extr;
%
for j = 1:I*J
for i = 1:T
    if phase_EN(i) ~= 7
        ev_extr_MJO_7(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(phase_EN(:) ~= 7); b = find(amp_EN(:) < 1); 
%
c = [a ; b]; c = unique(c);
%
for i = c(:)
ev_extr_MJO_7(c,:) = [];
end
%
%MJO phase 8-------------------
%
ev_extr_MJO_8 = ev_extr;
%
for j = 1:I*J
for i = 1:T
    if phase_EN(i) ~= 8
        ev_extr_MJO_8(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(phase_EN(:) ~= 8); b = find(amp_EN(:) < 1); 
%
c = [a ; b]; c = unique(c);
%
for i = c(:)
ev_extr_MJO_8(c,:) = [];
end
%
%Finding the number of days in each MJO phase in DJF:
%
[T1,n] = size(ev_extr_MJO_1); [T2,n] = size(ev_extr_MJO_2);
%
[T3,n] = size(ev_extr_MJO_3); [T4,n] = size(ev_extr_MJO_4);
%
[T5,n] = size(ev_extr_MJO_5); [T6,n] = size(ev_extr_MJO_6);
%
[T7,n] = size(ev_extr_MJO_7); [T8,n] = size(ev_extr_MJO_8);
%
n = [T1 T2 T3 T4 T5 T6 T7 T8];
%
cd G:/Lais_jasmin/u-bs743/data/ev_ext
%
save('ev_extr_MJO_EN','ev_extr_MJO_1','ev_extr_MJO_2',...
    'ev_extr_MJO_3','ev_extr_MJO_4','ev_extr_MJO_5','ev_extr_MJO_6',...
    'ev_extr_MJO_7','ev_extr_MJO_8','n','lat','lon','lat1','lon1',...
    'I','J','T','K','L','-v7.3');
%
cd G:/Lais_jasmin/u-bs743/scripts/prec
%