% Compute root mean square amplitude of AMO/PMO/NMO

clear

% CMIP5 All

load /crane/s1/smiller/CMIP5/tos/GISS-E2-R/HMO_out_2 AMO_NA_reg_lp_m AMO_det_lp PMO_NP_reg_lp_m PMO_det_lp

AMO_NA_reg_rms = rms(AMO_NA_reg_lp_m);
AMO_det_rms = rms(AMO_det_lp);

PMO_NP_reg_rms = rms(PMO_NP_reg_lp_m);
PMO_det_rms = rms(PMO_det_lp);

load /crane/s1/smiller/CMIP5/tos/GISS-E2-R/HMO_out_2 HMO_NH_reg_lp_m HMO_det_lp

HMO_NH_reg_rms = rms(HMO_NH_reg_lp_m);
HMO_det_rms = rms(HMO_det_lp);

% CMIP5 GISS

load /crane/s1/smiller/CMIP5/tos/GISS-E2-R/HMO_out_GISS AMO_NA_reg_lp PMO_NP_reg_lp

AMO_NA_reg_rms_GISS = rms(AMO_NA_reg_lp);

PMO_NP_reg_rms_GISS = rms(PMO_NP_reg_lp);

load /crane/s1/smiller/CMIP5/tos/GISS-E2-R/HMO_out_GISS HMO_NH_reg_lp

HMO_NH_reg_rms_GISS = rms(HMO_NH_reg_lp);

% CMIP AIE

load /crane/s1/smiller/CMIP5/tos/GISS-E2-R/HMO_out_2 AMO_NA_reg_aelp PMO_NP_reg_aelp

AMO_NA_reg_rms_ae = rms(AMO_NA_reg_aelp);

PMO_NP_reg_rms_ae = rms(PMO_NP_reg_aelp);

load /crane/s1/smiller/CMIP5/tos/GISS-E2-R/HMO_out_2 HMO_NH_reg_aelp 

HMO_NH_reg_rms_ae = rms(HMO_NH_reg_aelp);

% Compute bivariate regression correlations

load /crane/s1/smiller/CMIP5/tos/GISS-E2-R/HMO_out_2 HMO_NH_approx HMO_NH_reg_lp_m

HMO_corr = corr(HMO_NH_approx,HMO_NH_reg_lp_m);
HMO_corr = HMO_corr^2;

load /crane/s1/smiller/CMIP5/tos/GISS-E2-R/HMO_out_GISS HMO_NH_approx HMO_NH_reg_lp

HMO_corr_GISS = corr(HMO_NH_approx,HMO_NH_reg_lp);
HMO_corr_GISS = HMO_corr_GISS^2;

load /crane/s1/smiller/CMIP5/tos/GISS-E2-R/HMO_out_2 HMO_NH_approx_ae HMO_NH_reg_aelp

HMO_corr_ae = corr(HMO_NH_approx_ae,HMO_NH_reg_aelp);
HMO_corr_ae = HMO_corr_ae^2;

% Compute AMO-PMO correlation

AMO_PMO_det_corr = crosscorr(AMO_det_lp,PMO_det_lp,20);
AMO_PMO_det_corr = AMO_PMO_det_corr.^2;

AMO_PMO_reg_corr = crosscorr(AMO_NA_reg_lp_m,PMO_NP_reg_lp_m,20);
AMO_PMO_reg_corr = AMO_PMO_reg_corr.^2;

clear AMO_det_lp HMO_NH_reg_aelp HMO_NH_reg_lp HMO_NH_reg_lp_m HMO_det_lp PMO_NP_reg_aelp PMO_NP_reg_lp PMO_reg_lp_m PMO_det_lp
clear PMO_NP_reg_lp_m AMO_NA_reg_aelp AMO_NA_reg_lp AMO_NA_reg_lp_m HMO_NH_approx HMO_NH_approx_ae

save /crane/s1/smiller/CMIP5/tos/GISS-E2-R/AMO_PMO_rms