preca_MJO_be362.m

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%             Calculating the daily rainfall anomalies over SA            %
%                        in each MJO phase for DJF                        %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
path(path,'C:/Users/Laís/Desktop/tese/scripts/Plot');
path(path,'C:/Users/Laís/Desktop/tese/scripts/toolbox_geral');
path(path,'C:/Users/Laís/Desktop/tese/scripts/toolbox_ocean');
%
clc
clear
%
cd G:/for_lais;
%
ncdisp('be362.jan-dec_dmeans_ts.years1-30.rmm_indices.nc')
%360-day calendar. The year is made of 12 months, each of which has 30 days  
%10800 days in 30 years (30 * 360 = 10800)
%
amp = ncread('be362.jan-dec_dmeans_ts.years1-30.rmm_indices.nc',...
    'amplitude'); 
%
phase = ncread('be362.jan-dec_dmeans_ts.years1-30.rmm_indices.nc',...
    'phase'); 
%
phase(1:120) = NaN(120,1);% problem first 120 values = 2
%
cd G:/Lais_jasmin/u-be362/data
%
load preca
%
M = 360;%days in a year
%
amp = reshape(amp,M,K); phase = reshape(phase,M,K); 
%
%Reducing amp and phase for DJF:
%
amp = amp'; phase = phase';
%
amp = cat(2,amp(:,1:60),amp(:,331:end));
%
phase = cat(2,phase(:,1:60),phase(:,331:end));
%
amp = amp'; phase = phase';
%
amp = reshape(amp,T,1); phase = reshape(phase,T,1); 
%
preca = reshape(preca,I*J,T);
%
preca = preca';
%
preca_DJF = preca;%sample for all days in DJF for student t-test
%
save('preca_DJF','preca_DJF','-v7.3'); clear preca_DJF;
%
%Putting NaN in preca when the MJO is inactive (amplitude < 1)
%
for j = 1:I*J
for i = 1:T
if amp(i) < 1
preca(i,j) = NaN;
end
end
end
%
%Separating preca according to the MJO phases:
%
%MJO phase 1-------------------
%
preca_MJO_1 = preca;
%
for j = 1:I*J
for i = 1:T
    if phase(i) ~= 1
        preca_MJO_1(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(isnan(preca_MJO_1(:,1)));
%
for i = a(:)
preca_MJO_1(a,:) = [];
end
%
%MJO phase 2-------------------
%
preca_MJO_2 = preca;
%
for j = 1:I*J
for i = 1:T
    if phase(i) ~= 2
        preca_MJO_2(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(isnan(preca_MJO_2(:,1)));
%
for i = a(:)
preca_MJO_2(a,:) = [];
end
%
%MJO phase 3-------------------
%
preca_MJO_3 = preca;
%
for j = 1:I*J
for i = 1:T
    if phase(i) ~= 3
        preca_MJO_3(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(isnan(preca_MJO_3(:,1)));
%
for i = a(:)
preca_MJO_3(a,:) = [];
end
%
%MJO phase 4-------------------
%
preca_MJO_4 = preca;
%
for j = 1:I*J
for i = 1:T
    if phase(i) ~= 4
        preca_MJO_4(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(isnan(preca_MJO_4(:,1)));
%
for i = a(:)
preca_MJO_4(a,:) = [];
end
%
%MJO phase 5-------------------
%
preca_MJO_5 = preca;
%
for j = 1:I*J
for i = 1:T
    if phase(i) ~= 5
        preca_MJO_5(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(isnan(preca_MJO_5(:,1)));
%
for i = a(:)
preca_MJO_5(a,:) = [];
end
%
%MJO phase 6-------------------
%
preca_MJO_6 = preca;
%
for j = 1:I*J
for i = 1:T
    if phase(i) ~= 6
        preca_MJO_6(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(isnan(preca_MJO_6(:,1)));
%
for i = a(:)
preca_MJO_6(a,:) = [];
end
%
%MJO phase 7-------------------
%
preca_MJO_7 = preca;
%
for j = 1:I*J
for i = 1:T
    if phase(i) ~= 7
        preca_MJO_7(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(isnan(preca_MJO_7(:,1)));
%
for i = a(:)
preca_MJO_7(a,:) = [];
end
%
%MJO phase 8-------------------
%
preca_MJO_8 = preca;
%
for j = 1:I*J
for i = 1:T
    if phase(i) ~= 8
        preca_MJO_8(i,j) = NaN;
    end
end
end
%
%Eliminating NaN:
%
a = find(isnan(preca_MJO_8(:,1)));
%
for i = a(:)
preca_MJO_8(a,:) = [];
end
%
%Finding the number of days in each MJO phase in DJF:
%
[T1,n] = size(preca_MJO_1); [T2,n] = size(preca_MJO_2);
%
[T3,n] = size(preca_MJO_3); [T4,n] = size(preca_MJO_4);
%
[T5,n] = size(preca_MJO_5); [T6,n] = size(preca_MJO_6);
%
[T7,n] = size(preca_MJO_7); [T8,n] = size(preca_MJO_8);
%
n = [T1 T2 T3 T4 T5 T6 T7 T8];
%
save('preca_MJO','preca_MJO_1','preca_MJO_2','preca_MJO_3',...
    'preca_MJO_4','preca_MJO_5','preca_MJO_6','preca_MJO_7','preca_MJO_8',...
    'n','lat','lon','lat1','lon1','I','J','T','K','L','-v7.3');
%
cd G:/Lais_jasmin/u-be362/scripts